Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pair distribution functions”

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 631 records · Page 35

System and method for reasoning about the optimality of a configuration parameter of a distributed system

Embodiments provide a system and method for reasoning about the optimality of a configuration parameter of a distributed system. During operation, the system obtains a multi-layer graph for a system with a plurality of components, wherein the multi-layer graph comprises a configuration subgraph, a vulnerability subgraph, and a dependency subgraph. The system determines, based on the multi-layer graph, constraint relationships associated with configuration parameters for the components, wherein the constraint relationships include security constraints and functionality constraints. The system computes an unsatisfiable core which comprises a set of mutually incompatible constraints. The system resolves, based on a strategy and over multiple iterations, the unsatisfiable core by analyzing one pair of mutually incompatible constraints per a respective iteration, to obtain a new unsatisfiable core which comprises a smaller number of mutually incompatible constraints than the computed unsatisfiable core or a previously computed unsatisfiable core from a most recent iteration.

Soroush, Hamed↗

Unimodal Imaging of Monovalent Metal-Chelator Complexes and Lipids by MALDI Imaging Mass Spectrometry

Careful regulation of monovalent metal ions (M + ) is necessary to maintain a functional cellular system. Of these ions, appropriate sodium (Na + ) and potassium (K + ) concentrations are particularly integral for electrochemical signaling, as well as the secondary transport of nutrients and waste. Dysregulation of M + homeostasis can disrupt these mechanisms, potentially influencing the metabolism of downstream biomolecules such as lipids. Thus, the relationship between M + abundances and related biomolecular distributions must be elucidated to better understand the physiology of healthy and disordered tissues. Traditional techniques for imaging biological metal distributions include SIMS, LA-ICP-MS, and XRF; however, these capabilities are limited to elemental analysis or the analysis of molecular fragments and must be paired with other modalities to visualize distributions of more complex biomolecules within the same or similar samples. Conversely, matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) is a powerful tool often used for mapping such biomolecular distributions, but current methods are unable to detect metals within tissue. This study illustrates a novel methodology that adds metal detection to the MALDI IMS repertoire through which the simultaneous detection of M + metals and lipids is achievable. Using a robotic sprayer for homogeneous application, on-tissue deposition of the chelator deferiprone (DEF) enables subsequent detection of the ionizable metal-chelator complex by MALDI without hindering lipid detection. Our work provides proof-of-concept data for the simultaneous detection of K + , Na + , and intact lipids using MALDI IMS.

59 BASIC BIOLOGICAL SCIENCES↗

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of the $c$-jet mistagging efficiency in $t\bar{t}$ events using $pp$ collision data at $\sqrt{s}=13$ $\text {TeV}$ collected with the ATLAS detector

A technique is presented to measure the efficiency with which c-jets are mistagged as b-jets (mistagging efficiency) using tt¯ events, where one of the W bosons decays into an electron or muon and a neutrino and the other decays into a quark–antiquark pair. The measurement utilises the relatively large and known W → cs branching ratio, which allows a measurement to be made in an inclusive c-jet sample. The data sample used was collected by the ATLAS detector at √s = 13 TeV and corresponds to an integrated luminosity of 139 fb –1 . Events are reconstructed using a kinematic likelihood technique which selects the mapping between jets and tt¯ decay products that yields the highest likelihood value. The distribution of the b-tagging discriminant for jets from the hadronic W decays in data is compared with that in simulation to extract the mistagging efficiency as a function of jet transverse momentum. The total uncertainties are in the range 3–17%. The measurements generally agree with those in simulation but there are some differences in the region corresponding to the most stringent b-jet tagging requirement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Radiation Mapping for an Unmanned Aerial Vehicle: Development and Simulated Testing of Algorithms for Source Mapping and Navigation Path Generation

