Showing 1–50 of 171 results for author: Ball, D

ASKAP$-$EMU Discovery of "Raspberry": a new Galactic SNR Candidate G308.73+1.38

Authors: Sanja Lazarević, Miroslav D. Filipović, Bärbel S. Koribalski, Zachary J. Smeaton, Andrew M. Hopkins, Rami Z. E. Alsaberi, Velibor Velović, Brianna D. Ball, Roland Kothes, Denis Leahy, Adriano Ingallinera

Abstract: We report the ASKAP discovery of a new Galactic supernova remnant (SNR) candidate G308.73+1.38, which we name Raspberry. This new SNR candidate has an angular size of 20.7 arcmin $\times$ 16.7 arcmin, and we measure a total integrated flux of 407$\pm$50 mJy. We estimate Raspberry's most likely diameter of 10$-$30 pc which would place it at a distance of 3$-$5 kpc, on the near side of the Milky Way… ▽ More We report the ASKAP discovery of a new Galactic supernova remnant (SNR) candidate G308.73+1.38, which we name Raspberry. This new SNR candidate has an angular size of 20.7 arcmin $\times$ 16.7 arcmin, and we measure a total integrated flux of 407$\pm$50 mJy. We estimate Raspberry's most likely diameter of 10$-$30 pc which would place it at a distance of 3$-$5 kpc, on the near side of the Milky Way's Scutum$-$Centaurus Arm. We also find a Stokes$-$V point source close to the centre of Raspberry with a $\sim$5$σ$ significance. This point source may be the remaining compact source, a neutron star, or possibly a pulsar, formed during the initial supernova event. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 3 pages, 1 figure

Journal ref: Research Notes of the American Astronomical Society, 8, 107 (2024)

The Path to N$^3$LO Parton Distributions

Authors: The NNPDF Collaboration, Richard D. Ball, Andrea Barontini, Alessandro Candido, Stefano Carrazza, Juan Cruz-Martinez, Luigi Del Debbio, Stefano Forte, Tommaso Giani, Felix Hekhorn, Zahari Kassabov, Niccolò Laurenti, Giacomo Magni, Emanuele R. Nocera, Tanjona R. Rabemananjara, Juan Rojo, Christopher Schwan, Roy Stegeman, Maria Ubiali

Abstract: We extend the existing leading (LO), next-to-leading (NLO), and next-to-next-to-leading order (NNLO) NNPDF4.0 sets of parton distribution functions (PDFs) to approximate next-to-next-to-next-to-leading order (aN$^3$LO). We construct an approximation to the N$^3$LO splitting functions that includes all available partial information from both fixed-order computations and from small and large $x$ res… ▽ More We extend the existing leading (LO), next-to-leading (NLO), and next-to-next-to-leading order (NNLO) NNPDF4.0 sets of parton distribution functions (PDFs) to approximate next-to-next-to-next-to-leading order (aN$^3$LO). We construct an approximation to the N$^3$LO splitting functions that includes all available partial information from both fixed-order computations and from small and large $x$ resummation, and estimate the uncertainty on this approximation by varying the set of basis functions used to construct the approximation. We include known N$^3$LO corrections to deep-inelastic scattering structure functions and extend the FONLL general-mass scheme to $\mathcal{O}\left( α_s^3\right)$ accuracy. We determine a set of aN$^3$LO PDFs by accounting both for the uncertainty on splitting functions due to the incomplete knowledge of N$^3$LO terms, and to the uncertainty related to missing higher corrections (MHOU), estimated by scale variation, through a theory covariance matrix formalism. We assess the perturbative stability of the resulting PDFs, we study the impact of MHOUs on them, and we compare our results to the aN$^3$LO PDFs from the MSHT group. We examine the phenomenological impact of aN$^3$LO corrections on parton luminosities at the LHC, and give a first assessment of the impact of aN$^3$LO PDFs on the Higgs and Drell-Yan total production cross-sections. We find that the aN$^3$LO NNPDF4.0 PDFs are consistent within uncertainties with their NNLO counterparts, that they improve the description of the global dataset and the perturbative convergence of Higgs and Drell-Yan cross-sections, and that MHOUs on PDFs decrease substantially with the increase of perturbative order. △ Less

Submitted 5 July, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: 63 pages, 41 figures. Final version, as published in EPJC. Various clarifications added. Fig 2.9 and Tab. 4.6, 5.1 and 5.2 added. In particular, comparison to approximate YR4 estimate of PDF uncertainties added (Tab 4.1, 4.2 and associate discussion in sect. 5)

Report number: Nikhef-2023-020, TIF-UNIMI-2023-23, Edinburgh 2023/29

Deep Neural Network Initialization with Sparsity Inducing Activations

Authors: Ilan Price, Nicholas Daultry Ball, Samuel C. H. Lam, Adam C. Jones, Jared Tanner

Abstract: Inducing and leveraging sparse activations during training and inference is a promising avenue for improving the computational efficiency of deep networks, which is increasingly important as network sizes continue to grow and their application becomes more widespread. Here we use the large width Gaussian process limit to analyze the behaviour, at random initialization, of nonlinear activations tha… ▽ More Inducing and leveraging sparse activations during training and inference is a promising avenue for improving the computational efficiency of deep networks, which is increasingly important as network sizes continue to grow and their application becomes more widespread. Here we use the large width Gaussian process limit to analyze the behaviour, at random initialization, of nonlinear activations that induce sparsity in the hidden outputs. A previously unreported form of training instability is proven for arguably two of the most natural candidates for hidden layer sparsification; those being a shifted ReLU ($φ(x)=\max(0, x-τ)$ for $τ\ge 0$) and soft thresholding ($φ(x)=0$ for $|x|\leτ$ and $x-\text{sign}(x)τ$ for $|x|>τ$). We show that this instability is overcome by clipping the nonlinear activation magnitude, at a level prescribed by the shape of the associated Gaussian process variance map. Numerical experiments verify the theory and show that the proposed magnitude clipped sparsifying activations can be trained with training and test fractional sparsity as high as 85\% while retaining close to full accuracy. △ Less

Submitted 25 February, 2024; originally announced February 2024.

Comments: Published in the International Conference on Learning Representations (ICLR) 2024

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Determination of the theory uncertainties from missing higher orders on NNLO parton distributions with percent accuracy

Authors: The NNPDF Collaboration, Richard D. Ball, Andrea Barontini, Alessandro Candido, Stefano Carrazza, Juan Cruz-Martinez, Luigi Del Debbio, Stefano Forte, Tommaso Giani, Felix Hekhorn, Zahari Kassabov, Niccolò Laurenti, Giacomo Magni, Emanuele R. Nocera, Tanjona R. Rabemananjara, Juan Rojo, Christopher Schwan, Roy Stegeman, Maria Ubiali

Abstract: We include uncertainties due to missing higher order corrections to QCD computations (MHOU) used in the determination of parton distributions (PDFs) in the recent NNPDF4.0 set of PDFs. We use our previously published methodology, based on the treatment of MHOUs and their full correlations through a theory covariance matrix determined by scale variation, now fully incorporated in the new NNPDF theo… ▽ More We include uncertainties due to missing higher order corrections to QCD computations (MHOU) used in the determination of parton distributions (PDFs) in the recent NNPDF4.0 set of PDFs. We use our previously published methodology, based on the treatment of MHOUs and their full correlations through a theory covariance matrix determined by scale variation, now fully incorporated in the new NNPDF theory pipeline. We assess the impact of the inclusion of MHOUs on the NNPDF4.0 central values and uncertainties, and specifically show that they lead to improved consistency of the PDF determination with an ensuing moderate reduction of PDF uncertainties at NNLO. △ Less

Submitted 4 July, 2024; v1 submitted 18 January, 2024; originally announced January 2024.

Comments: 23 pages, 6 figures. Final version, as published in EPJC. Several small improvements; expanded discussion of the phi indicator in sect 4.2

Report number: TIF-UNIMI-2023-22, Edinburgh 2023/34, CERN-TH-2024-009

Photons in the proton: implications for the LHC

Authors: The NNPDF Collaboration, Richard D. Ball, Andrea Barontini, Alessandro Candido, Stefano Carrazza, Juan Cruz-Martinez, Luigi Del Debbio, Stefano Forte, Tommaso Giani, Felix Hekhorn, Zahari Kassabov, Niccolò Laurenti, Giacomo Magni, Emanuele R. Nocera, Tanjona R. Rabemananjara, Juan Rojo, Christopher Schwan, Roy Stegeman, Maria Ubiali

Abstract: We construct a set of parton distribution functions (PDFs), based on the recent NNPDF4.0 PDF set, that also include a photon PDF. The photon PDF is constructed using the LuxQED formalism, while QED evolution accounting for O(alpha), O(alpha alphas) and O(alpha^2) corrections is implemented and benchmarked by means of the EKO code. We investigate the impact of QED effects on NNPDF4.0, and compare o… ▽ More We construct a set of parton distribution functions (PDFs), based on the recent NNPDF4.0 PDF set, that also include a photon PDF. The photon PDF is constructed using the LuxQED formalism, while QED evolution accounting for O(alpha), O(alpha alphas) and O(alpha^2) corrections is implemented and benchmarked by means of the EKO code. We investigate the impact of QED effects on NNPDF4.0, and compare our results both to our previous NNPDF3.1QED PDF set and to other recent PDF sets that include the photon. We assess the impact of photon-initiated processes and electroweak corrections on a variety of representative LHC processes, and find that they can reach the 5% level in vector boson pair production at large invariant mass. △ Less

