Maximum-likelihood estimators in physics-informed neural networks for high-dimensional inverse problems
Not provided.
Computer Science↗
SEARCH · Search NASA
Results for “maximum-likelihood”
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.
Measurement of the top-quark mass using decays with a J/ψ meson at s=13 TeV with the ATLAS detector
The top-quark mass is measured using top-quark decays producing an isolated lepton and J/ψ meson reconstructed in its μ+μ− decay mode. The data sample was recorded with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of s=13$$ \sqrt{s}=13 $$ TeV during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 140 fb−1. The measurement is based on the invariant mass m(ℓμ+μ−) of the system made of the isolated lepton ℓ from the W boson decay and the non-isolated μ+μ− pair from a J/ψ decay of a b-hadron, exploiting its sensitivity to the top-quark mass. An unbinned maximum-likelihood fit to the m(ℓμ+μ−) distribution is performed to extract the top-quark mass. The top-quark mass is measured to be mtop = 172.17 ± 0.80 (stat) ± 0.81 (syst) ± 1.07 (recoil) GeV, with a total uncertainty of 1.56 GeV. The third uncertainty arises from changing the dipole parton shower gluon-recoil scheme used in top-quark decays.
Aad, G↗
Measurements of W + W − production cross-sections in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector
Measurements of W + W − → e ± νμ ∓ ν production cross-sections are presented, providing a test of the predictions of perturbative quantum chromodynamics and the electroweak theory. The measurements are based on data from pp collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS detector at the Large Hadron Collider in 2015–2018, corresponding to an integrated luminosity of 140 fb −1 . The number of events due to top-quark pair production, the largest background, is reduced by rejecting events containing jets with b-hadron decays. An improved methodology for estimating the remaining top-quark background enables a precise measurement of W + W − cross-sections with no additional requirements on jets. The fiducial W + W − cross-section is determined in a maximum-likelihood fit with an uncertainty of 3.1%. The measurement is extrapolated to the full phase space, resulting in a total W + W − cross-section of 127 ± 4 pb. Differential cross-sections are measured as a function of twelve observables that comprehensively describe the kinematics of W + W − events. The measurements are compared with state-of-the-art theory calculations and excellent agreement with predictions is observed. A charge asymmetry in the lepton rapidity is observed as a function of the dilepton invariant mass, in agreement with the Standard Model expectation. A CP-odd observable is measured to be consistent with no CP violation. Limits on Standard Model effective field theory Wilson coefficients in the Warsaw basis are obtained from the differential cross-sections.
Accelerator Physics↗
A broad specificity β-propeller enzyme from Rhodopseudomonas palustris that hydrolyzes many lactones including γ-valerolactone
Lactones are prevalent in biological and industrial settings, yet there is a lack of information regarding enzymes used to metabolize these compounds. One compound, γ-valerolactone (GVL), is used as a solvent to dissolve plant cell walls into sugars and aromatic molecules for subsequent microbial conversion to fuels and chemicals. Despite the promise of GVL as a renewable solvent for biomass deconstruction, residual GVL can be toxic to microbial fermentation. Here, we identified a Ca 2+ -dependent enzyme from Rhodopseudomonas palustris (Rpa3624) and showed that it can hydrolyze aliphatic and aromatic lactones and esters, including GVL. Maximum-likelihood phylogenetic analysis of other related lactonases with experimentally determined substrate preferences shows that Rpa3624 separates by sequence motifs into a subclade with preference for hydrophobic substrates. Additionally, we solved crystal structures of this β-propeller enzyme separately with either phosphate, an inhibitor, or a mixture of GVL and products to define an active site where calcium-bound water and calcium-bound aspartic and glutamic acid residues make close contact with substrate and product. Our kinetic characterization of WT and mutant enzymes combined with structural insights inform a reaction mechanism that centers around activation of a calcium-bound water molecule promoted by general base catalysis and close contacts with substrate and a potential intermediate. Similarity of Rpa3624 with other β-propeller lactonases suggests this mechanism may be relevant for other members of this emerging class of versatile catalysts.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Development of a compact fast-neutron spectrometer for nuclear emergency response applications
We have developed a Compact Fast Neutron Spectrometer (CFNS) for passive assay of special nuclear material (SNM) through the observation of fast neutrons. The CFNS consists of eight organic glass scintillators (OGS) coupled to silicon photomultipliers and a waveform digitizer, which are integrated within a human-portable box. The CFNS determines the neutron energy profile by spectrum unfolding using the Maximum-Likelihood Expectation Maximization method. The detector acquisition system was optimized to have a dynamic range of up to 10 MeV neutron energy. Bulk special nuclear material (SNM) measurements from the National Criticality Experiments Research Center were analyzed for SNM validation/examination. Additionally, the results show that the CFNS can be used to distinguish between fission and (α, n) neutron emitters, regardless of intervening material type (Cu and polyethylene) and thickness, by taking the ratio of neutron counts at different regions in the unfolded energy spectrum. Additionally, by fitting an exponential curve to the unfolded energy spectrum of PuO 2 and Pu neutron emitters, the CFNS showed the ability of distinguishing between pure Pu oxide, pure Pu metal and mixed oxide-metal configurations.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Community detection in hypergraphs via mutual information maximization
Abstract The hypergraph community detection problem seeks to identify groups of related vertices in hypergraph data. We propose an information-theoretic hypergraph community detection algorithm which compresses the observed data in terms of community labels and community-edge intersections. This algorithm can also be viewed as maximum-likelihood inference in a degree-corrected microcanonical stochastic blockmodel. We perform the compression/inference step via simulated annealing. Unlike several recent algorithms based on canonical models, our microcanonical algorithm does not require inference of statistical parameters such as vertex degrees or pairwise group connection rates. Through synthetic experiments, we find that our algorithm succeeds down to recently-conjectured thresholds for sparse random hypergraphs. We also find competitive performance in cluster recovery tasks on several hypergraph data sets.
97 MATHEMATICS AND COMPUTING↗
Sub-10 second fly-scan nano-tomography using machine learning
Abstract X-ray computed tomography is a versatile technique for 3D structure characterization. However, conventional reconstruction algorithms require that the sample not change throughout the scan, and the timescale of sample dynamics must be longer than the data acquisition time to fulfill the stable sample requirement. Meanwhile, concerns about X-ray-induced parasite reaction and sample damage have driven research efforts to reduce beam dosage. Here, we report a machine-learning-based image processing method that can significantly reduce data acquisition time and X-ray dose, outperforming conventional approaches like Filtered-Back Projection, maximum-likelihood, and model-based maximum-a-posteriori probability. Applying machine learning, we achieve ultrafast nano-tomography with sub-10 second data acquisition time and sub-50 nm pixel resolution in a transmission X-ray microscope. We apply our algorithm to study dynamic morphology changes in a lithium-ion battery cathode under a heating rate of 50 o C min −1 , revealing crack self-healing during thermal annealing. The proposed method can be applied to various tomography modalities.
36 MATERIALS SCIENCE↗
BEYONDPLANCK II. CMB mapmaking through Gibbs sampling
We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.
79 ASTRONOMY AND ASTROPHYSICS↗
Early-time γ -ray constraints on cosmic-ray acceleration in the core-collapse SN 2023ixf with the Fermi Large Area Telescope
Context. While supernova remnants (SNRs) have been considered the most relevant Galactic cosmic ray (CR) accelerators for decades, core-collapse supernovae (CCSNe) could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV γ -rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. Aims. The multi-wavelength observation of SN 2023ixf, a Type II supernova (SN) in the nearby galaxy M 101 (at a distance of 6.85 Mpc), opens the possibility to constrain CR acceleration within a few days after the collapse of the red super-giant stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. Methods. We performed a maximum-likelihood analysis of γ -ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy, non-thermal emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. Results. We do not find significant γ -ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Conclusions. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. Consequently, a more sophisticated, inhomogeneous characterisation of the shock and the progenitor’s environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times.
79 ASTRONOMY AND ASTROPHYSICS↗
A knock-on deuteron imager for measurements of fuel and hotspot asymmetry in direct-drive inertial confinement fusion implosions (invited)
A knock-on deuteron imager (KoDI) has been implemented to measure the fuel and hotspot asymmetry of cryogenic inertial confinement fusion implosions on OMEGA. Energetic neutrons produced by D–T fusion elastically scatter (“knock on”) deuterons from the fuel layer with a probability that depends on ρR. Deuterons above 10 MeV are produced by near-forward scattering, and imaging them is equivalent to time-integrated neutron imaging of the hotspot. Deuterons below 6 MeV are produced by a combination of side scattering and ranging in the fuel, and encode information about the spatial distribution of the dense fuel. Here, the KoDI instrument consists of a multi-penumbral aperture positioned 10–20 cm from the implosion using a ten-inch manipulator and a detector pack at 350 cm from the implosion to record penumbral images with magnification of up to 35×. Range filters and the intrinsic properties of CR-39 are used to distinguish different charged-particle images by energy along the same line of sight. Image plates fielded behind the CR-39 record a 10 keV x-ray image using the same aperture. A maximum-likelihood reconstruction algorithm has been implemented to infer the source from the projected penumbral images. The effects of scattering and aperture charging on the instrument point-spread function are assessed. Synthetic data are used to validate the reconstruction algorithm and assess an appropriate termination criterion. Significant aperture charging has been observed in the initial experimental dataset, and increases with aperture distance from the implosion, consistent with a simple model of charging by laser-driven EMP.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Characterization of a novel HIV-1 circulating recombinant form, CRF91_cpx, comprising CRF02_AG, G, J, and U, mostly among men who have sex with men
Prospective molecular studies of HIV-1 pol region (2253–5250 in HXB2 genome) sequences from sequenced samples of 269 HIV-1-infected patients in Cyprus (2017–2021) revealed a transmission cluster of 14 unknown HIV-1 recombinants that were not classified as previously established CRFs. The earliest recombinant was collected in September 2017, and the transmission cluster continued to grow until November 2020. Near full-length HIV-1 genome sequences of the 11 of the 14 recombinants were successfully obtained (790–8795 in HXB2 genome) and aligned against a reference dataset of HIV-1 subtypes and CRFs. We employed MEGAX for maximum-likelihood tree construction (GTR model, 1000 bootstrap replicates), Cluster-Picker for phylogenetic clustering analysis (genetic distance ≤0.045, bootstrap support value ≥70%), and REGA-3.0 for subtype determination. Bootscan and similarity plot analyses (sliding window of 400 nucleotides overlapped by 40 nucleotides) were conducted using SimPlot-v3.5.1, and subregion confirmatory neighbour-joining tree analyses were conducted using MEGAX (Kimura two-parameter model, 1000 bootstrap replicates, ≥70% bootstrap-support value). Exclusive clustering of the HIV-1 recombinants revealed their uniqueness. The recombination analyses illustrated the same unique mosaic pattern with six putative intersubtype recombination breakpoints, seven fragments of subtypes CRF02_AG, G, J and an unclassified fragment. We conclusively characterized the mosaic structure of the novel HIV-1 CRF, named CRF91_cpx, by the Los Alamos HIV Sequence Database. Additionally, we identified a URF of CRF91_cpx with two additional recombination sites, generated by a recombination event between subtype B and CRF91_cpx. Since the identification of CRF91_cpx, two additional patient samples have been entered into the CRF91_cpx transmission cluster, demonstrating active growth.
60 APPLIED LIFE SCIENCES↗
Monte Carlo method for constructing confidence intervals with unconstrained and constrained nuisance parameters in the NOvA experiment
Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. Wilks' theorem provides a simple way to construct confidence intervals on model parameters, but it only applies under certain conditions. These conditions, such as nested hypotheses and unbounded parameters, are often violated in neutrino oscillation measurements and other experimental scenarios. Monte Carlo methods can address these issues, albeit at increased computational cost. In the presence of nuisance parameters, however, the best way to implement a Monte Carlo method is ambiguous. Furthermore, this paper documents the method selected by the NOvA experiment, the profile construction. It presents the toy studies that informed the choice of method, details of its implementation, and tests performed to validate it. It also includes some practical considerations which may be of use to others choosing to use the profile construction.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Fast decay of classification error in variational quantum circuits
Variational quantum circuits (VQCs) have shown great potential in near-term applications. However, the discriminative power of a VQC, in connection to its circuit architecture and depth, is not understood. To unleash the genuine discriminative power of a VQC, we propose a VQC system with the optimal classical post-processing—maximum-likelihood estimation on measuring all VQC output qubits. Via extensive numerical simulations, we find that the error of VQC quantum data classification typically decays exponentially with the circuit depth, when the VQC architecture is extensive—the number of gates does not shrink with the circuit depth. This fast error suppression ends at the saturation towards the ultimate Helstrom limit of quantum state discrimination. On the other hand, non-extensive VQCs such as quantum convolutional neural networks are sub-optimal and fail to achieve the Helstrom limit, demonstrating a trade-off between ansatz complexity and classification performance in general. To achieve the best performance for a given VQC, the optimal classical post-processing is crucial even for a binary classification problem. To simplify VQCs for near-term implementations, we find that utilizing the symmetry of the input properly can improve the performance, while oversimplification can lead to degradation.
Zhang, Bingzhi↗
Joint inference of multiplicative and additive systematics in galaxy density fluctuations and clustering measurements
Galaxy clustering measurements are a key probe of the matter density field in the Universe. With the era of precision cosmology upon us, surveys rely on precise measurements of the clustering signal for meaningful cosmological analysis. However, the presence of systematic contaminants can bias the observed galaxy number density, and thereby bias the galaxy two-point statistics. As the statistical uncertainties get smaller, correcting for these systematic contaminants becomes increasingly important for unbiased cosmological analysis. We present and validate a new method for understanding and mitigating both additive and multiplicative systematics in galaxy clustering measurements (two-point function) by joint inference of contaminants in the galaxy overdensity field (one-point function) using a maximum-likelihood estimator (MLE). We test this methodology with Kilo-Degree Survey-like mock galaxy catalogues and synthetic systematic template maps. We estimate the cosmological impact of such mitigation by quantifying uncertainties and possible biases in the inferred relationship between the observed and the true galaxy clustering signal. Our method robustly corrects the clustering signal to the sub-percent level and reduces numerous additive and multiplicative systematics from 1.5σ to less than 0.1σ for the scenarios we tested. In addition, we provide an empirical approach to identifying the functional form (additive, multiplicative, or other) by which specific systematics contaminate the galaxy number density. Even though this approach is tested and geared towards systematics contaminating the galaxy number density, the methods can be extended to systematics mitigation for other two-point correlation measurements.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
The emergence and transmission dynamics of HIV-1 CRF07_BC in Mainland China
A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.
60 APPLIED LIFE SCIENCES↗
Seeing through noise in power laws
Despite widespread claims of power laws across the natural and social sciences, evidence in data is often equivocal. Modern data and statistical methods reject even classic power laws such as Pareto’s law of wealth and the Gutenberg–Richter law for earthquake magnitudes. We show that the maximum-likelihood estimators and Kolmogorov–Smirnov (K-S) statistics in widespread use are unexpectedly sensitive to ubiquitous errors in data such as measurement noise, quantization noise, heaping and censorship of small values. This sensitivity causes spurious rejection of power laws and biases parameter estimates even in arbitrarily large samples, which explains inconsistencies between theory and data. We show that logarithmic binning by powers of λ > 1 attenuates these errors in a manner analogous to noise averaging in normal statistics and that λ thereby tunes a trade-off between accuracy and precision in estimation. Binning also removes potentially misleading within-scale information while preserving information about the shape of a distribution over powers of λ, and we show that some amount of binning can improve sensitivity and specificity of K-S tests without any cost, while more extreme binning tunes a trade-off between sensitivity and specificity. We therefore advocate logarithmic binning as a simple essential step in power-law inference.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
Improved Heralded Single-Photon Source with a Photon-Number-Resolving Superconducting Nanowire Detector
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.
Davis, Samantha I.↗
Dark matter trigger for early dark energy coincidence
Current cosmological measurements present a persistent tension in the value of the current cosmic expansion rate, the Hubble constant, as inferred from cosmic microwave background (CMB) and large-scale structure (LSS) data compared to that inferred from the classical distance ladder. Early dark energy (EDE), whose cosmological role is localized in time around the epoch of matter-radiation equality just prior to the release of the CMB photons, is designed to resolve this "Hubble tension". However, the model introduces a new coincidence problem: Why should the EDE dynamics occur near matter-radiation equality if EDE is decoupled from both matter and radiation? The resolution of this problem may lie in an early dark sector (EDS), wherein the dark matter mass is dependent on the EDE scalar field. In this work, we construct such an EDS model and show that it naturally resolves the EDE coincidence problem at the background level without any fine-tuning of the coupling to dark matter or of the initial conditions. When fitting to current cosmological data, including that from the local distance ladder, CMB, and LSS, our EDS maximum-likelihood model performs comparably to EDE for resolving the Hubble tension. However, fitting the Planck CMB data requires a specific range of initial field positions to balance the scalar field fluctuations that drive acoustic oscillations, providing testable differences with other EDE models and a platform for future model-building.
79 ASTRONOMY AND ASTROPHYSICS↗