Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bayesian Analysis”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Do joint CMB and HST data support a scale invariant spectrum?

We combine current measurements of the local expansion rate, H {sub 0}, and Big Bang Nucleosynthesis (BBN) estimates of helium abundance with the latest cosmic microwave background (CMB) data from the Planck Collaboration to discuss the observational viability of the scale invariant Harrison-Zeldovch-Peebles (HZP) spectrum. We also analyze some of its extensions, namely, HZP + Y {sub P} and HZP + N {sub eff}, where Y {sub P} is the primordial helium mass fraction and N {sub eff} is the effective number of relativistic degrees of freedom. We perform a Bayesian analysis and show that the latter model is favored with respect to the standard cosmology for values of N {sub eff} lying in the interval 3.70 ± 0.13 (1σ), which is currently allowed by some independent analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Constraining the break of spatial diffeomorphism invariance with Planck data

The current most accepted paradigm for the early universe cosmology, the inflationary scenario, shows a good agreement with the recent Cosmic Microwave Background (CMB) and polarization data. However, when the inflation consistency relation is relaxed, these observational data exclude a larger range of red tensor tilt values, prevailing the blue ones which are not predicted by the minimal inflationary models. Recently, it has been shown that the assumption of spatial diffeomorphism invariance breaking (SDB) in the context of an effective field theory of inflation leads to interesting observational consequences. Among them, the possibility of generating a blue tensor spectrum, which can recover the specific consistency relation of the String Gas Cosmology, for a certain choice of parameters. We use the most recent CMB data to constrain the SDB model and test its observational viability through a Bayesian analysis assuming as reference an extended ΛCDM+tensor perturbation model, which considers a power-law tensor spectrum parametrized in terms of the tensor-to-scalar ratio, r , and the tensor spectral index, n {sub t} . If the inflation consistency relation is imposed, r =−8 n {sub t} , we obtain a strong evidence in favor of the reference model whereas if such relation is relaxed, a weak evidence in favor of the model with diffeomorphism breaking is found. We also use the same CMB data set to make an observational comparison between the SDB model, standard inflation and String Gas Cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

INSTANTANEOUS STARBURST OF THE MASSIVE CLUSTERS WESTERLUND 1 AND NGC 3603 YC

We present a new method to determine the age spread of resolved stellar populations in a starburst cluster. The method relies on a two-step process. In the first step, kinematic members of the cluster are identified based on multi-epoch astrometric monitoring. In the second step, a Bayesian analysis is carried out, comparing the observed photometric sequence of cluster members with sets of theoretical isochrones. When applying this methodology to optical and near-infrared high angular resolution Hubble Space Telescope (HST) and adaptive optics observations of the {approx}5 Myr old starburst cluster Westerlund 1 and {approx}2 Myr old starburst cluster NGC 3603 YC, we derive upper limits for the age spreads of 0.4 and 0.1 Myr, respectively. The results strongly suggest that star formation in these starburst clusters happened almost instantaneously.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE RISE AND FALL OF THE STAR FORMATION HISTORIES OF BLUE GALAXIES AT REDSHIFTS 0.2 < z < 1.4

Popular cosmological scenarios predict that galaxies form hierarchically from the merger of many progenitors, each with their own unique star formation history (SFH). We use a sophisticated approach to constrain the SFHs of 4517 blue (presumably star-forming) galaxies with spectroscopic redshifts in the range 0.2 < z < 1.4 from the All-Wavelength Extended Groth Strip International Survey. This consists in the Bayesian analysis of the observed galaxy spectral energy distributions with a comprehensive library of synthetic spectra assembled using realistic, hierarchical star formation, and chemical enrichment histories from cosmological simulations. We constrain the SFH of each galaxy in our sample by comparing the observed fluxes in the B, R, I, and K{sub s} bands and rest-frame optical emission-line luminosities with those of one million model spectral energy distributions. We explore the dependence of the resulting SFHs on galaxy stellar mass and redshift. We find that the average SFHs of high-mass galaxies rise and fall in a roughly symmetric bell-shaped manner, while those of low-mass galaxies rise progressively in time, consistent with the typically stronger activity of star formation in low-mass compared to high-mass galaxies. For galaxies of all masses, the star formation activity rises more rapidly at high than at low redshift. These findings imply that the standard approximation of exponentially declining SFHs widely used to interpret observed galaxy spectral energy distributions may not be appropriate to constrain the physical parameters of star-forming galaxies at intermediate redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Copacabana: a probabilistic membership assignment method for galaxy clusters

