Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tracking error”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-G-JeP1-08: Dual Modality Verification for Respiratory Gating Using New Real- Time Tumor Tracking Radiotherapy System

Purpose: The respirato ry gating system combined the TrueBeam and a new real-time tumor-tracking radiotherapy system (RTRT) was installed. The RTRT system consists of two x-ray tubes and color image intensifiers. Using fluoroscopic images, the fiducial marker which was implanted near the tumor was tracked and was used as the internal surrogate for respiratory gating. The purposes of this study was to develop the verification technique of the respiratory gating with the new RTRT using cine electronic portal image device images (EPIDs) of TrueBeam and log files of the RTRT. Methods: A patient who underwent respiratory gated SBRT of the lung using the RTRT were enrolled in this study. For a patient, the log files of three-dimensional coordinate of fiducial marker used as an internal surrogate were acquired using the RTRT. Simultaneously, the cine EPIDs were acquired during respiratory gated radiotherapy. The data acquisition was performed for one field at five sessions during the course of SBRT. The residual motion errors were calculated using the log files (E{sub log}). The fiducial marker used as an internal surrogate into the cine EPIDs was automatically extracted by in-house software based on the template-matching algorithm. The differences between the the marker positions of cine EPIDs and digitally reconstructed radiograph were calculated (E{sub EPID}). Results: Marker detection on EPID using in-house software was influenced by low image contrast. For one field during the course of SBRT, the respiratory gating using the RTRT showed the mean ± S.D. of 95{sup th} percentile E{sub EPID} were 1.3 ± 0.3 mm,1.1 ± 0.5 mm,and those of E{sub log} were 1.5 ± 0.2 mm, 1.1 ± 0.2 mm in LR and SI directions, respectively. Conclusion: We have developed the verification method of respiratory gating combined TrueBeam and new real-time tumor-tracking radiotherapy system using EPIDs and log files.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Report on Concepts & Approaches for SSBD for eCHEM

The verification of special nuclear material (SNM) in spent fuel pyroprocessing is an important safeguards challenge. The detection of spontaneous fission (SF) neutrons from curium is an accepted, non-destructive technique that has been applied to verify special nuclear material (SNM) content in used fuel and other materials in the fuel cycle. The nuclear material accounting (NMA) technique at the Korea Atomic Energy Research Institute’s Reference Engineering-scale Pyroprocessing Facility (REPF) is based on the Cm balance technique. Several publications have demonstrated the safeguards benefit from using process monitoring (PM) on nuclear facilities as a complementary measure to NMA. More recently, this concept was expanded and preliminarily demonstrated for pyroprocessing. The concept of Signature Based Safeguards (SBS) is part of this expansion, and is built around the interpretation of input from various sensors in a declared facility coupled with complementary NMA methods to increase confidence and lower standard error inventory differences (SEID). The SBS methodology was conceptually developed and relies on near real time analysis of process monitoring data to detect material diversion complemented by robust containment and surveillance (C/S) measures. This work demonstrates one example of how the SBS framework can be used in the electrorefiner. In this SBS application, a combination of cyclic voltammetry (CV) and neutron counting is applied to track and monitor Pu mass balance. The main purpose of this experiment is to determine if meaningful information can be gained from CV measurements with regard to the Mg/Gd ratio. This data will be coupled with ICP-MS to verify Gd concentrations and analyzed for statistical significance. It is expected the CV data will register a significant change under the off-normal operating conditions. Knowing how to identify and interpret those changes may help inform how to target more traditional neutron counting methods, which could support a more efficient safeguards system. The experimental results will be compared with theoretical calculations and the ERAD simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Deterministic High-Fidelity Neutronics Simulation of Pebble Bed Reactors Using Pebble Tracking Transport

The pebble tracking transport (PTT) algorithm offers a high-fidelity deterministic approach for neutron transport for pebble bed reactors (PBRs). This approach requires the mesh for the active-core region to consist exclusively of tetrahedral elements, where each node in the pebble-packing region represents a pebble centroid. This paper investigates the application of PTT for full-scale PBRs, considering both the isothermal and the temperature-dependent core conditions. Macroscopic cross sections are generated using Serpent 2 full-core eigenvalue simulations where pebbles are grouped into disjoint subsets using machine learning. To minimize the need for individual cross-section sets for each pebble in the core, K-means clustering is used to group pebbles by temperature and neutronic environment parameters. Here, we compare the multiplication factor and power rate distributions between PTT simulations using the Griffin reactor physics software and reference solutions from Serpent 2. Our analysis shows that a full-core, high-fidelity PTT calculation produces accurate results with minimal local (pebblewise) errors. Additionally, timing results indicate that PTT simulations converge rapidly on modern supercomputing platforms.

