Low Photon Statistics Research Articles

Dose-efficient automatic differentiation for ptychographic reconstruction

Ptychography, as a powerful lensless imaging method, has become a popular member of the coherent diffractive imaging family over decades of development. The ability to utilize low-dose X-rays and/or fast scans offers a big advantage in a ptychographic measurement (for example, when measuring radiation-sensitive samples), but results in low-photon statistics, making the subsequent phase retrieval challenging. Here, we demonstrate a dose-efficient automatic differentiation framework for ptychographic reconstruction (DAP) at low-photon statistics and low overlap ratio. As no reciprocal space constraint is required in this DAP framework, the framework, based on various forward models, shows superior performance under these conditions. It effectively suppresses potential artifacts in the reconstructed images, especially for the inherent periodic artifact in a raster scan. We validate the effectiveness and robustness of this method using both simulated and measured datasets.

Very High-energy (>50 GeV) Gamma-Ray Flux Variability of Bright Fermi Blazars

Understanding the high-energy emission processes and variability patterns are two of the most challenging research problems associated with relativistic jets. In particular, the long-term (months to years) flux variability at very high energies (VHE >50 GeV) has remained an unexplored domain so far. This is possibly due to the decreased sensitivity of the Fermi Large Area Telescope (LAT) above a few GeV, hence low photon statistics, and observing constraints associated with the ground-based Cherenkov telescopes. This paper reports the results obtained from the 0.05−2 TeV Fermi-LAT data analysis of a sample of 29 blazars with the primary objective to explore their months-to-year-long very high-energy (VHE) flux variability behavior. This systematic search has led to, for the first time, the detection of significant flux variations in five blazars at the >99% confidence level, whereas eight of them exhibit variability, albeit at a lower confidence level (∼95%–99%). A comparison of the 0.05–2 TeV flux variations with that observed at 0.1–50 GeV band has revealed similar variability behavior for most of the sources. However, complex variability patterns that are not reflected contemporaneously in both energy bands were also detected, thereby providing tantalizing clues about the underlying radiative mechanisms. These results open up a new dimension to unravel the VHE emission processes operating in relativistic jets, hence sowing the seeds for their future observations with the upcoming Cherenkov Telescope Array.

Pushing the limit of BGO-based dual-ended Cherenkov PET detectors through photon transit time correction

Objective. The high production cost of commonly used lutetium-based fast scintillators and the development of silicon photomultipliers technology have made bismuth germanate (BGO) a promising candidate for time-of-flight positron emission tomography (TOF PET) detectors owing to its generation of prompt Cherenkov photons. However, using BGO as a hybrid scintillator is disadvantageous owing to its low photon statistics and distribution that does not conform well to a single Gaussian. To mitigate this, a proposal was made to increase the likelihood of detecting the first Cherenkov photons by positioning two photosensors in opposition at the entrance and exit faces of the scintillator and subsequently selectively picking an earlier timestamp. Nonetheless, the timing variation arising from the photon transit time remains affected by the entire length of the crystal, thereby presenting a possibility for further enhancement. Approach. In this study, we aimed to improve the timing performance of the dual-ended BGO Cherenkov TOF PET detector by capitalizing on the synergistic advantages of applying depth-of-interaction (DOI) information and crystal surface finishes or reflector properties. A dual-ended BGO detector was implemented using a 3 × 3 × 15 mm3 BGO crystal. Coincidence events were acquired against a 3 × 3 × 3 mm3 LYSO:Ce:Mg reference detector. The timing performance of the dual-ended BGO detectors was analyzed using conventionally proposed timestamp methods before and after DOI correction. Results. Through a DOI-based correction of photon transit time spread, we demonstrated a further improvement in the timing resolution of the BGO-based Cherenkov TOF PET detector utilizing a dual-ended detector configuration and adaptive arrival time pickoff. We achieved further improvements in timing resolution by correcting the offset spread induced by the fluctuation of timing signal rise time in the dual-ended detector. Significance. Although polishing the crystal surface was still favorable in terms of full-width-half-maximum value, incorporating DOI information from the unpolished crystal to compensate for photon travel time facilitated additional enhancement in the overall timing performance, thereby surpassing that achieved with the polished crystal.