Cosmological analyses using galaxy clusters in optical/near-infrared photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for all galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89 per cent membership accuracy with a mild dependence on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass–cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities (⁠μ * ⁠) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation at low redshift for a Legacy Survey of Space and Time type survey. These results indicate the potential of Copacabana and μ * to be used in cosmological analyses of optically selected clusters in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Effect of recoils on soft-drop-groomed observables in 𝛾-tagged jets in a multistage approach

We investigate medium-induced modifications to jet substructure observables that characterize hard components in central Pb-Pb collisions at $\sqrt{𝑠_{𝑁⁢𝑁}}$ = 5.02 TeV. Using a multistage Monte Carlo simulation of in-medium jet shower evolution, we explore flavor-dependent medium effects through simulations of inclusive and 𝛾-tagged jets. The results show that quark jets undergo a nonmonotonic modification compared with gluon jets in observables such as the Pb-Pb to 𝑝−𝑝 ratio of the soft drop prong angle 𝑟 𝑔 , the relative prong transverse momentum 𝑘 𝑇,𝑔 , and the groomed mass 𝑚 𝑔 distributions. Due to this nonmonotonic modification, 𝛾-tagged jets, enriched in quark jets, provide surprisingly clear signals of medium-induced structural modifications, distinct from effects dominated by selection bias. Further systematic studies demonstrate that these effects are dominated by recoil medium response. This work highlights the potential of hard substructures in 𝛾-tagged jets as powerful tools for probing the jet-medium interactions in high-energy heavy-ion collisions, in particular by enabling detailed investigations of jet-medium parton scatterings via their associated medium response. All simulations for 𝛾-tagged jet analyses carried out in this paper used triggered events containing at least one hard photon, which highlights the utility of these observables for future Bayesian analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Dark matter clumps as sources of gravitational-wave glitches in LIGO-Virgo-KAGRA data

We consider the hypothetical possibility that nonstationary glitch features in the noise of ground-based gravitational-wave detectors could be produced by small dark matter clumps that pass through the earth in the vicinity of gravitational-wave detectors. Here, we first derive the gravitational-wave strain that would be generated by the passage of such a dark matter clump. We find that the strain is primarily sourced by the Newtonian gravitational acceleration of the mirrors toward the clump and by the Shapiro time delay of the photons in the laser beams as they pass through the gravitational potential created by the dark matter clump. We also find that the Newtonian acceleration effect dominates the gravitational-wave strain for both ground and space-based interferometers. We then compare our dark matter clump, gravitational-wave strain model to 84 Koi-Fish glitches detected during the second observing run of the LIGO/Virgo/KAGRA Collaboration through a Markov chain Monte Carlo Bayesian analysis. We find that all glitches but nine can be confidently rejected as having originated from dark matter clumps. For the remaining glitches, the dark matter hypothesis cannot be excluded, and the maximum a posteriori parameters yield minimum densities of about 10 −7 g⁡/cm 3 , within the model. These results allow us to place the first direct upper limits with gravitational-wave detectors on the local over-density of dark matter in the form of clumps in the local neighborhood of Earth, namely 𝜌 DM ⁢clumps ≲ 10 −15 g⁡/cm −3 .

Astronomy and AstroPhysics↗
DOE OSTI2011

Criteria for reliable entanglement quantification with finite data

We propose one and a half criteria for determining how many measurements are needed to quantify entanglement reliably. We base these criteria on Bayesian analysis of measurement results and apply our methods to four-qubit entanglement, but generalizations to more qubits are straightforward.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Landau-Ginzburg theory of charge density wave formation accompanying lattice and electronic long-range ordering