Griffin↗
DOE OSTI2018

A new method to calculate unsteady particle kinematics and drag coefficient in a subsonic post-shock flow

In this paper, we introduce a new method (piecewise integrated dynamics equation fit, PIDEF) that uses the particle dynamics equation to determine unsteady kinematics and drag coefficient (C D ) for a particle in subsonic post-shock flow. The uncertainty of this method is assessed based on simulated trajectories for both quasi-steady and unsteady flow conditions. Traditional piecewise polynomial fitting (PPF) shows high sensitivity to measurement error and the function used to describe C D , creating high levels of relative error (>>1) when applied to unsteady shock-accelerated flows. The PIDEF method provides reduced uncertainty in calculations of unsteady acceleration and drag coefficient for both quasi-steady and unsteady flows. This makes PIDEF a preferable method over PPF for complex flows where the temporal response of C D is unknown. Finally, we apply PIDEF to experimental measurements of particle trajectories from 8-pulse particle tracking and determine the effect of incident Mach number on relaxation kinematics and drag coefficient of micron-sized particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Fast dynamic aperture optimization with forward-reversal integration

A fast dynamic aperture (DA) optimization method for storage rings has been developed through the use of reversal integration. Even if dynamical systems have an exact reversal symmetry, a numerical forward integration differs from its reversal. For a chaotic trajectory, cumulative round-off errors are scaled, which results in an exponential growth on the difference. The exponential effect is a generic chaos indicator which represents the sensitivity of the chaotic motion to its initial condition. The chaos indicator of the charged particle motion can be obtained by comparing the forward integrations of particle trajectories with corresponding reversals, a.k.a. “backward integrations.” The indicator is observable even through short-term particle tracking simulations. Therefore, adopting it as an objective function could speed up optimization. Finally, the DA of the National Synchrotron Light Source II storage ring, and another test diffraction-limited light source ring, were optimized using this method for the purpose of demonstration.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

TH-AB-202-12: The First Clinical Implementation of a Real-Time Six Degree of Freedom Tracking System During Radiation Therapy

Purpose: In current practice, imaging is typically performed prior to treatment; the cancer target motion during treatment is unknown. We present the first clinical implementation of real time Kilovoltage Intrafraction Monitoring (KIM) system which tracks the cancer target translational and rotational motions during treatment. Methods: KIM technology: KIM estimates the 3D position of the target tumour based on segmented 2D positions of the three implanted fiducials in each of the kV images (125 kV, 10 mA at 11 fps) taken continuously during the treatment arc. The 2D-3D target estimation is based on a probability distribution function, obtained during pre-treatment CBCT. Rotations about each axis with the centroid of the markers as the pivot were calculated using the iterative closest point algorithm in real time. Patient: A patient with prostate adenocarcinoma undergoing stereotactic body radiotherapy (SBRT) with 36.25 Gy delivered in 5 fractions (Varian TrueBeam, 6X, VMAT) was enrolled in the study. The trial complies with Australian ethical and regulatory standards. Results: Of the 5 fractions of treatment the patient received, KIM was utilised successfully in 4 fractions with 3 couch shifts due to large persistent prostate movements (>2mm for more than 5 seconds). KIM translational accuracy and precision in comparison with post treatment kV-MV triangulation are 0.28±0.59 mm, −0.19±0.25 mm and 0.23±0.69 mm in the Left-Right, Superior-Inferior and Anterior-Posterior directions, respectively. KIM rotational accuracy as compared with triangulation is: 0.429°±2.22°, −0.44°±4.7° and 0.06°±1.08° in the roll, pitch and yaw direction, respectively. Conclusion: The first six degree of freedom KIM system was successfully implemented clinically. The presented KIM system has sub-millimeters accuracy and precision in all three translational axes, and less than 1° of mean error in all three rotational axes. Acknowledgement: This work is supported by Cancer Australia grants APPXXX, APPYYY.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Dynamic Aperture Evaluation for the Hadron Storage Ring in the Electron-Ion Collider

