SEARCH · Search NASA

Results for “Bayesian statistical modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Peak Ground Displacement Saturates Exactly When Expected: Implications for Earthquake Early Warning

The scaling of rupture properties with magnitude is of critical importance to earthquake early warning (EEW) systems that rely on source characterization using limited snapshots of waveform data. ShakeAlert, a prototype EEW system that is being developed for the western United States, provides real-time estimates of earthquake magnitude based on P-wave peak ground displacements calculated at stations triggered by the event. The algorithms used in ShakeAlert assume that the displacement measurements at each station are statistically independent and that there exists a linear and time-independent relation between log peak ground displacement and earthquake magnitude. Here we challenge this basic assumption using the largest dataset assembled for this purpose to date: a comprehensive database of more than 140,000 vertical-component waveforms from M4.5–M9 earthquakes occurring near Japan from 1997 through 2018 and recorded by the K-NET and KiK-net strong-motion networks. By analyzing the time-evolution of P-wave peak ground displacements for these earthquakes, we show that there is a break, or saturation, in the magnitude-displacement scaling that depends on the length of the measurement time window. We demonstrate that the magnitude at which this saturation occurs is well-explained by a simple and non-deterministic model of earthquake rupture growth. We then utilize the predictions of this saturation model to develop a Bayesian framework for estimating posterior uncertainties in real-time magnitude estimates.

58 GEOSCIENCES↗

DOE OSTI2014

The diskmass survey. VIII. On the relationship between disk stability and star formation

We study the relationship between the stability level of late-type galaxy disks and their star-formation activity using integral-field gaseous and stellar kinematic data. Specifically, we compare the two-component (gas+stars) stability parameter from Romeo and Wiegert (Q {sub RW}), incorporating stellar kinematic data for the first time, and the star-formation rate estimated from 21 cm continuum emission. We determine the stability level of each disk probabilistically using a Bayesian analysis of our data and a simple dynamical model. Our method incorporates the shape of the stellar velocity ellipsoid (SVE) and yields robust SVE measurements for over 90% of our sample. Averaging over this subsample, we find a meridional shape of σ{sub z}/σ{sub R}=0.51{sub −0.25}{sup +0.36} for the SVE and, at 1.5 disk scale lengths, a stability parameter of Q {sub RW} = 2.0 ± 0.9. We also find that the disk-averaged star-formation-rate surface density ( Σ-dot {sub e,∗}) is correlated with the disk-averaged gas and stellar mass surface densities (Σ {sub e,} {sub g} and Σ {sub e,} {sub *}) and anti-correlated with Q {sub RW}. We show that an anti-correlation between Σ-dot {sub e,∗} and Q {sub RW} can be predicted using empirical scaling relations, such that this outcome is consistent with well-established statistical properties of star-forming galaxies. Interestingly, Σ-dot {sub e,∗} is not correlated with the gas-only or star-only Toomre parameters, demonstrating the merit of calculating a multi-component stability parameter when comparing to star-formation activity. Finally, our results are consistent with the Ostriker et al. model of self-regulated star-formation, which predicts Σ-dot {sub e,∗}/Σ{sub e,g}∝Σ{sub e,∗}{sup 1/2}. Based on this and other theoretical expectations, we discuss the possibility of a physical link between disk stability level and star-formation rate in light of our empirical results.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Stopping criteria for ending autonomous, single detector radiological source searches

