Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability generating function”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

2D k -th nearest neighbour statistics: a highly informative probe of galaxy clustering

ABSTRACT Beyond standard summary statistics are necessary to summarize the rich information on non-linear scales in the era of precision galaxy clustering measurements. For the first time, we introduce the 2D k-th nearest neighbour (kNN) statistics as a summary statistic for discrete galaxy fields. This is a direct generalization of the standard 1D kNN by disentangling the projected galaxy distribution from the redshift-space distortion signature along the line-of-sight. We further introduce two different flavours of 2D kNNs that trace different aspects of the galaxy field: the standard flavour which tabulates the distances between galaxies and random query points, and a ‘DD’ flavour that tabulates the distances between galaxies and galaxies. We showcase the 2D kNNs’ strong constraining power both through theoretical arguments and by testing on realistic galaxy mocks. Theoretically, we show that 2D kNNs are computationally efficient and directly generate other statistics such as the popular two-point correlation function (2PCF), voids probability function, and counts-in-cell statistics. In a more practical test, we apply the 2D kNN statistics to simulated galaxy mocks that fold in a large range of observational realism and recover parameters of the underlying extended halo occupation distribution (HOD) model that includes velocity bias and galaxy assembly bias. We find unbiased and significantly tighter constraints on all aspects of the HOD model with the 2D kNNs, both compared to the standard 1D kNN, and the classical redshift-space 2PCF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Tunable Geometries in Sparse Clifford Circuits

We investigate the emergence of different effective geometries in stochastic Clifford circuits with sparse coupling. By changing the probability distribution for choosing two-site gates as a function of distance, we generate sparse interactions that either decay or grow with distance as a function of a single tunable parameter. Tuning this parameter reveals three distinct regimes of geometry for the spreading of correlations and growth of entanglement in the system. We observe linear geometry for short-range interactions, treelike geometry on a sparse coupling graph for long-range interactions, and an intermediate fast scrambling regime at the crossover point between the linear and treelike geometries. This transition in geometry is revealed in calculations of the subsystem entanglement entropy and tripartite mutual information. We also study emergent lightcones that govern these effective geometries by teleporting a single qubit of information from an input qubit to an output qubit. These tools help to analyze distinct geometries arising in dynamics and correlation spreading in quantum many-body systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
DOE OSTI2022

Constraining dark matter microphysics with the annihilation signal from subhalos

In the cold dark matter scenario, galactic dark matter halos are populated with a large number of smaller subhalos. Previous work has shown that dark matter annihilations in subhalos can generate a distinctive, non-Poisson signal in the gamma-ray photon counts probability distribution function (PDF). Here we show that the gamma-ray PDF also carries information about the velocity dependence of the dark matter annihilation cross section. After calculating the PDF assuming s-wave and Sommerfeld-enhanced annihilation, we perform a mock data analysis to illustrate how current and future observations can constrain the microphysics of the dark matter annihilation. We find that, with current Fermi data, and assuming a dark matter annihilation cross section roughly at the limit of current bounds from annihilation in dwarf spheroidal galaxies, one can potentially distinguish the non-Poissonian fluctuations expected from dark matter annihilation in subhalos from Poisson sources, as well as from dark matter models with an incorrect velocity-dependence. Here we explore how robust these results are to assumptions about the modeling of the galactic gamma-ray background, but further work is needed to determine the impact of realistic astrophysical source populations on our results. We also point out a four-parameter degeneracy between the velocity dependence of the dark matter annihilation, the minimum subhalo mass, the power law index of the subhalo mass function, and the normalization of the dark matter signal. This degeneracy can be broken with priors from N-body simulations or from observational constraints on the subhalo mass function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

How to Obtain the Redshift Distribution from Probabilistic Redshift Estimates