The Electron-Ion Collider (EIC) is aiming at a design luminosity of 1e34 cm⁻²s⁻¹. To maintain such a high luminosity, both beams in the EIC need an acceptable beam lifetime in the presence of the beam-beam interaction. For this purpose, we carried out weak-strong element-by-element particle tracking to evaluate the long-term dynamic aperture for the hadron ring lattice design. We improved our simulation code SimTrack to treat some new lattice design features, such as radially offset on-momentum orbits, coordinate transformations in the interaction region, etc. In this article, we will present the preliminary dynamic aperture calculation results with β^{*}- function scan, radial orbit shift, crossing angle collision, and magnetic field errors.

Luo, Y.↗
DOE OSTI2019

Evaluation of MOPITT Version 7 joint TIR–NIR X CO retrievals with TCCON

Observations of carbon monoxide (CO) from the Measurements Of Pollution In The Troposphere (MOPITT) instrument aboard the Terra spacecraft were expected to have an accuracy of 10 % prior to the launch in 1999. Here we evaluate MOPITT Version 7 joint (V7J) thermal-infrared and near-infrared (TIR–NIR) retrieval accuracy and precision and suggest ways to further improve the accuracy of the observations. We take five steps involving filtering or bias corrections to reduce scatter and bias in the data relative to other MOPITT soundings and ground-based measurements. (1) We apply a preliminary filtering scheme in which measurements over snow and ice are removed. (2) We find a systematic pairwise bias among the four MOPITT along-track detectors (pixels) on the order of 3–4 ppb with a small temporal trend, which we remove on a global scale using a temporally trended bias correction. (3) Using a small-region approximation (SRA), a new filtering scheme is developed and applied based on additional quality indicators such as the signal-to-noise ratio (SNR). After applying these new filters, the root-mean-squared error computed using the local median from the SRA over 16 years of global observations decreases from 3.84 to 2.55 ppb. (4) We also use the SRA to find variability in MOPITT retrieval anomalies that relates to retrieval parameters. We apply a bias correction to one parameter from this analysis. (5) After applying the previous bias corrections and filtering, we compare the MOPITT results with the GGG2014 ground-based Total Carbon Column Observing Network (TCCON) observations to obtain an overall global bias correction. These comparisons show that MOPITT V7J is biased high by about 6 %–8 %, which is similar to past studies using independent validation datasets on V6J. When using TCCON spectrometric column retrievals without the standard airmass correction or scaling to aircraft (WMO scale), the ground- and satellite-based observations overall agree to better than 0.5 %. GEOS-Chem data assimilations are used to estimate the influence of filtering and scaling to TCCON on global CO and tend to pull concentrations away from the prior fluxes and closer to the truth. We conclude with suggestions for further improving the MOPITT data products.

Hedelius, Jacob K. (ORCID:0000000320257519)↗
DOE OSTI2021

A Computational Information Criterion for Particle-Tracking with Sparse or Noisy Data

Traditional probabilistic methods for the simulation of advection-diffusion equations (ADEs) often overlook the entropic contribution of the discretization, e.g., the number of particles, within associated numerical methods. Many times, the gain in accuracy of a highly discretized numerical model is outweighed by its associated computational costs or the noise within the data. Herein, we address the question of how many particles are needed in a simulation to best approximate and estimate parameters in one-dimensional advective-diffusive transport. To do so, we use the well-known Akaike Information Criterion (AIC) and a recently-developed correction called the Computational Information Criterion (COMIC) to guide the model selection process. Random-walk and mass-transfer particle tracking methods are employed to solve the model equations at various levels of discretization. Numerical results demonstrate that the COMIC provides an optimal number of particles that can describe a more efficient model in terms of parameter estimation and model prediction compared to the model selected by the AIC even when the data is sparse or noisy, the sampling volume is not uniform throughout the physical domain, or the error distribution of the data is non-IID Gaussian.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Harmonizing solar induced fluorescence across spatial scales, instruments, and extraction methods using proximal and airborne remote sensing: A multi-scale study in a soybean field

