Instrument Response Function Research Articles

Single-photon smFRET. III. Application to pulsed illumination

Förster resonance energy transfer (FRET) using pulsed illumination has been pivotal in leveraging lifetime information in FRET analysis. However, there remain major challenges in quantitative single-photon, single-molecule FRET (smFRET) data analysis under pulsed illumination including 1) simultaneously deducing kinetics and number of system states; 2) providing uncertainties over estimates, particularly uncertainty over the number of system states; and 3) taking into account detector noise sources such as cross talk and the instrument response function contributing to uncertainty; in addition to 4) other experimental noise sources such as background. Here, we implement the Bayesian nonparametric framework described in the first companion article that addresses all aforementioned issues in smFRET data analysis specialized for the case of pulsed illumination. Furthermore, we apply our method to both synthetic as well as experimental data acquired using Holliday junctions.

Construction and study of instrument response functions for analysis of the National Ignition Facility (NIF) neutron time-of-flight detectors.

The analysis of the National Ignition Facility (NIF) neutron time-of-flight (nToF) detectors uses a forward-fit routine that depends critically on the instrument response functions (IRFs) of the diagnostics. The details of the IRFs used can have large impacts on measurements such as ion temperature and down-scattered ratio (DSR). Here, we report on the recent steps taken to construct and validate nToF IRFs at the NIF to an increased degree of accuracy, as well as remove the need for fixed DSR baseline offsets. The IRF is treated in two parts: a "core," measured experimentally with an x-ray impulse source, and a "tail" that occurs later in time and has limited experimental data. The tail region is calibrated with the data from indirect drive exploding pusher shots, which have little neutron scattering and are traditionally assumed to have zero DSR. Using analytic modeling estimates, the non-zero DSR for these shots is estimated. The impact of varying IRF tail components on DSR is investigated with a systematic parameter study, and good agreement is found with the non-zero DSR estimates. These approaches will be used to improve the precision and uncertainty of NIF nToF DSR measurements.

Conceptual Design of a Scintillator-Based Fast-Ion Loss Detector for the Wendelstein 7-X Stellarator

A conceptual design of a scintillator-based fast-ion loss detector (FILD) has been developed for the Wendelstein 7-X stellarator (W7-X). Simulations using the Monte Carlo codes ASCOT5 and BEAMS3D have been performed to calculate the expected flux of neutral beam injection (NBI)-generated fast hydrogen ions onto the conceptual detector probe head. These fast-ion loss fluxes have been calculated for several magnetic field configurations as well as probe insertion positions. At the maximum insertion position, both co- and counter-going losses with high incident pitch angles are observed; however, at retracted positions, only co-going fast ions reach the probe head. The FILDSIM code has been used to optimize the geometry of the detector entrance and collimating elements to achieve a wide velocity space coverage as well as a high velocity–space resolution. A synthetic FILD signal is calculated for the expected loss distribution via forward modeling using the instrument response function. The synthetic signal is found to largely retain the velocity space features of the loss distribution.

Measurements of Rayleigh–Brillouin scattering spectra of Ar, N[formula omitted], and CH[formula omitted] at pressures up to 20 atm

This paper describes the spectral measurement of Rayleigh–Brillouin scattering (RBS) in Ar, N2, and CH4 at 300 K and pressures up to 20 atm, which is relevant for the application of filtered Rayleigh scattering (FRS) diagnostics at elevated pressure. The measurement scheme for RBS spectra was relatively novel, based on an air-spaced virtually imaged phased array (VIPA), capable of rapid acquisition at high resolution. However, the VIPA’s nonuniform dispersion, irregular instrument response function, and sensitivity to local ambient conditions challenged the reconstruction of high-fidelity lineshapes in frequency space. Nonuniform dispersion was addressed using the VIPA dispersion function, and sensitivity of the air-spaced optic to the refractive index of air was minimized by tracking local temperature and pressure during measurements. Irregular instrument response was remedied with a computational model of the VIPA’s optical field, and by deconvolution. A comparison of experimental lineshapes with model predictions by Tenti S6 yielded good agreement. In the kinetic and nominal transition regimes, local residuals were within 3% for Ar and N2, and 5–10% for CH4. Beyond the transition regime, experimental lineshapes continued to display good qualitative agreement, with frequency shifts of the Brillouin peaks representing the main deviation from model prediction. Artifacts in the experimental data from uncertainty and the deconvolution algorithm were identified and discussed. This work establishes an experimental foundation for future measurements of RBS spectra, which will assist in the development of FRS diagnostics of reacting flows at high pressure.