Quantitative and Qualitative Improvement of Low-Count [68Ga]Citrate and [90Y]Microspheres PET Image Reconstructions Using Block Sequential Regularized Expectation Maximization Algorithm.

There are several important positron emission tomography (PET) imaging scenarios that require imaging with very low photon statistics, for which both quantitative accuracy and visual quality should not be neglected. For example, PET imaging with the low photon statistics is closely related to active efforts to significantly reduce radiation exposure from radiopharmaceuticals. We investigated two examples of low-count PET imaging: (a) imaging [90Y]microsphere radioembolization that suffers the very small positron emission fraction of Y-90's decay processes, and (b) cancer imaging with [68Ga]citrate with uptake time of 3-4 half-lives, necessary for visualizing tumors. In particular, we investigated a type of penalized likelihood reconstruction algorithm, block sequential regularized expectation maximization (BSREM), for improving both image quality and quantitative accuracy of these low-count PET imaging cases. The NEMA/IEC Body phantom filled with aqueous solution of Y-90 or Ga-68 was scanned to mimic the low-count scenarios of corresponding patient data acquisitions on a time-of-flight (TOF) PET/magnetic resonance imaging system. Contrast recovery, background variation, and signal-to-noise ratio were evaluated in different sets of count densities using both conventional TOF ordered subset expectation (TOF-OSEM) and TOF-BSREM algorithms. The regularization parameter, beta, in BSREM that controls the tradeoff between image noise and resolution was evaluated to find a value for improved confidence in image interpretation. Visual quality assessment of the images obtained from patients administered with [68Ga]citrate (n = 6) was performed. We also made preliminary visual image quality assessment for one patient with [90Y]microspheres. In Y-90 imaging, the effect of 511-keV energy window selection for minimizing the number of random events was also evaluated. Quantitatively, phantom images reconstructed with TOF-BSREM showed improved contrast recovery, background variation, and signal-to-noise ratio values over images reconstructed with TOF-OSEM. Both phantom and patient studies of delayed imaging of [68Ga]citrate show that TOF-BSREM with beta = 500 gives the best tradeoff between image noise and image resolution based on visual assessment by the readers. The NEMA-IQ phantom study with [90Y]microspheres shows that the narrow energy window (460-562keV) recovers activity concentrations in small spheres better than the regular energy window (425-650keV) with the beta value of 2000 using the TOF-BSREM algorithm. For the images obtained from patients with [68Ga]citrate using TOF-BSREM with beta = 500, the visual analogue scale (VAS) was improved by 17% and the Likert score was increased by 1 point on average, both in comparison to corresponding scores for images reconstructed using TOF-OSEM. Our investigation shows that the TOF-BSREM algorithm improves the image quality and quantitative accuracy in low-count PET imaging scenarios. However, the beta value in this algorithm needed to be adjusted for each radiopharmaceutical and counting statistics at the time of scans.

Two-shot X-ray dark-field imaging.

In this article, we report on a novel acquisition scheme for time- and dose-saving retrieval of dark-field data in grating-based phase-contrast imaging. In comparison to currently available techniques, the proposed approach only requires two phase steps. More importantly, our method is capable of accurately retrieving the dark-field signal where conventional approaches fail, for instance in the case of very low photon statistics. Finally, we successfully extend two-shot dark-field imaging to tomographic investigations, by implementing an iterative reconstruction with appropriate weights. Our results indicate an important progression towards the clinical feasibility of dark-field tomography.

A concise quantum efficiency measurement system for gaseous photomultipliers

We have been developing gaseous photomultiplier tubes (PMTs) with a bialkali photocathode combined with micropattern gas detectors (MPGDs). The sensitivity of a PMT is described in terms of its quantum efficiency (QE). The QE is very important in many experiments, particularly when dealing with low photon statistics. A concise QE measurement system has been developed to evaluate the characteristics of a sealed gaseous PMT with a bialkali photocathode. This QE measurement system consists of a Xe arc lamp source, reflective optics and filters, a monochromator, and a reference Si photodiode detector. Using the system, we evaluated the QE, gain, and long-term stability of the gaseous PMT. Here, we report the results of our evaluation.