Abstract A reliable estimate of the redshift distribution n ( z ) is crucial for using weak gravitational lensing and large-scale structures of galaxy catalogs to study cosmology. Spectroscopic redshifts for the dim and numerous galaxies of next-generation weak-lensing surveys are expected to be unavailable, making photometric redshift (photo- z ) probability density functions (PDFs) the next best alternative for comprehensively encapsulating the nontrivial systematics affecting photo- z point estimation. The established stacked estimator of n ( z ) avoids reducing photo- z PDFs to point estimates but yields a systematically biased estimate of n ( z ) that worsens with a decreasing signal-to-noise ratio, the very regime where photo- z PDFs are most necessary. We introduce Cosmological Hierarchical Inference with Probabilistic Photometric Redshifts ( CHIPPR ), a statistically rigorous probabilistic graphical model of redshift-dependent photometry that correctly propagates the redshift uncertainty information beyond the best-fit estimator of n ( z ) produced by traditional procedures and is provably the only self-consistent way to recover n ( z ) from photo- z PDFs. We present the chippr prototype code, noting that the mathematically justifiable approach incurs computational cost. The CHIPPR approach is applicable to any one-point statistic of any random variable, provided the prior probability density used to produce the posteriors is explicitly known; if the prior is implicit, as may be the case for popular photo- z techniques, then the resulting posterior PDFs cannot be used for scientific inference. We therefore recommend that the photo- z community focus on developing methodologies that enable the recovery of photo- z likelihoods with support over all redshifts, either directly or via a known prior probability density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Further development of finite-temperature density functional theory (Technical Report on FG02-08ER46496)

The final year of funding was spent supporting research into extending conditional probability theory to generate the thermal dependence of PBE. The necessary background work led to a paper being published on the uniform electron gas (Dennis Perchak, Ryan J. McCarty, and Kieron Burke, Phys. Rev. B 105, 165143 (2022).). The graduate student John Kozlowski also contributed to a mathematical paper about DFT (Steven Crisostomo, Ryan Pederson, John Kozlowski, Bhupalee Kalita, Antonio C. Cancio, Kiril Datchev, Adam Wasserman, Suhwan Song, and Kieron Burke, Letters in Mathematical Physics 113, 42 (2023).).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Estimating Cosmological Constraints from Galaxy Cluster Abundance using Simulation-Based Inference

Inferring the values and uncertainties of cosmological parameters in a cosmology model is of paramount importance for modern cosmic observations. In this paper, we use the simulation-based inference (SBI) approach to estimate cosmological constraints from a simplified galaxy cluster observation analysis. Using data generated from the Quijote simulation suite and analytical models, we train a machine learning algorithm to learn the probability function between cosmological parameters and the possible galaxy cluster observables. The posterior distribution of the cosmological parameters at a given observation is then obtained by sampling the predictions from the trained algorithm. Our results show that the SBI method can successfully recover the truth values of the cosmological parameters within the 2σ limit for this simplified galaxy cluster analysis, and acquires similar posterior constraints obtained with a likelihood-based Markov Chain Monte Carlo method, the current state-of the-art method used in similar cosmological studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Practical Probabilistic Programming

Recent advances in probabilistic programming languages (PPLs) have provided the capability for exact inference: computing a closed-form probability distribution for a given probabilistic program. In particular, the new language Roulette uses a language oriented programming (LOP) approach, wherein analysts build new programming languages on top of a set of primitives provided by Roulette, which then translates these structures into a weighted model counting problem which can be solved by automated reasoning tools. However, because Roulette provides few convenience features, developing these new languages is challenging even for expert users. We developed a standard library of common probability functions for Roulette with the goal of improved usability. This included approximation of continuous probability density functions using discrete probability mass functions. We demonstrated this approach by modeling a cosmic ray striking a RAM controller. We found that Roulette provides a powerful interface for highly expressive probabilistic programs to be generated. In collaboration with the NNSA Advanced Simulation and Computing program, which resulted in development of a tool called Circulette, we were able to model complex circuits expressed in Verilog using probabilistic programs with an expressivity not previously possible. Our research question that motivated the development of a Roulette standard library was to determine whether non-experts could use a PPL to model relevant problems regarding radiation effects on microelectronics. This standard library improved the expressivity of Roulette by implementing common probability density functions, mathematical operators on distributions, and support for empirical distributions. While Roulette is a powerful modeling language, the untyped, LOP approach makes error messages difficult to understand and requires expert aid. We recommend further research on Roulette, especially with its error messages, to enable improved usability. At the same time, this project demonstrated that for users familiar with Roulette and the LOP approach, Roulette provides powerful new capabilities that can be integrated with other Sandia modeling capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Fidelity Comparison of Time-Bin and Fock State Encoding in Hybrid Quantum Systems Under Channel and Transduction Effects