Current and Future γ-Ray Searches for Dark Matter Annihilation Beyond the Unitarity Limit

For decades, searches for electroweak-scale dark matter (DM) have been performed without a definitive detection. This lack of success may hint that DM searches have focused on the wrong mass range. A proposed candidate beyond the canonical parameter space is ultraheavy DM (UHDM). In this work, we consider indirect UHDM annihilation searches for masses between 30 TeV and 30 PeV—extending well beyond the unitarity limit at ∼100 TeV—and discuss the basic requirements for DM models in this regime. We explore the feasibility of detecting the annihilation signature, and the expected reach for UHDM with current and future very-high-energy (VHE; >100 GeV) γ-ray observatories. Specifically, we focus on three reference instruments: two Imaging Atmospheric Cherenkov Telescope arrays, modeled on VERITAS and CTA-North, and one extended air shower array, motivated by HAWC. With reasonable assumptions on the instrument response functions and background rate, we find a set of UHDM parameters (mass and cross section) for which a γ-ray signature can be detected by the aforementioned observatories. We further compute the expected upper limits for each experiment. With realistic exposure times, the three instruments can probe DM across a wide mass range. At the lower end, it can still have a point-like cross section, while at higher masses the DM could have a geometric cross section, indicative of compositeness.

Fast-Gated 16 × 16 SPAD Array With 16 on-Chip 6 ps Time-to-Digital Converters for Non-Line-of-Sight Imaging

We present the design and characterization of a fully-integrated array of <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$16 \times 16$ </tex-math></inline-formula> Single-Photon Avalanche Diodes (SPADs) with fast-gating capabilities and 16 on-chip 6 ps time-to-digital converters, which has been embedded in a compact imaging module. Such sensor has been developed for Non-Line-Of-Sight imaging applications, which require: i) a narrow instrument response function, for a centimeteraccurate single-shot precision; ii) fast-gated SPADs, for time-filtering of directly reflected photons; iii) high photon detection probability, for acquiring faint signals undergoing multiple scattering events. Thanks to a novel multiple differential SPAD-SPAD sensing approach, SPAD detectors can be swiftly activated in less than 500 ps and the full-width at half maximum of the instrument response function is always less than 75 ps (60 ps on average). Temporal responses are consistently uniform throughout the gate window, showing just few picoseconds of time dispersion when 30 ns gate pulses are applied, while the differential non-linearity is as low as 250 fs. With a photon detection probability peak of 70% at 490 nm, a fill-factor of 9.6% and up to <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$1.6 \cdot 10^{8}$ </tex-math></inline-formula> photon time-tagging measurements per second, such sensor fulfills the demand for fully-integrated imaging solutions optimized for non-line-of-sight imaging applications, enabling to cut exposure times while also optimizing size, weight, power and cost, thus paving the way for further scaled architectures.

A new neutron time-of-flight detector for yield and ion-temperature measurements at the OMEGA Laser Facility.