While the localization of radiological sources has traditionally been handled with statistical algorithms, such a task can be augmented with advanced machine learning methodologies. The combination of deep and reinforcement learning has provided learning-based navigation to autonomous, single-detector, mobile systems. However, these approaches lacked the capacity to terminate a surveying/search task without outside influence of an operator or perfect knowledge of source location (defeating the purpose of such a system). Two stopping criteria are investigated in this work for a machine learning navigated system: one based upon Bayesian and maximum likelihood estimation (MLE) strategies commonly used in source localization, and a second providing the navigational machine learning network with a “stop search” action. A convolutional neural network was trained via reinforcement learning in a 10 m × 10 m simulated environment to navigate a randomly placed detector-agent to a randomly placed source of varied strength (stopping with perfect knowledge during training). The network agent could move in one of four directions (up, down, left, right) after taking a 1 s count measurement at the current location. During testing, the stopping criteria for this navigational algorithm was based upon a Bayesian likelihood estimation technique of source presence, updating this likelihood after each step, and terminating once the confidence of the source being in a single location exceeded 0.9. A second network was trained and tested with similar architecture as the previous but which contained a fifth action: for self-stopping. The accuracy and speed of localization with set detector and source initializations were compared over 50 trials of MLE-Bayesian approach and 1000 trials of the CNN with self-stopping. The statistical stopping condition yielded a median localization error of ~1.41 m and median localization speed of 12 steps. The machine learning stopping condition yielded a median localization error of 0 m and median localization speed of 17 steps. This work demonstrated two stopping criteria available to a machine learning guided, source localization system.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2015

Constraining Data Mining with Physical Models: Voltage- and Oxygen Pressure-Dependent Transport in Multiferroic Nanostructures

Development of new generation electronic devices requires understanding and controlling the electronic transport in ferroic, magnetic, and optical materials, which is hampered by two factors. First, the complications of working at the nanoscale, where interfaces, grain boundaries, defects, and so forth, dictate the macroscopic characteristics. Second, the convolution of the response signals stemming from the fact that several physical processes may be activated simultaneously. Here, we present a method of solving these challenges via a combination of atomic force microscopy and data mining analysis techniques. Rational selection of the latter allows application of physical constraints and enables direct interpretation of the statistically significant behaviors in the framework of the chosen physical model, thus distilling physical meaning out of raw data. We demonstrate our approach with an example of deconvolution of complex transport behavior in a bismuth ferrite–cobalt ferrite nanocomposite in ambient and ultrahigh vacuum environments. Measured signal is apportioned into four electronic transport patterns, showing different dependence on partial oxygen and water vapor pressure. These patterns are described in terms of Ohmic conductance and Schottky emission models in the light of surface electrochemistry. Finally and furthermore, deep data analysis allows extraction of local dopant concentrations and barrier heights empowering our understanding of the underlying dynamic mechanisms of resistive switching.

36 MATERIALS SCIENCE↗

DOE OSTI2011

Search for the Standard Model Higgs boson produced in association with a W Boson in the isolated-track charged-lepton channel using the Collider Detector at Fermilab