Here we propose an analytical Landau-Ginzburg (LG) theory of the charge density waves coupled with lattice and electronic long-range order parameters. Examples of long-range order include the electronic wave function of superconducting Cooper pairs, structural distortions, electric polarization, and magnetization. We formulate the LG free energy density as a power expansion with respect to the charge density and other long-range order parameters as well as their spatial gradients and biquadratic coupling terms. We introduced a biquadratic coupling between the charge density gradient and long-range order parameters as well as nonlinear higher gradients of the long-range order parameters. The biquadratic gradient coupling is critical to the appearance of different spatially modulated phases in charge-ordered ferroics and high-temperature superconductors. We derived the thermodynamic conditions for the stability of the spatially modulated phases, which are the intertwined spatial waves of charge density and lattice/electronic long-range order. The analytical expressions for the energies of different phases, corresponding order parameters, charge density waves amplitudes, and modulation periods obtained in this paper can be employed to guide the comprehensive physical explanation, deconvolution, and Bayesian analysis of experimental data on quantum materials ranging from charge-ordered ferroics to high-temperature superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Bayesian inference of nuclear symmetry energy from measured and imagined neutron skin thickness in Sn 116 , 118 , 120 , 122 , 124 , 130 , 132 , Pb 208 , and Ca 48

The neutron skin thickness Δr np in heavy nuclei has been known as one of the most sensitive terrestrial probes of the nuclear symmetry energy E sym (ρ) around $\frac{2}{3}$ of the saturation density ρ 0 of nuclear matter. Existing neutron skin data mostly from hadronic observables suffer from large uncertainties and their extraction from experiments are often strongly model dependent. While waiting eagerly for the promised model-independent and high-precision neutron skin data for 208 Pb and 48 Ca from the parity-violating electron scattering experiments (PREX-II and CREX at JLab as well as MREX at MESA), within the Bayesian statistical framework using the Skyrme-Hartree-Fock model we infer the posterior probability distribution functions (PDFs) of the slope parameter L of the nuclear symmetry energy at ρ 0 from imagined Δr np ( 208 Pb)=0.15, 0.20, and 0.30 fm with a 1σ error bar of 0.02, 0.04, and 0.06 fm, respectively, as well as Δr np ( 48 Ca)=0.12, 0.15, and 0.25 fm, with different 1σ error bar of 0.01 and 0.02 fm, respectively. The results are compared with the PDFs of L inferred using the same approach from the available Δr np data for 116, 118, 120, 122, 124, 130, 132 Sn from hadronic probes. They are also compared with results from a recent Bayesian analysis of the radius and tidal deformability data of canonical neutron stars from GW170817 and NICER. The neutron skin data for Sn isotopes gives L = 45.5 $^{+ 26.5}_{-21.6}$ MeV surrounding its mean value or L = 53 . 4 $^{+ 18.6}_{ -29.5}$ MeV surrounding its maximum a posteriori value, respectively, with the latter smaller than but consistent with the L = 66 $^{+ 12}_{-20}$ MeV from the neutron star data within their 68% confidence intervals. We found that Δr np = 0.17 –0.18 fm in 208 Pb with an error bar of about 0.02 fm leads to a PDF of L compatible with that from analyzing the Sn data. To provide additionally useful information on L extracted from the Δr np of Sn isotopes, the experimental error bar of Δr np in 208 Pb should be at least smaller than 0.06 fm aimed by some current experiments. In addition, the Δr np ( 48 Ca) needs to be larger than 0.15 fm but smaller than 0.25 fm to be compatible with the Sn and/or neutron star results. To further improve our current knowledge about L and distinguish its PDFs in the examples considered, even higher precisions leading to significantly less than ±20 MeV error bars for L at 68% confidence level are necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Proton pairing in neutron stars from chiral effective field theory