A new neutron time-of-flight (nTOF) detector for deuterium-deuterium (DD)-fusion yield and ion-temperature measurements was designed, installed, and calibrated for the OMEGA Laser Facility. This detector provides an additional line of sight for DD neutron yield and ion-temperature measurements for yields exceeding 1 × 1010 with higher precision than existing detectors. The nTOF detector consists of a 90-mm-diam, 20-mm-thick BC-422 scintillator and a gated Photek photomultiplier tube (PMT240). The PMT collects scintillating light through the 20-mm side of the scintillator without the use of a light guide. There is no lead shielding from hard x rays in order to allow the x-ray instrument response function of the detector to be measured easily. Instead, hard x-ray signals generated in implosion experiments are gated out by the PMT. The design provides a place for glass neutral-density filters between the scintillator and the PMT to avoid PMT saturation at high yields. The nTOF detector is installed in the OMEGA Target Bay along the P8A sub-port line of sight at a distance of 5.3m from the target chamber center. In addition to DD measurements, the same detector can be used to measure the neutron yield and ion temperature from deuterium-tritium (DT) implosion targets in the 5 × 1010 to 2 × 1012 yield range. The design details and the calibration results of this nTOF detector for both D2 and DT implosions on OMEGA will be presented.

Deep Learning-Based Spectrum Reconstruction Method for Raman Spectroscopy

Raman spectroscopy, measured by a Raman spectrometer, is usually disturbed by the instrument response function and noise, which leads to certain measurement error and further affects the accuracy of substance identification. In this paper, we propose a spectral reconstruction method which combines the existing maximum a posteriori (MAP) method and deep learning (DL) to recover the degraded Raman spectrum. The proposed method first employs the MAP method to reconstruct the measured Raman spectra, so as to obtain preliminary estimated Raman spectra. Then, a convolutional neural network (CNN) is trained by using the preliminary estimated Raman spectra and the real Raman spectra to learn the mapping from the preliminary estimated Raman spectra to the real Raman spectra, so as to achieve a better spectral reconstruction effect than merely using the MAP method or a CNN. To prove the effectiveness of the proposed spectral reconstruction method, we employed the proposed method and some traditional spectral reconstruction methods to reconstruct the simulated and measured Raman spectra, respectively. The experimental results show that compared with traditional methods, the estimated Raman spectra reconstructed by the proposed method are closer to the real Raman spectra.

High‐Resolution Time‐Correlated Single‐Photon Counting Using Electro‐Optic Sampling

AbstractA simple, practical method based on electro‐optic gating is experimentally shown to improve the temporal resolution of single‐photon detection by more than 16 times. Delay times between ultrafast single photons and a reference clock are stretched by a desired programmable sampling gate factor, allowing reconstruction of delay histograms with ≈0.001 photons per pulse to within 60 ps. By transferring the bandwidth of RF electronics to single‐photon counting, complex single‐photon signals with large time‐bandwidth products &gt; 2000 are temporally stretched up so that they can be resolved by slow detectors, irrespective of the quality of the detector instrument response function. This method is also applied to biphotons, in order to reconstruct the 2D histogram of joint detection delays with 15 times better resolution. The phenomenon of nonlocal dispersion, which is not resolvable directly with the slow detectors, is then observable to within the 98 ps level. The proof‐of‐concept demonstration uses off‐the‐shelf commercial fiber‐integrated LiNbO3 modulators and RF electronics, and the method is readily integrable on‐chip and to speeds &gt;100 GHz, offering a practical solution to ultrafast time‐correlated single‐photon counting beyond the research laboratory.

Bootstrap estimation of the effect of instrument response function uncertainty on the reconstruction of fusion neutron sources

Neutron imagers are important diagnostics for the inertial confinement fusion implosions at the National Ignition Facility. They provide two- and three-dimensional reconstructions of the neutron source shape that are key indicators of the overall performance. To interpret the shape results properly, it is critical to estimate the uncertainty in those reconstructions. There are two main sources of uncertainties: limited neutron statistics, leading to random errors in the reconstructed images, and incomplete knowledge of the instrument response function (the pinhole-dependent point spread function). While the statistical errors dominate the uncertainty for lower yield deuterium-tritium (DT) shots, errors due to the instrument response function dominate the uncertainty for DT yields on the order of 1016 neutrons or higher. In this work, a bootstrapping method estimates the uncertainty in a reconstructed image due to the incomplete knowledge of the instrument response function. The main reconstruction is created from the fixed collection of pinhole images that are best aligned with the neutron source. Additional reconstructions are then built using subsets of that collection of images. Variations in the shapes of these additional reconstructions originate solely from uncertainties in the instrument response function, allowing us to use them to provide an additional systematic uncertainty estimate.