Submitted 4 July, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: 32 pages, 19 figures. Final; version, as published on EPJC. Pheno discussion in Sect 4.2 improved, a few typos corrected

Report number: TIF-UNIMI-2023-17, Edinburgh 2023/19, CERN-TH-2023-159

Grain-level effects on in-situ deformation-induced phase transformations in a complex-phase steel using 3DXRD and EBSD

Authors: James A. D. Ball, Claire Davis, Carl Slater, Himanshu Vashishtha, Mohammed Said, Louis Hébrard, Florian Steinhilber, Jonathan P. Wright, Thomas Connolley, Stefan Michalik, David M. Collins

Abstract: A novel complex-phase steel alloy is conceived with a deliberately unstable austenite, $γ$, phase that enables the deformation-induced martensitic transformations (DIMT) to be explored at low levels of plastic strain. The DIMT was thus explored, in-situ and non-destructively, using both far-field Three-Dimensional X-Ray Diffraction (3DXRD) and Electron Back-Scatter Diffraction (EBSD). Substantial… ▽ More A novel complex-phase steel alloy is conceived with a deliberately unstable austenite, $γ$, phase that enables the deformation-induced martensitic transformations (DIMT) to be explored at low levels of plastic strain. The DIMT was thus explored, in-situ and non-destructively, using both far-field Three-Dimensional X-Ray Diffraction (3DXRD) and Electron Back-Scatter Diffraction (EBSD). Substantial $α'$ martensite formation was observed under 10% applied strain with EBSD, and many $\varepsilon$ grain formation events were captured with 3DXRD, indicative of the indirect transformation of martensite via the reaction $γ\rightarrow \varepsilon \rightarrow α'$. Using $\varepsilon$ grain formation as a direct measurement of $γ$ grain stability, the influence of several microstructural properties, such as grain size, orientation and neighbourhood configuration, on $γ$ stability have been identified. Larger $γ$ grains were found to be less stable than smaller grains. Any $γ$ grains oriented with {100} parallel to the loading direction preferentially transformed with lower stresses. Parent $\varepsilon$-forming $γ$ grains possessed a neighbourhood with increased ferritic/martensitic volume fraction. This finding shows, unambiguously, that $α$/$α'$ promotes $\varepsilon$ formation in neighbouring grains. The minimum strain work criterion model for $\varepsilon$ variant prediction was also evaluated, which worked well for most grains. However, $\varepsilon$-forming grains with a lower stress were less well predicted by the model, indicating crystal-level behaviour must be considered for accurate $\varepsilon$ formation. The findings from this work are considered key for the future design of alloys where the deformation response can be controlled by tailoring microstructure and local or macroscopic crystal orientations. △ Less

Submitted 7 November, 2023; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: 20 pages, 3 supplementary pages, 14 figures, 3 supplementary figures. Preprint submitted to Acta Materialia. Updated to correct Figure 14

Journal ref: Acta Materialia 265 (2024) 119608

The intrinsic charm quark valence distribution of the proton

Authors: Richard D. Ball, Alessandro Candido, Juan Cruz-Martinez, Stefano Forte, Tommaso Giani, Felix Hekhorn, Giacomo Magni, Emanuele R. Nocera, Juan Rojo, Roy Stegeman

Abstract: We provide a first quantitative indication that the wave function of the proton contains unequal distributions of charm quarks and antiquarks, i.e. a nonvanishing intrinsic valence charm distribution. A significant nonvanishing valence component cannot be perturbatively generated, hence our results reinforce previous evidence that the proton contains an intrinsic (i.e., not radiatively generated)… ▽ More We provide a first quantitative indication that the wave function of the proton contains unequal distributions of charm quarks and antiquarks, i.e. a nonvanishing intrinsic valence charm distribution. A significant nonvanishing valence component cannot be perturbatively generated, hence our results reinforce previous evidence that the proton contains an intrinsic (i.e., not radiatively generated) charm quark component. We establish our result through a determination of the parton distribution functions (PDFs) of charm quarks and antiquarks in the proton. We propose two novel experimental probes of this intrinsic charm valence component: D-meson asymmetries in Z+c-jet production at the LHCb experiment, and flavor-tagged structure functions at the Electron-Ion Collider. △ Less

Submitted 15 March, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 11 pages, 11 figures. PDF files available for download in LHAPDF format from https://nnpdf.mi.infn.it/nnpdf4-0-charm-asymmetry-study/ Final version, to be published in PRD (letters). Figures 2, 6 and 7 added. Fit without imposing the charm valence sum rule added. Various small clarifications

Report number: Edinburgh 2023/22, TIF-UNIMI-2023-28, Nikhef 2023-012, CERN-TH-2023-196

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

A Catalogue of Radio Supernova Remnants and Candidate Supernova Remnants in the EMU/POSSUM Galactic Pilot Field

Authors: Brianna D. Ball, Roland Kothes, Erik Rosolowsky, Jennifer West, Werner Becker, Miroslav D. Filipović, B. M. Gaensler, Andrew M. Hopkins, Bärbel Koribalski, Tom Landecker, Denis Leahy, Joshua Marvil, Xiaohui Sun, Filomena Bufano, Ettore Carretti, Adriano Ingallinera, Cameron L. Van Eck, Tony Willis

Abstract: We use data from the pilot observations of the EMU/POSSUM surveys to study the "missing supernova remnant (SNR) problem", the discrepancy between the number of Galactic SNRs that have been observed and the number that are estimated to exist. The Evolutionary Map of the Universe (EMU) and the Polarization Sky Survey of the Universe's Magnetism (POSSUM) are radio sky surveys that are conducted using… ▽ More We use data from the pilot observations of the EMU/POSSUM surveys to study the "missing supernova remnant (SNR) problem", the discrepancy between the number of Galactic SNRs that have been observed and the number that are estimated to exist. The Evolutionary Map of the Universe (EMU) and the Polarization Sky Survey of the Universe's Magnetism (POSSUM) are radio sky surveys that are conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). We report on the properties of 7 known SNRs in the joint Galactic pilot field, with an approximate longitude and latitude of 323$^\circ\leq$ l $\leq$ 330$^\circ$ and -4$^\circ\leq$ b $\leq$ 2$^\circ$ respectively, and identify 21 SNR candidates. Of these, 4 have been previously identified as SNR candidates, 3 were previously listed as a single SNR, 13 have not been previously studied, and 1 has been studied in the infrared. These are the first discoveries of Galactic SNR candidates with EMU/POSSUM and, if confirmed, they will increase the SNR density in this field by a factor of 4. By comparing our SNR candidates to the known Galactic SNR population, we demonstrate that many of these sources were likely missed in previous surveys due to their small angular size and/or low surface brightness. We suspect that there are SNRs in this field that remain undetected due to limitations set by the local background and confusion with other radio sources. The results of this paper demonstrate the potential of the full EMU/POSSUM surveys to uncover more of the missing Galactic SNR population. △ Less

Submitted 4 July, 2023; originally announced July 2023.

Registration between DCT and EBSD datasets for multiphase microstructures

Authors: James A. D. Ball, Jette Oddershede, Claire Davis, Carl Slater, Mohammed Said, Himanshu Vashishtha, Stefan Michalik, David M. Collins

Abstract: The ability to characterise the three-dimensional microstructure of multiphase materials is essential for understanding the interaction between phases and associated materials properties. Here, laboratory-based diffraction-contrast tomography (DCT), a recently-established materials characterization technique that can determine grain phases, morphologies, positions and orientations in a voxel-based… ▽ More The ability to characterise the three-dimensional microstructure of multiphase materials is essential for understanding the interaction between phases and associated materials properties. Here, laboratory-based diffraction-contrast tomography (DCT), a recently-established materials characterization technique that can determine grain phases, morphologies, positions and orientations in a voxel-based reconstruction method, was used to map part of a dual-phase steel alloy sample. To assess the resulting microstructures that were produced by the DCT technique, an EBSD map was collected within the same sample volume. To identify the 2D slice of the 3D DCT reconstruction that best corresponded to the EBSD map, a novel registration technique based solely on grain-averaged orientations was developed -- this registration technique requires very little a priori knowledge of dataset alignment and can be extended to other techniques that only recover grain-averaged orientation data such as far-field 3D X-ray diffraction microscopy. Once the corresponding 2D slice was identified in the DCT dataset, comparisons of phase balance, grain size, shape and texture were performed between DCT and EBSD techniques. More complicated aspects of the microstructural morphology such as grain boundary shape and grains less than a critical size were poorly reproduced by the DCT reconstruction, primarily due to the difference in resolutions of the technique compared with EBSD. However, lab-based DCT is shown to accurately determine the centre-of-mass position, orientation, and size of the large grains for each phase present, austenite and martensitic ferrite. The results reveals a complex ferrite grain network of similar crystal orientations that are absent from the EBSD dataset. Such detail demonstrates that lab-based DCT, as a technique, shows great promise in the field of multi-phase material characterization. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 15 pages, 11 figures. Preprint submitted to Materials Characterization

Journal ref: Materials Characterization, October 2023, Volume 204, Page 113228

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Per-grain and neighbourhood stress interactions during deformation of a ferritic steel obtained using three-dimensional X-ray diffraction

Authors: James A. D. Ball, Anna Kareer, Oxana V. Magdysyuk, Stefan Michalik, Thomas Connolley, David M. Collins