Image reconstruction algorithms were developed for radiation source mapping and used for generating the search path of a moving radiation detector, such as one onboard an unmanned aerial vehicle. Simulations consisted of first assuming radioactive sources of varying complexity and estimating the radiation fields that would then be produced by that source distribution. Next, the "measurements" that would result from a pair of adjacent spatial locations were computed. A crude estimate of the source distribution likely to have produced such "measurements" was reconstructed based upon the limited measurements. Location of the next "measurement" was then determined as halfway between the location of the estimated source and the current "measurement." With each additional sample, improved source distribution reconstructions were made and used to inform the immediate direction of detector motion. Source reconstruction or mapping was formulated as an inverse problem solved with either maximum a posteriori or least squares (LS) regression deconvolution methods. Different amounts of noise were added to the simulated "measurements," allowing evaluation of the methods' performances as functions of signal-to-noise ratio of the measured map. As expected, methods that promote sparsity were better suited in reconstructing point sources. Reliable prior information of the source distribution also improved the reconstruction results, especially with distributed sources. With a non-negative least square algorithm and the suggested paths it generated, location of sources was successfully estimated to an accuracy of 0.014 m within nine iterations in a single-source scenario and 12 iterations in a two-source scenario, given a 10% error on the integrated counts and a Poisson distribution of the noise associated with the measured counts.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Measurements of Beam Spin Asymmetries in p+p0 and p´p0 Dihadron Production at CLAS12

Semi-Inclusive Deep Inelastic Scattering (SIDIS) is a powerful experimental tool for studying the internal structure and dynamics of the proton, revealing how quarks and gluons are distributed and interact within it. SIDIS describes a process where an elec tron scatters off one of the constituent quarks within the proton, causing it to undergo hadronization, creating multiple hadrons in the final state. Through factorization, the full process can be split into probabilistic components: one which describes the internal structure of the proton using Parton Distribution Functions (PDFs), and another which describes the hadronization process using Fragmentation Functions (FFs). These functions are non-perturbative quantities of Quantum Chromodynamics (QCD), meaning they cannot be calculated directly from first principles and must instead be extracted from experimental measurements. Acommon approach for accessing PDFs and FFs using SIDIS is to measure asymmetries. In this context, asymmetries correspond to subtle differences in the angular distribution of outgoing particles that arise when the spin orientation of the incoming beam or target is reversed. Because many of these effects only appear when spin is involved, they isolate specific, nuanced properties of the proton’s spin-structure that are otherwise hidden in spin averaged measurements. In practice, they show up as specific azimuthal modulations (e.g., sin ¿R, sin(¿h ´ ¿R)), whose amplitudes isolate convolutions of PDFs and FFs at leading and subleading twist. Non-zero asymmetries of these angular distributions can be traced back to unique combinations of PDFs and FFs, offering a way to probe them directly. In this work, we measure SIDIS by analyzing high energy electron-proton scattering events using the CLAS12 detector at Jefferson Lab. This study focuses on subset of SIDIS referred to as dihadron SIDIS, where pairs of hadrons — here p+p0 and p´p0 — are observed. We analyzed these dihadrons using detector data collected during Fall 2018 and Spring 2019, where longitudinally polarized electrons from the CEBAF accelerator were incident on a liquid hydrogen target. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained using Monte Carlo simulations to reduce the amount of iv false combinatorial background p0’s. When deployed on experimental data, the model in creases our dihadron statistics by up to five-fold compared to previous CLAS12 p0 analyses. This work reports the first measurements of beam spin asymmetries for p+p0 and p´p0 dihadron production in SIDIS. The measured asymmetries offer new insights to the spin-dependent structure and dynamics within the proton, as well as the spin-dependent properties of quark fragmentation. Non-zero twist-3 sin¿R amplitudes are observed, pro viding sensitivity to the subleading twist PDF e(x). The PDF e(x) encodes quark-gluon correlations within the proton — a property that is otherwise inaccessible at leading twist. Additionally, this work measured significant twist-2 modulations carried by sin(¿h ´ ¿R) and sin(2¿h ´2¿R), providing experimental access to the helicity dihadron fragmentation function (DiFF) GK 1 . Because there is no equivalent quark helicity-dependent FF in single pion SIDIS, the DiFF GK 1 offers a unique lens into novel spin-dependent fragmentation. For instance, the twist-2 modulations observed in this study are enhanced by vector mesons created during fragmentation — a behavior predicted by phenomenological models. This study broadens our understanding of dihadron fragmentation, revealing new details about the flavor and charge dependence of hadronization.

