SEARCH · Engineering Papers

Results for “FIT”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Did we hear the sound of the Universe boiling? Analysis using the full fluid velocity profiles and NANOGrav 15-year data

Abstract In this paper, we analyse sound waves arising from a cosmic phase transition where the full velocity profile is taken into account as an explanation for the gravitational wave spectrum observed by multiple pulsar timing array groups. Unlike the broken power law used in the literature, in this scenario the power law after the peak depends on the macroscopic properties of the phase transition, allowing for a better fit with pulsar timing array (PTA) data. We compare the best fit with that obtained using the usual broken power law and, unsurprisingly, find a better fit with the gravitational wave (GW) spectrum that utilizes the full velocity profile. Even more importantly, the thermal parameters that produce the best fit are quite different. We then discuss models that can produce the best-fit point and complementary probes using CMB experiments and searches for light particles in DUNE, IceCUBE-Gen2, neutrinoless doubleβ-decay, and forward physics facilities (FPF) at the LHC like FASERν, etc.

Astronomy & Astrophysics↗

DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

SDSS-IV MaStar: [α/Fe] for the MaNGA Stellar Library from synthetic model spectra

ABSTRACT We calculate the α-enhancement ratio [α/Fe] for the Mapping Nearby Galaxies at APO (MaNGA) Stellar Library (MaStar) while also fitting for the fundamental atmospheric parameters effective temperature, surface gravity, and metallicity – Teff, log g, [Fe/H]. This approach builds upon a previous catalogue of stellar parameters, whereby only the fundamental atmospheric parameters are fit with solar-scaled models. Here, we use the same Markov Chain Monte Carlo method with the additional free parameter [α/Fe]. Using the full spectral fitting code pPXF, we are able to fit multiple lines sensitive to [α/Fe] for a more robust measurement. Quality flags based on the convergence of the sampler, errors in [α/Fe] and a cut in the χ2 of the model fit are used to clean the final catalogue, returning 17 214 spectra and values in the range of −0.25 < [α/Fe] < 0.48. Comparing our calculated [α/Fe] with literature values reveals a degeneracy in cool stars with log g ≥ ∼4; this comparison is then used to create an alternative and calibrated parameter set. We also plot the final catalogue in an [Fe/H] versus [α/Fe] diagram and recover the expected result of increasing [α/Fe] with decreasing [Fe/H] for Milky Way disc-halo stars. We apply our method to a subsample of spectra of uniform resolution and higher signal to noise that finds that our results are independent of this higher signal to noise. In the context of stellar population models, we are able to cover a parameter space for the creation of intermediate to old age models at solar-scaled [α/Fe], high [Fe/H] and enhanced [α/Fe], low [Fe/H].

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Phase space distribution functions and energy distributions of dark matter particles in haloes

ABSTRACT For a spherical dark matter halo with isotropic velocity distribution, the phase space distribution function (DF), the energy distribution, and the density profile form a set of self-consistent description of its equilibrium state, and knowing one is sufficient to determine the other two. The Navarro–Frenk–White density profile (NFW profile) is known to be a good approximation to the spherically averaged density distribution in simulated haloes. The DARKexp energy distribution is also known to compare well with the simulated energy distribution. We present a quantitative assessment of the NFW and DARKexp fits to the simulated DF and energy distribution for a wide range of haloes in a dark-matter-only simulation from the IllustrisTNG Project. As expected, we find that the NFW fits work well except at low energy when the density at small radii deviates from the NFW profile. Further, the NFW and DARKexp fits have comparable accuracy in the region where both fit well, but the DARKexp fits are better at low energy because they require matching of the central gravitational potential. We also find an approximate relation between the energy scale parametrizing the DARKexp energy distribution and that defined by the characteristic density and radius of the NFW profile. This relation may be linked to the relaxation process during halo formation.

Astronomy & Astrophysics↗

DOE OSTI2025

Constraining the phase shift of relativistic species in DESI BAOs