Abstract: Three-dimensional X-ray diffraction (3DXRD) has been used to measure, in-situ, the evolution of $\sim 1800$ grains in a single phase low carbon ferritic steel sample during uniaxial deformation. The distribution of initial residual grain stresses in the material was observed to prevail as plasticity builds, though became less pronounced, and therefore less influential as strain increased. The init… ▽ More Three-dimensional X-ray diffraction (3DXRD) has been used to measure, in-situ, the evolution of $\sim 1800$ grains in a single phase low carbon ferritic steel sample during uniaxial deformation. The distribution of initial residual grain stresses in the material was observed to prevail as plasticity builds, though became less pronounced, and therefore less influential as strain increased. The initial Schmid factor of a grain was found to be strongly correlated to the intergranular stress change and the range of stresses that are permissible; a grain well aligned for easy slip is more likely to exhibit a range of stresses than those orientated poorly for dislocation motion. The orientation path of a grain, however, is not only dependent on its initial orientation, but hypothesised to be influenced by its stress state and the stress state of its grain environment. A grain neighbourhood effect is observed: the Schmid factor of serial adjoining grains influences the stress state of a grain of interest, whereas parallel neighbours are much less influential. This phenomenon is strongest at low plastic strains only, with the effect diminishing as plasticity builds. The influence of initial residual stresses becomes less evident, and grains rotate to eliminate any orientation dependent load shedding. The ability of the BCC ferrite to exhaust such neighbourhood interactions, which would otherwise be detrimental in crystal structures with lower symmetric and fewer slip systems, is considered key to the high ductility possessed by these materials. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 17 pages (+2 supplementary pages), 13 figures (+3 supplementary figures). Preprint submitted to Acta Materialia

Journal ref: Commun Mater 5, 27 (2024)

Response to "Parton distributions need representative sampling"

Authors: The NNPDF Collaboration, Richard D. Ball, Juan Cruz-Martinez, Luigi Del Debbio, Stefano Forte, Zahari Kassabov, Emanuele R. Nocera, Juan Rojo, Roy Stegeman, Maria Ubiali

Abstract: We respond to the criticism raised by Courtoy et al., in which the faithfulness of the NNPDF4.0 sampling is questioned and an under-estimate of the NNPDF4.0 PDF uncertainties is implied. We list, correct, and clarify in detail a number of inaccurate or misleading claims that are made in this Reference. Specifically, we explain and explicitly demonstrate why the central value of the PDF distributio… ▽ More We respond to the criticism raised by Courtoy et al., in which the faithfulness of the NNPDF4.0 sampling is questioned and an under-estimate of the NNPDF4.0 PDF uncertainties is implied. We list, correct, and clarify in detail a number of inaccurate or misleading claims that are made in this Reference. Specifically, we explain and explicitly demonstrate why the central value of the PDF distribution does not generally coincide with the absolute minimum of the $χ^{2}$ to the data. We examine some PDFs that have been constructed in the above study and claimed to be "good solutions": we show that similar PDFs are found with the NNPDF methodology, but with very low probability. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 18 pages, 12 figures

Report number: Edinburgh 2022/19, TIF-UNIMI-2022-21

Parton Distributions and New Physics Searches: the Drell-Yan Forward-Backward Asymmetry as a Case Study

Authors: Richard D. Ball, Alessandro Candido, Stefano Forte, Felix Hekhorn, Emanuele R. Nocera, Juan Rojo, Christopher Schwan

Abstract: We discuss the sensitivity of theoretical predictions of observables used in searches for new physics to parton distributions (PDFs) at large momentum fraction $x$. Specifically, we consider the neutral-current Drell-Yan production of gauge bosons with invariant masses in the TeV range, for which the forward-backward asymmetry of charged leptons from the decay of the gauge boson in its rest frame… ▽ More We discuss the sensitivity of theoretical predictions of observables used in searches for new physics to parton distributions (PDFs) at large momentum fraction $x$. Specifically, we consider the neutral-current Drell-Yan production of gauge bosons with invariant masses in the TeV range, for which the forward-backward asymmetry of charged leptons from the decay of the gauge boson in its rest frame is a traditional probe of new physics. We show that the qualitative behaviour of the asymmetry depends strongly on the assumptions made in determining the underlying PDFs. We discuss and compare the large-$x$ behaviour of various different PDF sets, and find that they differ significantly. Consequently, the shape of the asymmetry observed at lower dilepton invariant masses, where all PDF sets are in reasonable agreement because of the presence of experimental constraints, is not necessarily reproduced at large masses where the PDFs are mostly unconstrained by data. It follows that the shape of the asymmetry at high masses may depend on assumptions made in the PDF parametrization, and thus deviations from the traditionally expected behaviour cannot be taken as a reliable indication of new physics. We demonstrate that forward-backward asymmetry measurements could help in constraining PDFs at large $x$ and discuss the accuracy that would be required to disentangle the effects of new physics from uncertainties in the PDFs in this region. △ Less

Submitted 21 December, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 27 pages, 20 figures. v2: version accepted for publication in EPJC

Report number: Edinburgh 2022/14, Nikhef-2022-013, TIF-UNIMI-2022-13

Evidence for intrinsic charm quarks in the proton

Authors: Richard D. Ball, Alessandro Candido, Juan Cruz-Martinez, Stefano Forte, Tommaso Giani, Felix Hekhorn, Kirill Kudashkin, Giacomo Magni, Juan Rojo

Abstract: The theory of the strong force, quantum chromodynamics, describes the proton in terms of quarks and gluons. The proton is a state of two up quarks and one down quark bound by gluons, but quantum theory predicts that in addition there is an infinite number of quark-antiquark pairs. Both light and heavy quarks, whose mass is respectively smaller or bigger than the mass of the proton, are revealed in… ▽ More The theory of the strong force, quantum chromodynamics, describes the proton in terms of quarks and gluons. The proton is a state of two up quarks and one down quark bound by gluons, but quantum theory predicts that in addition there is an infinite number of quark-antiquark pairs. Both light and heavy quarks, whose mass is respectively smaller or bigger than the mass of the proton, are revealed inside the proton in high-energy collisions. However, it is unclear whether heavy quarks also exist as a part of the proton wavefunction, which is determined by non-perturbative dynamics and accordingly unknown: so-called intrinsic heavy quarks. It has been argued for a long time that the proton could have a sizable intrinsic component of the lightest heavy quark, the charm quark. Innumerable efforts to establish intrinsic charm in the proton have remained inconclusive. Here we provide evidence for intrinsic charm by exploiting a high-precision determination of the quark-gluon content of the nucleon based on machine learning and a large experimental dataset. We disentangle the intrinsic charm component from charm-anticharm pairs arising from high-energy radiation. We establish the existence of intrinsic charm at the 3-standard-deviation level, with a momentum distribution in remarkable agreement with model predictions. We confirm these findings by comparing to very recent data on Z-boson production with charm jets from the LHCb experiment. △ Less

Submitted 15 November, 2022; v1 submitted 17 August, 2022; originally announced August 2022.

Comments: 10 pages, 4 figures + Supplementary Information, published version available in https://www.nature.com/articles/s41586-022-04998-2. v2: updated comparison with the MBC model in Fig. 1 right (thanks to Tim Hobbs for the suggestion)

Report number: Nikhef 2021-032, Edinburgh 2021/28, TIF-UNIMI-2021-21

Journal ref: Nature 608, 483-487 (2022)

Snowmass 2021 whitepaper: Proton structure at the precision frontier

Authors: S. Amoroso, A. Apyan, N. Armesto, R. D. Ball, V. Bertone, C. Bissolotti, J. Bluemlein, R. Boughezal, G. Bozzi, D. Britzger, A. Buckley, A. Candido, S. Carrazza, F. G. Celiberto, S. Cerci, G. Chachamis, A. M. Cooper-Sarkar, A. Courtoy, T. Cridge, J. M. Cruz-Martinez, F. Giuli, M. Guzzi, C. Gwenlan, L. A. Harland-Lang, F. Hekhorn , et al. (32 additional authors not shown)

Abstract: An overwhelming number of theoretical predictions for hadron colliders require parton distribution functions (PDFs), which are an important ingredient of theory infrastructure for the next generation of high-energy experiments. This whitepaper summarizes the status and future prospects for determination of high-precision PDFs applicable in a wide range of energies and experiments, in particular in… ▽ More An overwhelming number of theoretical predictions for hadron colliders require parton distribution functions (PDFs), which are an important ingredient of theory infrastructure for the next generation of high-energy experiments. This whitepaper summarizes the status and future prospects for determination of high-precision PDFs applicable in a wide range of energies and experiments, in particular in precision tests of the Standard Model and in new physics searches at the high-luminosity Large Hadron Collider and Electron-Ion Collider. We discuss the envisioned advancements in experimental measurements, QCD theory, global analysis methodology, and computing that are necessary to bring unpolarized PDFs in the nucleon to the N2LO and N3LO accuracy in the QCD coupling strength. Special attention is given to the new tasks that emerge in the era of the precision PDF analysis, such as those focusing on the robust control of systematic factors both in experimental measurements and theoretical computations. Various synergies between experimental and theoretical studies of the hadron structure are explored, including opportunities for studying PDFs for nuclear and meson targets, PDFs with electroweak contributions or dependence on the transverse momentum, for incisive comparisons between phenomenological models for the PDFs and computations on discrete lattice, and for cross-fertilization with machine learning/AI approaches. [Submitted to the US Community Study on the Future of Particle Physics (Snowmass 2021).] △ Less

Submitted 5 April, 2023; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: 83 pages, 27 figures, contribution to Snowmass 2021; v.3: journal version

