Search NASA⌕ Search

SEARCH · Search NASA

Results for “Marginal 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.

At least 73 records · Page 4

Cosmological parameter forecasts for a CMB-HD survey

We present forecasts on cosmological parameters for a CMB-HD survey. For a $\Lambda$CDM + $N_{eff}$ + $\sum m_\nu$ model, we find $\sigma(n_s) = 0.0013$ and $\sigma(N_{eff}) = 0.014$ using CMB and CMB lensing multipoles in the range of $\ell \in [30, 20000]$, after adding anticipated residual foregrounds, delensing the acoustic peaks, and adding DESI BAO data. This is about a factor of two improvement in ability to probe inflation via $n_s$ compared to precursor CMB surveys. The $N_{eff}$ constraint can rule out light thermal particles back to the end of inflation with 95% CL; for example, it can rule out the QCD axion in a model-independent way assuming the Universe's reheating temperature was high enough that the axion thermalized. We find that delensing the acoustic peaks and adding DESI BAO tightens parameter constraints. We also find that baryonic effects can bias parameters if not marginalized over, and that uncertainties in baryonic effects can increase parameter error bars; however, the latter can be mitigated by including information about baryonic effects from kinetic and thermal Sunyaev-Zel'dovich measurements by CMB-HD. The CMB-HD likelihood and Fisher estimation codes used here are publicly available; the likelihood is integrated with Cobaya to facilitate parameter forecasting.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cosmological parameter estimation with a joint-likelihood analysis of the cosmic microwave background and big bang nucleosynthesis

Here, we present a joint-likelihood analysis of big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) data, consistently combining likelihoods and taking into account uncertainties in nuclear reaction rates for the first time. Bayesian inference is performed on the baryon abundance and the effective number of neutrino species, 𝑁 eff , using a CMB Boltzmann solver in combination with LINX , a new flexible and efficient BBN code. We marginalize over Planck nuisance parameters and nuclear rates to find 𝑁 eff =3.0⁢8$^{+0.15}_{−0.14}$, 2.9⁢4$^{+0.16}_{−0.15}$, or 2.96$^{+0.13}_{−0.14}$, for three separate reaction networks. This framework enables robust testing of the lambda cold dark matter paradigm and its variants with CMB and BBN data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Use of Thermoregulatory Models to Enhance Space Shuttle and Space Station operations and Review of Human Thermoregulatory Control

Thermoregulation in the space environment is critical for survival, especially in off- nominal operations. In such cases, mathematical models of thermoregulation are frequently employed to evaluate safety-of-flight issues in various human mission scenarious. In this study, the 225-node Wissler model and the 41-Node Metabolic Man model are employed to evaluate the effects of such a scenario. Metabolic loads on astronauts wearing the advanced crew escape suit (ACES) and liquid cooled ventilation garment (LCVG) are imposed on astronauts exposed to elevated cabin temperatures resulting from a systems failure. The study indicates that the performance of the ACES/LCVG cooling system is marginal. Increases in workload and or cabin temperature above nominal will increase rectal temperature, stored heat load, heart rate, and sweating, which could lead to deficits in the performance of cognitive and motor tasks. This is of concern as the ACES/LCVG is employed during Shuttle decent when the likelihood of a safe landing may be compromised. The study indicates that the most effective mitigation strategy would be to decrease the LCVG inlet temperature.

Pisacane, V. L.↗

AutoClass: A Bayesian Approach to Classification

We describe a Bayesian approach to the untutored discovery of classes in a set of cases, sometimes called finite mixture separation or clustering. The main difference between clustering and our approach is that we search for the "best" set of class descriptions rather than grouping the cases themselves. We describe our classes in terms of a probability distribution or density function, and the locally maximal posterior probability valued function parameters. We rate our classifications with an approximate joint probability of the data and functional form, marginalizing over the parameters. Approximation is necessitated by the computational complexity of the joint probability. Thus, we marginalize w.r.t. local maxima in the parameter space. We discuss the rationale behind our approach to classification. We give the mathematical development for the basic mixture model and describe the approximations needed for computational tractability. We instantiate the basic model with the discrete Dirichlet distribution and multivariant Gaussian density likelihoods. Then we show some results for both constructed and actual data.