Solar-induced chlorophyll fluorescence (SIF) has been widely used to track vegetation photosynthesis at different scales ranging from in-situ measurements to satellite products. Airborne platforms sample SIF data at a spatial scale intermediate between in-situ and satellite, matching that of ground measurement (e.g. flux tower footprints and other field sampling), enabling us to explore causes of SIF variation and validate satellite-based SIF products. However, harmonizing SIF across sensors and platforms (correcting for systematic errors to yield a consistent, comparable SIF product) is challenging because SIF can be retrieved in different absorption windows, with different instruments and methods complicating the comparison between different observational levels (i.e., ground, airborne, satellites) and between sites equipped with different instruments with varying optical properties (spectral resolution and sampling intervals, spatial resolution). Additionally, the spatial and temporal variability of atmospheric properties can influence the retrieval of the weak SIF signal. Because of these complications, direct comparisons of airborne and ground SIF across scales are rarely attempted. Here, in this study, we combined airborne SIF data with simultaneous ‘ground truth’ data collected by stationary and mobile platforms in a soybean field in Nebraska, USA. In this effort, we tested several SIF extraction methods, including Fraunhofer Line Discrimination (FLD), improved Fraunhofer Line Discrimination (iFLD), Spectral Fitting Method (SFM), SpecFit, and a Singular Vector Decomposition (SVD) method. The SpecFit method was sensitive to the 715–740 nm water bands and removing the water bands in the fitting process yielded better agreement between the airborne and ground SIF spectra. Accurate estimation of the ground level downwelling irradiance obtained by ground measurements over a calibration target improved agreement between airborne and ground SIF retrievals at the O 2 A band, and allowed us to derive a SIF dataset with improved agreement across platforms and sampling scales. This experimental approach provided a method for generating comparable SIF signals across instruments, methods and platforms, which is critical to understanding the SIF-GPP relationship at different scales and to cross-validate the diversity of platforms used for satellite products calibration and validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Experimental Examination of Additively Manufactured Patterns on Structural Nuclear Materials for Digital Image Correlation Strain Measurements

Abstract Background There are a limited number of commercially available sensors for monitoring the deformation of materials in-situ during harsh environment applications, such as those found in the nuclear and aerospace industries. Such sensing devices, including weldable strain gauges, extensometers, and linear variable differential transformers, can be destructive to material surfaces being investigated and typically require relatively large surface areas to attach (> 10 mm in length). Digital image correlation (DIC) is a viable, non-contact alternative to in-situ strain deformation. However, it often requires implementing artificial patterns using splattering techniques, which are difficult to reproduce. Objective Additive manufacturing capabilities offer consistent patterns using programmable fabrication methods. Methods In this work, a variety of small-scale periodic patterns with different geometries were printed directly on structural nuclear materials (i.e., stainless steel and aluminum tensile specimens) using an aerosol jet printer (AJP). Unlike other additive manufacturing techniques, AJP offers the advantage of materials selection. DIC was used to track and correlate strain to alternative measurement methods during cyclic loading, and tensile tests (up to 1100 µɛ) at room temperature. Results The results confirmed AJP has better control of pattern parameters for small fields of view and facilitate the ability of DIC algorithms to adequately process patterns with periodicity. More specifically, the printed 100 μm spaced dot and 150 μm spaced line patterns provided accurate measurements with a maximum error of less than 2% and 4% on aluminum samples when compared to an extensometer and commercially available strain gauges. Conclusion Our results highlight a new pattern fabrication technique that is form factor friendly for digital image correlation in nuclear applications.

Novich, K. A. (ORCID:0000000204466022)↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2015

WE-AB-303-01: FEATURED PRESENTATION: A Dual-Detector Phase-Matched Digital Tomosynthesis (DTS) Imaging Scheme Using Aggregated KV and MV Projections for Intra-Treatment Lung Tumor Tracking