Future quantum networks are expected to integrate a heterogeneous combination of quantum systems, owing to the diverse advantages offered by different physical platforms in terms of scalability, coherence times, and interfacing capabilities. Within the context of this emerging quantum communication paradigm, this paper presents an analytical comparison of two photonic encoding schemes—time-bin and Fock state (single rail)—in hybrid quantum systems where flying qubits are entangled with stationary matter qubits. We evaluate their resilience against fiber channel and quantum transduction effects by calculating the fidelity of the final states relative to their ideal forms. Employing the characteristic function approach, we derive analytical fidelity expressions and investigate their dependence on parameters such as transmissivity, noise levels, fiber length, and source generation success probability. Additionally, we simulate the scenario with a dedicated QuTiP software implementation to verify the validity of the theoretical models. Our findings reveal that due to its inherent single-mode structure, the Fock state encoding consistently outperforms time-bin encoding in fidelity, as this structure significantly minimizes susceptibility to losses compared to the two-mode nature of the time-bin scheme. This analysis offers valuable insights for future hybrid quantum communication and information processing applications.

Fiorini, Francesco [Pisa U.]↗
DOE OSTI2021

A consistent, Bayesian, approach to the cross section probability distribution in the unresolved resonance region

The cross sections of neutron-induced reactions can be divided into three energy ranges: the resolved resonance region (RRR), the unresolved resonance region (URR), and the fast region. In general, the cross sections in the URR show significant fluctuations that cannot be predicted and cannot be experimentally resolved, thus, it is commonly assumed that the cross section at a specific energy is given by a probability distribution function (PDF) over a range of values that can span several orders of magnitude. The current methodology used to describe such behavior is to construct the PDF by stochastically generating resonance ladders and numerically measuring the PDF. The resonance ladders are sampled using known resonance statistical properties and average resonance widths and spacings extrapolated from the RRR. Although this is a standard and widely used technique, it is computationally very expensive, therefore, an alternative, analytical, approach would be preferable due to the considerable speed up of the computational time in real life applications. Moreover, the current methodology does not take into account existing experimental data, such for total and capture cross sections, that are available for many nuclei. Finally, this approach was developed to be used in reactor-scale applications and it is not suited for use in single-event applications. In this work we will rethink the entire approach to the PDF construction using a Bayesian mindset. This will allow us to provide a different definition of the PDF that allows a much faster calculation of the higher-temperature PDFs and a proper combination of theoretical and experimental PDFs following the probability theory. We will also show that our definition is well suited for single-event applications and we will make an explicit connection between our method and the standard approach. We do this by showing that the central limit theorem applies and our method leads to the same PDF obtained with the standard methodology, for a large number of events per history.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Muon-induced fission as a probe of the underlying dynamics in nuclear fission

Muon-induced fission could be utilized as a probe to study the underlying dynamics of nuclear fission. Here, the probability of muon attachment to the light asymmetric fission fragment is sensitive to fission dynamics, such as the timescale and friction of the fission event, charge asymmetry, and possibly the shape of the fission fragments. We focus on muonic atoms that are formed with actinide nuclei. A relativistic approach is employed, solving the Dirac equation for the muonic wave function in the presence of a time-dependent electromagnetic field generated by the fissioning nucleus. Computations are carried out on a three-dimensional Cartesian lattice with no symmetry assumptions. The results show a strong dependence of the attachment probability on the fission charge asymmetry and a more modest dependence on friction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Turbine Damage Probability

Damage probability maps for offshore wind turbines exposed to tropical cyclones (TCs) under both historical and future climate scenarios along the U.S. Atlantic and Gulf Coasts are presented in this dataset. TCs are generated using the Risk Analysis Framework for Tropical Cyclones (RAFT), forced by CMIP6 historical and future global climate simulations. Maximum wind speeds for 20- and 50-year TCs are processed through a fragility function specific to offshore wind (OSW) turbines in order to estimate the probability of damage – specifically yielding and buckling – based on wind speed intensity.

17 WIND ENERGY↗
DOE OSTI2025

Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains

Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.

Bounded domains↗
DOE OSTI2022

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.↗
DOE OSTI2022

Stable High‐Conductivity Ethylenedioxythiophene Polymers via Borane‐Adduct Doping