Here, we study the 1 S 0 proton pairing gap in beta-equilibrated neutron star matter within the framework of chiral effective field theory. We focus on the role of three-body forces, which strongly modify the effective proton-proton spin-singlet interaction in dense matter. We find that three-body forces generically reduce both the size of the pairing gap and the maximum density at which proton pairing may occur. The pairing gap is computed within BCS theory using a single-particle dispersion relation calculated up to second order in perturbation theory. Model uncertainties are estimated by varying the nuclear potential (its order in the chiral expansion and high-momentum cutoff) and the choice of single-particle spectrum in the gap equation. We find that a second-order perturbative treatment of the single-particle spectrum suppresses the proton 1 S 0 pairing gap relative to the use of a free spectrum. We estimate the critical temperature for the onset of proton superconductivity to be T c = (3:2–5:1) 10 9 K, which is consistent with previous theoretical results in the literature and marginally within the range deduced from a recent Bayesian analysis of neutron star cooling observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Limiting masses and radii of neutron stars and their implications

We combine the equation of state of dense matter up to twice nuclear saturation density 𝑛 sat obtained using chiral effective field theory (𝜒⁢EFT) and recent observations of neutron stars to gain insights about the high-density matter encountered in their cores. A key element in our study is the recent Bayesian analysis of correlated EFT truncation errors based on order-by-order calculations up to next-to-next-to-next-to-leading order in the 𝜒⁢EFT expansion. We refine the bounds on the maximum mass imposed by causality at high densities and provide stringent limits on the maximum and minimum radii of ∼1.4M ⊙ and ∼2.0M ⊙ stars. Including 𝜒⁢EFT predictions from 𝑛sat to 2𝑛 sat reduces the permitted ranges of the radius of a ~1.4M ⊙ star, 𝑅 1.4 , by ∼3.5km. If observations indicate 𝑅 1.4 < 11.2km, then our study implies that either the squared speed of sound 𝑐$^{2}_{𝑠}$ > 1/2 for densities above 2𝑛 sat or that 𝜒⁢EFT breaks down below 2𝑛 sat . Here, we also comment on the nature of the secondary compact object in GW190814 with mass ≃ 2.6M ⊙ and discuss the implications of massive neutron stars > 2.1⁢M ⊙ ⁡(2.6M ⊙ ) in future radio and gravitational-wave searches. Some form of strongly interacting matter with 𝑐$^{2}_{𝑠}$ > 0.35⁢(0.55) must be realized in the cores of such massive neutron stars. In the absence of phase transitions below 2𝑛 sat , the small tidal deformability inferred from GW170817 lends support for the relatively small pressure predicted by 𝜒⁢EFT for the baryon density 𝑛 B in the range 1–2𝑛 sat . Together they imply that the rapid stiffening required to support a high maximum mass should occur only when 𝑛 B ≳ 1.5–1.8 𝑛 sat .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Effective field theory for dilute Fermi systems at fourth order

We discuss high-order calculations in perturbative effective field theory for fermions at low energy scales. The Fermi-momentum or k F a s as expansion for the ground-state energy of the dilute Fermi gas is calculated to fourth order, both in cutoff regularization and in dimensional regularization. For the case of spin one-half fermions we find from a Bayesian analysis that the expansion is well converged at this order for |k F a s |≲ 0.5. Furthermore, we show that Padé-Borel resummations can improve the convergence for |k F a s |≲ 1. Here, our results provide important constraints for nonperturbative calculations of ultracold atoms and dilute neutron matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

O 16 O 16 collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider comparing α clustering versus substructure

Collisions of . light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed 16 O 16 O run at the CERN Large Hadron Collider (LHC) and at the BNL Relativistic Heavy Ion Collider (RHIC) we study the potential for finding α clustering in 16 O. Here we use the state-of-the-art iEBE-VISHNU package with 16 O nucleonic configurations from ab initio nuclear lattice simulations. This setup was tuned using a Bayesian analysis on p Pb and PbPb systems. We find that the 16 O 16 O system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between α-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Constraints on partial half-lives of 136 Ce and 138 Ce double electron captures