Simple Double Hetero[5]helicenes Realize Highly Efficient and Narrowband Circularly Polarized Organic Light-Emitting Diodes

Simple Double Hetero[5]helicenes Realize Highly Efficient and Narrowband Circularly Polarized Organic Light-Emitting Diodes

Instrumental Tip-of-the-iceberg Effects on the Prompt Emission of Swift/BAT Gamma-ray Bursts

Abstract The observed durations of prompt gamma-ray emission from gamma-ray bursts (GRBs) are often used to infer the progenitors and energetics of the sources. Inaccurate duration measurements will have a significant impact on constraining the processes powering the bursts. The “tip-of-the-iceberg” effect describes how the observed signal is lost into background noise; lower instrument sensitivity leads to higher measurement bias. In this study, we investigate how observing conditions, such as the number of enabled detectors, background level, and incident angle of the source relative to the detector plane, affect the measured duration of GRB prompt emission observed with the Burst Alert Telescope on board the Neil Gehrels Swift Observatory (Swift/BAT). We generate “simple-pulse” light curves from an analytical fast rise exponential decay function and from a sample of eight real GRB light curves. We fold these through the Swift/BAT instrument response function to simulate light curves Swift/BAT would have observed for specific observing conditions. We find duration measurements are highly sensitive to observing conditions and the incident angle of the source has the highest impact on measurement bias. In most cases duration measurements of synthetic light curves are significantly shorter than the true burst duration. For the majority of our sample, the percentage of duration measurements consistent with the true duration is as low as ∼25%–45%. In this article, we provide quantification of the tip-of-the-iceberg effect on GRB light curves due to Swift/BAT instrumental effects for several unique light curves.

Development of a Monte Carlo-wave model to simulate time domain diffuse correlation spectroscopy measurements from first principles.

.Significance: Diffuse correlation spectroscopy (DCS) is an optical technique that measures blood flow non-invasively and continuously. The time-domain (TD) variant of DCS, namely, TD-DCS has demonstrated a potential to improve brain depth sensitivity and to distinguish superficial from deeper blood flow by utilizing pulsed laser sources and a gating strategy to select photons with different pathlengths within the scattering tissue using a single source–detector separation. A quantitative tool to predict the performance of TD-DCS that can be compared with traditional continuous wave DCS (CW-DCS) currently does not exist but is crucial to provide guidance for the continued development and application of these DCS systems.Aims: We aim to establish a model to simulate TD-DCS measurements from first principles, which enables analysis of the impact of measurement noise that can be utilized to quantify the performance for any particular TD-DCS system and measurement geometry.Approach: We have integrated the Monte Carlo simulation describing photon scattering in biological tissue with the wave model that calculates the speckle intensity fluctuations due to tissue dynamics to simulate TD-DCS measurements from first principles.Results: Our model is capable of simulating photon counts received at the detector as a function of time for both CW-DCS and TD-DCS measurements. The effects of the laser coherence, instrument response function, detector gate delay, gate width, intrinsic noise arising from speckle statistics, and shot noise are incorporated in the model. We have demonstrated the ability of our model to simulate TD-DCS measurements under different conditions, and the use of our model to compare the performance of TD-DCS and CW-DCS under a few typical measurement conditions.Conclusion: We have established a Monte Carlo-Wave model that is capable of simulating CW-DCS and TD-DCS measurements from first principles. In our exploration of the parameter space, we could not find realistic measurement conditions under which TD-DCS outperformed CW-DCS. However, the parameter space for the optimization of the contrast to noise ratio of TD-DCS is large and complex, so our results do not imply that TD-DCS cannot indeed outperform CW-DCS under different conditions. We made our code available publicly for others in the field to find use cases favorable to TD-DCS. TD-DCS also provides a promising way to measure deep brain tissue dynamics using a short source–detector separation, which will benefit the development of technologies including high density DCS systems and image reconstruction using a limited number of source–detector pairs.