Efficient doping of polymer semiconductors is required for high conductivity and efficient thermoelectric performance. Lewis acids, e.g., B(C 6 F 5 ) 3 , have been widely employed as dopants, but the mechanism is not fully understood. 1:1 “Wheland type” or zwitterionic complexes of B(C 6 F 5 ) 3 are created with small conjugated molecules 3,6-bis(5-(7-(5-methylthiophen-2-yl)-2,3-dihydrothieno[3,4-b][1,4]dioxin-5-yl)thiophen-2-yl)-2,5-dioctyl-2,5-dihydropyrrolo[3,4-c]pyrrole-1,4-dione [oligo_DPP(EDOT) 2 ] and 3,6-bis(5''-methyl-[2,2':5',2''-terthiophen]-5-yl)-2,5-dioctyl-2,5-dihydropyrrolo[3,4-c]pyrrole-1,4-dione [oligo_DPP(Th) 2 ]. Using a wide variety of experimental and computational approaches, the doping ability of these Wheland Complexes with B(C 6 F 5 ) 3 are characterized for five novel diketopyrrolopyrrole-ethylenedioxythiophene (DPP-EDOT)-based conjugated polymers. The electrical properties are a strong function of the specific conjugated molecule constituting the adduct, rather than acidic protons generated via hydrolysis of B(C 6 F 5 ) 3 , serving as the oxidant. It is highly probable that certain repeat units/segments form adduct structures in p -type conjugated polymers which act as intermediates for conjugated polymer doping. Electronic and optical properties are consistent with the increase in hole-donating ability of polymers with their cumulative donor strengths. The doped film of polymer (DPP(EDOT) 2 -(EDOT) 2 ) exhibits exceptionally good thermal and air-storage stability. The highest conductivities, ≈300 and ≈200 S cm -1 , are achieved for DPP(EDOT) 2 -(EDOT) 2 doped with B(C 6 F 5 ) 3 and its Wheland complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Retrieval of the molecular orientation distribution from atom-pair angular distributions

Imaging laser-induced rotational dynamics is an important and active field due to its applications in capturing reactions in the molecular frame and in molecular imaging. Experimental measurement of the molecular orientation distribution, as a function of the Euler angles, has only been demonstrated for special cases when the detectable signal is generated along the molecular symmetry axis. Here we developed the general theory that maps the probability density distribution of the molecular orientation to the atom-pair angular distributions for nonlinear molecules. With the theory, the molecular orientation distribution can be retrieved from the measured atom-pair angular distribution, which we demonstrate experimentally using ultrafast electron diffractive imaging of impulsively aligned trifluoro-iodomethane molecules. Here, the retrieved molecular orientation distribution is in good agreement with direct numerical simulations of the time-dependent Schrodinger equation using the experimental conditions. Unlike the existing retrieval methods, the retrieval method does not require solving Schrodinger equation, works for any alignment method, and is in principle applicable to asymmetric top molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Correlation energy of the uniform electron gas determined by ground-state conditional probability density functional theory

Conditional-probability density functional theory (CP-DFT) is a formally exact method for finding correlation energies from Kohn-Sham DFT without evaluating an explicit energy functional. We present details on how to generate accurate exchange-correlation energies for the ground-state uniform gas. We also use the exchange hole in a CP antiparallel spin calculation to extract the high-density limit. We give a highly accurate analytic solution to the Thomas-Fermi model for this problem, showing its performance relative to Kohn-Sham and may be useful at high temperatures. We explore several approximations to the CP potential. Furthermore, results are compared to accurate parameterizations for both exchange-correlation energies and holes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Refining fast calorimeter simulations with a Schrödinger Bridge

Machine learning-based simulations, especially calorimeter simulations, are promising tools for approximating the precision of classical high energy physics simulations with a fraction of the generation time. Nearly all methods proposed so far learn neural networks that map a random variable with a known probability density, like a Gaussian, to realistic-looking events. In many cases, physics events are not close to Gaussian and so these neural networks have to learn a highly complex function. We study an alternative approach: Schrödinger bridge Quality Improvement via Refinement of Existing Lightweight Simulations (SQuIRELS). SQuIRELS leverages the power of diffusion-based neural networks and Schrödinger bridges to map between samples where the probability density is not known explicitly. We apply SQuIRELS to the task of refining a classical fast simulation to approximate a full classical simulation. On simulated calorimeter events, we find that SQuIRELS is able to reproduce highly non-trivial features of the full simulation with a fraction of the generation time.

Calorimeter methods↗