Stutz, John↗

Baseband processor development for the Advanced Communications Satellite Program

An onboard-baseband-processor concept for a satellite-switched time-division-multiple-access (SS-TDMA) communication system was developed for NASA Lewis Research Center. The baseband processor routes and controls traffic on an individual message basis while providing significant advantages in improved link margins and system flexibility. Key technology developments required to prove the flight readiness of the baseband-processor design are being verified in a baseband-processor proof-of-concept model. These technology developments include serial MSK modems, Clos-type baseband routing switch, a single-chip CMOS maximum-likelihood convolutional decoder, and custom LSL implementation of high-speed, low-power ECL building blocks.

Moat, D.↗

High dynamic GPS receiver using maximum likelihood estimation and frequency tracking

A new high dynamic global positioning system (GPS) receiver is presented and its performance characterized by analysis, simulation, and demonstration. The demonstration receiver is a breadboard model capable of tracking a single simulated satellite signal in pseudorange and range rate. Pseudorange and range rate estimates are made once every 20 ms, using a maximum likelihood estimator, and are tracked by means of a third-order fading memory filter in a feedback configuration. The receiver tracks pseudorange with rms errors of under 1 m when subjected to simulated 50 g, 40 g/s circular trajectories. The tracking threshold is approximately 28 dB-Hz, which provides 12 dB margin relative to the minimum specified signal strength, assuming 3.5 dB system noise figure and 0 dBi antenna gain.

Hurd, W. J.↗

Dark Energy Survey Year 6 results: Clustering redshifts and importance sampling of self-organized-maps 𝑛⁡(𝑧) realizations for 3 × 2 ⁢pt samples

This work is part of a series establishing the redshift framework for the 3 × 2 ⁢pt analysis of the Dark Energy Survey Year 6 (DES Y6). For DES Y6, photometric redshift distributions are estimated using self-organizing maps (SOMs), calibrated with spectroscopic and many-band photometric data. To overcome limitations from color-redshift degeneracies and incomplete spectroscopic coverage, we enhance this approach by incorporating clustering-based redshift constraints (clustering-z, or WZ) from angular cross-correlations with BOSS and eBOSS galaxies and eBOSS quasar samples. We define a WZ likelihood and apply importance sampling to a large ensemble of SOM-derived 𝑛⁡(𝑧) realizations, selecting those consistent with the clustering measurements to produce a posterior sample for each lens and source bin. The analysis uses angular scales corresponding to 1.5–5 Mpc to optimize signal-to-noise ratio while mitigating modeling uncertainties and marginalizes over redshift-dependent galaxy bias and other systematics informed by the N-body simulation CARDINAL . While a sparser spectroscopic reference sample limits WZ constraining power at 𝑧 >1.1, particularly for source bins, we demonstrate that combining SOM with WZ improves redshift accuracy and enhances the overall cosmological constraining power of DES Y6. As a result, we estimate an improvement in 𝑆 8 of approximately 10% for cosmic shear and 3 ×2⁢pt analysis, primarily due to the WZ calibration of the source samples.

Cosmological parameters↗

Compressing the Cosmological Information in One-dimensional Correlations of the Lyman-$\alpha$ Forest

Observations of the Lyman-α forest from spectroscopic surveys such as the Baryon Oscillation Spectroscopic Survey or its extension, eBOSS, or the ongoing Dark Energy Spectroscopic Instrument (DESI) survey offer a unique window to study the growth of structure on megaparsec scales. Interpretation of these measurements is a complicated task, requiring hydrodynamical simulations to model and marginalize over the thermal and ionization state of the intergalactic medium. This complexity has limited the use of Lyα clustering measurements in joint cosmological analyses. In this work we show that the cosmological information content of the one-dimensional power spectrum (P 1D ) of the Lyman-α forest can be compressed into a simple two-parameter likelihood without any significant loss of constraining power. We simulate P 1D measurements from DESI using hydrodynamical simulations and show that the compressed likelihood is model independent and lossless, recovering unbiased results even in the presence of massive neutrinos or running of the primordial power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