The Higgs boson is the only elementary particle predicted by the Standard Model (SM) that has not yet been observed experimentally. If it exists, it explains the spontaneous electroweak symmetry breaking and the origin of mass for gauge bosons and fermions. We test the validity of the SM by performing a search for the associated production of a Higgs boson and a W boson in the channel where the Higgs boson decays to a bottom-antibottom quark pair and the W boson decays to a charged lepton and a neutrino (the WH channel). We study a dataset of proton-antiproton collisions at a centre-of-mass energy √s = 1.96 TeV provided by the Tevatron accelerator, corresponding to an integrated luminosity of 5.7 fb -1 , and recorded using the Collider Detector at Fermilab (CDF).We select events consistent with the signature of exactly one charged lepton (electron or muon), missing transverse energy due to the undetected neutrino (MET) and two collimated streams of particles (jets), at least one of which is required to be identified as originating from a bottom quark. We improve the discrimination of Higgs signal from backgrounds through the use of an artificial neural network. Using a Bayesian statistical inference approach, we set for each hypothetical Higgs boson mass in the range 100-150 GeV/c 2 with 5 GeV/c 2 increments a 95% credibility level (CL) upper limit on the ratio between the Higgs production cross section times branching fraction and the SM prediction. Our main original contributions are the addition of a novel charged lepton reconstruction algorithm with looser requirements (ISOTRK) with respect the electron or muon tight criteria (TIGHT), as well as the introduction of a novel trigger-combination method that allows to maximize the event yield while avoiding trigger correlations and that is used for the ISOTRK category. The ISOTRK candidate is a high-transverse-momentum good-quality track isolated from other activity in the tracking system and not required to match a calorimeter cluster, as for a tight electron candidate, or an energy deposit in the muon detector, as for a tight muon candidate. The ISOTRK category recovers real charged leptons that otherwise would be lost in the non-instrumented regions of the detector. This allows the reconstruction of more W boson candidates, which in turn increases the number of reconstructed WH signal candidate events, and therefore improves the sensitivity of the WH search. For the TIGHT charged lepton categories, we employ charged-lepton-dedicated triggers to improve the rate of WH signal acceptance during data taking. Since there is no ISOTRK-dedicated trigger at CDF, for the ISOTRK charged lepton category we employ three MET-plus-jets-based triggers. For each trigger we first identify the jet selection where the trigger efficiency is flat with respect to jet information (transverse energy and direction of motion in the transverse plane for the two jets in the event) and then we parametrize the trigger efficiency as a function of trigger MET. On an event-by-event basis, for each trigger we compute a trigger efficiency as a function of trigger parametrization, trigger MET, jet information, trigger prescale and information about whether the trigger is defined or not. For the ISOTRK category we combine the three triggers using a novel method, which allows the combination of any number of triggers in order to maximize the event yield while avoiding trigger correlations. On an event-by-event basis, only the trigger with the largest efficiency is used. By avoiding a logical 'OR' between triggers, the loss in the yield of events accepted by the trigger combination is compensated by a smaller and easier-to-compute corresponding systematic uncertainty. The addition of the ISOTRK charged lepton category to the TIGHT category produces an increase of 33% in the WH signal yield and a decrease of 15.5% to 19.0% in the median expected 95% CL cross-section upper limits across the entire studied Higgs mass interval. The improvement in analysis sensitivity is smaller than the improvement in signal yield because the ISOTRK category has a smaller signal over background ratio than the TIGHT category, due to the looser ISOTRK reconstruction criteria. The observed (median expected) 95% CL SM Higgs upper limits on cross section times branching ratio vary between 2.39 x SM (2.73 x SM) for a Higgs mass of 100 GeV/c 2 to 31.1 x SM (31.2 x SM) for a Higgs mass of 150 GeV/c 2 , while the value for a 115 GeV/c 2 Higgs boson is that of 5.08 x SM (3.79 x SM). The novel trigger combination method is already in use by several CDF analyses. It is applicable to any analysis that uses triggers based on MET and jets, such as supersymmetry searches at the ATLAS and CMS experiments at the Large Hadron Collider. In its most general form, the method can be used by any analysis that combines any number of different triggers.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Constraining the Milky Way Mass Profile with Phase-space Distribution of Satellite Galaxies

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M={1.23}_{-0.18}^{+0.21}\times {10}^{12}{M}_{\odot }$ and a concentration of $c={9.4}_{-2.1}^{+2.8}$ with the prior based on the M–c relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of the kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. Finally, by varying the disruption rate of such satellites in this simulation, we estimate a ~5% scatter in the inferred MW halo mass among hydrodynamics-based simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Bayesian estimation of Karhunen–Loève expansions; A random subspace approach

One of the most widely-used statistical procedures for dimensionality reduction of high dimensional random fields is Principal Component Analysis (PCA), which is based on the Karhunen-Lo eve expansion (KLE) of a stochastic process with finite variance. The KLE is analogous to a Fourier series expansion for a random process, where the goal is to find an orthogonal transformation for the data such that the projection of the data onto this orthogonal subspace is optimal in the L 2 sense, i.e, which minimizes the mean square error. In practice, this orthogonal transformation is determined by performing an SVD (Singular Value Decomposition) on the sample covariance matrix or on the data matrix itself. Sampling error is typically ignored when quantifying the principal components, or, equivalently, basis functions of the KLE. Furthermore, it is exacerbated when the sample size is much smaller than the dimension of the random field. In this paper, we introduce a Bayesian KLE procedure, allowing one to obtain a probabilistic model on the principal components, which can account for inaccuracies due to limited sample size. The probabilistic model is built via Bayesian inference, from which the posterior becomes the matrix Bingham density over the space of orthonormal matrices. We use a modified Gibbs sampling procedure to sample on this space and then build a probabilistic Karhunen-Lo eve expansions over random subspaces to obtain a set of low-dimensional surrogates of the stochastic process. We illustrate this probabilistic procedure with a finite dimensional stochastic process inspired by Brownian motion.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Bayesian-based response expansion and uncertainty quantification using sparse measurement sets