Report number: Edinburgh 2022/08, FERMILAB-PUB-22-222-QIS-SCD-T, MPP-2022-32, SLAC-PUB-17652, SMU-HEP-22-02, TIF-UNIMI-2022-6

Journal ref: Acta Phys.Polon.B 53 (2022) 12, A1

The PDF4LHC21 combination of global PDF fits for the LHC Run III

Authors: Richard D. Ball, Jon Butterworth, Amanda M. Cooper-Sarkar, Aurore Courtoy, Thomas Cridge, Albert De Roeck, Joel Feltesse, Stefano Forte, Francesco Giuli, Claire Gwenlan, Lucian A. Harland-Lang, T. J. Hobbs, Tie-Jiun Hou, Joey Huston, Ronan McNulty, Pavel M. Nadolsky, Emanuele R. Nocera, Tanjona R. Rabemananjara, Juan Rojo, Robert S. Thorne, Keping Xie, C. -P. Yuan

Abstract: A precise knowledge of the quark and gluon structure of the proton, encoded by the parton distribution functions (PDFs), is of paramount importance for the interpretation of high-energy processes at present and future lepton-hadron and hadron-hadron colliders. Motivated by recent progress in the PDF determinations carried out by the CT, MSHT, and NNPDF groups, we present an updated combination of… ▽ More A precise knowledge of the quark and gluon structure of the proton, encoded by the parton distribution functions (PDFs), is of paramount importance for the interpretation of high-energy processes at present and future lepton-hadron and hadron-hadron colliders. Motivated by recent progress in the PDF determinations carried out by the CT, MSHT, and NNPDF groups, we present an updated combination of global PDF fits: PDF4LHC21. It is based on the Monte Carlo combination of the CT18, MSHT20, and NNPDF3.1 sets followed by either its Hessian reduction or its replica compression. Extensive benchmark studies are carried out in order to disentangle the origin of the differences between the three global PDF sets. In particular, dedicated fits based on almost identical theory settings and input datasets are performed by the three groups, highlighting the role played by the respective fitting methodologies. We compare the new PDF4LHC21 combination with its predecessor, PDF4LHC15, demonstrating their good overall consistency and a modest reduction of PDF uncertainties for key LHC processes such as electroweak gauge boson production and Higgs boson production in gluon fusion. We study the phenomenological implications of PDF4LHC21 for a representative selection of inclusive, fiducial, and differential cross sections at the LHC. The PDF4LHC21 combination is made available via the LHAPDF library and provides a robust, user-friendly, and efficient method to estimate the PDF uncertainties associated to theoretical calculations for the upcoming Run III of the LHC and beyond. △ Less

Submitted 23 March, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: 87 pages, 44 figures, 11 tables

Report number: Edinburgh 2021/31, FERMILAB-PUB-22-121-QIS-SCD-T, MSUHEP-22-010, Nikhef 2021-033, SMU-HEP-22-01

Implementing and evaluating far-field 3D X-ray diffraction at the I12 JEEP beamline, Diamond Light Source

Authors: James A. D. Ball, Anna Kareer, Oxana V. Magdysyuk, Stefan Michalik, Anastasia Vrettou, Neal Parkes, Thomas Connolley, David M. Collins

Abstract: Three-dimensional X-ray diffraction (3DXRD) is shown to be feasible at the I12 Joint Engineering, Environmental and Processing (JEEP) beamline of Diamond Light Source. As a demonstration, a microstructually simple low-carbon ferritic steel was studied in a highly textured and annealed state. A processing pipeline suited to this beamline was created, using software already established in the 3DXRD… ▽ More Three-dimensional X-ray diffraction (3DXRD) is shown to be feasible at the I12 Joint Engineering, Environmental and Processing (JEEP) beamline of Diamond Light Source. As a demonstration, a microstructually simple low-carbon ferritic steel was studied in a highly textured and annealed state. A processing pipeline suited to this beamline was created, using software already established in the 3DXRD user community, enabling grain centre-of-mass positions, orientations and strain tensor elements to be determined. Orientations, with texture measurements independently validated from electron backscatter diffraction (EBSD) data, possessed a $\sim 0.1^\circ$ uncertainty, comparable with other 3DXRD instruments. The spatial resolution was limited by the far-field detector pixel size; the average of the grain centre of mass position errors was determined as $\pm \sim 80 μ$m. An average per-grain error of $\sim 1 \times 10^{-3}$ for the elastic strains was also measured; this could be reduced in future experiments by improving sample preparation, data collection and analysis techniques. Application of 3DXRD onto I12 shows great potential, where its implementation is highly desirable due the flexible, open architecture of the beamline. User-owned or designed sample environments can be used, thus 3DXRD could be applied to previously unexplored scientific areas. △ Less

Submitted 20 February, 2023; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 13 pages, 10 figures. Published in Journal of Synchrotron Radiation

Journal ref: Journal of Synchrotron Radiation, July 2022, Volume 29, Pages 1043-1053

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

Correlation and Combination of Sets of Parton Distributions

Authors: Richard D. Ball, Stefano Forte, Roy Stegeman

Abstract: We study the correlation between different sets of parton distributions (PDFs). Specifically, viewing different PDF sets as distinct determinations, generally correlated, of the same underlying physical quantity, we examine the extent to which the correlation between them is due to the underlying data. We do this both for pairs of PDF sets determined using a given fixed methodology, and between se… ▽ More We study the correlation between different sets of parton distributions (PDFs). Specifically, viewing different PDF sets as distinct determinations, generally correlated, of the same underlying physical quantity, we examine the extent to which the correlation between them is due to the underlying data. We do this both for pairs of PDF sets determined using a given fixed methodology, and between sets determined using different methodologies. We show that correlations have a sizable component that is not due to the underlying data, because the data do not determine the PDF uniquely. We show that the data-driven correlations can be used to assess the efficiency of methodologies used for PDF determination. We also show that the use of data-driven correlations for the combination of different PDFs into a joint set can lead to inconsistent results, and thus that the statistical combination used in constructing the widely used PDF4LHC15 PDF set remains the most reliable method. △ Less

Submitted 8 December, 2021; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: 21 pages, 9 figures. Final version published in EPJC: clarifications on data replicas and functional uncertainties added, typos corrected, references updated

Report number: Edinburgh 2021/19, TIF-UNIMI-2021-5

An open-source machine learning framework for global analyses of parton distributions

Authors: Richard D. Ball, Stefano Carrazza, Juan Cruz-Martinez, Luigi Del Debbio, Stefano Forte, Tommaso Giani, Shayan Iranipour, Zahari Kassabov, Jose I. Latorre, Emanuele R. Nocera, Rosalyn L. Pearson, Juan Rojo, Roy Stegeman, Christopher Schwan, Maria Ubiali, Cameron Voisey, Michael Wilson

Abstract: We present the software framework underlying the NNPDF4.0 global determination of parton distribution functions (PDFs). The code is released under an open source licence and is accompanied by extensive documentation and examples. The code base is composed by a PDF fitting package, tools to handle experimental data and to efficiently compare it to theoretical predictions, and a versatile analysis f… ▽ More We present the software framework underlying the NNPDF4.0 global determination of parton distribution functions (PDFs). The code is released under an open source licence and is accompanied by extensive documentation and examples. The code base is composed by a PDF fitting package, tools to handle experimental data and to efficiently compare it to theoretical predictions, and a versatile analysis framework. In addition to ensuring the reproducibility of the NNPDF4.0 (and subsequent) determination, the public release of the NNPDF fitting framework enables a number of phenomenological applications and the production of PDF fits under user-defined data and theory assumptions. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: 14 pages, 4 figures. A companion paper describes the NNPDF4.0 PDFs

Report number: Edinburgh 2021/13, Nikhef-2021-020, TIF-UNIMI-2021-12

The Path to Proton Structure at One-Percent Accuracy

Authors: Richard D. Ball, Stefano Carrazza, Juan Cruz-Martinez, Luigi Del Debbio, Stefano Forte, Tommaso Giani, Shayan Iranipour, Zahari Kassabov, Jose I. Latorre, Emanuele R. Nocera, Rosalyn L. Pearson, Juan Rojo, Roy Stegeman, Christopher Schwan, Maria Ubiali, Cameron Voisey, Michael Wilson

Abstract: We present a new set of parton distribution functions (PDFs) based on a fully global dataset and machine learning techniques: NNPDF4.0. We expand the NNPDF3.1 determination with 44 new datasets, mostly from the LHC. We derive a novel methodology through hyperparameter optimisation, leading to an efficient fitting algorithm built upon stochastic gradient descent. We use NNLO QCD calculations and ac… ▽ More We present a new set of parton distribution functions (PDFs) based on a fully global dataset and machine learning techniques: NNPDF4.0. We expand the NNPDF3.1 determination with 44 new datasets, mostly from the LHC. We derive a novel methodology through hyperparameter optimisation, leading to an efficient fitting algorithm built upon stochastic gradient descent. We use NNLO QCD calculations and account for NLO electroweak corrections and nuclear uncertainties. Theoretical improvements in the PDF description include a systematic implementation of positivity constraints and integrability of sum rules. We validate our methodology by means of closure tests and "future tests" (i.e. tests of backward and forward data compatibility), and assess its stability, specifically upon changes of PDF parametrization basis. We study the internal compatibility of our dataset, and investigate the dependence of results both upon the choice of input dataset and of fitting methodology. We perform a first study of the phenomenological implications of NNPDF4.0 on representative LHC processes. The software framework used to produce NNPDF4.0 is made available as an open-source package together with documentation and examples. △ Less

