skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Bit Blasting Probabilistic Programs
Probabilistic programming languages (PPLs) are an expressive means for creating and reasoning about probabilistic models. Unfortunatelyhybridprobabilistic programs that involve both continuous and discrete structures are not well supported by today’s PPLs. In this paper we develop a new approximate inference algorithm for hybrid probabilistic programs that first discretizes the continuous distributions and then performs discrete inference on the resulting program. The key novelty is a form of discretization that we callbit blasting, which uses a binary representation of numbers such that a domain of 2 b discretized points can be succinctly represented as a discrete probabilistic program overpoly b Boolean random variables. Surprisingly, we prove that many common continuous distributions can be bit blasted in a manner that incurs no loss of accuracy over an explicit discretization and supports efficient probabilistic inference. We have built a probabilistic programming system for hybrid programs calledHyBit, which employs bit blasting followed by discrete probabilistic inference. We empirically demonstrate the benefits of our approach over existing sampling-based and symbolic inference approaches  more » « less
Award ID(s):
2220408
PAR ID:
10612728
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Association for Computing Machinery (ACM)
Date Published:
Journal Name:
Proceedings of the ACM on Programming Languages
Volume:
8
Issue:
PLDI
ISSN:
2475-1421
Format(s):
Medium: X Size: p. 865-888
Size(s):
p. 865-888
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We studyℓnorms ofℓ2-normalized eigenfunctions of quantum cat maps. For maps with short quantum periods (constructed by Bonechi and de Biévre in F Bonechi and S De Bièvre (2000,Communications in Mathematical Physics,211, 659–686)) we show that there exists a sequence of eigenfunctionsuwith u log N 1 / 2 . For general eigenfunctions we show the upper bound u log N 1 / 2 . Here the semiclassical parameter is h = 2 π N 1 . Our upper bound is analogous to the one proved by Bérard in P Bérard (1977,Mathematische Zeitschrift,155, 249-276) for compact Riemannian manifolds without conjugate points. 
    more » « less
  2. Abstract Objective.Combined proton–photon treatments, where most fractions are delivered with photons and only a few are delivered with protons, may represent a practical approach to optimally use limited proton resources. It has been shown that, when organs at risk (OARs) are located within or near the tumor, the optimal multi-modality treatment uses protons to hypofractionate parts of the target volume and photons to achieve near-uniform fractionation in dose-limiting healthy tissues, thus exploiting the fractionation effect. These plans may be sensitive to range and setup errors, especially misalignments between proton and photon doses. Thus, we developed a novel stochastic optimization method to directly incorporate these uncertainties into the biologically effective dose (BED)-based simultaneous optimization of proton and photon plans.Approach.The method considers the expected value E b and standard deviation σ b of the cumulative BED b in every voxel of a structure. For the target, a piecewise quadratic penalty function of the form b min E b 2 σ b + 2 is minimized, aiming for plans in which the expected BED minus two times the standard deviation exceeds the prescribed BED b min . Analogously, E b + 2 σ b b max + 2 is considered for OARs.Main results.Using a spinal metastasis case and a liver cancer patient, it is demonstrated that the novel stochastic optimization method yields robust combined treatment plans. Tumor coverage and a good sparing of the main OARs are maintained despite range and setup errors, and especially misalignments between proton and photon doses. This is achieved without explicitly considering all combinations of proton and photon error scenarios.Significance.Concerns about range and setup errors for safe clinical implementation of optimized proton–photon radiotherapy can be addressed through an appropriate stochastic planning method. 
    more » « less
  3. Abstract We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise ratio integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central ∼150 pc out to ∼2.5 effective radii. Observations from Hubble Space Telescope WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6D galaxy model parameter space to determine NGC 2693's SMBH mass (MBH), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We find M BH = 1.7 ± 0.4 × 10 9 M and a triaxial intrinsic shape with axis ratiosp=b/a= 0.902 ± 0.009 and q = c / a = 0.721 0.010 + 0.011 , triaxiality parameterT= 0.39 ± 0.04. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 prefer M BH = 2.4 ± 0.6 × 10 9 M and M BH = 2.9 ± 0.3 × 10 9 M , respectively. Neither model can account for the non-axisymmetric stellar velocity features present in the IFS data. 
    more » « less
  4. MXenes have demonstrated potential for various applications owing to their tunable surface chemistry and metallic conductivity. However, high temperatures can accelerate MXene film oxidation in air. Understanding the mechanisms of MXene oxidation at elevated temperatures, which is still limited, is critical in improving their thermal stability for high-temperature applications. Here, we demonstrate that Ti 3 C 2 T x MXene monoflakes have exceptional thermal stability at temperatures up to 600 ° C in air, while multiflakes readily oxidize in air at 300 ° C. Density functional theory calculations indicate that confined water between Ti 3 C 2 T x flakes has higher removal energy than surface water and can thus persist to higher temperatures, leading to oxidation. We demonstrate that the amount of confined water correlates with the degree of oxidation in stacked flakes. Confined water can be fully removed by vacuum annealing Ti 3 C 2 T x films at 600 ° C, resulting in substantial stability improvement in multiflake films (can withstand 600 ° C in air). These findings provide fundamental insights into the kinetics of confined water and its role in Ti 3 C 2 T x oxidation. This work enables the use of stable monoflake MXenes in high-temperature applications and provides guidelines for proper vacuum annealing of multiflake films to enhance their stability. 
    more » « less
  5. Abstract Ultra-high-energy cosmic rays (UHECRs), particles characterized by energies exceeding 1018eV, are generally believed to be accelerated electromagnetically in high-energy astrophysical sources. One promising mechanism of UHECR acceleration is magnetized turbulence. We demonstrate from first principles, using fully kinetic particle-in-cell simulations, that magnetically dominated turbulence accelerates particles on a short timescale, producing a power-law energy distribution with a rigidity-dependent, sharply defined cutoff well approximated by the form f cut E , E cut = sech ( E / E cut ) 2 . Particle escape from the turbulent accelerating region is energy dependent, withtesc∝E−δandδ∼ 1/3. The resulting particle flux from the accelerator follows dN / dEdt E s sech ( E / E cut ) 2 , withs∼ 2.1. We fit the Pierre Auger Observatory’s spectrum and composition measurements, taking into account particle interactions between acceleration and detection, and show that the turbulence-associated energy cutoff is well supported by the data, with the best-fitting spectral index being s = 2.1 0.13 + 0.06 . Our first-principles results indicate that particle acceleration by magnetically dominated turbulence may constitute the physical mechanism responsible for UHECR acceleration. 
    more » « less