Systems subjected to dynamic loads often require monitoring of their vibrational response, but limitations on the total number and placement of the measurement sensors can hinder the data-collection process. Here, we present an indirect approach to estimate a system’s full-field dynamic response, including all uninstrumented locations, using response measurements from sensors sparsely located on the system. This approach relies on Bayesian inference that utilizes a system model to estimate the full-field response and quantify the uncertainty in these estimates. By casting the estimation problem in the frequency domain, this approach utilizes the modal frequency response functions as a natural, frequency-dependent weighting scheme for the system mode shapes to perform the expansion. This frequency-dependent weighting scheme enables an accurate expansion, even with highly correlated mode shapes that may arise from spatial aliasing due to the limited number of sensors, provided these correlated modes do not have natural frequencies that are closely spaced. Furthermore, the inherent regularization mechanism that arises in this Bayesian-based procedure enables the utilization of the full set of system mode shapes for the expansion, rather than any reduced subset. This approach can produce estimates when considering a single realization of the measured responses, and with some modification, it can also produce estimates for power spectral density matrices measured from many realizations of the responses from statistically stationary random processes. A simply supported beam provides an initial numerical validation, and a cylindrical test article excited by acoustic loads in a reverberation chamber provides experimental validation.

42 ENGINEERING↗

DOE OSTI2015

Joint Planck and WMAP assessment of low CMB multipoles

The remarkable progress in cosmic microwave background (CMB) studies over past decade has led to the era of precision cosmology in striking agreement with the ΛCDM model. However, the lack of power in the CMB temperature anisotropies at large angular scales (low-ℓ), as has been confirmed by the recent Planck data also (up to 0ℓ=4), although statistically not very strong (less than 3σ), is still an open problem. One can avoid to seek an explanation for this problem by attributing the lack of power to cosmic variance or can look for explanations i.e., different inflationary potentials or initial conditions for inflation to begin with, non-trivial topology, ISW effect etc. Features in the primordial power spectrum (PPS) motivated by the early universe physics has been the most common solution to address this problem. In the present work we also follow this approach and consider a set of PPS which have features and constrain the parameters of those using WMAP 9 year and Planck data employing Markov-Chain Monte Carlo (MCMC) analysis. The prominent feature of all the models of PPS that we consider is an infra-red cut off which leads to suppression of power at large angular scales. We consider models of PPS with maximum three extra parameters and use Akaike information criterion (AIC) and Bayesian information criterion (BIC) of model selection to compare the models. For most models, we find good constraints for the cut off scale k{sub c}, however, for other parameters our constraints are not that good. We find that sharp cut off model gives best likelihood value for the WMAP 9 year data, but is as good as power law model according to AIC. For the joint WMAP 9 + Planck data set, Starobinsky model is slightly preferred by AIC which is also able to produce CMB power suppression up to 0ℓ≤3 to some extent. However, using BIC criteria, one finds model(s) with least number of parameters (power law model) are always preferred.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

Integrating biological knowledge into variable selection: an empirical Bayes approach with an application in cancer biology