Ridge-trench collision in Archean and Post-Archean crustal growth: Evidence from southern Chile

The growth of continental crust at convergent plate margins involves both continuous and episodic processes. Ridge-trench collision is one episodic process that can cause significant magmatic and tectonic effects on convergent plate margins. Because the sites of ridge collision (ridge-trench triple junctions) generally migrate along convergent plate boundaries, the effects of ridge collision will be highly diachronous in Andean-type orogenic belts and may not be adequately recognized in the geologic record. The Chile margin triple junction (CMTJ, 46 deg S), where the actively spreading Chile rise is colliding with the sediment-filled Peru-Chile trench, is geometrically and kinematically the simplest modern example of ridge collision. The south Chile margin illustrates the importance of the ridge-collision tectonic setting in crustal evolution at convergent margins. Similarities between ridge-collision features in southern Chile and features of Archean greenstone belts raise the question of the importance of ridge collision in Archean crustal growth. Archean plate tectonic processes were probably different than today; these differences may have affected the nature and importance of ridge collision during Archean crustal growth. In conclusion, it is suggested that smaller plates, greater ridge length, and/or faster spreading all point to the likelihood that ridge collision played a greater role in crustal growth and development of the greenstone-granite terranes during the Archean. However, the effects of modern ridge collision, and the processes involved, are not well enough known to develop specific models for the Archean ridge collison.

Nelson, E. P.↗

Performance and operational considerations in the design of vehicle antennas for mobile satellite communications

This paper examines the vehicle antenna requirements for mobile satellite systems. The antenna parameters are discussed in the light of the requirements and the limitations in performance imposed by the physical constraints of antenna and by vehicle geometries. Measurements of diffraction and antenna noise temperature in an operational environment are examined, as well as their effects on system margins. Mechanical versus electronic designs are compared with regards to performance, cost, reliability, and design complexity. Comparisons between open-loop and close-loop tracking systems are made and the effects of bandwidth, sidelobe levels, operational constraints, vehicle angular velocity, and acceleration are discussed. Some consideration is given to the use of hybrid systems employing both open and closed-loop tracking. Changes to antenna/terminal specifications are recommended which will provide greater design flexibility and increase the likelihood of meeting the performance and operational requirements.

Milne, R.↗

Planck constraints on the tensor-to-scalar ratio

In this work, we present constraints on the tensor-to-scalar ratio r using Planck data. We use the latest release of Planck maps (PR4), processed with the NPIPE code, which produces calibrated frequency maps in temperature and polarization for all Planck channels from 30 GHz to 857 GHz using the same pipeline. We computed constraints on r using the BB angular power spectrum, and we also discuss constraints coming from the TT spectrum. Given Planck's noise level, the TT spectrum gives constraints on r that are cosmic-variance limited (with $\sigma$ r = 0.093), but we show that the marginalized posterior peaks towards negative values of r at about the 1.2$\sigma$ level. We derived Planck constraints using the BB power spectrum at both large angular scales (the 'reionization bump') and intermediate angular scales (the 'recombination bump') from $\ell$=2 to 150, and find a stronger constraint than that from TT, with $\sigma$ r = 0.069. The Planck BB spectrum shows no systematic bias, and is compatible with zero, given both the statistical noise and the systematic uncertainties. The likelihood analysis using B modes yields the constraint r < 0.158 at 95% confidence using more than 50% of the sky. This upper limit tightens to r < 0.069 when Planck EE, BB, and EB power spectra are combined consistently, and it tightens further to r < 0.056 when the Planck TT power spectrum is included in the combination. Finally, combining Planck with BICEP2/Keck 2015 data yields an upper limit of r < 0.044.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmic background radiation anisotropy at degree angular scales - Further results from the South Pole