In the early Universe, neutrinos decouple quickly from the primordial plasma and propagate without further interactions. The impact of free-streaming neutrinos is to create a temporal shift in the gravitational potential that impacts the acoustic waves known as baryon acoustic oscillations (BAOs), resulting in a non-linear spatial shift in the Fourier-space BAO signal. In this work, we make use of and extend upon an existing methodology to measure the phase shift amplitude $\beta _{\phi }$ and apply it to the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) BAOs with an anisotropic BAO fitting pipeline. We validate the fitting methodology by testing the pipeline with two publicly available fitting codes applied to highly precise cubic box simulations and realistic simulations representative of the DESI DR1 data. We find further study towards the methods used in fitting the BAO signal will be necessary to ensure accurate constraints on $\beta _{\phi }$ in future DESI data releases. Using DESI DR1, we present individual measurements of the anisotropic BAO distortion parameters and the $\beta _{\phi }$ for the different tracers, and additionally a combined fit to $\beta _{\phi }$ resulting in $\beta _{\phi } = 2.7 \pm 1.7$. After including a prior on the distortion parameters from constraints using Planck we find $\beta _{\phi } = 2.7^{+0.60}_{-0.67}$ suggesting $\beta _{\phi } > 0$ at 4.3$\sigma$ significance. This result may hint at a phase shift that is not purely sourced from the standard model expectation for $N_{\rm {eff}}$ or could be a upwards statistical fluctuation in the measured $\beta _{\phi }$; this result relaxes in models with additional freedom beyond Lambda-cold dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

A Detection of the Baryon Acoustic Oscillation features in the SDSS BOSS DR12 Galaxy Bispectrum

We present the first high-significance detection (4.1σ) of the Baryon Acoustic Oscillations (BAO) feature in the galaxy bispectrum of the twelfth data release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS) constant mass sample (0.43 ≤ z ≤ 0.7). We measured the scale dilation parameter, α, using the power spectrum, bispectrum, and both simultaneously for DR12, plus 2048 MultiDark-PATCHY mocks in the North and South Galactic Caps (NGC and SGC, respectively), and the volume-weighted averages of those two samples (N+SGC). The fitting to the mocks validated our analysis pipeline, yielding values consistent with the mock cosmology. By fitting to the power spectrum and bispectrum separately, we tested the robustness of our results, finding consistent values from the NGC, SGC, and N+SGC in all cases. We found DV = 2032 ± 24(stat.) ± 15(sys.) Mpc, DV = 2038 ± 55(stat.) ± 15(sys.) Mpc, and DV= 2031 ± 22(stat.) ± 10(sys.) Mpc from the N+SGC power spectrum, bispectrum, and simultaneous fitting, respectively. Our bispectrum measurement precision was mainly limited by the size of the covariance matrix. Based on the fits to the mocks, we showed that if a less noisy estimator of the covariance were available, from either a theoretical computation or a larger suite of mocks, the constraints from the bispectrum and simultaneous fits would improve to 1.1 per cent (1.3 per cent with systematics) and 0.7 per cent (0.9 per cent with systematics), respectively, with the latter being slightly more precise than the power spectrum only constraints from the reconstructed field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Toward shell model interactions with credible uncertainties

Background: The nuclear shell model is a powerful framework for predicting nuclear structure observables, but relies on interaction matrix elements fit to experimental data as its inputs. Extending the shell model's applicability, particularly toward dripline nuclei, requires efficient fitting methods and credible uncertainty quantification. Traditional approaches face computational challenges and may underestimate uncertainties. Purpose: We develop and test a framework combining eigenvector continuation and Markov chain Monte Carlo to efficiently fit shell model interaction matrix elements and quantify their uncertainties. Methods: Eigenvector continuation is used to emulate shell model calculations, reducing computational costs. The emulator enables Markov chain Monte Carlo sampling to optimize interaction matrix elements and rigorously assess parametric uncertainties. Here, the framework is benchmarked using the USDB interaction in the 𝑠⁢𝑑 shell. Results: The emulator reproduces the USDB interaction with negligible error, validating its use in shell model fitting applications. However, we find that to obtain credible predictive intervals, the model defect of the shell model itself, rather than experimental or emulator error, must be taken into account in order to obtain credible uncertainties. Conclusions: The proposed framework provides an efficient and rigorous approach for fitting shell model interactions and quantifying uncertainties. Further, the normality assumption used in the past appears sufficient to describe the distribution of interaction matrix elements. However, it is crucial to account for model correlations to avoid underestimating uncertainties.

Nuclear forces↗

DOE OSTI2019

Helical antiferromagnetic ordering in EuNi 1.95 As 2 single crystals