Submitted 31 May, 2022; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 139 pages, more than 60 figures. A companion paper describes the open-source NNPDF code. Figs 9.1 and 9.3 replaced after correction of a plotting bug (thanks to J.Huston and G.Salam for spotting it). A few more typos corrected and a reference updated. Final published version. References updated, several typos corrected, 72 issues raised by the referee addressed

Next generation proton PDFs with deuteron and nuclear uncertainties

Authors: Rosalyn L. Pearson, Richard D. Ball, Emanuele R. Nocera

Abstract: As data become more precise, estimating theoretical uncertainties in global PDF determinations is likely to become increasingly necessary to obtain correspondingly precise PDFs. Here we present a next generation of global proton PDFs (NNPDF4.0) that include theoretical uncertainties due to the use of heavy nuclear and deuteron data in the fit. We estimate these uncertainties by comparing the value… ▽ More As data become more precise, estimating theoretical uncertainties in global PDF determinations is likely to become increasingly necessary to obtain correspondingly precise PDFs. Here we present a next generation of global proton PDFs (NNPDF4.0) that include theoretical uncertainties due to the use of heavy nuclear and deuteron data in the fit. We estimate these uncertainties by comparing the values of the nuclear observables computed with the nuclear PDFs against those computed with proton PDFs. For heavy nuclear PDFs we use the nuclear nNNPDF2.0 set, while for deuteron PDFs we develop an iterative procedure to determine proton and deuteron PDFs simultaneously, each including the uncertainties in the other. Accounting for nuclear uncertainties resolves some of the tensions in the global fit of the proton PDFs, especially those between the nuclear data and the extended LHC data set used in NNPDF4.0. △ Less

Submitted 2 August, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: Revised submission to SciPost

Constraints on black-hole charges with the 2017 EHT observations of M87*

Authors: Prashant Kocherlakota, Luciano Rezzolla, Heino Falcke, Christian M. Fromm, Michael Kramer, Yosuke Mizuno, Antonios Nathanail, Hector Olivares, Ziri Younsi, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell, Wilfred Boland , et al. (212 additional authors not shown)

Abstract: Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87*… ▽ More Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 15 pages, 4 figures, published in PRD on May 19

Correlation of Theoretical Uncertainties in PDF Fits and Theoretical Uncertainties in Predictions

Authors: Richard D. Ball, Rosalyn L. Pearson

Abstract: We show how to account for correlations between theoretical uncertainties incorporated in parton distribution function (PDF) fits, and the theoretical uncertainties in the predictions made using these PDFs. We demonstrate by explicit calculations, both analytical and numerical, that these correlations can lead to corrections to the central values of the predictions, and reductions in both the PDF… ▽ More We show how to account for correlations between theoretical uncertainties incorporated in parton distribution function (PDF) fits, and the theoretical uncertainties in the predictions made using these PDFs. We demonstrate by explicit calculations, both analytical and numerical, that these correlations can lead to corrections to the central values of the predictions, and reductions in both the PDF uncertainties and the theoretical uncertainties in the prediction. We illustrate our results with predictions for top production rapidity distributions and the Higgs total cross-section at the LHC, using the NLO NNPDF3.1 PDF set which incorporates missing higher order uncertainties. We conclude that the inclusion of correlations can increase both the accuracy and precision of predictions involving PDFs, particularly for processes with data already included in the PDF fit. △ Less

Submitted 8 August, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: 35 pages, 14 figures: final published version

Report number: Edinburgh 2019/17

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

Forecasts of the trend in global-mean temperature to 2100 arising from the scenarios of first-difference CO2 and peak fossil fuel

Authors: L. Mark W. Leggett, David. A. Ball

Abstract: Two future scenarios that are not explicitly in the range of scenarios (the Representative Concentration Pathway scenarios) utilised by the IPCC. These two scenarios are the emissions trend under peak fossil fuel (for example, Mohr et al., 2015); and the climate sensitivity determinable from the relationship between first-difference CO2 and temperature recently shown by Leggett and Ball (2015). Th… ▽ More Two future scenarios that are not explicitly in the range of scenarios (the Representative Concentration Pathway scenarios) utilised by the IPCC. These two scenarios are the emissions trend under peak fossil fuel (for example, Mohr et al., 2015); and the climate sensitivity determinable from the relationship between first-difference CO2 and temperature recently shown by Leggett and Ball (2015). This paper provides forecasts of a global surface temperature trajectory to 2100 resulting from the effect of these two scenarios. The time-series models developed both displayed high statistical significance and converged in their forecasts, so adding to the potential robustness of the findings. Under the effect of the combination of the peak fossil fuel and first-difference CO2 scenarios, we found that temperature is forecast to continue to rise, but only gently, until around 2023 where it reaches a level slightly higher than at present; and from then to display an also gentle steady decrease. It is shown that this trajectory is markedly lower than that generated by the IPCC business-as-usual level-of-CO2 (RCP8.5) model (Riahi et al. 2011). These lower results are evidence that the climate problem may require less future preventative action than is presently being considered necessary. If so, it is stressed that the same evidence is support for the case that the peak fossil fuel problem would require ameliorative action. This globally required action is the same as it would have been for climate (as embodied in the Paris Agreement on Climate Change of the 21st Conference of the Parties of the UNFCCC in Paris adopted in 2015): the rapid transition to a predominantly renewable global energy system. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Empirical evidence for a global atmospheric temperature control system: physical structure

Authors: Mark Leggett, David Ball

Abstract: There is evidence that a natural control system influences global atmospheric surface temperature (Leggett and Ball, 2020). The present paper sets up and tests a hypothesis concerning the physical makeup of the sequential elements of the control system and its outcomes. The final outcome the control system influences is defined as global atmospheric surface temperature. The terms used for the cont… ▽ More There is evidence that a natural control system influences global atmospheric surface temperature (Leggett and Ball, 2020). The present paper sets up and tests a hypothesis concerning the physical makeup of the sequential elements of the control system and its outcomes. The final outcome the control system influences is defined as global atmospheric surface temperature. The terms used for the control system element types in the hypothesis are, in sequence, leading element, controller and actuator. Actuators are hypothesised to affect in turn the final outcome either directly or via penultimate outcomes. The existence of the control system is evidenced by demonstration of statistically significant one way Granger causality across each step of the hypothesised control system sequence. Evidence is presented that the leading element of the control system, represented by the Normalized Difference Vegetation Index, is the global biosphere. The biosphere as a control system has previously been referred to as Gaia (Lovelock and Margulis, 1974). A fourth (second derivative) term is found to enhance the Proportional, Integral and Derivative process terms of the control system shown in Leggett and Ball (2020). The main actuators of the control system found are shown to be wind speed and cloud cover. Cloud cover is shown to influence the final outcome, global surface temperature, directly. It and wind speed also influence the penultimate outcomes found, those of enhanced ocean heat uptake and enhanced outgoing longwave radiation. These together lead to control system output to the final outcome, global atmospheric temperature. Overall, evidence for the activity of the control system is present across many major physical dimensions of the global atmosphere. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Deuteron Uncertainties in the Determination of Proton PDFs

Authors: Richard D. Ball, Emanuele R. Nocera, Rosalyn L. Pearson

Abstract: We evaluate the uncertainties due to nuclear effects in global fits of proton parton distribution functions (PDFs) that utilise deep-inelastic scattering and Drell-Yan data on deuterium targets. To do this we use an iterative procedure to determine proton and deuteron PDFs simultaneously, each including the uncertainties in the other. We apply this procedure to determine the nuclear uncertainties… ▽ More We evaluate the uncertainties due to nuclear effects in global fits of proton parton distribution functions (PDFs) that utilise deep-inelastic scattering and Drell-Yan data on deuterium targets. To do this we use an iterative procedure to determine proton and deuteron PDFs simultaneously, each including the uncertainties in the other. We apply this procedure to determine the nuclear uncertainties in the SLAC, BCDMS, NMC and DYE866/NuSea fixed target deuteron data included in the NNPDF3.1 global fit. We show that the effect of the nuclear uncertainty on the proton PDFs is small, and that the increase in overall uncertainties is insignificant once we correct for nuclear effects. △ Less

Submitted 13 January, 2021; v1 submitted 30 October, 2020; originally announced November 2020.

Comments: Version accepted for publication in Eur. Phys. J. C; typo in Table 3 corrected; discussion about previous NNPDF studies added; various minor points clarified

Report number: Edinburgh 2020/15, Nikhef/2020-033

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

A Plasmoid Model for the Sgr A* Flares Observed with GRAVITY and Chandra

Authors: David Ball, Feryal Özel, Pierre Christian, Chi-Kwan Chan, Dimitrios Psaltis

Abstract: The Galactic Center black hole Sgr A* shows significant variability and flares in the submillimeter, infrared, and X-ray wavelengths. Owing to its exquisite resolution in the IR bands, the GRAVITY experiment for the first time spatially resolved the locations of three flares and showed that a bright region moves in ellipse-like trajectories close to but offset from the black hole over the course o… ▽ More The Galactic Center black hole Sgr A* shows significant variability and flares in the submillimeter, infrared, and X-ray wavelengths. Owing to its exquisite resolution in the IR bands, the GRAVITY experiment for the first time spatially resolved the locations of three flares and showed that a bright region moves in ellipse-like trajectories close to but offset from the black hole over the course of each event. We present a model for plasmoids that form during reconnection events and orbit in the coronal region around a black hole to explain these observations. We utilize general-relativistic radiative transfer calculations that include effects from finite light travel time, plasmoid motion, particle acceleration, and synchrotron cooling and obtain a rich structure in the flare lightcurves. This model can naturally account for the observed motion of the bright regions observed by the GRAVITY experiment and the offset between the center of the centroid motion and the position of the black hole. It also explains why some flares may be double-peaked while others have only a single peak and uncovers a correlation between the structure in the lightcurve and the location of the flare. Finally, we make predictions for future observations of flares from the inner accretion flow of Sgr A* that will provide a test of this model. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: 10 pages, 9 figures

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabili… ▽ More Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