Background: An important question in the analysis of biochemical data is that of identifying subsets of molecular variables that may jointly influence a biological response. Statistical variable selection methods have been widely used for this purpose. In many settings, it may be important to incorporate ancillary biological information concerning the variables of interest. Pathway and network maps are one example of a source of such information. However, although ancillary information is increasingly available, it is not always clear how it should be used nor how it should be weighted in relation to primary data. Results: We put forward an approach in which biological knowledge is incorporated using informative prior distributions over variable subsets, with prior information selected and weighted in an automated, objective manner using an empirical Bayes formulation. We employ continuous, linear models with interaction terms and exploit biochemically-motivated sparsity constraints to permit exact inference. We show an example of priors for path way and network-based information and illustrate our proposed method on both synthetic response data and by an application to cancer drug response data. Comparisons are also made to alternative Bayesian and frequentist penalised-likelihood methods for incorporating network-based information. Conclusions: The empirical Bayes method proposed here can aid prior elicitation for Bayesian variable selection studies and help to guard against mis-specification of priors. Empirical Bayes, together with the proposed pathway-based priors, results in an approach with a competitive variable selection performance. In addition, the overall procedure is fast, deterministic, and has very few user-set parameters, yet is capable of capturing interplay between molecular players. The approach presented is general and readily applicable in any setting with multiple sources of biological prior knowledge.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression. Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to previous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Exploring physics of ferroelectric domain walls via Bayesian analysis of atomically resolved STEM data

The physics of ferroelectric domain walls is explored using the Bayesian inference analysis of atomically resolved STEM data. We demonstrate that domain wall profile shapes are ultimately sensitive to the nature of the order parameter in the material, including the functional form of Ginzburg-Landau-Devonshire expansion, and numerical value of the corresponding parameters. The preexisting materials knowledge naturally folds in the Bayesian framework in the form of prior distributions, with the different order parameters forming competing (or hierarchical) models. Here, we explore the physics of the ferroelectric domain walls in BiFeO 3 using this method, and derive the posterior estimates of relevant parameters. More generally, this inference approach both allows learning materials physics from experimental data with associated uncertainty quantification, and establishing guidelines for instrumental development answering questions on what resolution and information limits are necessary for reliable observation of specific physical mechanisms of interest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

Bayesian Design of Experiments for Logistic Regression to Evaluate Multiple Forensic Algorithms

When evaluating the performance of several forensic classification algorithms, it is desirable to construct a design that considers a variety of performance levels for each of the algorithms. We describe a strategy to use Bayesian design of experiments with multiple prior estimates to capture anticipated performance. Our goal is to characterize results from the different classification algorithms as a function of multiple explanatory variables and use this to choose a design about which units to test. Bayesian design of experiments has been successful for generalized linear models, including logistic regression models. Here, we develop methodology for the case where there are several potentially nonoverlapping priors for anticipated performance under consideration. The weighted priors method performs well for a broad range of true underlying model parameter choices and is more robust when compared to other candidate design choices. Additionally, we show how this can be applied in the multivariate input case and provide some useful summary measures. The shared information plot is used to evaluate design point allocation, and the D-value difference plot allows for the comparison of design performance across multiple potential parameter values in higher dimensions. Here, we illustrate the methods with several examples.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Statistical sensitivity on right-handed currents in presence of eV scale sterile neutrinos with KATRIN

The KATRIN experiment aims to determine the absolute neutrino mass by measuring the endpoint region of the tritium β-spectrum. As a large-scale experiment with a sharp energy resolution, high source luminosity and low background it may also be capable of testing certain theories of neutrino interactions beyond the standard model (SM). An example of a non-SM interaction are right-handed currents mediated by right-handed W bosons in the left-right symmetric model (LRSM). In this extension of the SM, an additional SU(2){sub R} symmetry in the high-energy limit is introduced, which naturally includes sterile neutrinos and predicts the seesaw mechanism. In tritium β decay, this leads to an additional term from interference between left- and right-handed interactions, which enhances or suppresses certain regions near the endpoint of the beta spectrum. In this work, the sensitivity of KATRIN to right-handed currents is estimated for the scenario of a light sterile neutrino with a mass of some eV. This analysis has been performed with a Bayesian analysis using Markov Chain Monte Carlo (MCMC). The simulations show that, in principle, KATRIN will be able to set sterile neutrino mass-dependent limits on the interference strength. The sensitivity is significantly increased if the Q value of the β decay can be sufficiently constrained. However, the sensitivity is not high enough to improve current upper limits from right-handed W boson searches at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