Purpose: To develop a dual-detector phase-matched DTS technique for continuous and fast intra-treatment lung tumor localization. Methods: Tumor localization accuracy of limited-angle DTS imaging is affected by low inter-slice resolution. The dual-detector DTS technique aims to overcome this limitation through combining orthogonally acquired beam’s eye view MV projections and kV projections for intra-treatment DTS reconstruction and localization. To aggregate the kV and MV projections for reconstruction, the MV projections were linearly converted to synthesize corresponding kV projections. To further address the lung motion induced localization errors, this technique uses respiratory phase-matching to match the motion information between on-board DTS and reference DTS to offset the adverse effects of motion blurriness in tumor localization.A study was performed using the CIRS008A lung phantom to simulate different on-board target variation scenarios for localization. The intra-treatment kV and MV acquisition was achieved through the Varian TrueBeam Developer Mode. Four methods were compared for their localization accuracy: 1. the proposed dual-detector phase-matched DTS technique; 2. the single-detector phase-matched DTS technique; 3. the dual-detector 3D-DTS technique without phase-matching; and 4. the single-detector 3D-DTS technique without phase-matching. Results: For scan angles of 2.5°, 5°, 10°, 20° and 30°, the dual-detector phase-matched DTS technique localized the tumor with average(±standard deviations) errors of 0.4±0.3 mm, 0.5±0.3 mm, 0.6±0.2 mm, 0.9±0.4 mm and 1.0±0.3 mm, respectively. The corresponding values of single-detector phase-matched DTS technique were 4.0±2.5 mm, 2.7±1.1 mm, 1.7±1.2 mm, 2.2±0.9 mm and 1.5±0.8 mm, respectively. The values of dual-detector 3D-DTS technique were 6.2±1.7 mm, 6.3±1.2 mm, 5.3±1.3 mm, 2.0±2.2 mm and 1.5±0.5 mm, respectively. And the values of single-detector 3D-DTS technique were 9.7±8.9 mm, 9.8±8.8 mm, 10.0±9.7 mm, 3.9±2.7 mm and 2.2±1.3 mm, respectively. Conclusion: The dual-detector phase-matched DTS technique substantially improves the tumor localization accuracy, which can be applied to real-time intra-treatment lung tumor localization. The research was funded by the National Institutes of Health Grant No. R01-CA184173 and a grant from Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

S2M: A Stochastic Simulation Model of Poliovirus Genetic State Transition

Modeling the molecular mechanisms that govern genetic variation can be useful in understanding the dynamics that drive genetic state transition in quasispecies viruses. For example, there is considerable interest in understanding how the relatively benign vaccine strains of poliovirus eventually revert to forms that confer neurovirulence and cause disease (ie, vaccine-derived poliovirus). This report describes a stochastic simulation model, S2M, which can be used to generate hypothetical outcomes based on known mechanisms of genetic diversity. S2M begins with predefined genotypes based on the Sabin-1 and Mahoney wild-type sequences, constructs a set of independent cell-based populations, and performs in-cell replication and cell-to-cell infection cycles while quantifying genetic changes that track the transition from Sabin-1 toward Mahoney. Realism is incorporated into the model by assigning defaults for variables that constrain mechanisms of genetic variability based roughly on metrics reported in the literature, yet these values can be modified at the command line in order to generate hypothetical outcomes driven by these parameters. To demonstrate the utility of S2M, simulations were performed to examine the effects of the rates of replication error and recombination and the presence or absence of defective interfering particles, upon reaching the end states of Mahoney resemblance (semblance of a vaccine-derived state), neurovirulence, genome fitness, and cloud diversity. Simulations provide insight into how modeled biological features may drive hypothetical outcomes, independently or in combination, in ways that are not always intuitively obvious.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Learning turbulent flows with generative models for super resolution and sparse flow reconstruction

Neural operators are promising surrogates for dynamical systems but when trained with standard L 2 losses they tend to oversmooth fine-scale turbulent structures. Here, we show that combining operator learning with generative modeling overcomes this limitation. We consider three practical turbulent-flow challenges where conventional neural operators fail: spatio-temporal super-resolution, forecasting, and sparse flow reconstruction. For Schlieren jet super-resolution, an adversarially trained neural operator (adv-NO) reduces the energy-spectrum error by 15 × while preserving sharp gradients at neural operator-like inference cost. For 3D homogeneous isotropic turbulence, adv-NO trained on only 160 timesteps from a single trajectory forecasts accurately for five eddy-turnover times and offers 114 × wall-clock speed-up at inference than the baseline diffusion-based forecasters, enabling near-real-time rollouts. For reconstructing cylinder wake flows from highly sparse Particle Tracking Velocimetry-like inputs, a conditional generative model infers full 3D velocity and pressure fields with correct phase alignment and statistics. These advances enable accurate reconstruction and forecasting at low compute cost, bringing near-real-time analysis and control within reach in experimental and computational fluid mechanics.