Reconstruction of optical coefficients in turbid media using time-resolved reflectance and calibration-free instrument response functions.

Measurements of time-resolved reflectance from a homogenous turbid medium can be employed to retrieve the absolute values of its optical transport coefficients. However, the uncertainty in the temporal shift of the experimentally determined instrument response function (IRF) with respect to the real system response can lead to errors in optical property reconstructions. Instrument noise and measurement of the IRF in a reflectance geometry can exacerbate these errors. Here, we examine three reconstruction approaches that avoid requiring direct measurements of photon launch times. They work by (a) fitting relative shapes of the reflectance profile with a pre-determined constraint on the scattering coefficient, (b) calibrating launch-time differences via a reference sample, and (c) freely fitting for the launch-time difference within the inverse problem. Analysis methods that can place a tight bound on the scattering coefficient can produce errors within 5-15% for both absorption and scattering at source-detector separations of 10 and 15 mm. Including the time-shift in the fitting procedure also recovered optical coefficients to under 20% but showed large crosstalk between extracted scattering and absorption coefficients. We find that the uncertainty in the temporal shift greatly impacts the reconstructed reduced scattering coefficient compared to absorption.

Signal recovery of a Fabry–Pérot interferometric x-ray pulse detector based on the RadOptic effect

The signal recovery of a Fabry–Pérot interferometric x-ray pulse detector based on the RadOptic effect in the non-limiting case was investigated in this research. A Fe-doped InP with an invariant excess carrier recombination mechanism was used as the interference cavity material to achieve a constant temporal instrumental response function (tIRF). A linear and time-invariant detection system described by the convolution of the time-varying x-ray pulse and the constant tIRF was established based on the transient refractive index variation model determined by the three effects of band filling, band shrinkage, and free-carrier absorption. For the non-limiting case, the accumulation of excess carriers enhanced the sensitivity but altered the fluctuations of the real x-ray pulse. To realistically reconstruct the x-ray pulse, two-photon absorption of the infrared ultrashort pulse was used to simulate the ultrashort x-ray excitation to obtain the tIRF. Finally, using the conjugate gradient method, the original signal recorded by the detection system was deconvoluted to recover the signal. The success of signal recovery in the non-limiting case provided the basis for the development of detectors with adjustable sensitivity controlled by carrier lifetime.

High Resolution Fluorescence Lifetime Maps from Minimal Photon Counts.

Fluorescence lifetime imaging microscopy (FLIM) may reveal subcellular spatial lifetime maps of key molecular species. Yet, such a quantitative picture of life necessarily demands high photon budgets at every pixel under the current analysis paradigm, thereby increasing acquisition time and photodamage to the sample. Motivated by recent developments in computational statistics, we provide a direct means to update our knowledge of the lifetime maps of species of different lifetimes from direct photon arrivals, while accounting for experimental features such as arbitrary forms of the instrument response function (IRF) and exploiting information from empty laser pulses not resulting in photon detection. Our ability to construct lifetime maps holds for arbitrary lifetimes, from short lifetimes (comparable to the IRF) to lifetimes exceeding interpulse times. As our method is highly data efficient, for the same amount of data normally used to determine lifetimes and photon ratios, working within the Bayesian paradigm, we report direct blind unmixing of lifetimes with subnanosecond resolution and subpixel spatial resolution using standard raster scan FLIM images. We demonstrate our method using a wide range of simulated and experimental data.

The throughput calibration of the VERITAS telescopes

