QCD corrections in Standard Model EFT fits to W Z and W W production
Not Available
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
SEARCH · Engineering Papers
Results for “Model fitting”
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.
Observations of the Ultraviolet-bright Star Y453 in the Globular Cluster M4 (NGC 6121)
We present a spectral analysis of the UV-bright star Y453 in M4. Model fits to the star’s optical spectrum yield T{sub eff} ∼ 56,000 K. Fits to the star’s FUV spectrum, obtained with the Cosmic Origins Spectrograph on board the Hubble Space Telescope, reveal it to be considerably hotter, with T{sub eff} ∼ 72,000 K. We adopt T{sub eff} = 72,000 ± 2000 K and logg = 5.7 ± 0.2 as our best-fit parameters. Scaling the model spectrum to match the star’s optical and near-infrared magnitudes, we derive a mass M{sub ∗}=0.53±0.24 M{sub ⊙} and luminosity logL/L{sub ⊙}=2.84±0.05, consistent with the values expected of an evolved star in a globular cluster. Comparing the star with post-horizontal-branch evolutionary tracks, we conclude that it most likely evolved from the blue horizontal branch, departing the asymptotic giant branch before third dredge-up. It should thus exhibit the abundance pattern (O-poor and Na-rich) characteristic of the second-generation (SG) stars in M4. We derive the star’s photospheric abundances of He, C, N, O, Si, S, Ti, Cr, Fe, and Ni. CNO abundances are roughly 0.25 dex greater than those of the cluster’s SG stars, while the Si and S abundances match the cluster values. Abundances of the iron-peak elements (except for iron itself) are enhanced by 1–3 dex. Rather than revealing the star’s origin and evolution, this pattern reflects the combined effects of diffusive and mechanical processes in the stellar atmosphere.
79 ASTRONOMY AND ASTROPHYSICS↗
Line-shape code comparison through modeling and fitting of experimental spectra of the C II 723-nm line emitted by the ablation cloud of a carbon pellet.
Abstract not provided.
Hansen, Stephanie B.↗
Modeling CH 4 and CO 2 cycling using porewater stable isotopes in a thermokarst bog in Interior Alaska: results from three conceptual reaction networks
Quantifying rates of microbial carbon transformation in peatlands is essential for gaining mechanistic understanding of the factors that influence methane emissions from these systems, and for predicting how emissions will respond to climate change and other disturbances. In this study, we used porewater stable isotopes collected from both the edge and center of a thermokarst bog in Interior Alaska to estimate in situ microbial reaction rates. We expected that near the edge of the thaw feature, actively thawing permafrost and greater abundance of sedges would increase carbon, oxygen and nutrient availability, enabling faster microbial rates relative to the center of the thaw feature. We developed three different conceptual reaction networks that explained the temporal change in porewater CO2, CH4, δ13C-CO2 and δ13C-CH4. All three reaction-network models included methane production, methane oxidation and CO2 production, and two of the models included homoacetogenesis — a reaction not previously included in isotope-based porewater models. All three models fit the data equally well, but rates resulting from the models differed. Most notably, inclusion of homoacetogenesis altered the modeled pathways of methane production when the reaction was directly coupled to methanogenesis, and it decreased gross methane production rates by up to a factor of five when it remained decoupled from methanogenesis. The ability of all three conceptual reaction networks to successfully match the measured data indicate that this technique for estimating in-situ reaction rates requires other data and information from the site to confirm the considered set of microbial reactions. Despite these differences, all models indicated that, as expected, rates were greater at the edge than in the center of the thaw bog, that rates at the edge increased more during the growing season than did rates in the center, and that the ratio of acetoclastic to hydrogenotrophic methanogenesis was greater at the edge than in the center. In both locations, modeled rates (excluding methane oxidation) increased with depth. A puzzling outcome from the effort was that none of the models could fit the porewater dataset without generating “fugitive” carbon (i.e., methane or acetate generated by the models but not detected at the field site), indicating that either our conceptualization of the reactions occurring at the site remains incomplete or our site measurements are missing important carbon transformations and/or carbon fluxes. This model–data discrepancy will motivate and inform future research efforts focused on improving our understanding of carbon cycling in permafrost wetlands.
13CH4↗
Photoionization Emission Models for the Cyg X-3 X-Ray Spectrum
We present model fits to the X-ray line spectrum of the well-known high-mass X-ray binary Cyg X-3. The primary observational data set is a spectrum taken with the Chandra X-ray Observatory High Energy Transmission Grating in 2006, though we compare it to all the other observations of this source taken so far by this instrument. We show that the density must be ≥10 12 cm –3 in the region responsible for most of the emission. We discuss the influence of the dust scattering halo on the broadband spectrum, and we argue that dust scattering and extinction is not the most likely origin for the narrow feature seen near the Si K edge. We identify the features of a wind in the profiles of the strong resonance lines and show that the wind is more apparent in the lines from the lighter elements. We argue that this wind is most likely associated with the companion star. We show that the intensities of most lines can be fitted, crudely, by a single-component photoionized model. Yet, the iron K lines do not fit with this model. We show that the iron K line variability as a function of orbital phase is different from the lower-energy lines, which indicates that the lines arise in physically distinct regions. We discuss the interpretation of these results in the context of what is known about the system and similar systems.
79 ASTRONOMY AND ASTROPHYSICS↗
THE NANOGRAV NINE-YEAR DATA SET: OBSERVATIONS, ARRIVAL TIME MEASUREMENTS, AND ANALYSIS OF 37 MILLISECOND PULSARS
We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or “red,” timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.
79 ASTRONOMY AND ASTROPHYSICS↗
Response surface data for sensitivity study of industrial spray injected fluidized bed reactor
An industrial fluidized bed reactor was designed to convert an aqueous solid laden stream into a consistent granular product. CFD simulations were run using the MFiX two-fluid model for this fluidizing bed operating at 650 °C. A set of simulations were run over a Latin-hypercube sample of five model parameters – bed particle size, bed particle density, coal particle size, spray feed flow rate, and fluidizing gas flow rate. Data from the simulations were collected on three quantities of interest – bed differential temperature, low solids velocity, and bed voidage. The data presented here is the full span of response surfaces generated using the process Gaussian response surface model in the Dakota code, as well as the table of data for coefficients of the fitted model. The fits to the five-dimensional Gaussian Process models were 0.7797, 0.8664, and 0.9440 for the temperature, velocity, and solids packing, respectively.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
ASP (Adaptive Splines for Prediction) [SWR-22-66]
This software package provides functionality to fit piecewise cubic splines to model the relationship between a univariate X (independent variable) and a univariate Y (dependent variable). The output is akin to linear regression but has higher "capacity" in that it can learn more complicated relationships despite leveraging only a single predictor. It thus can learn very reasonable relationships between X and Y with minimal effort. The canonical intended use case is for predicting electric demand (Y) from temperature (X). There can be complicated, non-linear relationships between these two variables, but these can be reasonably approximated by the sum of: 1) an "aggregate" relationship between daily average temperature and daily average demand, and 2) intra-day patterns represented as hourly offsets from the daily average load. The spline functionality is used for learning a robust "aggregate" relationship without the inconvenience of manual feature engineering while still ensuring robust results. The model fitting workhorse is R's built-in smooth.spline. This package provides a suite of options for controlling how smooth.spline fits a spline to the data. In particular, it is designed to return a model with a single critical point (i.e., a single location along the domain/support where the first derivative is zero). This is to ensure that predicted changes in electric demand are always positive as temperatures become more extreme. This helps to prevent overfitting, particularly when the input data has few observations or has significant uncertainty from other sources.
Murphy, Sinnott↗
Near-IR polarized scattered light imagery of the DoAr 28 transitional disk
We present the first spatially resolved polarized scattered light H-band detection of the DoAr 28 transitional disk. Our two epochs of imagery detect the scattered light disk from our effective inner working angle of 0.″10 (13 AU) out to 0.″50 (65 AU). This inner working angle is interior to the location of the system’s gap inferred by previous studies using spectral energy distribution modeling (15 AU). We detected a candidate point source companion 1.″08 northwest of the system; however, our second epoch of imagery strongly suggests that this object is a background star. We constructed a grid of Monte Carlo Radiative Transfer models of the system, and our best fit models utilize a modestly inclined (50°), 0.01 M{sub ⊙} disk that has a partially depleted inner gap from the dust sublimation radius out to ∼8 AU. Subtracting this best fit, axi-symmetric model from our polarized intensity data reveals evidence for two small asymmetries in the disk, which could be attributable to a variety of mechanisms.
79 ASTRONOMY AND ASTROPHYSICS↗
Elucidation of Fe-N-C electrocatalyst active site functionality via in-situ X-ray absorption and operando determination of oxygen reduction reaction kinetics in a PEFC
In the past decade the notable effort placed on improving intrinsic electrochemical kinetics of platinum group metal (PGM)-free electrocatalysts for the oxygen reduction reaction (ORR) has led to a significant improvement in both performance and understanding of this class of electrocatalysts. However, a limited amount of this development and understanding has been undertaken using operando electrochemical diagnostics at the membrane electrode assembly (MEA) level. In this work, the operando ORR kinetics on an atomically dispersed iron-nitrogen-carbon ((AD)Fe-N-C) PGM-free electrocatalyst have been examined to extract the reaction order and the activation energy of the ORR. The experiments were carefully designed to ensure the stability/predictability of the electrocatalyst during the data collection process and thus validate the relevance of the values obtained for the aforementioned parameters. A kinetic model that considers a potential-dependent availability of active sites (θ) is proposed. Active site availability is shown to be a function of both the change in the oxidation state (nR) and the redox potential at which the metal center transitions from a higher oxidation state to a lower one (). The resulting model fitting parameters for nR and (0.71 and 0.788 V, respectively) obtained from the analysis of operando data correlate well with those from in situ X-ray absorption near edge structure measurements (nR=0.57) and in situ cyclic voltammetry measurements (0.75 V < <0.8 V) in the MEA environment. The resulting model provides an excellent fit of MEA performance across the range of pressures, temperatures, and potentials under which the data were collected.
Activation energy↗
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↗
Optimization and supervised machine learning methods for fitting numerical physics models without derivatives
Here, we address the calibration of a computationally expensive nuclear physics model for which derivative information with respect to the fit parameters is not readily available. Of particular interest is the performance of optimization-based training algorithms when dozens, rather than millions or more, of training data are available and when the expense of the model places limitations on the number of concurrent model evaluations that can be performed. As a case study, we consider the Fayans energy density functional model, which has characteristics similar to many model fitting and calibration problems in nuclear physics. We analyze hyperparameter tuning considerations and variability associated with stochastic optimization algorithms and illustrate considerations for tuning in different computational settings.
97 MATHEMATICS AND COMPUTING↗
Survival of Immature Gopher Tortoises Recruited into a Translocated Population
Population manipulations such as translocation and head-starting are increasingly used as recovery tools for chelonians. However, evaluating success of individual projects can require decades of monitoring to detect population trends in these long-lived species. Furthermore, there are often few benchmarks from stable, unmanipulated populations against which to compare demographic rates, particularly for the immature stages. Here, we used 8 years of mark-recapture data to estimate apparent survival of immature gopher tortoises recruited into an introduced population of gopher tortoises (Gopherus polyphemus) first established on St. Catherines Island, Georgia, USA in 1987. During 2006 29 -2013, we conducted targeted trapping of immature gopher tortoises and compared survival of the hatchling, juvenile and subadult stages among treatments: 1) individuals released shortly after hatching from eggs obtained from gravid female founders (‘direct releases’); 2) individuals reared in captivity for 6-9 months following hatching (‘head-starts’); and 3) individuals first encountered as free-ranging, wild-recruited offspring (‘wild recruits’). Among the candidate models we examined, the best fit model included additive effects of tortoise stage and treatment, however, overlapping 95% credible intervals among treatments (CI) suggested that survival did not vary significantly among treatments. Annual apparent survival increased over the immature period, highlighting the importance of calculating separate estimates for the different immature stages. Across all treatments, the additive model estimated annual apparent survival probability to be 0.37 (CI: 0.25 40 – 0.48) for hatchlings, 0.71 (CI: 0.61 – 0.81) for juveniles, and 0.83 (CI: 0.74 – 0.94) for subadults. Our study, in combination with previous monitoring efforts at St. Catherines Island, provides strong evidence that the translocation and subsequent population augmentation efforts have been successful in establishing a robust population of gopher tortoises. Additionally, our results provide estimates of demographic rates for life stages that are poorly understood but critical to understanding population dynamics of this imperiled species.
59 BASIC BIOLOGICAL SCIENCES↗
Differentiating Structural Changes of Glycoproteins in Solution using Small Angle Scattering Analysis
Alpha 1-Antitrypsin (A1AT) and Fibrinogen are plasma glycoproteins with different, but specific biological functions. A1AT (which forms part of the serpin family) has been shown to have protective roles of lung cells against emphysema, a disease characterized by lung tissue destruction, while fibrinogen is a major factor in the blood clotting process. While biologically distinct, the two proteins present similarities in biophysical studies due to their size and polymerization tendency. Previous approximations, using radius of gyration analysis (shape independent), have shown that there is a phase transition for A1AT under osmotic stress. The goals of this study are to better understand: the stability of A1AT's native conformation when its environment deviates from physiological conditions; and the effects on fibrin's final structure as fibrinogen undergoes polymerization under different conditions. Environmental variables include osmotic stress, crowding, and presence of other molecules (e.g. lipids, catalytic enzyme thrombin, or cholesterol) found in plasma. The method used to assess the stability of the native conformation was to observe changes in the structural envelope of the protein. Static and dynamic small angle x-ray scattering (SAXS) and small angle neutron scattering (SANS) experiments of the proteins in solution were conducted for this purpose. Scattering data was fitted using various shape dependent models in SASView. Pair distance functions (PDF) were used to validate any changes identified by the model fits. Multi-variable phase diagrams were constructed from the fitting results in order to quantify structural and conformational changes of the proteins.
Segally, Taylor N.↗
Inverse problems-based maximum likelihood estimation of ground reflectivity for selected regions of interest from stripmap SAR data [Regularized maximum likelihood estimation of ground reflectivity from stripmap SAR data]
In this study, we derive a comprehensive forward model for the data collected by stripmap synthetic aperture radar (SAR) that is linear in the ground reflectivity parameters. It is also shown that if the noise model is additive, then the forward model fits into the linear statistical model framework, and the ground reflectivity parameters can be estimated by statistical methods. We derive the maximum likelihood (ML) estimates for the ground reflectivity parameters in the case of additive white Gaussian noise. Furthermore, we show that obtaining the ML estimates of the ground reflectivity requires two steps. The first step amounts to a cross-correlation of the data with a model of the data acquisition parameters, and it is shown that this step has essentially the same processing as the so-called convolution back-projection algorithm. The second step is a complete system inversion that is capable of mitigating the sidelobes of the spatially variant impulse responses remaining after the correlation processing. We also state the Cramer-Rao lower bound (CRLB) for the ML ground reflectivity estimates.We show that the CRLB is linked to the SAR system parameters, the flight path of the SAR sensor, and the image reconstruction grid.We demonstrate the ML image formation and the CRLB bound for synthetically generated data.
47 OTHER INSTRUMENTATION↗
A Bayesian nonparametric analysis for zero-inflated multivariate count data with application to microbiome study
High-throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. Here, we develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community-level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Search for Sterile-Neutrino-Based Muon Neutrino Disappearance Using the MicroBooNe Deep Learning Analysis
We describe a search for νµ disappearance using the MicroBooNE Deep Learning analysis 1µ1p selection. Presently, the unexplained MiniBooNE and LSND anomalies could be explained by a sterile neutrino impacting neutrino oscillations. Our analysis searches for the allowed parameter space that could describe such a sterile neutrino. We determine the allowed and excluded region of a 3+1 sterile-based muon neutrino disappearance model in MicroBooNE at 90% confidence. Our allowed region includes both the null model, and current global best fit model. Context for the underlying Deep Learning analysis is provided and several validation studies surrounding both the disappearance search, and originating 1µ1p selection are performed to strengthen confidence in the result. In addition, a next-generation deep learning tool for cosmic-ray-muon event discrimination is proposed and evaluated, demonstrating a removal of 70% of the remaining event background when added to current methods, under the cut criteria used.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