An official website of the United States government
Energy-resolving photon-counting detectors (PCDs) separate photons from a polychromatic X-ray source into a number of separate energy bins. This spectral information from PCDs would allow advancements in X-ray imaging, such as improving image contrast, quantitative imaging, and material identification and characterization. However, aspects like detector spectral distortions and scattered photons from the object can impede these advantages if left unaccounted for. Scattered X-ray photons act as noise in an image and reduce image contrast, thereby significantly hindering PCD utility. In this paper, we explore and outline several important characteristics of spectral X-ray scatter with examples of soft-material imaging (such as cancer imaging in mammography or explosives detection in airport security). Our results showed critical spectral signatures of scattered photons that depend on a few adjustable experimental factors. Additionally, energy bins over a large portion of the spectrum exhibit lower scatter-to-primary ratio in comparison to what would be expected when using a conventional energy-integrating detector. These important findings allow flexible choice of scatter-correction methods and energy-bin utilization when using PCDs. Our findings also propel the development of efficient spectral X-ray scatter correction methods for a wide range of PCD-based applications.
more »
« less
An accurate treatment of scattering and diffusion in piecewise power-law models for cosmic ray and radiation/neutrino transport
ABSTRACT A popular numerical method to model the dynamics of a ‘full spectrum’ of cosmic rays (CRs), also applicable to radiation/neutrino hydrodynamics, is to discretize the spectrum at each location/cell as a piecewise power law in ‘bins’ of momentum (or frequency) space. This gives rise to a pair of conserved quantities (e.g. CR number and energy) that are exchanged between cells or bins, which in turn give the update to the normalization and slope of the spectrum in each bin. While these methods can be evolved exactly in momentum-space (e.g. considering injection, absorption, continuous losses/gains), numerical challenges arise dealing with spatial fluxes, if the scattering rates depend on momentum. This has often been treated either by neglecting variation of those rates ‘within the bin,’ or sacrificing conservation – introducing significant errors. Here, we derive a rigorous treatment of these terms, and show that the variation within the bin can be accounted for accurately with a simple set of scalar correction coefficients that can be written entirely in terms of other, explicitly evolved ‘bin-integrated’ quantities. This eliminates the relevant errors without added computational cost, has no effect on the numerical stability of the method, and retains manifest conservation. We derive correction terms both for methods that explicitly integrate flux variables (e.g. two-moment or M1-like) methods, as well as single-moment (advection-diffusion, FLD-like) methods, and approximate corrections valid in various limits.
more »
« less
- PAR ID:
- 10385408
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 518
- Issue:
- 4
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 5882-5892
- Size(s):
- p. 5882-5892
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Unfolding is an ill-posed inverse problem in particle physics aiming to infer a true particle-level spectrum from smeared detector-level data. For computational and practical reasons, these spaces are typically discretized using histograms, and the smearing is modeled through a response matrix corresponding to a discretized smearing kernel of the particle detector. This response matrix depends on the unknown shape of the true spectrum, leading to a fundamental systematic uncertainty in the unfolding problem. To handle the ill-posed nature of the problem, common approaches regularize the problem either directly via methods such as Tikhonov regularization, or implicitly by using wide-bins in the true space that match the resolution of the detector. Unfortunately, both of these methods lead to a non-trivial bias in the unfolded estimator, thereby hampering frequentist coverage guarantees for confidence intervals constructed from these methods. We propose two new approaches to addressing the bias in the wide-bin setting through methods called One-at-a-time Strict Bounds (OSB) and Prior-Optimized (PO) intervals. The OSB intervals are a bin-wise modification of an existing guaranteed-coverage procedure, while the PO intervals are based on a decision-theoretic view of the problem. Importantly, both approaches provide well-calibrated frequentist confidence intervals even in constrained and rank-deficient settings. These methods are built upon a more general answer to the wide-bin bias problem, involving unfolding with fine bins first, followed by constructing confidence intervals for linear functionals of the fine-bin counts. We test and compare these methods to other available methodologies in a wide-bin deconvolution example and a realistic particle physics simulation of unfolding a steeply falling particle spectrum.more » « less
-
Background: When phenotypic characters are described in the literature, they may be constrained or clarified with additional information such as the location or degree of expression, these terms are called “modifiers”. With effort underway to convert narrative character descriptions to computable data, ontologies for such modifiers are needed. Such ontologies can also be used to guide term usage in future publications. Spatial and method modifiers are the subjects of ontologies that already have been developed or are under development. In this work, frequency (e.g., rarely, usually), certainty (e.g., probably, definitely), degree (e.g., slightly, extremely), and coverage modifiers (e.g., sparsely, entirely) are collected, reviewed, and used to create two modifier ontologies with different design considerations. The basic goal is to express the sequential relationships within a type of modifiers, for example, usually is more frequent than rarely, in order to allow data annotated with ontology terms to be classified accordingly. Method: Two designs are proposed for the ontology, both using the list pattern: a closed ordered list (i.e., five-bin design) and an open ordered list design. The five-bin design puts the modifier terms into a set of 5 fixed bins with interval object properties, for example, one_level_more/less_frequently_than, where new terms can only be added as synonyms to existing classes. The open list approach starts with 5 bins, but supports the extensibility of the list via ordinal properties, for example, more/less_frequently_than, allowing new terms to be inserted as a new class anywhere in the list. The consequences of the different design decisions are discussed in the paper. CharaParser was used to extract modifiers from plant, ant, and other taxonomic descriptions. After a manual screening, 130 modifier words were selected as the candidate terms for the modifier ontologies. Four curators/experts (three biologists and one information scientist specialized in biosemantics) reviewed and categorized the terms into 20 bins using the Ontology Term Organizer (OTO) (http://biosemantics.arizona.edu/OTO). Inter-curator variations were reviewed and expressed in the final ontologies. Results: Frequency, certainty, degree, and coverage terms with complete agreement among all curators were used as class labels or exact synonyms. Terms with different interpretations were either excluded or included using “broader synonym” or “not recommended” annotation properties. These annotations explicitly allow for the user to be aware of the semantic ambiguity associated with the terms and whether they should be used with caution or avoided. Expert categorization results showed that 16 out of 20 bins contained terms with full agreements, suggesting differentiating the modifiers into 5 levels/bins balances the need to differentiate modifiers and the need for the ontology to reflect user consensus. Two ontologies, developed using the Protege ontology editor, are made available as OWL files and can be downloaded from https://github.com/biosemantics/ontologies. Contribution: We built the first two modifier ontologies following a consensus-based approach with terms commonly used in taxonomic literature. The five-bin ontology has been used in the Explorer of Taxon Concepts web toolkit to compute the similarity between characters extracted from literature to facilitate taxon concepts alignments. The two ontologies will also be used in an ontology-informed authoring tool for taxonomists to facilitate consistency in modifier term usage.more » « less
-
ABSTRACT We derive a consistent set of moment equations for cosmic ray (CR)-magnetohydrodynamics, assuming a gyrotropic distribution function (DF). Unlike previous efforts, we derive a closure, akin to the M1 closure in radiation hydrodynamics (RHD), that is valid in both the nearly isotropic DF and/or strong-scattering regimes, and the arbitrarily anisotropic DF or free-streaming regimes, as well as allowing for anisotropic scattering and transport/magnetic field structure. We present the appropriate two-moment closure and equations for various choices of evolved variables, including the CR phase space DF f, number density n, total energy e, kinetic energy ϵ, and their fluxes or higher moments, and the appropriate coupling terms to the gas. We show that this naturally includes and generalizes a variety of terms including convection/fluid motion, anisotropic CR pressure, streaming, diffusion, gyro-resonant/streaming losses, and re-acceleration. We discuss how this extends previous treatments of CR transport including diffusion and moment methods and popular forms of the Fokker–Planck equation, as well as how this differs from the analogous M1-RHD equations. We also present two different methods for incorporating a reduced speed of light (RSOL) to reduce time-step limitations: In both, we carefully address where the RSOL (versus true c) must appear for the correct behaviour to be recovered in all interesting limits, and show how current implementations of CRs with an RSOL neglect some additional terms.more » « less
-
In this work, a field representation of the conservation law of linear momentum is derived from the atomistic, using the theory of distributions as the mathematical tool, and expressed in terms of temperature field by defining temperature as a derived quantity as that in molecular kinetic theory and atomistic simulations. The formulation leads to a unified atomistic and continuum description of temperature and a new linear momentum equation that, supplemented by an interatomic potential, completely governs thermal and mechanical processes across scales from the atomic to the continuum. The conservation equation can be used to solve atomistic trajectories for systems at finite temperatures, as well as the evolution of field quantities in space and time, with atomic or multiscale resolution. Four sets of numerical examples are presented to demonstrate the efficacy of the formulation in capturing the effect of temperature and thermal fluctuations, including phonon density of states, thermally activated dislocation motion, dislocation formation during epitaxial processes, and attenuation of longitudinal acoustic waves as a result of their interaction with thermal phonons.more » « less