Matousek, Gregory [Duke Univ., Durham, NC (United ↗

Drell-Yan angular lepton distributions at small x from TMD factorization.

The Drell-Yan process is studied in the framework of TMD factorization in the Sudakov region s » Q 2 » \( {q}_{\perp}^2 \) corresponding to recent LHC experiments with Q 2 of order of mass of Z-boson and transverse momentum of DY pair ~ few tens GeV. The DY hadronic tensors are expressed in terms of quark and quark-gluon TMDs with \( \frac{1}{Q^2} \) and \( \frac{1}{N_c^2} \) accuracy. It is demonstrated that in the leading order in N c the higher-twist quark-quark-gluon TMDs reduce to leading-twist TMDs due to QCD equation of motion. The resulting hadronic tensors depend on two leading-twist TMDs: f 1 responsible for total DY cross section, and Boer-Mulders function \( {h}_1^{\perp } \) . The corresponding qualitative and semi-quantitative predictions seem to agree with LHC data on five angular coefficients A 0 – A 4 of DY pair production. The remaining three coefficients A 5 – A 7 are determined by quark-quark-gluon TMDs multiplied by extra \( \frac{1}{N_c} \) so they appear to be relatively small in accordance with LHC results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

$t\bar{t}b\bar{b}$ at the LHC: on the size of corrections and b -jet definitions

We report on the calculation of the next-to-leading order QCD corrections to the production of a $t\bar t$ pair in association with two heavy-flavour jets. We concentrate on the di-lepton $t\bar t$ decay channel at the LHC with $\sqrt{s}$ = 13 TeV. The computation is based on pp → $e^+ν_eμ^-\bar{ν}_μb\bar{b}b\bar{b}$ matrix elements and includes all resonant and non-resonant diagrams, interferences and off-shell effects of the top quark and the W gauge boson. As it is customary for such studies, results are presented in the form of inclusive and differential fiducial cross sections. We extensively investigate the dependence of our results upon variation of renormalisation and factorisation scales and parton distribution functions in the quest for an accurate estimate of the theoretical uncertainties. We additionally study the impact of the contributions induced by the bottom-quark parton density. Results presented here are particularly relevant for measurements of $t\bar{t}H(H → b\bar{b})$ and the determination of the Higgs coupling to the top quark. In addition, they might be used for precise measurements of the top-quark fiducial cross sections and to investigate top-quark decay modelling at the LHC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of the top quark pair production cross section in PbPb collisions at $\sqrt{s_{\mathrm{NN}}}=5.36$ TeV

The inclusive cross section for top quark pair ($\mathrm{t}\overline{\mathrm{t}}$) production in lead-lead (PbPb) collisions is reported for the first time at a center-of-mass energy per nucleon pair of 5.36 TeV. The analysis uses data corresponding to an integrated luminosity of 1.58 nb −1 collected by the CMS experiment at the CERN LHC in 2023. The $\mathrm{t}\overline{\mathrm{t}}$ production cross section, ${\sigma}_{\mathrm{t}\overline{\mathrm{t}}}={3.42}_{-0.51}^{+0.54}{\left(\mathrm{stat}\right)}_{-0.43}^{+0.50}\left(\mathrm{syst}\right)$ μb, is measured in dilepton final states using a fit to a multivariate discriminator that combines the decay electron and muon kinematic properties with the multiplicity of bottom quark jets. The result is consistent with perturbative quantum chromodynamics calculations at next-to-next-to-leading order (NNLO) accuracy employing several nuclear parton distribution functions. In addition, the Drell–Yan production cross section (σ DY ) for dilepton masses above 10 GeV and the ratio of $\mathrm{t}\overline{\mathrm{t}}$ to DY cross sections $\left({R}_{\mathrm{t}\overline{\mathrm{t}}/\mathrm{DY}}\right)$ are found to be compatible with the NNLO predictions. The observables ${\sigma}_{\mathrm{t}\overline{\mathrm{t}}}$, σ DY , and ${R}_{\mathrm{t}\overline{\mathrm{t}}/\mathrm{DY}}$ are measured separately for central and semicentral PbPb collisions to investigate for the first time the dependence of top quark production on the collision impact parameter.

Heavy Ion Experiments↗

α-cluster structure of 18 Ne

In this work, we investigate α-clustering in 18 Ne and compare it with what is known about clustering in the mirror nucleus 18 O. The excitation function for resonant elastic scattering of 14 O + α was measured in inverse kinematics using the TexAT active-target detector system, covering excitation energies from 8 to 17 MeV. The data were analyzed using a multi-channel R-matrix approach. Detailed spectroscopic information—including excitation energies, spin and parity assignments, and partial α and total widths—was extracted from the R-matrix analysis and compared with theoretical models and previously available data. Correspondence between levels in 18 O and 18 Ne is established. An extensive shell-model analysis of the 18 O– 18 Ne mirror pair shows very good agreement between theory and experiment, especially in the distribution of clustering strength. The comparison reveals that certain states, particularly at higher excitation energies, exhibit significantly stronger α-cluster character than predicted. This suggests that these states have a collective structure aligned with the corresponding α-reaction channel.

6 ≤ A ≤ 19↗

Extraction and validation of a new set of CMS pythia8 tunes from underlying-event measurements

New sets of CMS underlying-event parameters ("tunes") are presented for the PYTHIA8 event generator. These tunes use the NNPDF3.1 parton distribution functions (PDFs) at leading (LO), next-to-leading (NLO), or next-to-next-to-leading (NNLO) orders in perturbative quantum chromodynamics, and the strong coupling evolution at LO or NLO. Measurements of charged-particle multiplicity and transverse momentum densities at various hadron collision energies are fit simultaneously to determine the parameters of the tunes. Comparisons of the predictions of the new tunes are provided for observables sensitive to the event shapes at LEP, global underlying event, soft multiparton interactions, and double-parton scattering contributions. In addition, comparisons are made for observables measured in various specific processes, such as multijet, Drell-Yan, and top quark-antiquark pair production including jet substructure observables. The simulation of the underlying event provided by the new tunes is interfaced to a higher-order matrix-element calculation. For the first time, predictions from PYTHIA8 obtained with tunes based on NLO or NNLO PDFs are shown to reliably describe minimum-bias and underlying-event data with a similar level of agreement to predictions from tunes using LO PDF sets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Combination of inclusive top-quark pair production cross-section measurements using ATLAS and CMS data at $ \sqrt{s} $ = 7 and 8 TeV

A combination of measurements of the inclusive top-quark pair production cross-section performed by ATLAS and CMS in proton–proton collisions at centre-of-mass energies of 7 and 8 TeV at the LHC is presented. The cross-sections are obtained using top-quark pair decays with an opposite-charge electron–muon pair in the final state and with data corresponding to an integrated luminosity of about 5 fb$^{−1}$ at $ \sqrt{s} $ = 7 TeV and about 20 fb$^{−1}$ at $ \sqrt{s} $ = 8 TeV for each experiment. The combined cross-sections are determined to be 178.5 ± 4.7 pb at $ \sqrt{s} $ = 7 TeV and $ {243.3}_{-5.9}^{+6.0} $ pb at $ \sqrt{s} $ = 8 TeV with a correlation of 0.41, using a reference top-quark mass value of 172.5 GeV. The ratio of the combined cross-sections is determined to be R$_{8/7}$ = 1.363 ± 0.032. The combined measured cross-sections and their ratio agree well with theory calculations using several parton distribution function (PDF) sets. The values of the top-quark pole mass (with the strong coupling fixed at 0.118) and the strong coupling (with the top-quark pole mass fixed at 172.5 GeV) are extracted from the combined results by fitting a next-to-next-to-leading-order plus next-to-next-to-leading-log QCD prediction to the measurements. Using a version of the NNPDF3.1 PDF set containing no top-quark measurements, the results obtained are $ {m}_t^{\textrm{pole}}={173.4}_{-2.0}^{+1.8} $ GeV and $ {\alpha}_{\textrm{s}}\left({m}_Z\right)={0.1170}_{-0.0018}^{+0.0021} $.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Atomic-Scale Imaging of Condensed Counterions

Here, the functioning of a wide variety of charged macromolecules, from DNA to fuel cell membranes, is dependent on how the counterions surrounding them are arranged. In order to decrease Coulombic repulsion, some of the fixed charges on these molecules are neutralized by a fraction of the counterions-this phenomenon is called counterion condensation. The nature of counterion condensation can be only be inferred indirectly from traditional experiments such as X-ray scattering and modern experiments such as single molecule electrometry. The prevalent conclusion in the literature, based on both theory and experiment, is that the distribution of counterions is peaked right next to the macromolecule, i.e., condensation results in the formation of contact ion pairs. In this study, cryogenic electron microscopy (cryo-EM) was used to study the arrangement of condensed halide counterions near a positively charged polypeptoid nanofiber. The locations of both condensed and fixed charges were determined directly from atomic-scale images. Our experimentally determined counterion distributions were peaked at distances of about 5 Å away from the fixed positive charge, indicating the presence of a layer of water molecules between condensed ion pairs. We posit that this distribution is driven by the entropy of the condensed ions.

36 MATERIALS SCIENCE↗

Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in i-band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established redMaGiC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016)MNRAA40035-871110.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies i<4zphot+18 and has ∼30% wider redshift distributions but ∼3.5 times more galaxies than redMaGiC. Assuming a wCDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of Ωm, w, and σ8, and gains of 16% in σ8, 10% in Ωm, and 12% in w. Similarly, in ΛCDM, we find an improvement of 19% and 27% on σ8 and Ωm, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i<22.2, has ∼7 times more galaxies and ∼20% wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the EMCCD point-source response using x-rays

Electron Multiplying Charge Coupled Devices, EMCCD are used as x-ray detectors. The NSLS-II Soft Inelastic x-ray Scattering (SIX) beam line has two EMCCDs for x-ray detection in the spectrometer arm. The spectrometer with high resolving power disperses x-rays vertically. The x-ray vertical position on the sensor plane is related to its energy. This allows for very accurate x-ray energy measurements through x-ray coordinates. X-rays interact with silicon and create a number of electron-hole pairs proportional to the x-ray energy. Electrons drift and diffuse toward pixel gates and are collected there. The diffused electrons form a charge cloud distributed over several neighboring pixels. This charge sharing enables coordinate measurements with accuracy better than the pixel pitch. The charge distribution shape has to be taken into account to achieve ultimate accuracy in coordinate measurements. In this paper, we present a method of the charge distribution shape analysis and demonstrate its applications. The drift and diffusion of electrons ftom the point of generation to pixel gates results in the bell-shaped electron cloud usually approximated by Gaussian shape. The number of electrons collected under a pixel is proportional to the shape function integral. These electron packets get transferred to the sense node of the output amplifier. The transfer process could introduce distortions to the original charge distribution. For example, during transfers, electrons in the packet could be exposed to traps if they are present in the sensor. The trapping and later the release processes distort the apparent shape of the charge distribution. Therefore, deviations of the charge distribution shape ftom the originally symmetrical form can indicate the presence of trap centers in the sensor and can be used for sensor diagnostics.

47 OTHER INSTRUMENTATION↗

Generative Ensemble Regression: Learning Particle Dynamics from Observations of Ensembles with Physics-informed Deep Generative Models

Here, we propose a new method for inferring the governing stochastic ordinary differential equations (SODEs) by observing particle ensembles at discrete and sparse time instants, i.e., multiple “snapshots.” Particle coordinates at a single time instant, possibly noisy or truncated, are recorded in each snapshot but are unpaired across the snapshots. By training a physics-informed generative model that generates “fake” sample paths, we aim to fit the observed particle ensemble distributions with a curve in the probability measure space, which is induced from the inferred particle dynamics. We employ different metrics to quantify the differences between distributions, e.g., the sliced Wasserstein distances and the adversarial losses in generative adversarial networks. We refer to this method as generative “ensemble-regression” (GER), in analogy to the classic “point-regression,” where we infer the dynamics by performing regression in the Euclidean space. We illustrate the GER by learning the drift and diffusion terms of particle ensembles governed by SODEs with Brownian motions and Lévy processes up to 100 dimensions. We also discuss how to treat cases with noisy or truncated observations. Apart from systems consisting of independent particles, we also tackle nonlocal interacting particle systems with unknown interaction potential parameters by constructing a physics-informed loss function. Finally, we investigate scenarios of paired observations and discuss how to reduce the dimensionality in such cases by proving a convergence theorem that provides theoretical support.

97 MATHEMATICS AND COMPUTING↗

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