The Eu + 2 spins-7/2 in EuNi 2 As 2 with the body-centered tetragonal ThCr 2 Si 2 structure order antiferromagnetically below the Néel temperature T N = 15 K into a helical antiferromagnetic (AFM) structure with the helix axis aligned along the tetragonal c axis and the Eu ordered moments aligned ferromagnetically within the ab plane as previously reported from neutron diffraction measurements [T. Jin et al., Phys. Rev. B 99, 014425 (2019)]. Here we study the crystallographic, magnetic, thermal, and electronic transport properties of Bi-flux-grown single crystals using single-crystal x-ray diffraction, anisotropic magnetic susceptibility χ, isothermal magnetization M, heat capacity Cp, and electrical resistivity ρ measurements versus applied magnetic field H and temperature T. Vacancies are found on the Ni sites corresponding to the composition EuNi 1.95 (1) As 2. A good fit of the ρ (T) data by the Bloch-Grüneisen theory for metals was obtained. The χab (T) data below T N are fitted well by molecular field theory (MFT), and the helix turn angle kd and the Eu-Eu Heisenberg exchange constants are extracted from the fit parameters. The kd value is in good agreement with the neutron-diffraction result. The magnetic contribution to the zero-field heat capacity below T N is also fitted by MFT. The isothermal in-plane magnetization M ab exhibits two metamagnetic transitions versus H, whereas M c (T = 2 K ) is nearly linear up to H = 14 T, both behaviors being consistent with MFT. The M c (H, T) , ρ(H c , T) , and Cp (H c , T) data yielded a H c - T phase diagram separating the AFM and paramagnetic phases in good agreement with MFT. Anisotropic χ(T) literature data for the ThCr 2 Si 2 -type helical antiferromagnet EuRh 2 As 2 were also fitted well by MFT. A comparison is made between the crystallographic and magnetic properties of ThCr 2 Si 2 -type Eu M 2 Pn 2 compounds with M = Fe , Co, Ni, Cu, or Rh, and Pn = P or As, where only ferromagnetic and c -axis helical AFM structures are found.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Investigation of secondary γ -ray angular distributions using the N 15 ( p , α 1 γ ) C * 12 reaction

The observation of secondary γ-rays provides an alternative method of measuring cross sections that populate excited final states in nuclear reactions. The angular distributions of these γ-rays also provide information on the underlying reaction mechanism. Despite a large amount of data of this type in the literature, publicly available R-matrix codes do not have the ability to calculate these types of angular distributions. In this report, the mathematical formalism derived in Brune and deBoer is implemented in the R-matrix code AZURE2 and calculations are compared with previous data from the literature for the 15 N( p, α1γ ) 12 C* reaction. In addition, new measurements, made at the University of Notre Dame Nuclear Science Laboratory using the Hybrid Array of Gamma Ray Detectors (HAGRiD), are reported that span an energy range from E p = 0.88 to 4.0 MeV. Excellent agreement between data and the phenomenological fit is obtained up to the limit of the previous fit at E p = 2.0 MeV and the R-matrix fit is extended from E x ≈ 13.5 up to E x ≈ 15.3 MeV, where 15 N+p and 12 C+α reactions are fit simultaneously for the first time. An excellent reproduction of the 15 N( p, α1γ ) 12 C* and 12 C(α, α) 12 C data is achieved, but inconsistencies and difficulty in fitting other data is encountered and discussed.

6 ≤ A ≤ 19↗

DOE OSTI2024

Single-channel and single-energy partial-wave analysis with continuity improved through minimal phase constraints

Single-energy partial-wave analysis has often been applied as a way to fit data with minimal model dependence. However, remaining unconstrained, partial waves at neighboring energies will vary discontinuously because the overall amplitude phase cannot be determined through single-channel measurements. This problem can be mitigated through the use of a constraining penalty function based on an associated energy-dependent fit. However, the weight given to this constraint results in a biased fit to the data. In this paper, for the first time, we explore a constraining function which does not influence the fit to data. The constraint comes from the overall phase found in multichannel fits which, in the present study, are the Bonn-Gatchina and Jülich-Bonn multichannel analyses. The data are well reproduced and weighting of the penalty function does not influence the result. The method is applied to K⁢Λ photoproduction data and all observables can be maximally well reproduced. While the employed multichannel analyses display very different multipole amplitudes, we show that the major difference between two sets of multipoles can be related to the different overall phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Effects of renormalon scheme and perturbative scale choices on determinations of the strong coupling from e + e − event shapes