A Bayesian framework for cosmic string searches in CMB maps

There exists various proposals to detect cosmic strings from Cosmic Microwave Background (CMB) or 21 cm temperature maps. Current proposals do not aim to find the location of strings on sky maps, all of these approaches can be thought of as a statistic on a sky map. We propose a Bayesian interpretation of cosmic string detection and within that framework, we derive a connection between estimates of cosmic string locations and cosmic string tension G μ. We use this Bayesian framework to develop a machine learning framework for detecting strings from sky maps and outline how to implement this framework with neural networks. The neural network we trained was able to detect and locate cosmic strings on noiseless CMB temperature map down to a string tension of G μ=5 ×10{sup −9} and when analyzing a CMB temperature map that does not contain strings, the neural network gives a 0.95 probability that G μ≤2.3×10{sup −9}.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2013

A HOT URANUS ORBITING THE SUPER METAL-RICH STAR HD 77338 AND THE METALLICITY-MASS CONNECTION

We announce the discovery of a low-mass planet orbiting the super metal-rich K0V star HD 77338 as part of our ongoing Calan-Hertfordshire Extrasolar Planet Search. The best-fit planet solution has an orbital period of 5.7361 {+-} 0.0015 days and with a radial velocity semi-amplitude of only 5.96 {+-} 1.74 ms{sup -1}, we find a minimum mass of 15.9{sup +4.7}{sub -5.3} M{sub Circled-Plus }. The best-fit eccentricity from this solution is 0.09{sup +0.25}{sub -0.09}, and we find agreement for this data set using a Bayesian analysis and a periodogram analysis. We measure a metallicity for the star of +0.35 {+-} 0.06 dex, whereas another recent work finds +0.47 {+-} 0.05 dex. Thus HD 77338b is one of the most metal-rich planet-host stars known and the most metal-rich star hosting a sub-Neptune-mass planet. We searched for a transit signature of HD 77338b but none was detected. We also highlight an emerging trend where metallicity and mass seem to correlate at very low masses, a discovery that would be in agreement with the core accretion model of planet formation. The trend appears to show that for Neptune-mass planets and below, higher masses are preferred when the host star is more metal-rich. Also a lower boundary is apparent in the super metal-rich regime where there are no very low mass planets yet discovered in comparison to the sub-solar metallicity regime. A Monte Carlo analysis shows that this low-mass planet desert is statistically significant with the current sample of 36 planets at the {approx}4.5{sigma} level. In addition, results from Kepler strengthen the claim for this paucity of the lowest-mass planets in super metal-rich systems. Finally, this discovery adds to the growing population of low-mass planets around low-mass and metal-rich stars and shows that very low mass planets can now be discovered with a relatively small number of data points using stable instrumentation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

A Statistical Framework for Operational lnfrasound Monitoring [Slides]

Historically, much infrasound research has been event-based, with known events providing the ground-truth required to inform the search for associated signals. For monitoring nuclear explosions in the atmosphere or underground, reliable techniques are required for operational infrasound data processing that effectively minimize both false alarms and missed detections. A difficulty in working with infrasound data is to adequately account for dynamic variations in the atmosphere, which occur over a variety of temporal and spatial scales. However, contrary to one's intuition, it is not necessary to have a high fidelity atmospheric model in order to perform robust detection, association and location. We present detection, association and location algorithms that assume no a priori knowledge of the atmosphere; however, where knowledge is available we show how it can be used to reduce the associated uncertainties. Rather than exporting seismic techniques to the problem of infrasound monitoring, our techniques were developed from the ground-up to properly account for the unique considerations of atmospheric acoustics. In contrast with other routinely used infrasound detectors, we utilize a contextual detection hypothesis that adaptively accounts for temporally variable correlated ambient noise. Our location scheme allows for the location of infrasound events without any meteorological data through a Bayesian framework. We demonstrate the application of the techniques presented here for a series of ground-truth events, and for routine data processing of large amounts of regional infrasound data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