The Mechanism of Electron Injection and Acceleration in Trans-Relativistic Reconnection

Authors: David Ball, Lorenzo Sironi, Feryal Özel

Abstract: Electron acceleration during magnetic reconnection is thought to play a key role in time-variable high-energy emission from astrophysical systems. By means of particle-in-cell simulations of trans-relativistic reconnection, we investigate electron injection and acceleration mechanisms in low-$β$ electron-proton plasmas. We set up a diversity of density and field structures (e.g., X-points and plas… ▽ More Electron acceleration during magnetic reconnection is thought to play a key role in time-variable high-energy emission from astrophysical systems. By means of particle-in-cell simulations of trans-relativistic reconnection, we investigate electron injection and acceleration mechanisms in low-$β$ electron-proton plasmas. We set up a diversity of density and field structures (e.g., X-points and plasmoids) by varying the guide field strength and choosing whether to trigger reconnection or let it spontaneously evolve. We show that the number of X-points and plasmoids controls the efficiency of electron acceleration, with more X-points leading to a higher efficiency. Using on-the-fly acceleration diagnostics, we also show that the non-ideal electric fields associated with X-points play a critical role in the first stages of electron acceleration. As a further diagnostic, we include two populations of test particles that selectively experience only certain components of electric fields. We find that the out-of-plane component of the parallel electric field determines the hardness of the high-energy tail of the electron energy distribution. These results further our understanding of electron acceleration in this regime of magnetic reconnection and have implications for realistic models of black hole accretion flows. △ Less

Submitted 16 August, 2019; originally announced August 2019.

Comments: 21 pages, 20 figures

Parton Distributions with Theory Uncertainties: General Formalism and First Phenomenological Studies

Authors: The NNPDF Collaboration, Rabah Abdul Khalek, Richard D. Ball, Stefano Carrazza, Stefano Forte, Tommaso Giani, Zahari Kassabov, Rosalyn L. Pearson, Emanuele R. Nocera, Juan Rojo, Luca Rottoli, Maria Ubiali, Cameron Voisey, Michael Wilson

Abstract: We formulate a general approach to the inclusion of theoretical uncertainties, specifically those related to the missing higher order uncertainty (MHOU), in the determination of parton distribution functions (PDFs). We demonstrate how, under quite generic assumptions, theory uncertainties can be included as an extra contribution to the covariance matrix when determining PDFs from data. We then rev… ▽ More We formulate a general approach to the inclusion of theoretical uncertainties, specifically those related to the missing higher order uncertainty (MHOU), in the determination of parton distribution functions (PDFs). We demonstrate how, under quite generic assumptions, theory uncertainties can be included as an extra contribution to the covariance matrix when determining PDFs from data. We then review, clarify, and systematize the use of renormalization and factorization scale variations as a means to estimate MHOUs consistently in deep inelastic and hadronic processes. We define a set of prescriptions for constructing a theory covariance matrix using scale variations, which can be used in global fits of data from a wide range of different processes, based on choosing a set of independent scale variations suitably correlated within and across processes. We set up an algebraic framework for the choice and validation of an optimal prescription by comparing the estimate of MHOU encoded in the next-to-leading order (NLO) theory covariance matrix to the observed shifts between NLO and NNLO predictions. We perform a NLO PDF determination which includes the MHOU, assess the impact of the inclusion of MHOUs on the PDF central values and uncertainties, and validate the results by comparison to the known shift between NLO and NNLO PDFs. We finally study the impact of the inclusion of MHOUs in a global PDF determination on LHC cross-sections, and provide guidelines for their use in precision phenomenology. In addition, we also compare the results based on the theory covariance matrix formalism to those obtained by performing PDF determinations based on different scale choices. △ Less

Submitted 17 October, 2019; v1 submitted 25 June, 2019; originally announced June 2019.

Comments: 73 pages, 27 figures. Final version, to be published in EPJC. Numerous changes and improvements including: discussion of phi added, Eq.s 6.3-6.5 and revised Tab 6.2 & 6.4; improved Figs 7.1-7.3; revised Sect. 81; revised Figs. B.1-B.4

Report number: Edinburgh 2019/9, Nikhef/2019-014, TIF-UNIMI-2019-9 DAMTP-2019-24, CAVENDISH-HEP-19-11

A First Determination of Parton Distributions with Theoretical Uncertainties

Authors: The NNPDF Collaboration, Rabah Abdul Khalek, Richard D. Ball, Stefano Carrazza, Stefano Forte, Tommaso Giani, Zahari Kassabov, Emanuele R. Nocera, Rosalyn L. Pearson, Juan Rojo, Luca Rottoli, Maria Ubiali, Cameron Voisey, Michael Wilson

Abstract: The parton distribution functions (PDFs) which characterize the structure of the proton are currently one of the dominant sources of uncertainty in the predictions for most processes measured at the Large Hadron Collider (LHC). Here we present the first extraction of the proton PDFs that accounts for the missing higher order uncertainty (MHOU) in the fixed-order QCD calculations used in PDF determ… ▽ More The parton distribution functions (PDFs) which characterize the structure of the proton are currently one of the dominant sources of uncertainty in the predictions for most processes measured at the Large Hadron Collider (LHC). Here we present the first extraction of the proton PDFs that accounts for the missing higher order uncertainty (MHOU) in the fixed-order QCD calculations used in PDF determinations. We demonstrate that the MHOU can be included as a contribution to the covariance matrix used for the PDF fit, and then introduce prescriptions for the computation of this covariance matrix using scale variations. We validate our results at next-to-leading order (NLO) by comparison to the known next order (NNLO) corrections. We then construct variants of the NNPDF3.1 NLO PDF set that include the effect of the MHOU, and assess their impact on the central values and uncertainties of the resulting PDFs. △ Less

Submitted 17 October, 2019; v1 submitted 10 May, 2019; originally announced May 2019.

Comments: 8 pages, 6 figures, PDF sets available from https://data.nnpdf.science/nnpdf31_th/. Final Version, to be published in EPJC letters. Appendix removed, several small improvements including brief discussion of usage

Report number: CAVENDISH-HEP-19-07, DAMTP-2019-20, Edinburgh 2019/07, Nikhef 2019-013, TIF-UNIMI-2019-4

Nuclear Uncertainties in the Determination of Proton PDFs

Authors: Richard D. Ball, Emanuele R. Nocera, Rosalyn L. Pearson

Abstract: We show how theoretical uncertainties due to nuclear effects may be incorporated into global fits of proton parton distribution functions (PDFs) that include deep-inelastic scattering and Drell-Yan data on nuclear targets. We specifically consider the CHORUS, NuTeV and E605 data included in the NNPDF3.1 fit, which used Pb, Fe and Cu targets, respectively. We show that the additional uncertainty in… ▽ More We show how theoretical uncertainties due to nuclear effects may be incorporated into global fits of proton parton distribution functions (PDFs) that include deep-inelastic scattering and Drell-Yan data on nuclear targets. We specifically consider the CHORUS, NuTeV and E605 data included in the NNPDF3.1 fit, which used Pb, Fe and Cu targets, respectively. We show that the additional uncertainty in the proton PDFs due to nuclear effects is small, as expected, and in particular that the effect on the $\bar{d}/\bar{u}$ ratio, the total strangeness $s+\bar{s}$, and the strange valence distribution $s-\bar{s}$ is negligible. △ Less

Submitted 23 March, 2019; v1 submitted 21 December, 2018; originally announced December 2018.

Comments: 34 pages, 15 figures, 4 tables - version accepted for publication in Eur. Phys. J. C

Report number: Edinburgh 2018/4; Nikhef 2018-062

Observational evidence that a Gaia-type feedback control system with proportional-integral-derivative characteristics is operating on atmospheric surface temperature at global scale

Authors: L. Mark W. Leggett, David A. Ball

Abstract: The Gaia hypothesis (Lovelock and Margulis, 1974) proposes that there is a control system operating at global level that regulates climate and chemistry at a habitable state for the biota. Here we provide statistically significant observational evidence that a feedback control system moderating atmospheric temperature is presently operating coherently at global scale, that is to say, observational… ▽ More The Gaia hypothesis (Lovelock and Margulis, 1974) proposes that there is a control system operating at global level that regulates climate and chemistry at a habitable state for the biota. Here we provide statistically significant observational evidence that a feedback control system moderating atmospheric temperature is presently operating coherently at global scale, that is to say, observational evidence for Gaia. Further, this control system is of a sophisticated type, involving the corrective feedback not only of a linear error term but also its derivative and its integral. This makes it of the same type as the most commonly used control system developed by humans, the proportional-integral-derivative (PID) control system. △ Less

Submitted 1 October, 2018; originally announced October 2018.

Electron and Proton Acceleration in Trans-Relativistic Magnetic Reconnection: Dependence on Plasma Beta and Magnetization

Authors: David Ball, Lorenzo Sironi, Feryal Özel

Abstract: Non-thermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection, where the magnetization $σ$, defined as the ratio of magnetic energy density to enthalpy density, is $\sim 1$, is frequently encountered in su… ▽ More Non-thermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection, where the magnetization $σ$, defined as the ratio of magnetic energy density to enthalpy density, is $\sim 1$, is frequently encountered in such flows. By means of a large suite of two-dimensional particle-in-cell simulations, we investigate electron and proton acceleration in the trans-relativistic regime. We focus on the dependence of the electron energy spectrum on $σ$ and the proton $β$ (i.e., the ratio of proton thermal pressure to magnetic pressure). We find that the electron spectrum in the reconnection region is non-thermal and can be generally modeled as a power law. At $β\lesssim 3 \times 10^{-3}$, the slope, $p$, is independent of $β$ and it hardens with increasing $σ$ as $p\simeq 1.8 +0.7/\sqrtσ$. Electrons are primarily accelerated by the non-ideal electric field at X-points, either in the initial current layer or in current sheets generated in between merging magnetic islands. At higher values of $β$, the electron power law steepens for all values of $σ$. At values of $β$ near $β_{\rm max}\approx1/4σ$, when both electrons and protons are relativistically hot prior to reconnection, the spectra of both species display an additional component at high energies, containing a few percent of particles. These particles are accelerated via a Fermi-like process by bouncing in between the reconnection outflow and a stationary magnetic island. We provide an empirical prescription for the dependence of the power-law slope and the acceleration efficiency on $β$ and $σ$, which can be used in global simulations of collisionless accretion disks. △ Less

Submitted 14 March, 2018; originally announced March 2018.

Comments: 19 pages, 23 figures, submitted to ApJ

Precision determination of the strong coupling constant within a global PDF analysis

Authors: Richard D. Ball, Stefano Carrazza, Luigi Del Debbio, Stefano Forte, Zahari Kassabov, Juan Rojo, Emma Slade, Maria Ubiali

Abstract: We present a determination of the strong coupling constant $α_s(m_Z)$ based on the NNPDF3.1 determination of parton distributions, which for the first time includes constraints from jet production, top-quark pair differential distributions, and the $Z$ $p_T$ distributions using exact NNLO theory. Our result is based on a novel extension of the NNPDF methodology - the correlated replica method - wh… ▽ More We present a determination of the strong coupling constant $α_s(m_Z)$ based on the NNPDF3.1 determination of parton distributions, which for the first time includes constraints from jet production, top-quark pair differential distributions, and the $Z$ $p_T$ distributions using exact NNLO theory. Our result is based on a novel extension of the NNPDF methodology - the correlated replica method - which allows for a simultaneous determination of $α_s$ and the PDFs with all correlations between them fully taken into account. We study in detail all relevant sources of experimental, methodological and theoretical uncertainty. At NNLO we find $α_s(m_Z) = 0.1185 \pm 0.0005^\text{(exp)}\pm 0.0001^\text{(meth)}$, showing that methodological uncertainties are negligible. We conservatively estimate the theoretical uncertainty due to missing higher order QCD corrections (N$^3$LO and beyond) from half the shift between the NLO and NNLO $α_s$ values, finding $Δα^{\rm th}_s =0.0011$. △ Less

Submitted 9 July, 2018; v1 submitted 9 February, 2018; originally announced February 2018.

Comments: 23 pages, 10 figures.Final version, published in EPJC. New plots added to figs 3.5 and 3.6 and improved discussion of the d'Agostini bias. Discussion of ABMP clarified. Some references added

Report number: CAVENDISH-HEP-18-02, CERN-TH-2018-024, Edinburgh 2017/19, Nikhef/2017-046, OUTP-17-14P, DAMTP-2018-3, TIF-UNIMI-2018-1

The Proton Spin, Semi-Inclusive processes, and a future Electron Ion Collider

Authors: Richard D. Ball, A. Deshpande

Abstract: We discuss spin physics, Guido Altarelli's contribution to it, and what we still have to learn. We set out in particular a programme for incorporating constraints from semi-inclusive data into global fits of polarized PDFs, and discuss the need for the EIC to increase the precision and kinematic coverage of current measurements. We discuss spin physics, Guido Altarelli's contribution to it, and what we still have to learn. We set out in particular a programme for incorporating constraints from semi-inclusive data into global fits of polarized PDFs, and discuss the need for the EIC to increase the precision and kinematic coverage of current measurements. △ Less

Submitted 23 January, 2018; v1 submitted 15 January, 2018; originally announced January 2018.

Comments: Contribution to the volume "From my Vast Repertoire: the Legacy of Guido Altarelli" : 23 pages, 6 figures

Report number: Edinburgh 2018/3, BNL-200028-2018-JAAM

Parton distributions with small-x resummation: evidence for BFKL dynamics in HERA data

Authors: Richard D. Ball, Valerio Bertone, Marco Bonvini, Simone Marzani, Juan Rojo, Luca Rottoli

Abstract: We present a determination of the parton distribution functions of the proton in which NLO and NNLO fixed-order calculations are supplemented by NLLx small-x resummation. Deep inelastic structure functions are computed consistently at NLO+NLLx or NNLO+NLLx, while for hadronic processes small-x resummation is included only in the PDF evolution, with kinematic cuts introduced to ensure the fitted da… ▽ More We present a determination of the parton distribution functions of the proton in which NLO and NNLO fixed-order calculations are supplemented by NLLx small-x resummation. Deep inelastic structure functions are computed consistently at NLO+NLLx or NNLO+NLLx, while for hadronic processes small-x resummation is included only in the PDF evolution, with kinematic cuts introduced to ensure the fitted data lie in a region where the fixed-order calculation of the hard cross-sections is reliable. In all other respects, the fits use the same methodology and are based on the same global dataset as the recent NNPDF3.1 analysis. We demonstrate that the inclusion of small-x resummation leads to a quantitative improvement in the perturbative description of the HERA inclusive and charm-production reduced cross-sections in the small x region. The impact of the resummation in our fits is greater at NNLO than at NLO, because fixed-order calculations have a perturbative instability at small x due to large logarithms that can be cured by resummation. We explore the phenomenological implications of PDF sets with small-x resummation for the longitudinal structure function $F_L$ at HERA, for parton luminosities and LHC benchmark cross-sections, for ultra-high energy neutrino-nucleus cross-sections, and for future high-energy lepton-proton colliders such as the LHeC. △ Less

Submitted 8 March, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 70 pages, many figures. Discussion on uncertainties due to subleading logarithmic contributions. Discussion on fits with pseudodata from future high-energy lepton-proton colliders. Updated references. Version to be published in EPJC

Report number: Edinburgh 2017/15, Nikhef/2017-027, OUTP-17-12P

Parton distributions from high-precision collider data

Authors: The NNPDF Collaboration, Richard D. Ball, Valerio Bertone, Stefano Carrazza, Luigi Del Debbio, Stefano Forte, Patrick Groth-Merrild, Alberto Guffanti, Nathan P. Hartland, Zahari Kassabov, José I. Latorre, Emanuele R. Nocera, Juan Rojo, Luca Rottoli, Emma Slade, Maria Ubiali

Abstract: We present a new set of parton distributions, NNPDF3.1, which updates NNPDF3.0, the first global set of PDFs determined using a methodology validated by a closure test. The update is motivated by recent progress in methodology and available data, and involves both. On the methodological side, we now parametrize and determine the charm PDF alongside the light quarks and gluon ones, thereby increasi… ▽ More We present a new set of parton distributions, NNPDF3.1, which updates NNPDF3.0, the first global set of PDFs determined using a methodology validated by a closure test. The update is motivated by recent progress in methodology and available data, and involves both. On the methodological side, we now parametrize and determine the charm PDF alongside the light quarks and gluon ones, thereby increasing from seven to eight the number of independent PDFs. On the data side, we now include the D0 electron and muon W asymmetries from the final Tevatron dataset, the complete LHCb measurements of W and Z production in the forward region at 7 and 8 TeV, and new ATLAS and CMS measurements of inclusive jet and electroweak boson production. We also include for the first time top-quark pair differential distributions and the transverse momentum of the Z bosons from ATLAS and CMS. We investigate the impact of parametrizing charm and provide evidence that the accuracy and stability of the PDFs are thereby improved. We study the impact of the new data by producing a variety of determinations based on reduced datasets. We find that both improvements have a significant impact on the PDFs, with some substantial reductions in uncertainties, but with the new PDFs generally in agreement with the previous set at the one sigma level. The most significant changes are seen in the light-quark flavor separation, and in increased precision in the determination of the gluon. We explore the implications of NNPDF3.1 for LHC phenomenology at Run II, compare with recent LHC measurements at 13 TeV, provide updated predictions for Higgs production cross-sections and discuss the strangeness and charm content of the proton in light of our improved dataset and methodology. The NNPDF3.1 PDFs are delivered for the first time both as Hessian sets, and as optimized Monte Carlo sets with a compressed number of replicas. △ Less

Submitted 22 September, 2017; v1 submitted 1 June, 2017; originally announced June 2017.

Comments: Revised version, to be publ.in EPJC. Typos corrected and points clarified. Figs 2.4, 4.21, 4.30 and 4.31 added. PDF comparison plots added. Fit with nuclear corrections added in Sect 4.11. Labeling of all distance plots corrected (u and d were interchanged). Points in charm cross-sections without cuts added to Tab 3.2

Report number: CAVENDISH-HEP-17-06, CERN-TH-2017-077, Edinburgh 2017/08, Nikhef/2017-006, OUTP-17-04P,TIF-UNIMI-2017-3

The Properties of Reconnection Current Sheets in GRMHD Simulations of Radiatively Inefficient Accretion Flows

Authors: David Ball, Feryal Özel, Dimitrios Psaltis, Chi-Kwan Chan, Lorenzo Sironi

Abstract: Non-ideal MHD effects may play a significant role in determining the dynamics, thermal properties, and observational signatures of radiatively inefficient accretion flows onto black holes. In particular, particle acceleration during magnetic reconnection events may influence black hole spectra and flaring properties. We use representative GRMHD simulations of black hole accretion flows to identify… ▽ More Non-ideal MHD effects may play a significant role in determining the dynamics, thermal properties, and observational signatures of radiatively inefficient accretion flows onto black holes. In particular, particle acceleration during magnetic reconnection events may influence black hole spectra and flaring properties. We use representative GRMHD simulations of black hole accretion flows to identify and explore the structures and properties of current sheets as potential sites of magnetic reconnection. In the case of standard and normal (SANE) disks, we find that, in the reconnection sites, the plasma beta ranges from $0.1$ to $1000$, the magnetization ranges from $10^{-4}$ to $1$, and the guide fields are weak compared to the reconnecting fields. In magnetically arrested (MAD) disks, we find typical values for plasma beta from $10^{-2}$ to $10^3$, magnetizations from $10^{-3}$ to $10$, and typically stronger guide fields, with strengths comparable to or greater than the reconnecting fields. These are critical parameters that govern the electron energy distribution resulting from magnetic reconnection and can be used in the context of plasma simulations to provide microphysics inputs to global simulations. We also find that ample magnetic energy is available in the reconnection regions to power the fluence of bright X-ray flares observed from the black hole in the center of the Milky Way. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: 8 pages, 8 figures, submitted to ApJ

Charm Production: Pole Mass or Running Mass?

Authors: Richard D. Ball

Abstract: We consider the role of the charm mass in PDF fits, and in particular the reliability of recent determinations of the MSbar running mass from fits to inclusive charm production data. We explain how the use of the pole mass, with a charm PDF determined empirically, can reduce PDF uncertainties from charm mass effects, and thereby lead to more reliable predictions for hard processes at the LHC. We consider the role of the charm mass in PDF fits, and in particular the reliability of recent determinations of the MSbar running mass from fits to inclusive charm production data. We explain how the use of the pole mass, with a charm PDF determined empirically, can reduce PDF uncertainties from charm mass effects, and thereby lead to more reliable predictions for hard processes at the LHC. △ Less

Submitted 12 December, 2016; originally announced December 2016.

Comments: 9 pages, no figures, proceedings article for Diffraction 2016

Report number: Edinburgh 2016/21

A Determination of the Charm Content of the Proton

Authors: The NNPDF Collaboration, Richard D. Ball, Valerio Bertone, Marco Bonvini, Stefano Carrazza, Stefano Forte, Alberto Guffanti, Nathan P. Hartland, Juan Rojo, Luca Rottoli

Abstract: We present an unbiased determination of the charm content of the proton, in which the charm parton distribution function (PDF) is parametrized on the same footing as the light quarks and the gluon in a global PDF analysis. This determination relies on the calculation of deep-inelastic structure functions in the FONLL scheme, generalized to account for massive charm-initiated contributions. In cont… ▽ More We present an unbiased determination of the charm content of the proton, in which the charm parton distribution function (PDF) is parametrized on the same footing as the light quarks and the gluon in a global PDF analysis. This determination relies on the calculation of deep-inelastic structure functions in the FONLL scheme, generalized to account for massive charm-initiated contributions. In contrast to the usual situation in which the charm PDF is assumed to be generated perturbatively by pair radiation off gluons and light quarks, vanishing at a scale set by the value of the charm mass m_c, we find that the fitted charm PDF vanishes within uncertainties at a scale Q~1.5 GeV for all x<~0.1, independent of the value of m_c used in the coefficient functions. We also find some evidence that the charm PDF at large x>~0.1 and low scales does not vanish, but rather has an "intrinsic" component, very weakly scale dependent and almost independent of the value of m_c, carrying about 1% of the total momentum of the proton. The uncertainties in all other PDFs are only slightly increased by the inclusion of fitted charm, while the dependence of these PDFs on m_c is significantly reduced. When the EMC charm structure function dataset is included, it is well described by the fit, and PDF uncertainties in the fitted charm PDF are significantly reduced, though we verify that excluding the EMC data does not qualitatively modify any of our findings. The increased stability with respect to m_c persists at high scales and is the main implication of our results for LHC phenomenology. Fitting the charm PDF modifies the predictions for processes such as high p_T and large rapidity charm pair production and Z+c production, and thus we expect that future LHC data will further constrain the charm content of the proton. △ Less

Submitted 23 November, 2016; v1 submitted 20 May, 2016; originally announced May 2016.

Comments: 46 pages, 32 figures. Final version, to be published in EPJC. Sections 3.2, 3.3, 3.4 and conclusions substantially rewritten and improved. Figs 11 and 14 improved, figs 10, 13, 16 added. Last typos fixed

Report number: Edinburgh 2016/06, CERN-TH-2016-087, OUTP-16-03P, TIF-UNIMI-2016-3

The asymptotic behaviour of parton distributions at small and large $x$

Authors: Richard D. Ball, Emanuele R. Nocera, Juan Rojo

Abstract: It has been argued from the earliest days of quantum chromodynamics (QCD) that at asymptotically small values of $x$ the parton distribution functions (PDFs) of the proton behave as $x^α$, where the values of $α$ can be deduced from Regge theory, while at asymptotically large values of $x$ the PDFs behave as $(1-x)^β$, where the values of $β$ can be deduced from the Brodsky-Farrar quark counting r… ▽ More It has been argued from the earliest days of quantum chromodynamics (QCD) that at asymptotically small values of $x$ the parton distribution functions (PDFs) of the proton behave as $x^α$, where the values of $α$ can be deduced from Regge theory, while at asymptotically large values of $x$ the PDFs behave as $(1-x)^β$, where the values of $β$ can be deduced from the Brodsky-Farrar quark counting rules. We critically examine these claims by extracting the exponents $α$ and $β$ from various global fits of parton distributions, analysing their scale dependence, and comparing their values to the naive expectations. We find that for valence distributions both Regge theory and counting rules are confirmed, at least within uncertainties, while for sea quarks and gluons the results are less conclusive. We also compare results from various PDF fits for the structure function ratio $F_2^n/F_2^p$ at large $x$, and caution against unrealistic uncertainty estimates due to overconstrained parametrisations. △ Less

Submitted 27 March, 2017; v1 submitted 31 March, 2016; originally announced April 2016.

Comments: 20 pages, 9 figures, this version matches the version accepted for publication in EPJC

Particle Acceleration and the Origin of X-ray Flares in GRMHD simulations of Sgr A*

Authors: David Ball, Feryal Ozel, Dimitrios Psaltis, Chi-kwan Chan

Abstract: Significant X-ray variability and flaring has been observed from Sgr A* but is poorly understood from a theoretical standpoint. We perform GRMHD simulations that take into account a population of non-thermal electrons with energy distributions and injection rates that are motivated by PIC simulations of magnetic reconnection. We explore the effects of including these non-thermal electrons on the p… ▽ More Significant X-ray variability and flaring has been observed from Sgr A* but is poorly understood from a theoretical standpoint. We perform GRMHD simulations that take into account a population of non-thermal electrons with energy distributions and injection rates that are motivated by PIC simulations of magnetic reconnection. We explore the effects of including these non-thermal electrons on the predicted broadband variability of Sgr A* and find that X-ray variability is a generic result of localizing non-thermal electrons to highly magnetized regions, where particles are likely to be accelerated via magnetic reconnection. The proximity of these high-field regions to the event horizon forms a natural connection between IR and X-ray variability and accounts for the rapid timescales associated with the X-ray flares. The qualitative nature of this variability is consistent with observations, producing X-ray flares that are always coincident with IR flares, but not vice versa, i.e., there are a number of IR flares without X-ray counterparts. △ Less

Submitted 23 August, 2016; v1 submitted 18 February, 2016; originally announced February 2016.

Journal ref: 2016 ApJ 826 77

Charm in Deep-Inelastic Scattering

Authors: Richard D. Ball, Marco Bonvini, Luca Rottoli

Abstract: We show how to extend systematically the FONLL scheme for inclusion of heavy quark mass effects in DIS to account for the possible effects of an intrinsic charm component in the nucleon. We show that when there is no intrinsic charm, FONLL is equivalent to S-ACOT to any order in perturbation theory, while when an intrinsic charm component is included FONLL is identical to ACOT, again to all orders… ▽ More We show how to extend systematically the FONLL scheme for inclusion of heavy quark mass effects in DIS to account for the possible effects of an intrinsic charm component in the nucleon. We show that when there is no intrinsic charm, FONLL is equivalent to S-ACOT to any order in perturbation theory, while when an intrinsic charm component is included FONLL is identical to ACOT, again to all orders in perturbation theory. We discuss in detail the inclusion of top and bottom quarks to construct a variable flavour number scheme, and give explicit expressions for the construction of the structure functions $F^c_2$, $F^c_L$ and $F^c_3$ to NNLO. △ Less

Submitted 20 November, 2015; v1 submitted 8 October, 2015; originally announced October 2015.

Comments: 43 pages, 1 figure. Final JHEP version

Report number: Edinburgh 2015/06, CERN-PH-TH/2015-118, OUTP-15-25P