Fluid dynamics↗
DOE OSTI2015

Systematic non-LTE study of the −2.6 ≤ [Fe/H] ≤ 0.2 F and G dwarfs in the solar neighborhood. I. Stellar atmosphere parameters

We present atmospheric parameters for 51 nearby F and G dwarf and subgiant stars uniformly distributed over the −2.60<[Fe/H]<+0.20 metallicity range that is suitable for the Galactic chemical evolution research. Lines of iron in the two ionization stages, Fe i and Fe ii, were used to derive a homogeneous set of effective temperatures, surface gravities, iron abundances, and microturbulence velocities. Our spectroscopic analyses took advantage of employing high-resolution (R ≥ 60,000) Shane/Hamilton and Canada–France–Hawaii Telescope/ESPaDOnS observed spectra and non-LTE (NLTE) line formation for Fe i and Fe ii in the classical one-dimensional model atmospheres. The spectroscopic method was tested in advance with the 20 benchmark stars, for which there are multiple measurements of the infrared flux method effective temperature and their Hipparcos parallax error is less than 10%. We found NLTE abundances from lines of Fe i and Fe ii to be consistent within 0.06 dex for every benchmark star, when applying a scaling factor of S{sub H} = 0.5 to the Drawinian rates of inelastic Fe+H collisions. The obtained atmospheric parameters were checked for each program star by comparing its position in the log g–T{sub eff} plane with the theoretical evolutionary track of given metallicity and α-enhancement in the Yi et al. grid. Our final effective temperatures lie exactly in between the T{sub IRFM} scales of Alonso et al. and Casagrande et al., with a mean difference of +46 and −51 K, respectively. NLTE leads to higher surface gravity compared with that for LTE. The shift in log g is smaller than 0.1 dex for stars with [Fe/H] ⩾−0.75, T{sub eff} ≤ 5750 K, or log g ≥ 4.20. NLTE analysis is crucial for the very metal-poor turnoff and subgiant stars, for which the shift in log g between NLTE and LTE can be up to 0.5 dex. The obtained accurate atmospheric parameters will be used in the forthcoming papers to determine NLTE abundances of important astrophysical elements from lithium to europium and to improve observational constraints on the chemodynamical models of the Galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Modelling Curvilinear Beamline Effects on Beam- and Spin-Dynamics in the Fermilab Muon g-2 Storage Ring.

A model of the Fermilab Muon g-2 Experiment storage ring is developed with the BMAD simulation library for beam- and spin-dynamics studies. Particular focus is given to the electrostatic quadrupoles, the plates of which follow the curvature of the particle design trajectory of the magnetically uniform storage ring. To model these, standard electric multipoles and a novel method based on cylindrical harmonics are investigated and found to accurately match the electric field maps produced for the bulk and fringe regions of the curved quadrupoles respectively, despite these methods being specifically for straight beamlines. The toroidal multipole description is then developed, to include curvature and longitudinal field variations (to model the fringe field regions). Various combinations of the different descriptions (including use of the field maps) are used for the bulk and fringe fields of the quadrupoles, and the toroidal multipoles are found to provide the best agreement with the measured tunes of the storage ring. Two integrators are also developed, based on the method of Wu, Forest and Robin. The first of these describes motion in a uniform magnetic dipole field with arbitrary electric field, and produces results comparable to those of the adaptive-step fourth-order Runge–Kutta integrator implemented in BMAD, suitable for tracking over thousands of turns. The second is a generalisation of the first to arbitrary magnetic vector potentials superimposed on that of the uniform dipole, with arbitrary electric fields, and was also found to agree with the Runge–Kutta integrator. Simulations of antimuon storage for 2000 turns of the storage ring were performed, where losses were studied along with the evolution of beam polarisation. The storage losses, important for reducing the systematic error in the measurement of the muon g-2, were found to be largely due to collimator shape. The correlation between momentum and spin precession frequency reflect the effects expected from the electric field correction in the T-BMT equation, on the order of 10 -6 . Particle injection was also studied, for various injection angles, quadrupole strengths and strengths of injection kick, and in the cases where experimental data was provided, the model was found to be consistent with the experiment.

43 PARTICLE ACCELERATORS↗