SEARCH · Engineering Papers
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.
This paper reexamines the notion of closed loop carrier phase synchronization motivated by the theory of maximum a posterior (MAP) phase estimation with emphasis on the development of new structures based on both maximum likelihood (ML) and average-likelihood (AL) functions.
This article reexamines the notion of closed-loop carrier phase synchronization motivated by the theory of maximum a posteriori phase estimation with emphasis on the development of new structures based on both maximum-likelihood and average-likelihood functions. The criterion of performance used for comparison of all the closed-loop structures discussed is the mean-squared phase error for a fixed-loop bandwidth.
The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems, The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA I pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models. within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle. maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius. under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry
In many cosmological models, the large angular scale anisotropy in the cosmic microwave background is parameterized by a spectral index, n, and a quadrupolar amplitude, Q. For a Harrison-Peebles-Zel'dovich spectrum, n = 1. Using data from the Far Infrared Survey (FIRS) and a new statistical measure, a contour plot of the likelihood for cosmological models for which -1 less than n less than 3 and 0 equal to or less than Q equal to or less than 50 micro K is obtained. Depending upon the details of the analysis, the maximum likelihood occurs at n between 0.8 and 1.4 and Q between 18 and 21 micro K. Regardless of Q, the likelihood is always less than half its maximum for n less than -0.4 and for n greater than 2.2, as it is for Q less than 8 micro K and Q greater than 44 micro K.
Here this manuscript details measurements of the anisotropic light-output response of trans-stilbene, henceforth referred to as stilbene, along the three primary crystal planes and analyzes how this anisotropic response impacts neutron image reconstruction in a handheld dual particle imager (H2DPI). A prototype H2DPI composed of stilbene pillars (6 x 6 x 50 mm 3 ) coupled to silicon photomultipliers was built and is capable of imaging both fast neutrons and gamma rays from kilogram quantities of special nuclear material. Stilbene was chosen as the scintillating medium for the imager because of its pulse shape discrimination capability, time resolution and relatively high light output. A drawback to using stilbene, however, is the anisotropic response. A recoiling proton from a neutron elastic-scattering event in stilbene will yield different amounts of scintillation light depending on the direction of the recoiling proton with respect to the crystal lattice. This manuscript analyzes how this anisotropic response impacts neutron image reconstruction in 4π. The light output in the three primary crystal planes of stilbene were independently measured 6–7 times for proton recoil energies ranging from 0.5–5.0 MeV using a quasi-monoenergetic time-of-flight neutron source. The measured light-output data in the three crystal planes were fit with a semi-empirical function based on Birks’ formula. These fits were used to reconstruct the location of a 252 Cf source in front of the imager (0°, -0.8°), to the left of the imager (-90°, -0.8°) and directly above the imager (0°, 90°) to determine if it is necessary to apply the directionally-dependent response of stilbene to accurately reconstruct source locations. A bootstrapping technique was applied to the measured data sets to produce 1000 images composed of 2000 cone projections for each source location and each measured light-output curve. List-mode maximum likelihood expectation maximization was applied to each image; the highest pixel location and the full width at half maximum (FWHM) were then recorded. The average and standard deviations of these parameters were taken for each set of images. The maximum difference in the average azimuthal pixel location when varying between the appropriate and incorrect light-output curves for the source locations at (0°, -0.8°) and (-90°, -0.8°) were respectively found to be 1.85 ± 1.39° and 2.03 ± 1.78°. All other pixel locations and FWHMs were within a single standard deviation of uncertainty. Taken together, the anisotropic response of stilbene has negligible impact on the neutron image reconstruction capability of the H2DPI.
The advent of deep space small spacecraft, as exemplified by the Mars Cubesat One (MarCO), Lunar Trailblazer, Janus, the Escape and Plasma Acceleration and Dynamics Explorers (EscaPADE), and the thirteen Artemis 1 missions, opens the possibility that a much larger number of deep space spacecraft may be launched over the next 10 years and beyond. While scientifically exciting, the prospect of a (much) larger mission suite raises significant challenges for the current approach to ground stations and mission operations. We have been investigating an integrated approach for ground stations and missions operations to enable new modes of operation while maintaining the capabilities of the current operational techniques. This integrated approach is built around three core capabilities: (1) A queuing antenna that enables monitoring the status of a much larger number of spacecraft, and allows spacecraft to transmit requests for telemetry with NASA’s Deep Space Network (DSN); (2) a flexible scheduling system that expands the current DSN scheduling services to enable allocating time on DSN antennas in near real-time; and (3) a cloud-based ground data system that can be spun up and down according to how tracks are assigned by the flexible scheduling system. We shall show that an 18 meter DSN queuing antenna equipped with cyrogenic receivers would enable use of the DSN Demand Access Service for small spacecraft throughout the inner Solar System, thus providing service to a large mission suite. We first discuss the architecture of the queuing antenna and its supporting systems, including, for instance, the service required to generate the schedule for the queueing antenna (which dictates how it slews to monitor multiple spacecraft in a day of operations). Next, we describe the signaling scheme used to encode a request, which is inherited from the already operational DSN Beacon Tone Service, and describe two alternative ways to detect the incoming tone at the ground station, one based on maximum likelihood estimation (MLE), and another one based on Fast-Fourier Transfer (FFT) processing. We then use these results to estimate the maximum range at which a request can be reliably detected as a function of the spacecraft and ground station communication capabilities. Finally, the last part of this part of this paper briefly describes the prototyping effort undertaken at Morehead State University (MSU) and JPL to demonstrate the viability of this new DSN demand access. In particular, we describe the suite of tests conducted using MSU’s 21 meter ground station to validate its use a queuing antenna.
An unsupervised selection of polarimetric features useful for the segmentation and analysis of polarimetric synthetic aperture radar (SAR) data is presented. The technique is based on multidimensional clustering of the parameters composing the polarimetric covariance matrix of the data. Clustering is performed on the logarithm of these quantities. Once the polarimetric cluster centers have been determined, segmentation of the polarimetric data into regions is performed using a maximum likelihood polarimetric classifier. Segmentation maps are further improved using a Markov random field to describe the statistics of the regions and computing the maximum of the product of the local conditional densities. Examples with real polarimetric SAR imagery are given to illustrate the potential of this method.
Event reconstruction is a central step in many particle physics experiments, turning detector observables into parameter estimates; for example estimating the energy of an interaction given the sensor readout of a detector. A corresponding likelihood function is often intractable, and approximations need to be constructed. Here, in our work, we first show how the full likelihood for a many-sensor detector can be broken apart into smaller terms, and secondly how we can train neural networks to approximate all terms solely based on forward simulation. Our technique results in a fast, flexible, and close-to-optimal surrogate model proportional to the likelihood and can be used in conjunction with standard inference techniques allowing for a consistent treatment of uncertainties. We illustrate our technique for parameter inference in neutrino telescopes based on maximum likelihood and Bayesian posterior sampling. Given its great flexibility, we also showcase our method for geometry optimization enabling to learn optimal detector designs. Lastly, we apply our method to realistic simulation of a ton-scale water-based liquid scintillator detector.
A considerable experimental effort is currently under way to test the persistent hints for oscillations due to an eV-scale sterile neutrino in the data of various reactor neutrino experiments. The assessment of the statistical significance of these hints is usually based on Wilks’ theorem, whereby the assumption is made that the log-likelihood is χ 2 -distributed. However, it is well known that the preconditions for the validity of Wilks’ theorem are not fulfilled for neutrino oscillation experiments. In this work we derive a simple asymptotic form of the actual distribution of the log-likelihood based on reinterpreting the problem as fitting white Gaussian noise. From this formalism we show that, even in the absence of a sterile neutrino, the expectation value for the maximum likelihood estimate of the mixing angle remains non-zero with attendant large values of the log-likelihood. Our analytical results are then confirmed by numerical simulations of a toy reactor experiment. Finally, we apply this framework to the data of the Neutrino-4 experiment and show that the null hypothesis of no-oscillation is rejected at the 2.6 σ level, compared to 3.2 σ obtained under the assumption that Wilks’ theorem applies.
Deterministic generation of single photons is essential for many quantum information technologies. A bulk optical nonlinearity emitting a photon pair, where the measurement of one of the photons heralds the presence of the other, is commonly used with the caveat that the single-photon emission rate is constrained due to a trade-off between multiphoton events and pair emission rate. Using an efficient and low noise photon-number-resolving superconducting nanowire detector we herald, in real time, a single photon at telecommunication wavelength. We perform a second-order photon correlation g 2 ( 0 ) measurement of the signal mode conditioned on the measured photon number of the idler mode for various pump powers and demonstrate an improvement of a heralded single-photon source. We develop an analytical model using a phase-space formalism that encompasses all multiphoton effects and relevant imperfections, such as loss and multiple Schmidt modes. We perform a maximum-likelihood fit to test the agreement of the model to the data and extract the best-fit mean photon number μ of the pair source for each pump power. A maximum reduction of 0.118 ± 0.012 in the photon g 2 ( 0 ) correlation function at μ = 0.327 ± 0.007 is obtained, indicating a strong suppression of multiphoton emissions. For a fixed g 2 ( 0 ) = 7 × 10 − 3 , we increase the single pair generation probability by 25%. Our experiment, built using fiber-coupled and off-the-shelf components, delineates a path to engineering ideal sources of single photons.
For the recapture debugging design introduced by Nayak (1988) the problem of estimating the hitting rates of the faults remaining in the system is considered. In the context of a conditional likelihood, moment estimators are derived and are shown to be asymptotically normal and fully efficient. Fixed sample properties of the moment estimators are compared, through simulation, with those of the conditional maximum likelihood estimators. Properties of the conditional model are investigated such as the asymptotic distribution of linear functions of the fault hitting frequencies and a representation of the full data vector in terms of a sequence of independent random vectors. It is assumed that the residual hitting rates follow a log linear rate model and that the testing process is truncated when the gaps between the detection of new errors exceed a fixed amount of time.
Early results are reported from measurements of the diurnal variability of total and clear-sky regional radiative parameters by the ERBE instruments on one dedicated satellite and the polar-orbiting NOAA-9 satellite. Attention is focused on November 1984, the first complete data set. The scene is identified in terms of longwave and shortwave radiances (daytime) or longwave radiation (night) and maximum likelihood estimates carried out with the addition of Earth Radiation Budget data from Nimbus-7. Analysis of the first data set revealed significant differences between total and clear-sky albedo. The clear-sky and LRE both reach maximum around noon and minimum values at midnight.
An array feed combining system for the recovery of signal-to-noise ratio (SNR) loss due to antenna reflector deformation has been implemented and is currently being evaluated on the Jet Propulsion Laboratory 34-m DSS-13 antenna. The current signal-combining system operates under the assumption that the white Gaussian noise processes in the received signals from different array elements are mutually uncorrelated. However, experimental data at DSS 13 indicate that these noise processes are indeed mutually correlated. The objective of this work is to develop a signal-combining system optimized to account for the mutual correlations between these noise processes. The set of optimum combining weight coefficients that maximizes the combined signal SNR in the correlated noises environment is determined. These optimum weights depend on unknown signal and noise covariance parameters. A maximum-likelihood approach is developed to estimate these unknown parameters to obtain estimates of the optimum weight coefficients based on residual carrier signal samples. The actual combined signal SNR using the estimated weight coefficients is derived and shown to converge to the maximum achievable SNR as the number of signal samples increases. These results are also verified by simulation. A numerical example shows a significant improvement in SNR performance can be obtained, especially when the amount of correlation increases.
An array feed combining system for the recovery of signal-to-noise ratio (SNR) loss due to antenna reflector deformation has been implemented and is currently being evaluated on the Jet Propulsion Laboratory 34-m DSS-13 antenna. The current signal-combining system operates under the assumption that the white Gaussian noise processes in the received signals from different array elements are mutually uncorrelated. However, experimental data at DSS 13 indicate that these noise processes are indeed mutually correlated. The objective of this work is to develop a signal-combining system optimized to account for the mutual correlations between these noise processes. The set of optimum combining weight coefficients that maximizes the combined signal SNR in the correlated noises environment is determined. These optimum weights depend on unknown signal and noise covariance parameters. A maximum-likelihood approach is developed to estimate these unknown parameters to obtain estimates of the optimum weight coefficients based on residual carrier signal samples. The actual combined signal SNR using the estimated weight coefficients is derived and shown to converge to the maximum achievable SNR as the number of signal samples increases. These results are also verified by simulation. A numerical example shows a significant improvement in SNR performance can be obtained, especially when the amount of correlation increases.
We derive an asymptotic expansion for the log-likelihood of Gaussian mixture models (GMMs) with equal covariance matrices in the low signal-to-noise regime. The expansion reveals an intimate connection between two types of algorithms for parameter estimation: the method of moments and likelihood optimizing algorithms such as Expectation-Maximization (EM). We show that likelihood optimization in the low SNR regime reduces to a sequence of least squares optimization problems that match the moments of the estimate to the ground truth moments one by one. This connection is a stepping stone towards the analysis of EM and maximum likelihood estimation in a wide range of models. A motivating application for the study of low SNR mixture models is cryo-electron microscopy data, which can be modeled as a GMM with algebraic constraints imposed on the mixture centers. We discuss the application of our expansion to algebraically constrained GMMs, among other example models of interest. © 2022 The Authors. Communications on Pure and Applied Mathematics published by Wiley Periodicals LLC.
We study the occurrence frequency of solar flares as a function of the hard X-ray peak count rate, using observations of the Solar Maximum Mission. The size distributions are well represented by power-law distributions with negative indices. As a better alternative to the conventional method, we devise a maximum likelihood method of determining the power-law index of the size distribution. We find that the power-law index of the size distribution changes with time and with the phase of the 154-day periodicity. The size distribution is steeper during the maximum years of solar cycle 21 (1980 and 1981) than during the declining phase (1982-1984). The size distribution, however, is flatter during the maximum phase of the 154-day periodicity than during the minimum phase. The implications of these findings are discussed.
We propose a complete framework for Bayesian image reconstruction and uncertainty quantification based on a Gaussian process prior (GPP) to overcome limitations of maximum likelihood expectation maximization (ML-EM) image reconstruction algorithm. The prior distribution is constructed with a zero-mean Gaussian process (GP) with a choice of a covariance function, and a link function is used to map the Gaussian process to an image. Unlike many other maximum a posteriori approaches, our method offers highly interpretable hyperparamters that are selected automatically with the empirical Bayes method. Furthermore, the GP covariance function can be modified to incorporate a priori structural priors, enabling multi-modality imaging or contextual data fusion. Lastly, we illustrate that our approach lends itself to Bayesian uncertainty quantification techniques, such as the preconditioned Crank–Nicolson method and the Laplace approximation. The proposed framework is general and can be employed in most radiation image reconstruction problems, and we demonstrate it with simulated free-moving single detector radiation source imaging scenarios. We compare the reconstruction results from GPP and ML-EM, and show that the proposed method can significantly improve the image quality over ML-EM, all the while providing greater understanding of the source distribution via the uncertainty quantification capability. Furthermore, significant improvement of the image quality by incorporating a structural prior is illustrated.
The measurement of diffuse PeV gamma-ray emission from the Galactic plane would provide information about the energy spectrum and propagation of Galactic cosmic rays, and the detection of a pointlike source of PeV gamma-rays would be strong evidence for a Galactic source capable of accelerating cosmic rays up to at least a few PeV. This paper presents several unbinned maximum-likelihood searches for PeV gamma-rays in the Southern Hemisphere using 5 yr of data from the IceTop air shower surface detector and the in-ice array of the IceCube Observatory. The combination of both detectors takes advantage of the low muon content and deep shower maximum of gamma-ray air showers and provides excellent sensitivity to gamma-rays between ~0.6 and 100 PeV. Our measurements of pointlike and diffuse Galactic emission of PeV gamma-rays are consistent with the background, so we constrain the angle-integrated diffuse gamma-ray flux from the Galactic plane at 2 PeV to 2.61 × 10 -19 cm -2 s -1 TeV -1 at 90% confidence, assuming an E -3 spectrum, and we estimate 90% upper limits on pointlike emission at 2 PeV between 10 -21 and 10 -20 cm -2 s -1 TeV -1 for an E -2 spectrum, depending on decl. Furthermore, we exclude unbroken power-law emission up to 2 PeV for several TeV gamma-ray sources observed by the High Energy Spectroscopic System and calculate upper limits on the energy cutoffs of these sources at 90% confidence. We also find no PeV gamma-rays correlated with neutrinos from IceCube's high-energy starting event sample. Lastly, these are currently the strongest constraints on PeV gamma-ray emission.