We study the role of renormalon cancellation schemes and perturbative scale choices in extractions of the strong coupling constant α s ( m Z ) and the leading nonperturbative shift parameter Ω 1 from resummed predictions of the e + e − event shape thrust. We calculate the thrust distribution to N L 3 L ′ resummed accuracy in soft-collinear effective theory (SCET) matched to the fixed-order O ( α s 2 ) prediction, and perform a new high-statistics computation of the O ( α s 3 ) matching in , although we do not include the latter in our final α s fits due to some observed systematics that require further investigation. We are primarily interested in testing the phenomenological impact sourced from varying amongst three renormalon cancellation schemes and two sets of perturbative scale profile choices. We then perform a global fit to available data spanning center-of-mass energies between 35–207 GeV in each scenario. Relevant subsets of our results are consistent with prior SCET-based extractions of α s ( m Z ) , but we are also led to a number of novel observations. Notably, we find that the combined effect of altering the renormalon cancellation scheme and profile parameters can lead to few-percent-level impacts on the extracted values in the α s − Ω 1 plane, indicating a potentially important systematic theory uncertainty that should be accounted for. We also observe that fits performed over windows dominated by dijet events are typically of a higher quality than those that extend into the far tails of the distributions, possibly motivating future fits focused more heavily in this region. Finally, we discuss how different estimates of the three-loop soft matching coefficient c S ˜ 3 can also lead to measurable changes in the fitted { α s , Ω 1 } values. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Pati-Salam SUSY GUT with Yukawa unification

Prior studies of a three family Yukawa unified supersymmetric grand unified theory (SUSY GUT) with SO(10) or Pati-Salam (PS) gauge symmetry proposed by Raby and students show that this model is able to fit low energy and inflation observables. Yet, the fit to low energy observables is not great especially for $\sin2\beta$, and up and down quark masses. In this paper, we show that by choosing PS as the gauge group and modifying the Yukawa sector, the low energy fit improves significantly while other qualities of the model are maintained. In particular, the lightest SUSY particle is the neutralino with mass of order 300 - 500 GeV, the lightest stop and sbottom have mass of order 3 - 5 TeV and the CP odd Higgs mass is of order 5 - 6 TeV, so we are in the decoupling limit for the light Higgs. In addition, we reinterpret gluino simplified model analyses by the ATLAS and CMS collaborations and find that the most stringent gluino mass bound for our model is $M_{\tilde{g}}\sim1.9\,\text{TeV}$. The current best fit point, consistent with this bound, with gluino mass $M_{\tilde{g}}=1.9\,\text{TeV}$ has $\chi^2/\text{dof}\approx1.12$, compared to the best fit point of the previous model with $\chi^2/\text{dof}=1.90$. We find that predictions for the electric dipole moment of the electron, the branching ratio $BR(\mu \rightarrow e \gamma)$ and the CP violating angle in the lepton sector, $\sin\delta$, are affected significantly as compared to previous results. In conclusion, we are unable to rule out this model or place an upper bound on gluino mass accessible by this run of the LHC because the $\chi^2/\text{dof}$ of this model is well below $2\sigma$ even for a gluino mass as high as 2.7 TeV. On the bright side, this means that this model is still viable and we might find low energy SUSY particles in the near future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Physically interpretable approximations of many-body spectral functions

The rational function approximation provides a natural and interpretable representation of response functions such as the many-body spectral functions. We apply the vector fitting (VFIT) algorithm to fit a variety of spectral functions calculated from the Holstein model of electron-phonon interactions. We show that the resulting rational functions are highly efficient in their fitting of sharp features in the spectral functions, and could provide a means to infer physically relevant information from a spectral data set. The position of the peaks in the approximated spectral function are determined by the location of poles in the complex plane. Additionally, we developed a variant of VFIT that incorporates regularization to improve the quality of fits. With this procedure, we demonstrate it is possible to achieve accurate spectral function fits that vary smoothly as a function of physical conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

MapsTorch : automatic differentiation for X-ray fluorescence data analysis