The γ-ray emissions from a radiopure cerium-bromide crystal with a mass of 4381 g were measured for a total of 497.4 d by means of high-resolution γ-ray spectrometry in the HADES underground laboratory at a depth of 500 m.w.e. A search for 0/2νϵϵ and 0/2νϵβ + double beta decay transitions of 136 Ce and 13 8Ce was performed using Bayesian analysis techniques. No signals were observed for a total of 35 investigated decay modes. 90% credibility limits were set in the order of 10 18 -10 19 yr. Existing constraints from a cerium oxide powder measurement were tested with a different cerium compound and half-life limits could be improved for most of the decay modes. The most likely accessible decay mode of the 136 Ce 2νϵϵ transition into the $0_{1}^{+}$ state of 136 Ba results in a new best 90% credibility limit of 5.0×10 18 yr.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Sensitivity of Au + Au collisions to the symmetric nuclear matter equation of state at 2–5 nuclear saturation densities

We demonstrate that proton and pion flow measurements in heavy-ion collisions at incident energies ranging from 1 to 20 GeV per nucleon in the fixed target frame can be used for an accurate determination of the symmetric nuclear matter equation of state at baryon densities equal 2–4 times nuclear saturation density n 0 . We simulate Au + Au collisions at these energies using a hadronic transport model with an adjustable vector mean-field potential dependent on baryon density n B . Here, we show that the mean field can be parametrized to reproduce a given density dependence of the speed of sound at zero temperature $c$$^{2}_{s}$ (n B , T = 0), which we vary independently in multiple density intervals to probe the differential sensitivity of heavy-ion observables to the equation of state at these specific densities. Recent flow data from the STAR experiment at the center-of-mass energies √ s NN = {3.0, 4.5} GeV can be described by our model, and a Bayesian analysis of these data indicates a hard equation of state at n B ϵ (2, 3)n 0 and a possible phase transition at n B ϵ (3, 4)n 0 . More data at √ s NN = 2–5 GeV, as well as a more thorough analysis of the model systematic uncertainties will be necessary for a more precise conclusion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Linearized Boltzmann-Langevin model for heavy quark transport in hot and dense QCD matter

In relativistic heavy-ion collisions, the production of heavy quarks at large transverse momenta is strongly suppressed compared to proton-proton collisions. In addition, an unexpectedly large azimuthal anisotropy was observed for the emission of charmed hadrons in noncentral collisions. Both observations pose challenges to the theoretical understanding of the coupling between heavy quarks and the quark-gluon plasma produced in these collisions. Transport models for the evolution of heavy quarks in a QCD medium offer the opportunity to study these effects; two of the most successful approaches are based on the linearized Boltzmann transport equation and the Langevin equation. In this work, we develop a hybrid transport model that combines the strengths of both of these approaches: Heavy quarks scatter with medium partons using matrix-elements calculated in perturbative QCD, while between these discrete hard scatterings they evolve using a Langevin equation with empirical transport coefficients to capture the nonperturbative soft part of the interaction. With the hybrid transport model coupled to a state-of-the-art event-by-event bulk evolution model based on 2+1D relativistic viscous fluid dynamics, we study the azimuthal anisotropy and nuclear modification factor of heavy quarks in Pb+Pb collisions at √s = 5.02 TeV. The parameters related to heavy-flavor transport are calibrated using a Bayesian analysis comparing them to available D -meson and B -meson data at the Large Hadron Collider. Lastly, using the calibrated model, we study the implications on heavy-flavor transport properties and predict observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Astrometry in two-photon interferometry using an Earth rotation fringe scan

Optical interferometers may not require a phase-stable optical link between the stations if instead sources of quantum-mechanically entangled pairs could be provided to them, enabling long baselines. We developed a new variation of this idea, proposing that photons from two different astronomical sources could be interfered at two decoupled stations. Interference products can then be calculated in post-processing or requiring only a slow, classical connection between stations. In this work, we investigated practical feasibility of this approach. We developed a Bayesian analysis method for the earth rotation fringe scanning technique and showed that in the limit of high signal-to-noise ratio it reproduced the results from a simple Fisher matrix analysis. Here we identify candidate stair pairs in the northern hemisphere, where this technique could be applied. With two telescopes with an effective collecting area of ~ 2 m 2 , we could detect fringing and measure the astrometric separation of the sources at ~100 μas precision in a few hours of observations, in agreement with previous estimates.

47 OTHER INSTRUMENTATION↗