Single-Molecule Fluorescence Resonance Energy Transfer in Nanopipets: Improving Distance Resolution and Concentration Range

In recent years fluorescence resonance energy transfer (FRET) has widely been used to measure distances, binding, and distance dynamics at the single-molecule (sm) level. Some basic constraints of smFRET are the limited distance resolution owing to low photon statistics and the restriction to high affinity interactions. We demonstrate that by confining molecules in nanopipets with an inner diameter of approximately 100 nm at the tip, FRET can be measured with improved photon statistics and at up to 50-fold higher concentrations. The flow of the donor/acceptor (Cy3B/ATTO647N) labeled double-stranded DNA conjugates was established by electrokinetic forces. Because of the small inner diameter of the nanopipet, every molecule passing the tip is detected applying alternating laser excitation (ALEX). Thus, the technique offers the advantage to study interactions with smaller association constants (<10(9) M-1) using minute sample amounts (<5 microL). The improved photon statistics reduces shot-noise contributions and results in sharper FRET distributions. Experimental results are supported by Monte Carlo simulations which also explain the occurrence of two populations in burst size distributions measured in nanopipet experiments. Because of the confinement of the molecules in nanopipets, the widths of FRET histograms are reduced to a degree where shot-noise is not the only limiting factor but also conformational dynamics of the linkers used to attach the chromophores have to be considered. In addition, our experiments emphasize the influence of photoinduced dark states on both the mean energy transfer efficiency and the width of FRET histograms.

X-ray pulsar radiation from polar caps heated by back-flow bombardment

Abstract We consider the problem of the thermal X-ray radiation from the hot polar caps of radio pulsars that show evidence of E×B subpulse drift in radio band. In our recent Paper I, using the partially screened gap (PSG) model of inner acceleration region we derived a simple relationship between the drift rate of subpulses observed in a radio-band and the thermal X-ray luminosity from polar caps heated by the back-flow particle bombardment. This relationship can be tested for pulsars in which the so-called carousel rotation time P4, reflecting the E×B plasma drift, and the thermal X-ray luminosity Lx from the hot polar cap are known. To test the model we used only two available pulsars: PSRs B0943+10 and B1133+16. They both satisfied the model prediction, although due to low photon statistics the thermal component could not be firmly identified from the X-ray data. Nevertheless, these pulsars were at least consistent with PSG pulsar model. In this Letter we consider two more pulsars: PSRs B0656+14 and B0628−28, the data for which have recently become available. In PSR B0656+14 the thermal radiation from the hot polar cap was clearly detected, and PSR B0628−28 also seems to have such a component. In all cases for which both P4 and Lx are presently known, the PSG pulsar model seems to be fully confirmed. Other available models of inner acceleration region fail to explain the observed relationship between radio and X-ray data. The pure vacuum gap model predicts too high Lx and too low P4, while the space charge limited model predicts too low Lx and the origin of the subpulse drift has no natural explanation.

A motion-incorporated reconstruction method for gated PET studies

Cardiac and respiratory motion artefacts in PET imaging have been traditionally resolved by acquiring the data in gated mode. However, gated PET images are usually characterized by high noise content due to their low photon statistics. In this paper, we present a novel 4D model for the PET imaging system, which can incorporate motion information to generate a motion-free image with all acquired data. A computer simulation and a phantom study were conducted to test the performance of this approach. The computer simulation was based on a digital phantom that was continuously scaled during data acquisition. The phantom study, on the other hand, used two spheres in a tank of water, all of which were filled with 18F water. One of the spheres was stationary while the other moved in a sinusoidal fashion to simulate tumour motion in the thorax. Data were acquired using both 4D CT and gated PET. Motion information was derived from the 4D CT images and then used in the 4D PET model. Both studies showed that this 4D PET model had a good motion-compensating capability. In the phantom study, this approach reduced quantification error of the radioactivity concentration by 95% when compared to a corresponding static acquisition, while signal-to-noise ratio was improved by 210% when compared to a corresponding gated image.