We report further results from the University of California at Santa Barbara program to measure anisotropy in the cosmic background radiation at angular scales near 1 deg, an angular range corresponding to the largest scales where structure is observed. A 30 GHz high electron mobility transistor amplifier-based detector was coupled to the Advanced Cosmic Microwave Explorer, a 1 m off-axis Gregorian telescope. We present data that represent 64 of the total of 500 hr acquired with this system during the 1990-1991 season. The data have a statistical error of 13.5/micro-K/pixel. These are the smallest error bars of any data set of this type published to date. The data contain a significant signal with a maximum likelihood Delta T/T roughly 1 x 10 exp -5. The spectrum of the signal seen in slightly less than 2 sigma away from the thermal spectrum expected of primordial fluctuations in the cosmic background radiation. If the source of the fluctuations is primordial, then the data are consistent with cold dark matter scenarios when normalized to the large-scale anisotropy observed by COBE, while if the origin of the signal is foreground emission or another form of contaminant then the data are marginally inconsistent with standard cold dark matter models. In either case, the data are sufficiently sensitive to provide a crucial test of many models.

Schuster, Jeffrey↗

Parameter inference from event ensembles and the top-quark mass

One of the key tasks of any particle collider is measurement. In practice, this is often done by fitting data to a simulation, which depends on many parameters. Sometimes, when the effects of varying different parameters are highly correlated, a large ensemble of data may be needed to resolve parameter-space degeneracies. An important example is measuring the top-quark mass, where other physical and unphysical parameters in the simulation must be profiled when fitting the top-quark mass parameter. We compare four different methodologies for top-quark mass measurement: a classical histogram fit similar to one commonly used in experiment augmented by soft-drop jet grooming; a 2D profile likelihood fit with a nuisance parameter; a machine-learning method called DCTR; and a linear regression approach, either using a least-squares fit or with a dense linearly-activated neural network. Despite the fact that individual events are totally uncorrelated, we find that the linear regression methods work most effectively when we input an ensemble of events sorted by mass, rather than training them on individual events. Although all methods provide robust extraction of the top-quark mass parameter, the linear network does marginally best and is remarkably simple. For the top study, we conclude that the Monte-Carlo-based uncertainty on current extractions of the top-quark mass from LHC data can be reduced significantly (by perhaps a factor of 2) using networks trained on sorted event ensembles. More generally, machine learning from ensembles for parameter estimation has broad potential for collider physics measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bayesian quantification of observability and equation of state of twin stars

The possibility of discovering twin stars, two neutron stars (NSs) with the same mass but different radii, is usually studied in forward modelings by using a restricted number of NS matter equations of state (EOSs) encapsulating a first-order phase transition from hadronic to quark matter (QM). Informing our likelihood function with the NS radius data from GW170817 and using a metamodel with nine parameters capable of mimicking most NS EOSs available in the literature, we conduct a Bayesian quantification of the observability and underlying EOSs of twin stars. Of the accepted EOSs, between 12 and 18% yield twin stars, depending on the restrictions we place on the second branch. The possibility of twin stars remains robust even under recent observational constraints. Here, we show that many of these twin star scenarios are observable with currently available levels of accuracy in measuring NS radii. We also present the marginalized posterior probability density functions (PDFs) of every EOS parameter for each of four mass-radius correlation topologies. We find that the inferred EOS depends sensitively on not only whether twin stars are present, but also the category of twin stars, indicating that the observation of twin stars would provide a strong constraint on the underlying EOS. In particular, for two coexisting hybrid stars having QM cores at different densities, the PDF for QM speed of sound squared 𝑐$^2_{qm}$ has two peaks, one below and another above the conformal limit 𝑐$^2_{qm}$ = 1/3 predicted by perturbative QCD.

QCD in nuclear reactions↗

BeyondPlanck X. Bandpass and beam leakage corrections