Context. The response of imaging atmospheric Cherenkov telescopes to incident γ-ray-initiated showers in the atmosphere changes as the telescopes age due to exposure to light and weather. These aging processes affect the reconstructed energies of the events and γ-ray fluxes. Aims. This work discusses the implementation of signal calibration methods for the Very Energetic Radiation Imaging Telescope Array System (VERITAS) to account for changes in the optical throughput and detector performance over time. Methods. The total throughput of a Cherenkov telescope is the product of camera-dependent factors, such as the photomultiplier tube gains and their quantum efficiencies, and the mirror reflectivity and Winston cone response to incoming radiation. This document summarizes different methods to determine how the camera gains and mirror reflectivity have evolved over time and how we can calibrate this changing throughput in reconstruction pipelines for imaging atmospheric Cherenkov telescopes. The implementation is validated against seven years of observations with the VERITAS telescopes of the Crab Nebula, which is a reference object in very-high-energy astronomy. Results. Regular optical throughput monitoring and the corresponding signal calibrations are found to be critical for the reconstruction of extensive air shower images. The proposed implementation is applied as a correction to the signals of the photomultiplier tubes in the telescope simulation to produce fine-tuned instrument response functions. This method is shown to be effective for calibrating the acquired γ-ray data and for recovering the correct energy of the events and photon fluxes. At the same time, it keeps the computational effort of generating Monte Carlo simulations for instrument response functions affordably low.

Intermodal fiber interferometer based on broadband light source and optical spectrum analyzer for external perturbations measurement

An intermodal fiber interferometer based on broadband light source and spectral signal processing is proposed and experimentally investigated. The perturbation of the parameters of the spectrum analyzer (scanning range and width of the instrumental spectral response function) on the output signal is demonstrated. It is shown that the method of correlation signal processing makes it possible to obtain linear and stable response to external fiber perturbations. The proposed interferometer scheme can be prospective for development of real-time sensors of physical quantities. Keywords: intermodal fiber interferometer, fiber-optic sensors, fiber optics, multimode fiber, interferometers.

Radiative dynamics and delayed emission in pure and doped InP/ZnSe/ZnS quantum dots.

We have used time-correlated single photon counting to elucidate the radiative dynamics of InP/ZnSe/ZnS core/shell/shell quantum dots (QDs) that differ in the amount and distribution of excess indium. Stoichiometric QDs having an In:P atom ratio very near unity exhibit simple luminescence kinetics. The photoluminescence (PL) rises with the 40ps instrument response function and exhibits a decay that is close to a single exponential with a time constant that decreases from 32 to 28ns with increasing shell thickness. QDs having excess indium (In:P ratio of 1.15-1.63) show a significant component of a slower rise time assigned to transient population of indium-based hole traps in the ZnSe shell. They also have a slower PL decay, attributed to an equilibrium between these traps, which are optically dark, and the emissive valence-band state. This results in a radiative lifetime that increases from 32 to 48ns with increasing shell thickness. Different treatments of the InP cores prior to shell deposition result in different core/shell interfaces as indicated by resonance Raman spectroscopy, as well as differences in the amplitude and timescale of the slow PL rise and the PL decay time. These are interpreted in terms of different radial distributions of the indium-based hole traps, which can be related to differences in the interfacial lattice strain.

Lifetime-based analysis of binary fluorophores mixtures in the low photon count limit

SummarySingle biomolecule sensing often requires the quantification of multiple fluorescent species. Here, we theoretically and experimentally use time-resolved fluorescence via Time Correlated Single Photon Counting (TCSPC) to accurately quantify fluorescent species with similar chromatic signatures. A modified maximum likelihood estimator is introduced to include two fluorophore species, with compensation of the instrument response function. We apply this algorithm to simulated data of a simplified two-fluorescent species model, as well as to experimental data of fluorophores' mixtures and to a model protein, doubly labeled with different fluorophores' ratio. We show that 100 to 200 photons per fluorophore, in a 10-ms timescale, are sufficient to provide an accurate estimation of the dyes' ratio on the model protein. Our results provide estimation for the desired photon integration time toward implementation of TCSPC in systems with fast occurring events, such as translocation of biomolecules through nanopores or single-molecule burst analyses experiments.