X-ray fluorescence (XRF) is a popular spectroscopy technique for elemental analysis. Spectrum fitting and parameter tuning are at the core of XRF analysis and are conventionally manually intensive, especially for synchrotron experiments involving large amounts of diverse samples. This work introduces the automatic differentiation (AD) technique to XRF and an open-source package called MapsTorch. By transforming an analytical model of the XRF spectrum into a differentiable computation graph with AD, MapsTorch enables robust optimization of parameters and elemental intensities. We evaluate MapsTorch by conducting computational experiments on a large number of historical synchrotron XRF datasets and compare its performance with the currently practiced fitting tool NLopt. The results show that MapsTorch consistently achieves high-quality fits and often leads to better fitting quality than NLopt, particularly in tasks such as initial spectrum fitting and elemental intensity refinement. The robust performance of MapsTorch paves the way for developing automated and high-throughput XRF data analysis workflows to handle the increasing data volumes expected from next-generation synchrotron facilities.

X-ray fluorescence↗

DOE OSTI2017

Optimal Pulse Processing, Pile-Up Decomposition, and Applications of Silicon Drift Detectors at LCLS

Silicon drift detectors (SDDs) revolutionized spectroscopy in fields as diverse as geology and dentistry. For a subset of experiments at ultrafast, X-ray free-electron lasers (FELs), SDDs can make substantial contributions. Often the unknown spectrum is interesting, carrying science data, or the background measurement is useful to identify unexpected signals. Many measurements involve only several discrete photon energies known a priori, allowing single-event decomposition of pile-up and spectroscopic photon counting. We designed a pulse function and demonstrated that the signal amplitude (i.e., proportional to the detected energy and obtained from fitting with the pulse function), rise time, and pulse height are interrelated, and at short peaking times, the pulse height and pulse area are not optimal estimators for detected energy; instead, the signal amplitude and rise time are obtained for each pulse by fitting, thus removing the need for pulse shaping. By avoiding pulse shaping, rise times of tens of nanoseconds resulted in reduced pulse pile-up and allowed decomposition of remaining pulse pile-up at photon separation times down to hundreds of nanoseconds while yielding time-of-arrival information with the precision of 10 ns. Waveform fitting yields simultaneously high energy resolution and high counting rates (two orders of magnitude higher than current digital pulse processors). At pulsed sources or high photon rates, photon pile-up still occurs. We showed that pile-up spectrum fitting is relatively simple and preferable to pile-up spectrum deconvolution. We then developed a photon pile-up statistical model for constant intensity sources, extended it to variable intensity sources (typical for FELs), and used it to fit a complex pileup spectrum. We subsequently developed a Bayesian pile-up decomposition method that allows decomposing pile-up of single events with up to six photons from six monochromatic lines with 99% accuracy. The usefulness of SDDs will continue into the X-ray FEL era of science. Their successors, the ePixS hybrid pixel detectors, already offer hundreds of pixels, each with a similar performance to an SDD, in a compact, robust and affordable package.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2017

Physical Activity and School Absenteeism Due to Illness in Adolescents

ABSTRACT BACKGROUND Knowledge about the beneficial role of physical activity (PA) for health and school performance is growing. Studies investigating the link between PA and school absenteeism due to illness are lacking. Therefore, we investigated associations between habitual PA and school absenteeism due to illness in adolescents and explored whether mental health and cardiovascular fitness mediated this association. METHODS We studied 328 students in grades 7 and 9 (mean age 13.8 years; 49% boys). The PA was measured objectively by an ActivPAL3™ accelerometer attached on the thigh during 1 full week (24 hours/day). Depressive symptoms and self‐esteem were self‐reported by the Center for Epidemiologic Studies Depression Scale (CES‐D) and Rosenberg Self‐Esteem Scale, respectively, and included as a proxy for mental health in the analyses. Cardiovascular fitness was measured by the 20‐m shuttle‐run test. School absenteeism due to illness data was provided by the school administration. RESULTS The PA was not significantly associated with school absenteeism, though there was an indirect association between PA and school absenteeism by cardiovascular fitness. CONCLUSIONS Cardiovascular fitness mediates the association between PA and school absenteeism due to illness. Thus, cardiovascular fitness of students should be improved to reduce school absenteeism due to illness.

de Groot, Renate↗

DOE OSTI2014

Radiochromic film calibration for low-energy seed brachytherapy dose measurement