We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regularization. The main net impact of these modifications is an overall shift in the 70 GHz bandpass of +0.6 GHz; we argue that any analysis of LFI data products, either from Planck or BeyondPlanck, should use these new bandpasses. In addition, we fit a single free bandpass parameter for each radiometer of the form Δ i = Δ 0 + $δ$ i , where Δ 0 represents an absolute frequency shift per frequency band and $δ$i is a relative shift per detector. The absolute correction is only fitted at 30 GHz with a full $\chi$ 2 -based likelihood, resulting in a correction of Δ 30 = 0.24 ± 0.03 GHz. The relative corrections are fitted using a spurious map approach, fundamentally similar to the method pioneered by the WMAP team, but without introducing many additional degrees of freedom. All bandpass parameters are sampled using a standard Metropolis sampler within the main BeyondPlanck Gibbs chain, and bandpass uncertainties are thus propagated to all other data products in the analysis. In total, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopt the official Planck LFI beam estimates without additional degrees of freedom, and only marginalize over the underlying sky model. We note that this is the first-time leakage from beam mismatch has been included for Planck LFI maps.

79 ASTRONOMY AND ASTROPHYSICS↗

GPS high dynamic receiver tracking demonstration results

Demonstration results are presented for a high dynamic GPS receiver. The receiver tested is a breadboard unit capable of tracking one simulated satellite signal in pseudorange and range rate. The receiver makes approximate maximum likelihood estimates of pseudorange and range rate each 20 ms, and tracks these observables using a third order filter with a time constant of 0.14 s. Carrier phase is not tracked, which eliminates the typical failure mode of loss of carrier lock associated with PLLs at high dynamics. The receiver tracks with pseudorange lag errors of under 0.06 m when subjected to simulated 50 g turns with 40 g/s peak jerk. Pseudorange errors due to receiver noise alone are approximately 0.6 m rms at a carrier power to noise spectral density ratio of 34 dB-Hz. The tracking threshold SNR is approximately 28 dB-Hz, which provides 12 dB margin relative to the 40 dB-Hz that occurs with minimum specified satellite signal strength, 3.5 dB system noise figure, and 0 dBi antenna gain.

Hurd, W. J.↗

Simulation Based Inference with Domain Adaptation for Strong Gravitational Lensing

Simulation based inference leverages machine learning to carry out Bayesian inference for systems with intractable likelihoods. However, transitioning a network trained on simulated data to real data runs the risk of encountering domain shift, leading to performance losses. We attempt to implement domain adaptation into the sbi neural posterior estimation framework using the Maximum Mean Discrepancy as an additional network loss. We test two network architectures and use masked autoregressive flow for density estimation. We test the network on a set of 400,000 simulated strong gravitational lensing images generated using deeplenstronomy. The source domain is defined as low noise whereas the target domain has a noise profile sampled from experimentally derived DES survey conditions. We find that, while DA does lead to performance improvements, they are marginal at ~6% less inference error. We also find a similar marginal improvement in uncertainty calibration at around 8%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

CMB-S4: Foreground-cleaning Pipeline Comparison for Measuring Primordial Gravitational Waves

We compare multiple foreground-cleaning pipelines for estimating the tensor-to-scalar ratio, r, using simulated maps of the planned CMB-S4 experiment within the context of the South Pole Deep Patch. To evaluate robustness, we analyze bias and uncertainty on r across various foreground suites using map-based simulations. The foreground-cleaning methods include: a parametric maximum likelihood approach applied to auto- and cross-power spectra between frequency maps; a map-based parametric maximum-likelihood method; and a harmonic-space internal linear combination using frequency maps. We summarize the conceptual basis of each method to highlight their similarities and differences. To better probe the impact of foreground residuals, we implement an iterative internal delensing step, leveraging a map-based pipeline to generate a lensing B-mode template from the large aperture telescope frequency maps. Our results show that the performance of the three approaches is comparable for simple and intermediate-complexity foregrounds, with σ(r) ranging from 3–5 ×10 −4 . However, biases at the 1σ–2σ level appear when analyzing more complex forms of foreground emission. By extending the baseline pipelines to marginalize over foreground residuals, we demonstrate that contamination can be reduced to within statistical uncertainties, albeit with a pipeline-dependent impact on σ(r), which translates to a detection significance between 2σ and 4σ for an input value of r = 0.003. These findings suggest varying levels of maturity among the tested pipelines, with the auto- and cross-spectra-based approach demonstrating the best stability and overall performance. Moreover, given the extremely low noise levels, mutual validation of independent foreground-cleaning pipelines is essential to ensure the robustness of any potential detection.

astronomy data analysis↗