Purpose: Radiochromic film dosimetry is typically performed for high energy photons and moderate doses characterizing external beam radiotherapy (XRT). The purpose of this study was to investigate the accuracy of previously established film calibration procedures used in XRT when applied to low-energy, seed-based brachytherapy at higher doses, and to determine necessary modifications to achieve similar accuracy in absolute dose measurements. Methods: Gafchromic EBT3 film was used to measure radiation doses upwards of 35 Gy from 75 kVp, 200 kVp, 6 MV, and (∼28 keV) I-125 photon sources. For the latter irradiations a custom phantom was built to hold a single I-125 seed. Film pieces were scanned with an Epson 10000XL flatbed scanner and the resulting 48-bit RGB TIFF images were analyzed using both FilmQA Pro software andMATLAB. Calibration curves relating dose and optical density via a rational functional form for all three color channels at each irradiation energy were determined with and without the inclusion of uncertainties in the measured optical densities and dose values. The accuracy of calibration curve variations obtained using piecewise fitting, a reduced film measurement area for I-125 irradiation, and a reduced number of dose levels was also investigated. The energy dependence of the film lot used was also analyzed by calculating normalized optical density values. Results: Slight differences were found in the resulting calibration curves for the various fitting methods used. The accuracy of the calibration curves was found to improve at low doses and worsen at high doses when including uncertainties in optical densities and doses, which may better represent the variability that could be seen in film optical density measurements. When exposing the films to doses > 8 Gy, two-segment piecewise fitting was found to be necessary to achieve similar accuracies in absolute dose measurements as when using smaller dose ranges. When reducing the film measurement area for the I-125 irradiations, the accuracy of the calibration curve was degraded due to the presence of localized film heterogeneities. No degradation in the calibration curves was found when reducing the number of calibration points down to only 4, but with piecewise fitting, 6 calibration points as well as a blank film are required. Variations due to photon energy in film optical density of up to 3% were found above doses of 2 Gy. Conclusions: A modified procedure for performing EBT3 film calibration was established for use with low-energy brachytherapy seeds and high dose exposures. The energy dependence between 6 MV and I-125 photons is significant such that film calibrations should be done with an appropriately low-energy source when performing low-energy brachytherapy dose measurements. Two-segment piecewise fitting with the inclusion of errors in measured optical density as well as dose was found to result in the most accurate calibration curves. Above doses of 1 Gy, absolute dose measurements can be made with an accuracy of 1.6% for 6 MV beams and 5.7% for I-125 seed exposures if using the I-125 source for calibration, or 2.3% if using the 75 kVp photon beam for calibration.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2014

SU-C-9A-04: Alternative Analytic Solution to the Paralyzable Detector Model to Calculate Deadtime and Deadtime Loss

Purpose: Some common methods to solve for deadtime are (1) dual-source method, which assumes two equal activities; (2) model fitting, which requires multiple acquisitions as source decays; and (3) lossless model, which assumes no deadtime loss at low count rates. We propose a new analytic alternative solution to calculate deadtime for paralyzable gamma camera. Methods: Deadtime T can be calculated analytically from two distinct observed count rates M1 and M2 when the ratio of the true count rates alpha=N2/N1 is known. Alpha can be measured as a ratio of two measured activities using dose calibrators or via radioactive decay. Knowledge of alpha creates a system with 2 equations and 2 unknowns, i.e., T and N1. To verify the validity of the proposed method, projections of a non-uniform phantom (4GBq 99mTc) were acquired in using Siemens SymbiaS multiple times over 48 hours. Each projection has >100kcts. The deadtime for each projection was calculated by fitting the data to a paralyzable model and also by using the proposed 2-acquisition method. The two estimates of deadtime were compared using the Bland-Altmann method. In addition, the dependency of uncertainty in T on uncertainty in alpha was investigated for several imaging conditions. Results: The results strongly suggest that the 2-acquisition method is equivalent to the fitting method. The Bland-Altman analysis yielded mean difference in deadtime estimate of ∼0.076us (95%CI: -0.049us, 0.103us) between the 2-acquisition and model fitting methods. The 95% limits of agreement were calculated to be -0.104 to 0.256us. The uncertainty in deadtime calculated using the proposed method is highly dependent on the uncertainty in the ratio alpha. Conclusion: The 2-acquisition method was found to be equivalent to the parameter fitting method. The proposed method offers a simpler and more practical way to analytically solve for a paralyzable detector deadtime, especially during physics testing.

07 ISOTOPE AND RADIATION SOURCES↗