Arxiv:2108.13388V1

Eﬀect of Photon Counting Shot Noise on TIRM

Fan Cuia and David J. Pineab Total internal reflection microscopy (TIRM) measures changes in the distance between a colloidal particle and a transparent substrate by measuring the intensity of light scattered by the particle when it is illuminated by the evanescent field that is created from light totally internally reflected at the substrate interface. From these measurements, the height-dependent effective potential ϕ(z) between the colloidal particle and the substrate can be measured. The spatial resolution with which TIRM can resolve the height z and effective potential ϕ(z) is limited by the intrinsic shot noise of the photon counting process used to measure the scattered light intensity. We develop a model to determine the spatial resolution with which TIRM can measure ϕ(z) and verify its validity with simulations and experiments. We establish the critical role of shot noise and intensity integration time τ in TIRM measurements of colloidal interaction potentials, which is a trade-off between minimizing shot noise and preserving the spatial resolution of the potential. 1 Introduction evanescent wave, which is typically in the range of 70-200 nm. Total internal reflection microscopy 1 (TIRM) is a powerful Because the intensity I of the light scattered by the particle is pro- method for measuring the microscopic interactions of colloidal portional to the intensity Ie of the evanescent field that is incident particles in a liquid suspension. Since its development some 30 on the particle, the scattered intensity also depends exponentially years ago, 2 it has been used to measure various colloidal interac- on the height of the particle, 1,8 tions, including screened electrostatic repulsion, 3 steric repulsion −βz due to grafted or adsorbed polymers, 4,5 van der Waals attrac- I(z) = I0 e . (3) tion, 4 depletion attraction, 6 and critical Casimir interactions. 7 Thus, we see that the scattered intensity is related directly to the Spatial resolutions as small as 1 nm have been reported. 2 As such, height z of a particle. This gives us a way to determine the height TIRM has become an invaluable tool for understanding colloidal z of the particle. interactions at a microscopic scale. The probability P(I)dI that the scattered intensity is between I Figure 1 shows a schematic of a typical TIRM experiment and and I +dI is the equal to the probability p(z)dz that the particle is summarizes the basic experimental setup. between a height of z and z + dz, and thus are related by The TIRM technique is based on two simple ideas. First, the probability that a particle at equilibrium in a liquid suspension P(I)|dI| = p(z)|dz| . (4) is at a height z above the substrate is given by the Boltzmann distribution Solving for p(z) yields,

p(z) = Ae−ϕ(z)/kBT , (1) dI p(z) = P(I) = β P(I)I(z) , (5) dz where ϕ(z) is the effective interaction potential between the particle and the substrate. Inverting this equation, we obtain the where we have used eqn (3) to evaluate the derivative. Thus, the expression probability distribution of heights p(z) appearing in eqn (2) can

ϕ(z) − ϕ(z ) p(z ) r = ln r , (2) k T p(z) B photon counter where we have introduced a reference height zr to eliminate the normalization factor A that appears in eqn (1). Equation (2) tells us that if we can measure the distribution of particle heights p(z) arXiv:2109.11651v1 [cond-mat.soft] 23 Sep 2021 above the substrate, we can determine the effective interaction between the particle and the substrate. particle The second idea starts with the observation that the evanescent intensity Ie of light totally internally reﬂected from the z substrate decays exponentially with height above the substrate, laser beam −1 Ie(z) = Ie0 exp(−βz), where β is the penetration depth of the prism

Fig. 1 TIRM schematic. Laser light enters from the right and is totally internally reﬂected at the substrate, launching an exponentially-damped a Department of Physics, Center for Soft Matter Research, New York University, 726 evanescent wave towards the particle, which is a distance z above the sub- Broadway, New York, NY 10003, USA. E-mail: [email protected], [email protected] strate in a liquid suspension. Light scattered by the particle is collected b Department of Chemical & Biomolecular Engineering, Tandon School of Engineering, by a microscope objective lens and directed toward a photon counter. New York University, 6 MetroTech Center, Brooklyn, New York 11201, USA.

1–10 | 1 Reconstruction in ALICE and calibration of TPC space-charge distortions in Run 3

Ernst Hellbär푎,∗for the ALICE collaboration 푎GSI - Helmholtzzentrum für Schwerionenforschung GmbH E-mail: [email protected]

The ALICE experiment will run with continuous readout at interaction rates of up to 50 kHz in Pb–Pb collisions during Run 3 of the LHC. In order to achieve this goal, a new data processing scheme and software are developed. This scheme strongly relies on GPUs (Graphics Process- ing Unit) for fast online and offline calibration and reconstruction as well as on efficient data compression. On the hardware side, the Time Projection Chamber (TPC), among other detector systems, received major upgrades to its readout chambers and readout electronics. The multiwire proportional chambers were replaced by stacks of four Gas Electron Multiplier foils to allow for continuous readout while keeping the ion backflow below 1%, minimizing space-charge effects from amplification ions entering the drift volume. Nevertheless, significant space-point distortions due to space charge are expected at the highest interaction rates in Pb–Pb collisions. In addition, space-charge density fluctuations lead to distortion fluctuations which have to be corrected on time scales of 10 ms in order to preserve the intrinsic tracking resolution of the TPC. While the average space-charge distortions can be corrected using information from external detectors as a reference, data-driven machine learning algorithms and convolutional neural networks are foreseen to provide the correction for the distortion fluctuations. arXiv:2109.12000v1 [physics.ins-det] 24 Sep 2021 The Ninth Annual Conference on Large Hadron Collider Physics - LHCP2021 7-12 June 2021 Online

∗Speaker

Reconstruction in ALICE and calibration of TPC space-charge distortions in Run 3 Ernst Hellbär

1. Introduction

The ALICE experiment aims to record a Pb–Pb minimum-bias data sample in Run 3 and Run 4 of the LHC which is 50–100 times larger than that recorded in the first two running periods. This will be realized by fully exploiting the capabilities of the LHC, which will deliver Pb–Pb collisions at interaction rates of up to 50 kHz. Contrary to the event-based triggered readout which was employed before Run 3, ALICE will now continuously process and record the data provided by the detectors, i.e. data from every collision will be available. Continuous readout poses several challenges and new requirements for the global data processing strategy [1], including the calibration and reconstruction software and, for some of the detector systems, also requiring hardware upgrades. The data processing is performed in two stages, synchronous processing and asynchronous processing. The former takes place during data taking and includes a first calibration and tracking, while the final calibration and full reconstruction of the data is performed during the asynchronous processing at a later point in time. At both stages, the application of GPUs for a majority of the processing leads to a significant gain in computing time, which is required to continuously process the amounts of data expected to be delivered in Run 3 and Run 4.

The Time Projection Chamber (TPC) [2] is the main detector for tracking and particle identifi- cation in ALICE. During the long shutdown of the LHC before Run 3, the readout chambers and electronics were upgraded [3] in order to provide the capability to run in continuous readout mode. Instead of multiwire proportional chambers and a triggered gating grid, which effectively limited the readout rate to about 3 kHz in order to avoid space-charge effects from ion backflow (IBF), the upgraded readout chambers are based on Gas Electron Multiplier (GEM) foils. A combination of four GEM foils with different hole pitches and an optimized high-voltage configuration suppresses the IBF to below 1%, eliminating the need for a gating grid and allowing the operation of the detector with continuous readout. Taking into account the high interaction rates foreseen for Pb–Pb collisions in Run 3, significant space-point distortions due to space charge from IBF are still expected, requiring a sophisticated calibration approach to preserve the intrinsic tracking performance of the TPC.

2. Space-charge eﬀects in the TPC

The space-charge density depends on the amount of ions piling up within one ion drift time. It is driven by the number of ion pile-up events, which is the product of the collision rate and the ion drift time, as well as by the eﬀective gain and the IBF, which are combined to 휀 = gain × IBF. The expected space-charge density distribution at 50 kHz Pb–Pb is derived using the hit distributions from Monte Carlo simulations. As dedicated measurements of the ion mobility in Ne-CO2-N2 (90-10-5) [4] indicate an ion drift time around 200 ms for the full drift, ions from 10000 events pile up on average in the TPC drift volume at any given moment at 50 kHz interaction rate. The space- charge distortions d푟, d푟휑 and d푧 for any given space-charge density distribution are calculated using numerical algorithms to solve the Poisson equation and the Langevin equation. The average radial space-charge distortions for 50 kHz Pb–Pb collisions are shown in the left plot of Fig.1 as a function of the TPC radius and 푧-position.

2 ? GEANT4-based calibration of an organic liquid scintillator

N. Mauritzsona, K.G. Fissuma,∗, H. Perreya, J.R.M. Annandc, R.J.W. Frosta, R. Hall-Wiltonb,c,d, R. Al Jebalib,c, K. Kanakib, V. Maulerovaa,b, F. Messia,1, E. Roforsa

aDivision of Nuclear Physics, Lund University, SE-221 00 Lund, Sweden bDetector Group, European Spallation Source ERIC, SE-221 00 Lund, Sweden cSUPA School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, Scotland, UK dDipartimento di Fisica “G. Occhialini”, Università degli Studi di Milano-Bicocca, Piazza della Scienza 3, 20126 Milano, Italy

Abstract

A light-yield calibration of an NE 213A organic liquid scintillator detector has been performed using both monoenergetic and polyenergetic gamma-ray sources. Scintillation light was detected in a photomultiplier tube, and the corresponding pulses were subjected to waveform digitization on an event-by-event basis. The resulting Compton edges have

been analyzed using a GEANT4 simulation of the detector which models both the interactions of the ionizing radiation as well as the transport of scintillation photons. The simulation is calibrated and also compared to well-established prescriptions used to determine the Compton edges, resulting ultimately in light-yield calibration functions. In the process, the simulation-based method produced information on the gain and intrinsic pulse- height resolution of the detector. It also facilitated a previously inaccessible understanding of the systematic uncertainties associated with the calibration of the scintillation-light yield. The simulation-based method was also compared to well-established numerical prescriptions for locating the Compton edges. Ultimately, the simulation predicted as much as 17% lower light-yield calibrations than the prescriptions. These calibrations indicate that approximately 35% of the scintillation light associated with a given gamma-ray reaches the photocathode. It is remarkable how well two 50 year old prescriptions for calibrating scintillation-light yield in organic scintillators have stood the test of time. Keywords: scintillation light-yield calibration, organic liquid scintillator, NE 213A,

gamma-rays, Compton edge, GEANT4

?The data set doi:10.5281/zenodo.5524234 is available for download from https://zenodo.org/record/5524234. arXiv:2109.11963v1 [physics.ins-det] 24 Sep 2021 ∗Corresponding author. Telephone: +46 46 222 9677; Fax: +46 46 222 4709 1present address: DVel AB, Scheelevägen 32, SE-223 63 Lund, Sweden

Preprint submitted to Nuclear Instruments and Methods A September 27, 2021 Prepared for submission to JINST

FPGA-based low-cost synchronized ﬁber network for experimental setups in space

Tim Oberschulte푎 Thijs Wendrich푏 Holger Blume푎 푎Institute of Microelectronic Systems, Leibniz Universität Hannover, Appelstraße 4, 30167 Hannover, Germany 푏Institute of Quantum Optics, Leibniz Universität Hannover, Welfengarten 1, 30167 Hannover, Germany

E-mail: [email protected]

Abstract: Custom experiment setups in physics often require control electronics to execute actions and measurements on a small time scale. When further constraints limit the experiment’s environment, for example when the experiment is inside a sounding rocket, conventional network systems will not suffice those constraints because of weight, heat or budget limitations. This paper proposes a network architecture with a time resolution of less than 1 ns over a pair of plastic fibers while using low-cost commercial hardware. The plastic fibers in comparison to copper fibers have a low weight and additionally can isolate parts of the setup galvanically. Data rates of 40 Mbit s−1 enable the network to transfer large amounts of measurements and configuration data over the network. Proof-of-concept implementations of network endpoints and switches on small FPGAs are analyzed in terms of synchronicity, data rate and resource usage. Using commercial parts the resolution of 1 ns is reached with a standard deviation of less than 100 ps. Compared to a copper wire implementation the weight is reduced by about one order of magnitude. With its low weight at a low cost, the network is useful in space or laboratory setups which require high time resolution.

Keywords: Trigger concepts and systems, Space instrumentation, VLSI circuits arXiv:2109.11638v1 [physics.ins-det] 23 Sep 2021 Prepared for submission to JINST

22nd iWoRID 2021 27.06-01.07.2021 Ghent

GBTX emulator for development and special versions of GBT-based readout chains

Wojciech M. Zabołotny,푎 Adrian P. Byszuk,푎 Dmitrii Dementev,푏 David Emschermann,푐 Marek Gumiński,푎 Michał Kruszewski,푎 Piotr Miedzik,푎 Krzysztof Poźniak,푎 Ryszard Romaniuk,푎 Christian J. Schmidt푐 and Mikhail Shitenkov푏 푎Institute of Electronic Systems, Warsaw University of Technology, Nowowiejska 15/19, 00-665 Warszawa, Poland 푏Veksler and Baldin Laboratory of High Energy Physics, Joint Institute for Nuclear Research, Baldin str. 6, Dubna 141980, Russian Federation 푐GSI - Helmholtzzentrum für Schwerionenforschung GmbH, Planckstraße 1, 64291 Darmstadt, Germany

E-mail: [email protected]

Abstract: The GBTX ASIC is a standard solution for providing fast control and data readout for radiation detectors used in HEP experiments. However, it is subject to export control restrictions due to the usage of radiation-hard technology. An FPGA-based GBTX emulator (GBTxEMU) has been developed to enable the development of GBT-based readout chains in countries where the original GBTX cannot be imported. Thanks to utilizing a slightly modiﬁed GBT-FGPA core, it maintains basic compatibility with standard GBT-based systems. The GBTxEMU also may be an interesting solution for developing GBT-based readout chains for less demanding experiments.

Keywords: Data acquisition circuits, Front-end electronics for detector readout, Digital electronic circuits arXiv:2109.11591v1 [physics.ins-det] 23 Sep 2021

1Corresponding author. RIPTIDE: a novel recoil-proton track imaging detector for fast neutrons

A. Musumarra,a,b F. Leone,a,b C. Massimi,c,b,1 M.G. Pellegriti,b F. Romano,b R. Spighi,b and M. Villac,b a University of Catania, Department of Physics and Astronomy Via Santa Sofia 64, I-95123 Catania, Italy b Istituto Nazionale di Fisica Nucleare (INFN) Via Enrico Fermi 54, I-00044 Frascati (Rome) Italy c University of Bologna, Department of Physics and Astronomy I-40126 Bologna Via Irnerio 46, I-40126 Bologna, Italy E-mail: [email protected]

ABSTRACT: Neutron detectors are an essential tool for the development of many research fields, as nuclear, particle and astroparticle physics as well as radiotherapy and radiation safety. Since neutrons cannot directly ionize, their detection is only possible via nuclear reactions. Consequently, neutron-based experimental techniques are related to the detection of charged particle or electromagnetic radiation originating from neutron-induced reactions. The study of fast neutrons is often based on the neutron-proton elastic scattering reaction. In this case, the ionization induced by the recoil protons in a hydrogenous material constitutes the basic information for the design and development of neutron detectors. Although experimental techniques have continuously improved and refined, so far, proton-recoil track imaging is still weak in laboratory rate environments because of the extremely small detection efficiency. To address this deficiency, we propose a novel recoil-proton track imaging system in which the light deriving from a fast scintillation signal is used to perform a complete reconstruction in space and time of the event. In particular, we report the idea of RIPTIDE (RecoIl Proton Track Imaging DEtector): an innovative system which combines a plastic scintillator coupled to imaging devices, based on CMOS technology, or Micro Channel Plate sensors. The proposed apparatus aims at providing neutron spectrometry capability by stereoscopically imaging the recoil-protons tracks, correlating the spatial information with the time information. RIPTIDE intrinsically enable the online analysis of the ionization track, thus retrieving the neutron direction and energy, without spoiling the overall efficiency of the detection system. Finally, the spatial and topological event reconstruction enables particle discrimination – a crucial requirement for neutron detection - by deducing the specific energy loss along the track.

KEYWORDS: fast neutron detection; imaging; scintillators.

1 Corresponding author. The Radial Action from Probe Amplitudes to All Orders

Uri Kol,1, ∗ Donal O’Connell,2, † and Ofri Telem3, 4, ‡ 1Center for Cosmology and Particle Physics, Department of Physics New York University, 726 Broadway, New York, NY 10003, USA 2Higgs Centre for Theoretical Physics, School of Physics and Astronomy, Peter Guthrie Tate Road, The University of Edinburgh, EH9 3FD, Scotland 3Department of Physics, University of California, Berkeley, CA 94720, USA 4Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA We extract the relativistic classical radial action from scattering amplitudes, to all orders in perturbation theory, in the probe limit. Our sources include point charges and monopoles, as well as the Schwarzschild and pure-NUT gravitational backgrounds. A characteristic relativistic eﬀect, that scattering trajectories may wind around these sources any number of times, can be recovered when all-order amplitudes are available. We show that the amplitude for scattering a probe oﬀ a pure NUT is given by the solution of a transcendental equation involving continued fractions, and explain how to solve this equation to any desired loop order.

I. INTRODUCTION

Recent years have seen tremendous progress in the field of Post-Newtonian (PN) and Post-Minkowskian (PM) Effective Field Theory (EFT), relevant for the analytic calculation of compact object coalescence events. For compre- hensive reviews of PN methods in GR, see [1–3] and references within, as well as more recent progress up to 6-7PN [4–9]. This progress has been of fundamental importance for reliable analytic and semi-analytic determination of gravitational waveforms, especially in the context of Effective-One-Body Hamiltonians [10–23]. Besides the more conventional PN methods, there has also been a renewed interest in PM methods [24–61], which involve an expansion in Newton’s constant G but not in v/c. A central idea which underlies much of the PM approach (as well as some of the PN approach) is the matching of General Relativity (GR) to an effective field theory [62, 63]. This matching can be done at the level of classical observables such as the scattering angle [64, 65] or the impulse [66, 67]. One promising approach to PM dynamics involves an unexpected connection√ — to relativistic scattering amplitudes in quantum field theory [24, 28, 66–71]. Since gravitons couple with a G factor, the PM expansion goes over to the usual perturbative loop expansion of the corresponding scattering amplitude, with the (n − 1)-loop amplitude responsible for the n-PM EFT coefficients. The basic advantage of the use of scattering amplitudes is they can be evaluated using cutting edge tools such as generalized unitarity [72, 73], the double copy [74–76], as well as sophisticated loop integration methods [77]. Impressive concrete progress in this direction has been made [28, 78–82] by matching full GR to an effective 2-body theory at the level of quantum scattering amplitudes, or more precisely their ~ → 0 limit. Though very promising, the on-shell approach is not without its difficulties, especially at order G4 and above where it need to capture tail effects [83–91] which are non-local in time, e.g. an outgoing gravitational wave which is reflected in the far zone back into the inspiral. R r 0 0 One striking feature of the PM EFT approach has been the importance of the radial action Ir(r) ≡ pr(r )dr , rturn which arises in the classical Hamilton-Jacobi equation for the system [92]. This should not come as a surprise, as the radial action encapsulates all of the classical dynamics of the system. Inspired by the eikonal approximation, the arXiv:2109.12092v1 [hep-th] 24 Sep 2021 authors of [82] suggest the following amplitude-radial action relation Z iA ∼ eiIr − 1 . (1) j

m1m2 The aim of this paper is to derive this relation in the classical limit and in the probe limit κ ≡ 2 → 0, to (m1+m2) all orders in G/j. We do this by solving the QM scattering problem for a probe mass in curved space — the same calculation which leads to BH greybody factors [93–95] and quasinormal modes [96, 97]. However, we are interested in

∗Electronic address: [email protected] †Electronic address: [email protected] ‡Electronic address: [email protected] Cosmological Implications of Axion-Matter Couplings

, Daniel Green,F, Yi GuoF and Benjamin WallischF ♠

F Department of Physics, University of California, San Diego, La Jolla, CA 92093, USA

♠ School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, USA

Abstract Axions and other light particles appear ubiquitously in physics beyond the Standard Model, with a variety of possible couplings to ordinary matter. Cosmology oﬀers a unique probe of these particles as they can thermalize in the hot environment of the early universe for any such coupling. For sub-MeV particles, their entropy must leave a measurable cosmological signal,

usually via the eﬀective number of relativistic particles, Neﬀ . In this paper, we will revisit the cosmological constraints on the couplings of axions and other pseudo-Nambu-Goldstone bosons to Standard Model fermions from thermalization below the electroweak scale, where these couplings

are marginal and give contributions to the radiation density of ∆Neﬀ > 0.027. We update the calculation of the production rates to eliminate unnecessary approximations and ﬁnd that the cosmological bounds on these interactions are complementary to astrophysical constraints, e.g. from supernova SN 1987A. We additionally provide quantitative explanations for these bounds and their relationship. arXiv:2109.12088v1 [astro-ph.CO] 24 Sep 2021 Notes on symmetries in particle physics

Akash Jain Department of Physics & Astronomy, University of Victoria, PO Box 1700 STN CSC, Victoria, BC, V8W 2Y2, Canada

Abstract: These are introductory notes on symmetries in quantum field theory and how they apply to particle physics. The notes cover the fundamentals of group theory, their representations, Lie groups, and Lie algebras, along with an elaborate discussion of the representations of SU(N), Lorentz, and Poincaré groups and their respective algebras. We spend a lot of time on the realisation of these symmetry groups in quantum field theory, as both global and gauge symmetries, as well as their spontaneous breaking and the Higgs mechanism. In the end, we culminate all the lessons from the course to enumerate the symmetries and field content of the Standard Model of particle physics and write down the Standard Model Lagrangian. Special consideration is given to how the weak-force gauge bosons and the matter fields obtain their mass via the Higgs mechanism. These lecture notes were developed while teaching the graduate course PHYS 509 - Standard Model Phenomenology at the University of Victoria in the spring of 2021. They build upon and extend the hand-written notes of Prof. Adam Ritz who had developed and delivered the course in previous years. That prior work is acknowledged, and included here with permission.

Date: September 27, 2021

Comments: 86 + 1 pages, 1 ﬁgure, 1 table, LATEX arXiv:2109.12087v1 [hep-th] 24 Sep 2021

E-mail: [email protected] September 27, 2021 0:34 ws-procs961x669 WSPC Proceedings - 9.61in x 6.69in MG16˙JSP page 1

Running vacuum interacting with dark matter or with running gravitational coupling. Phenomenological implications

Joan Sol`aPeracaulaa,b aDepartament de F´ısica Qu`antica i Astrof´ısica, and b Institute of Cosmos Sciences, Universitat de Barcelona, Av. Diagonal 647, E-08028 Barcelona, Catalonia, Spain, ∗E-mail: [email protected]

The cosmological term, Λ, in Einstein’s equations is an essential ingredient of the ‘concordance’ ΛCDM model of cosmology. In this mini-review presentation, we assess the possibility that Λ can be a dynamical quantity, more specifically a ‘running quantity’ in quantum field theory in curved spacetime. A great deal of phenomenological works have shown in the last few years that this option (sometimes accompanied with a running gravitational coupling) may cure some of the tensions afflicting the ΛCDM. The ‘running vacuum models’ (RVM’s) are characterized by the vacuum energy density, ρvac, being a series of (even) powers of the Hubble rate and its time derivatives. Here we describe the technical quantum field theoretical origin of the RVM structure in FLRW spacetime, which goes well-beyond the original semi-qualitative renormalization group arguments. In particular, we compute the renormalized energy-momentum tensor using the adiabatic regularization procedure and show that it leads to the RVM form. In other words, we find that the renormalized vacuum energy density, ρvac(H) evolves as a (constant) additive term plus a leading dynamical components O(H2). There are also O(H4) contributions, which can be relevant for the early universe. Remarkably enough, the renormalized ρvac(H) does not exhibit dangerous terms proportional to the quartic power of the masses (∼ m4) of the fields. It is well-known that these terms have been the main source of trouble since they are responsible for the extreme fine tuning and ultimately for the cosmological constant problem. In this context, the current ρvac(H) is dominated by a constant term, as it should be, but it acquires a mild dynamical component ∼ νH2 (|ν| 1) which makes the RVM to mimic quintessence.

Keywords: Running Vacuum, Dynamical Dark Energy, Inﬂation, Dark Matter

1. Introduction The standard or concordance model of cosmology, the so-called ΛCDM, is a very successful theoretical framework for the description of the Universe 1. Crucial ingre-

arXiv:2109.12086v1 [gr-qc] 24 Sep 2021 dients of it are left, though, without direct observational evidence and/or a proper theoretical interpretation based on fundamental principles. Such is the situation with dark matter, but also with the cosmological term, Λ, which has been traditionally associated with the vacuum energy density (VED) of the Universe, ρvac. In fact, the notion of VED in cosmology is a most subtle concept, which is challenging theoretical physicists and cosmologists for many decades, specially with the advent of Quantum Theory in general and the more sophisticated machinery of Quantum Field Theory (QFT). The roots of the problem reside in the interpretation of the arXiv:2109.12084v1 [nucl-th] 24 Sep 2021 hr h aktretno field tensor three rank the where tensor qal ai arb en faPed-ag Transfor- another Pseudo-Gauge [22]: generate a (PGT) of always mation means can by one pair valid operators, equally of pair one edti hiei o nqein unique not is choice this deed n atclrfr fteqatmenergy-momentum quantum (EMT) the of tensor form the particular of choos- a impact by ing decompositions obtained the are different current also momentum terms, boost-angular hydrodynamics. and general relativistic [9–16] more in [17–21] tensor In potential spin spin a of of role the ine a equilibrium. of thermal could predictions local part at polarization spin fluid the spin and of the we orbital decomposition affect a side, significantly an of into theory choice momentum the the angular from that However, ignore exper- not 8]. and can [7, theory between data agreement imental 6] a the that [5, found restore coupling was could shear it Recently, spin-thermal review. overlooked a previously for [4] see polariza- [3], momentum-dependent global tion on fails data of the but measurements reproduce polarization, to the spin hydrody- momenta) explain all relativistic over to (integrated of able are Predictions in for namics physics possibility spin fluids. the of relativistic opened investigations relativistic has phenomenological in 2] new [1, formed collisions (QGP) ion Plasma heavy Gluon Quark the ∗ [email protected] o hsrao,rcn tde aebgnt exam- to begun have studies recent reason, this For in particles of polarization spin of measurement The Φ suoguedpnec ftesi oaiainado the of and polarization spin the of dependence Pseudo-gauge b λ,µν S b Z ′ T b b λ,µν ′ µν utstsytefloigsymmetries: following the satisfy must = = = − banta vnteailvria ffc odciiydepe conductivity effect express vortical D explicit axial non-interacting the the a even calculate of that also function obtain me We Wigner polarization transformations. the the t for PG that to predictions show shear the we thermal affects eral, of PG contribution a pot of the spin choice f affect the We choices with PG. vorticity PG Belinfante thermal the Other the in replacing obtained by one spi the written The from momentum. differs angular PG (PGs) spin ical and Pseudo-Gauges orbital different of in decomposition (LTE) Equilibrium Thermal S T Φ b b b eoti h ensi oaiainvco fsi / parti 1/2 spin of vector polarization spin mean the obtain We T λ,µν µν λ,νµ b .INTRODUCTION I. eateto hsc n srnm,Iw tt University State Iowa Astronomy, and Physics of Department µν + , − n ftesi tensor spin the of and 2 1 Φ ∇ b Z b λ,µν λ µν,λρ Φ b + λ,µν = ∂ ρ − special Z b − Z µν,λρ b n h akfour rank the and Φ b Φ b νµ,λρ µ,λν eaiiy Given relativity. = − − Φ b S ν,λµ Z b b λ,µν µν,ρλ In- . .Buzzegoli M. . rmtelcltemleulbimsaitcloperator: statistical equilibrium thermal local starting the obtained hyper- are from lambda collisions heavy-ion of in In- polarization produced ons spin reported. for been predictions yet deed, not in have polarization spin pseudo-gauges for different formulae explicit literature, the in see equilibrium, only thermal contribution global [27]. a of also out bring is might where system they the statistical [23–26], when that in here shown PGTs analyzed was was of it theory role field The dif- quantum to to results. lead opposed would physical as PGTs ferent the - therefore local quantities, constraining - global by equa- obtained hydrodynamic are relativistic tions ir- the angular physically However, as and regarded relevant. momentum generally energy, are of they total equation momentum, the the affect and not motions do transformations such Since where ihacrepnigvnsigsi esradwhere EMT and velocity of fluid tensor the spin form is vanishing Belinfante corresponding a the with choosing by obtained eutn pnplrzto fapril ihmomentum with particle a k of polarization spin resulting s[,2,29] 28, [5, is with n mass and S ept h atta hsisewsaraydiscussed already was issue this that fact the Despite S ξ µ ̟ µ ∗ ( ε ( k k k do h PG. the on nd = ) nili h sa oaiainformula. polarization usual the in ential ρ b = ) oso h enailcretadwe and current axial mean the of ions = LTE srdi ev-o olsos ngen- In collisions. heavy-ion in asured B oaiainotie ntecanon- the in obtained polarization n esi oaiain hrfr,the Therefore, polarization. spin he √ − rcfil tLEi ffce ythe by affected is LTE at field irac udta hsdffrnecnntbe not can difference this the that ound for choices different with i.e. , − m lsi eaiitcfli tLocal at fluid relativistic a in cles k = 4 8 2 1 m ms oa501 USA 50011, Iowa Ames, , 1 ̟ tfis re ftemdnmcgradients thermodynamic of order first at m + Z ξ 1 S ǫ µν µν ǫ µλστ B m µ exp µρστ ( u = = 2 k ) iie ytetemperature the by divided and 2 1 − k k ≃ − τ ε ( τ ∂ k k 1 2 S R Z µ ρ Σ ( ̟ µ β xa otcleffect vortical axial ∂ R dΣ dΣ ( ν µ Σ k β + dΣ + ) µ ν · R ∂ n k − Σ ν · T S b β R dΣ n k ∂ B ξ µ µν F Σ µ ν ( ) β (1 dΣ k β F · , µ ) ν n k (1 , ) − − · , n k − n F b j F µ n ) F ζ F ̟ ) ρσ t ˆ λ T ξ , The . ρσ (3) (2) (1) (4) T , b β B arXiv:2109.11831v1 [gr-qc] 24 Sep 2021 hsc,wihcni unb npicpedmntae nex in demonstrated principle in be turn under in si can deeper very which creates conceptually physics, analogy This hole [2]. sonic-to-black [1] BH the a plore of th analogue kinetic fluid a the on solely to hinges that methodology th novel discuss a and within computations the perform We determin geometry, propagation. effective and acousti an static to related a as of mechanism horizon the i the performed at investigation emitted the spontaneously on are report b we standard paper of present horizon the event the o of presence analogue the by hydrodynamic characterized the is fluid a of flow transonic The Introduction 1 1 asm Mannarelli Massimo NN aoaoiNzoaidlGa as,VaG Acitelli, G. Via Sasso, Gran del Nazionali Laboratori INFN, 2,3,4 2,3,4 iatmnod iia nvri` iPs,Pl Fibonacc Polo Pisa, Università di Fisica, di Dipartimento NNSzoed ia ooFbnci ag .Pneov 3, Pontecorvo B. Largo Fibonacci, Polo Pisa, di Sezione INFN hnngs ti elsie ob xeddt tnadblack standard to extended be to suited systems. statistica well out-of-equilibrium the is on it t and on gas, horizon depends phonon acoustic only the method of our properties Since cal a this holes. the doing black to proportional standard In is in hole gas. as acoustic phonon the of the entropy of the gain that cor entropy the loss entropy the obtain the with we equating horizon by hole, temperature acoustic Hawking the the of by sion emitted gas phonon the kinetic covariant Using distribu phonons. the generated with spontaneously horizon the hole of black properties standard geometrical of the analogue hydrodynamic the holes, tic hnneiso yaosi lc holes black acoustic by emission Phonon epeetanvlitrrtto fteHwigtemperatur Hawking the of interpretation novel a present We 1 ai Grasso Dario , Chiofalo Abstract Italy 1 2 ivaTrabucco Silvia , 4 2 -70 seg (AQ), Assergi I-67100 22, ,LroB otcro3 Pontecorvo B. Largo i, infnto fthe of function tion tnigo h horizon the of standing oeadw interpret we and hole c e ftehorizon, the of rea rpriso the of properties l n hnnproduction phonon ing hoyt describe to theory oyapoc applied approach eory 1 nhwphonons how on [1] n eiet urnl at currently periments 3 naosi horizon, acoustic an f ,b connecting by s, neligphysics underlying e n ai Luisa Maria and , akhls(H) In (BHs). holes lack pemne oex- to manner mple fteacoustic the of 62 ia Italy Pisa, 56127 egeometri- he oe n to and holes eassume we , etexpres- rect facous- of e Viable Intermediate Inflation in the Mimetic DBI Model

Narges Rashidia,b1 and Kourosh Nozaria,b2

a Department of Theoretical Physics, Faculty of Science, University of Mazandaran, P. O. Box 47416-95447, Babolsar, IRAN

b ICRANet-Mazandaran, University of Mazandaran, P. O. Box 47416-95447, Babolsar, IRAN

Abstract We study the intermediate inflation in the mimetic Dirac-Born-Infeld model. By considering β the scale factor as a = a0 exp(bt ), we show that in some ranges of the intermediate parameters b and β, the model is free of the ghost and gradient instabilities. We study the scalar spectral index, tensor spectral index, and the tensor-to-scalar ratio in this model and compare the results with Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data at 68% and 95% CL. In this regard, we find some constraints on the intermediate parameters that lead to the observationally viable values of the perturbation parameters. We also seek the non-gaussian features of the primordial perturbations in the equilateral configuration. By performing the numerical analysis on the nonlinearity parameter in this configuration, we show that the amplitude of the non- gaussianity in the intermediate mimetic DBI model is predicted to be in the range −16.7 < f equil < −12.5. We show that, with 0 < b ≤ 10 and 0.345 < β < 0.387, we have an instabilities- free intermediate mimetic DBI model that gives the observationally viable perturbation and non-gaussianity parameters. arXiv:2109.11829v1 [astro-ph.CO] 24 Sep 2021

[email protected] (Corresponding Author) [email protected]

1 arXiv:2109.11794v1 [nucl-th] 24 Sep 2021 a tt 71] nsc ae npicpe vnthe that even with principle, Indeed, Λ pp in the extracted. case, mind in a be aim fi- such can the In lengths in scattering momentum [7–12]. relative two state low the nal where with reactions emerge in baryons interaction their studying by n h orsodn w-oymmnu correlation momentum two-body corresponding the ing N Y did end, results. the conclusive at anticipated and, the statistics provide low not from suffered [14] tion Λ strength the spin-triplet re- determine the experiment to but of double-polarization efforts Unfortunately, reaction, even 16]. that or [15, single- in and/or a spin-triplet possible de- the quires is spin of the amplitude, separation of singlet the determination A i.e. pendence, 14]. Jülich [13, in ator htn rprbascnb rprdi re oper- the to on order information of addi- in essential source prepared An tional be experiments. so scattering can standard hyperons beams form the proper of no character that short-lived the of the sequence of dependence spin N Λ Y the the only about for specifically far and those formation, so established, all experimentally However, of are first properties decades. bulk many the for studied been e,eg,[]adtercn eiw 5 ] nfc,tein- the fact, In 6]. ( [5, reviews hyperons recent of objects, the teractions of and compact generation [4] such the e.g., of see, for mergers is particular in latter in waves of properties The gravitational and studies the stars, by understanding matter. neutron for driven of neutron relevance side of of state other world- especially of the facilities equation at on hypernuclei the and of [1–3] study the wide by side the one on triggered recently, interactions hyperon-nucleon in ∗ † [email protected] [email protected] oercnl,i a elzdta nomto nthe on information that realized was it recently, More Introduction: → cteiglntscnas eotie rmstudy- from obtained be also can lengths scattering ocs scmltl akn.Telte sacon- a is latter The lacking. completely is forces, K xlrn h Σ the Exploring + Λ p 2 rt h ige n h triplet the and singlet comb the collisions heavy-ion arate or proton-proton high-energy in nw.W hwta ne esnbeasmtos measure a assumptions, reasonable under that show We known. nttt o dacdSmlto,Isiu für Kernphys Institut Simulation, Advanced for Institute a enmaue tteCS acceler- COSY the at measured been has n ¨lc etrfrHdo hsc,FrcugznrmJ Forschungszentrum Physics, Hadron for Jülich Center and h pndpnetcmoet ftefnaetlhyperon- fundamental the of components spin-dependent The 1 emot ntttfu taln n enhskadBethe and Kernphysik und für Strahlen- Institut Helmholtz p hr a enarnwlo interest of renewal a been has There nlsaeitrcini h reaction the in interaction final-state p neato yteTFcollabora- TOF the by interaction Y Λ = N + ytm oedtie in- detailed More system. 1 , nttt o dacdSmlto,Isiu für Kernphys Institut Simulation, Advanced for Institute p Σ oshnsetu ¨lc,D545Juih Germany Jülich, D-52425 Jülich, Forschungszentrum , neato ymaueet ftecreainfunction correlation the of measurements by interaction )wt uloshave nucleons with Ξ) N Y 3 bls tt nvriy 16Tiii Georgia Tbilisi, 0186 University, State Tbilisi nvri¨tBn,D515Bn,Germany Bonn, D-53115 Bonn, Universität oc sprovided is force S wv contributions. -wave .G Meißner U.-G. .Haidenbauer J. eto aaa o nrishv enaalbefra for available isospin has been channel have This energies [22]. time low long at data section section. cross the elastic on the data on has and This one function disentan- provided section. correlation contributions, to cross spin one reaction the allows gle the interplay of al- different case they qualitatively in whereas up amplitude, scattering add the ways of depend sign they the because a corre- on cancel, the such partially to in will states repul- below, function spin lation is detail the spin more one of contributions in the other the discussed situation the of be in will one As that in sive. while interaction attractive is the specific without states in when - possible namely data is section func- cases, It cross correlation experiment. with spin-dependent the combination any of in - measurements tion from free feasible in is as weight 20]. statistical [17, same assumed the commonly sim- as contributions with scattering, triplet up and add singlet- can ply The dependence spin deduced. the on be information detailed no here ucina esrdi ev-o olsoso high- or collisions heavy-ion in energetic measured as function r aafo h TRClaoain[9 rmamea- a at from collisions [19] Au+Au Collaboration in STAR surement the from data are h LC olbrto in collaboration ALICE the 2]and [20] 2,2]aeasn.A o nris h bevbe are singlet) observables the (spin energies, the channels low by open At from dominated absent. complications are 24] i.e. [23, elastic, purely is 3 h neato nthe in interaction the hr sn opigt h Λ the to coupling no is there esrmnso ( of measurements ilwv fteΣ the of wave tial sso ee hfsadwdh fΣ of widths anal- and the shifts from level came of nature repulsive ysis its for indications Σ eol civdwe h neato nthe can in medium interaction the the in when potential achieved Σ only repulsive a be [26], Ref. in S − eeepiyti dafrteΣ the for idea this exemplify We states spin the of separation a that out point we Here, 1 1 fraini ev uli[5.A icse ndetail in discussed As [25]. nuclei heavy in -formation , ata ae.Teei nmiuu vdnethat evidence unambiguous is There waves. partial 2 k n ¨lc etrfrHdo Physics, Hadron for Jülich Center and ik, 1 ndwt rs eto aaalw osep- to allows data section cross with ined , ∗ 3 † √ pp eto h Σ the of ment s u ih -22 Jülich, Germany ülich, D-52425 ulo neatosaelreyun- largely are interactions nucleon etrfrTertclPhysics, Theoretical for Center olsos[7 8.Wt epc oΛ to respect With 18]. [17, collisions 3TV[1,rsetvl.Hwvr also However, respectively. [21], TeV 13 = + π p − hne srplie oerecently, More repulsive. is channel 3 K , + S p ik, 1 + orlto function correlation ata aei eusv.First repulsive. is wave partial nlsv pcr eae to related spectra inclusive ) pp N olsosat collisions ytmadtereaction the and system √ + 1 S p s 0 ytm hr cross where system, − 0 e n by and GeV 200 = n si triplet) (spin and tm n from and atoms I 3 = √ s / othat so 2 TeV 7 = 3 S p 1 there par- Microphysical manifestations of viscosity and consequences for anisotropies in the very early universe

Chandrima Ganguly1, 2, ∗ and Jerome Quintin3, † 1Department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom 2Wolfson College, Barton Road, Cambridge CB3 9BB, United Kingdom 3Max Planck Institute for Gravitational Physics (Albert Einstein Institute), D-14476 Potsdam, Germany It has been known that a non-perfect fluid that accounts for dissipative viscous effects can evade a highly anisotropic chaotic mixmaster approach to a singularity. Viscosity is often simply param- eterised in this context, so it remains unclear whether isotropisation can really occur in physically motivated contexts. We present a few examples of microphysical manifestations of viscosity in fluids that interact either gravitationally or, for a scalar field for instance, through a self-coupling term in the potential. In each case, we derive the viscosity coefficient and comment on the applicability of the approximations involved when dealing with dissipative non-perfect fluids. Upon embedding the fluids in a cosmological context, we then show the extent to which these models allow for isotropisation of the universe in the approach to a singularity. We first do this in the context of expansion anisotropy only, i.e., in the case of a Bianchi type-I universe. We then include anisotropic 3-curvature modelled by the Bianchi type-IX metric. It is found that a self-interacting scalar field at finite temperature allows for efficient isotropisation, whether in a Bianchi type-I or type-IX spacetime, although the model is not tractable all the way to a singularity. Mixmaster chaotic behaviour, which is well known to arise in anisotropic models including anisotropic 3-curvature, is found to be suppressed in the latter case as well. We find that the only model permitting an isotropic singularity is that of a dense gas of black holes.

I. INTRODUCTION some semi-classical higher-derivative theories of gravity have stable isotropic cosmological singularities (e.g., [11]) While the currently observable universe is isotropic to or can limit the growth of anisotropies [12, 13]. Other a very high degree, this is not a generic feature of space- theories of modified gravity can similarly bound shear times near singularities — rather the opposite. In fact, anisotropies [14] or screen them [15, 16]. Moreover, within general relativity and under certain assumptions in a gravitational ultraviolet-complete theory such as in the matter sector, the most generic approaches to sin- quadratic gravity [17], requiring semi-classical cosmolog- gularities (such as cosmological big bang or big crunch ical transition amplitudes from the big bang to today to singularities, but also black hole singularities) are highly be well defined and finite severely constrains anisotropic anisotropic. They display a chaotic mixmaster behaviour singularities [18, 19]. [1], with infinite chaotic oscillations on a finite time in- In the context of alternative very early universe scenar- terval. This behaviour proceeds in epochs, the beginning ios such as models of bouncing cosmology, the question and end of which are well approximated by the Kasner of the evolution of anisotropies also plays an important metric [2]. This is known as the Belinski-Khalatnikov- role, even well before the approach to the high-curvature Lifshitz (BKL) singularity [3] (see, e.g., [4–6] for reviews). big crunch/bounce singularity (or before a non-singular In the context of the current paradigm of very early bounce occurs). For instance, in a Bianchi type-I uni- universe cosmology, anisotropies from the big bang ini- verse, the contribution from shear anisotropies to the to- arXiv:2109.11701v1 [gr-qc] 24 Sep 2021 tial singularity would be quickly washed away by a pe- tal energy density budget is proportional to 1/a6, where riod of accelerated expansion (i.e. inflation). However, a is the spatially averaged scale factor. While this de- our knowledge of possible pre-inflationary physics is very cays very rapidly in an expanding universe, it conversely scarce, and the question of what happened near the grows much faster than for other known matter types big bang remains of fundamental interest. In particu- (e.g., pressureless dust and radiation), thus representing lar, semi-classical general relativity most likely does not an instability to standard isotropic contracting models. hold anymore at such high energy scales, and what was This is not an issue for ekpyrotic cosmology since the the nature of the initial big bang singularity (if there isotropic background scaling solution arises from a scalar was one) remains an open question, especially whether it field with energy density growing as 1/a2 with > 3, was of BKL type or rather isotropic. Within string the- thus effectively diluting anisotropies. As such, the ekpy- ory, chaotic anisotropies are expected (e.g., [7–10]), while rotic scenario has been shown to be very robust with respect to dynamically producing an isotropic universe, as demonstrated by analytic and numerical studies [20– ∗ [email protected] 25]. In fact, formally, ekpyrosis implies a no-hair theorem † [email protected] stating that the future big crunch is a stable isotropic Prepared for submission to JHEP

The ChPT: top-down and bottom-up

Karol Kampf Institute of Particle and Nuclear Physics, Charles University V Holešovičkách 2, 180 00 Prague 8, Czech Republic E-mail: [email protected]

Abstract: In this work, higher-derivative corrections of the non-linear sigma model of both even and odd intrinsic-parity sectors are systematically studied, focusing on ordered amplitudes of flavor scalars in massless limit. It should correspond to a theory known as Chiral perturbation theory (ChPT) without external sources and with only single-trace operators. We briefly overview its formal development and apply new S-matrix methods to its amplitude constructions. The bottom-up analysis of the tree-level amplitudes of different orders and multiplicities focuses on the formal structure of general ChPT. Possible theoretical simplifications based on the Kleiss-Kuijf and Bern-Carrasco-Johansson relations are presented. Finally, in the same context, the comparison with the so-called Z-function, which is connected with string theory, is also discussed. arXiv:2109.11574v1 [hep-th] 23 Sep 2021 Draft version September 27, 2021 Typeset using LATEX twocolumn style in AASTeX631

On the Accuracy of Underground Muon Intensity Calculations

A. Fedynitch,1, 2 W. Woodley,3 and M.-C. Piro3

1Institute for Cosmic Ray Research, the University of Tokyo, 5-1-5 Kashiwa-no-ha, Kashiwa, Chiba 277-8582, Japan 2Institute of Physics, Academia Sinica, Taipei, 11529, Taiwan 3Department of Physics, University of Alberta, Edmonton, Alberta T6G 2E1, Canada

(Dated: September 27, 2021)

ABSTRACT Cosmic ray muons detected by deep underground or underwater detectors have served as an information source on the high-energy cosmic ray spectrum and hadronic interactions in air showers for almost a century. The theoretical interest in underground muons has nearly faded because space-borne experiments probe the cosmic ray spectrum more directly, and accelerators provide precise measurements of hadron yields. However, underground muons probe unique hadron interaction energies and phase space, which are still inaccessible to present accelerator experiments. The cosmic ray nucleon energies reach the hundred-TeV and PeV ranges, which are barely accessible with space-borne experiments. Our new calculation combines two modern computational tools, mceq, for surface muon fluxes, and proposal, for underground transport. We demonstrate excellent agreement with measurements of cosmic ray muon intensities underground within estimated errors. Beyond that, the precision of historical data turns out to be significantly smaller than our error estimates. This result shows that the sources of high-energy atmospheric lepton flux uncertainties at the surface or underground can be significantly constrained without taking more data or building new detectors. The reduction of uncertainties can be expected to impact data analyses at large-volume neutrino telescopes and for the design of future ton-scale direct dark matter detectors.

1. INTRODUCTION characterise the prompt flux and to search for its sig- For more than 80 years, it has been known that muons nature in data. With the advent of large-volume neu- observed at the surface or by experiments underground trino telescopes, the searches have been extended from reflect the properties of high-energy cosmic rays inter- facilities underground to those underwater and under- acting with the Earth’s atmosphere (Barrett et al. 1952; ice (Gaisser et al. 1995; Gondolo et al. 1996; Sinegov- Bugaev et al. 1998; Gaisser et al. 2016; Zyla et al. 2020). skaya & Sinegovsky 2001). All attempts to detect the As the topic evolved over many decades, the theoreti- flux of prompt muons or neutrinos remain unsuccessful, cal activity focused on searches of prompt atmospheric and only a few extreme models have been excluded, e.g. muons from decays of charm or heavier flavor mesons by LVD (Aglietta et al. 1999) or recently by IceCube and the characterisation of the cosmic ray composition (Aartsen et al. 2016a,b). These current results are in at energies around the knee of the cosmic ray spectrum good agreement with the conclusions of BMN. Another field with an interest in accurate deep under- arXiv:2109.11559v1 [astro-ph.HE] 23 Sep 2021 (Elbert et al. 1983; Bugaev et al. 1989; Hoerandel 2003; Becker Tjus & Merten 2020). ground muon calculations is direct detection dark mat- One of the latest theoretical analyses of underground ter (DM) experiments. Cosmic ray muons are a signif- muon intensities was performed 20 years ago by Bugaev, icant background for DM searches (see, e.g. Formaggio Misaki, and Naumov (BMN) (Bugaev et al. 1998) to & Martoff(2004) for a review), producing secondary neutrons in interactions with rock or surrounding ma- terials. These muon-induced neutrons pose challenges [email protected] for the direct detection of DM nuclear recoils since they [email protected] can mimic the expected energy spectrum from WIMP- nucleus (Weakly Interacting Massive Particles) scatter- [email protected] ing. By installing DM experiments in deep underground UTTG-18-2021 NORDITA-2021-078 Natural Chain Inflation

Katherine Freese,1, 2, 3, ∗ Aliki Litsa,2, † and Martin Wolfgang Winkler1, 2, ‡ 1Department of Physics, The University of Texas at Austin, Austin, 78712 TX, USA 2Stockholm University and The Oskar Klein Centre for Cosmoparticle Physics, Alba Nova, 10691 Stockholm, Sweden 3Nordita, KTH Royal Institute of Technology and Stockholm University, Hannes Alfvénsväg12, 10691 Stockholm, Sweden In Chain Inflation the universe tunnels along a series of false vacua of ever-decreasing energy. The main goal of this paper is to embed Chain Inflation in high energy fundamental physics. We begin by illustrating a simple effective formalism for calculating Cosmic Microwave Background (CMB) observables in Chain Inflation. Density perturbations seeding the anisotropies emerge from the probabilistic nature of tunneling (rather than from quantum fluctuations of the inflation). To obtain the correct normalization of the scalar power spectrum and the scalar spectral index, we 1/4 12 find an upper limit on the scale of inflation at horizon crossing of CMB scales, V∗ < 10 GeV. We then provide an explicit realization of chain inflation, in which the inflaton is identified with an axion in supergravity. The axion enjoys a perturbative shift symmetry which is broken to a discrete remnant by instantons. The model, which we dub ‘natural chain inflation’ satisfies all cosmological constraints and can be embedded into a standard ΛCDM cosmology. Our work provides a major step towards the ultraviolet completion of chain inflation in string theory.

I. INTRODUCTION sity perturbations in the primordial plasma which later manifest as anisotropies in the CMB. 4 Cosmic inflation, an early accelerated epoch of the Uni- It was previously shown that for ∼10 phase tran- verse, was proposed to solve several conundrums of the sitions per e-fold of inflation, the correct amplitude Hot Big Bang [1], and in addition provides density fluctu- of CMB temperature fluctuations is obtained [10, 11]. ations that serve as seeds for the formation of large scale The large number of transitions automatically ensures structure. In recent years, the paradigm most frequently that vacuum bubbles from individual nucleation sites studied consists of a scalar field, the inflaton, which rolls quickly percolate, thus evading the “empty universe down its potential and drives the exponential expansion problem” [12] of old inflation. A nearly scale-invariant of space [2, 3]. Quantum fluctuations of the inflaton seed scalar power spectrum – as preferred by observation – is the temperature anisotropies we observe in the CMB [4]. obtained if the Hubble rate H and the tunneling rate Γ However, there exist less known, though equally suc- vary (at most) slowly from vacuum to vacuum [11]. cessful theories of inflation, in which the rolling of the In this work we provide simple analytic expressions scalar field is replaced by its tunneling from higher en- to determine the CMB observables in concrete models of ergy metastable vacua to lower energy minima. During chain inflation. The formalism is based on our previously the time spent in the higher energy minimum, the poten- derived approximation of tunneling rates which replaces tial energy drives accelerated expansion. Guth’s original the thin-wall approximation for generic quasiperiodic po- “old inflation”, with a single tunneling event, was the tentials [11]. original tunneling model, but it suffered from a failure to The main goal of this paper is to illustrate an explicit reheat after inflation: whereas pockets of the Universe do supergravity realization of chain inflation in which the undergo the required first order phase transition, most of inflaton is identified with an axion. Two instanton terms the Universe remains in the false vacuum and the bubbles break the axionic shift symmetry and induce a periodic of true vacuum never percolate or reheat. potential as in Eq. (24). Such a setting carries profound Several solutions to this reheating problem have been motivation from string theory in which axions with mul- proposed, including Double Field Inflation [5, 6] (a sec- tiple instanton terms are ubiquitous. The leading instan- ond rolling field gives rise to time dependence of the tun- ton term generates an overall cosine shape of the axion neling rate) and Chain Inflation, which is the subject of potential, while the subleading instanton induces a se- this paper. In Chain Inflation [7–9], the Universe under- ries of metastable minima by causing small wiggles in arXiv:2109.11556v1 [hep-th] 23 Sep 2021 goes a series of first order phase transitions. Instead of a the potential. We show that the model allows for a suc- rolling scalar field, the inflaton tunnels along a series of cessful regime of chain inflation in which all cosmological metastable minima in its potential. While quantum fluc- constraints are satisfied. tuations are suppressed due to the inflaton mass in each vacuum, the probabilistic nature of tunneling causes den- II. CMB OBSERVABLES

In this section we describe a quick method to derive the ∗ [email protected] CMB observables of chain inﬂation which solely relies on † [email protected] the determination of the extrema in the potential. We ‡ [email protected] consider a (quasi)periodic potential V (ϕ) containing a se- SciPost Physics Submission

Determination of the rapidity evolution kernel from Drell-Yan data at low transverse momenta

F. Hautmann1, I. Scimemi2 and A. Vladimirov3

1 University of Oxford / Universiteit Antwerpen 2 Universidad Complutense de Madrid 3 Universität Regensburg

Proceedings for the XXVIII International Workshop on Deep-Inelastic Scattering and Related Subjects, Stony Brook University, New York, USA, 12-16 April 2021 doi:10.21468/SciPostPhysProc.?

Abstract

We report recent results on the extraction of the rapidity evolution kernel for Sudakov processes from experimental measurements of Drell-Yan (DY) transverse momentum distributions. We discuss the role of high-precision Large Hadron Collider (LHC) data, and illustrate the interplay of this analysis with transverse momentum dependent and collinear parton density functions.

The experimental program of high-precision measurements in DY lepton-pair production at the LHC influences determinations of the rapidity evolution kernel for Sudakov processes in QCD, complementing the information which may be obtained from lower-energy DY experiments as well as semi-inclusive deep inelastic scattering. The kernel can be extracted from precise measurements of the DY transverse momentum distribution in the region of low transverse momenta. In this article we discuss recent results on this based on our work [1]. The rapidity evolution kernel was introduced by Collins and Soper [2,3 ] in the context of transverse momentum dependent (TMD) factorization [4] for the DY differential cross section at arXiv:2109.12051v1 [hep-ph] 24 Sep 2021 low transverse momenta [5–7]. It controls the evolution in rapidity of partonic TMD distribution functions [8] as a function of mass µ and transverse coordinate bT . Rapidity divergences are treated by cut-off methods [2,4 ] or by subtraction methods [9–12]. The kernel’s rate of change with mass is given by the perturbatively-calculable cusp anomalous dimension (currently known > 1 numerically to four loops [13, 14]). For large transverse distances bT Λ−QCD the kernel is nonperturbative. ∼ Information on the non-perturbative rapidity evolution kernel may be obtained by lattice QCD methods. Studies in this direction have recently started [15–20]. Alternatively, it may be obtained by comparison of theory with experimental measurements. This has a long history, from the pioneering ‘ResBos’ fits [21–24] of Tevatron data to the recent ‘artemide’ [25] and ‘NangaParbat’ [26] global fits to LHC, Tevatron, RHIC and fixed-target data. The large-distance behavior of the rapidity evolution kernel, currently subject of active in- 2 vestigation, has traditionally been modeled as a quadratic behavior bT . The fits [21–24], for ∝ 1 Anomalous diffusion in QCD matter

Paul Caucal1, ∗ and Yacine Mehtar-Tani1, 2, † 1Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 2RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA We study the effects of quantum corrections on transverse momentum broadening of a fast parton passing through dense QCD matter. We show that, at leading logarithmic accuracy the broadening distribution tends at late times or equivalently for large system sizes L to a universal distribution 2 2 that only depends on a single scaling variable k⊥/Qs where the typical transverse momentum scale ln Q2 ' (1 + 2β) ln L − 3 (1 + β) ln ln L increases with time as s√ 2 up to non-universal terms, with an anomalous dimension β ∼ αs. This property is analogous to geometric scaling of gluon distributions in the saturation regime and traveling waves solutions to reaction-diffusion processes. We note that since β > 0 the process is super-diffusive, which is also reﬂected at large transverse momentum 4−2β where the scaling distribution exhibits a heavy tail k⊥ akin to Levy´ random walks.

Transverse momentum broadening (TMB) of ener- extends over a large range of transverse momenta above getic quarks and gluons while traversing QCD matter the typical transverse momentum scale. plays a central role in a variety of processes studied at Levy´ flights are ubiquitous in nature and span a colliders to probe QCD, ranging from jet suppression in wide variety of stochastic processes in biological sys- heavy ion collisions [1,2] to transverse momentum de- tems [11, 12], molecular chemistry [13], optical lattice pendent gluon distribution function that encode infor- [14], turbulent diffusion and polymer transport theory mation on the 3D structure of the proton and nuclei in [15, 16]. Furthemore, heavy tailed distributions are also high energy collisions in particular at small Bjorken x observed in self-organized critical states [17, 18]. In this where gluon saturation is expected to take place [3–5]. work, we point out for the first time another occurrence High energy partons experience random kicks in hot of such random walks in the transport of eikonal par- or cold nuclear matter that cause their transverse mo- tons in dense QCD matter and we compute the anoma- mentum k⊥ w.r.t. to their direction of motion to increase lous exponents that characterize the deviation to stan- over time. The leading order elastic process is given by a dard diffusion. single gluon exchange via Coulomb scattering and leads to an approximate brownian motion in transverse momentum space, so long as the mean free path is larger QUANTUM CORRECTIONS TO TRANSVERSE than the in-medium correlation length, where the typ- MOMENTUM BROADENING IN QCD MEDIA ical transverse momentum square scales linearly with 2 system size L, namely k typ qˆ L, where qˆ is the dif- h ⊥i ∼ The TMB distribution is related to the forward scat- fusion coefficient. tering amplitude (x⊥) of an effective dipole in color S Moreover, radiative processes can also increase trans- representation R = A, F with transverse size x⊥(see e.g. verse momentum of the leading parton due to recoil ef- [7, 19, 20]) via a Fourier transform, fects. This question was recently addressed, and it has Z 2 −ik⊥·x⊥ been shown that such contributions, albeit suppressed (k⊥) = d x⊥e (x⊥) . (1) by the coupling constant αs, are enhanced by double P S

arXiv:2109.12041v1 [hep-ph] 24 Sep 2021 logarithms which must be resummed to all orders when 2 Considering the dipole formulation in position space αs ln L 1 [6–9]. ∼ allows for a straightforward resummation of multiple In this letter we go beyond this result by investigating scattering by exponentiating the single scattering cross- in more detail the consequences of the non-local nature section, so long as the interactions between the dipole of the quantum corrections on the TMB distribution. We and the medium are local and instantaneous. Thus, we 1 CR 2 2 ﬁnd in particular that the latter reaches a universal scal- may write (x⊥) = exp( qˆ(1/x )L x ) [21]. The S − 4 Nc ⊥ ⊥ ing solution at late times (large L) that we compute an- latter relation deﬁnes the quenching parameter in the alytically along with its sub-asymptotic deviations ex- adjoint representation which is assumed to be a slowly ploiting a formal analogy between the present problem varying function of x⊥. At tree-level it reads and traveling wave solutions to reaction-diffusion processes [10]. As a consequence of the self-similarity char- (0) 2 1 qˆ (1/x⊥) =q ˆ0 ln 2 2 , (2) acterizing the anomalous random walk, the TMB distri- x⊥µ bution is of Levy´ type. It is in particular associated with a heavy power law tail describing long rare steps which up to powers of xT µ suppressed terms. For a weakly coupled QGP the bare quenching parameter qˆ0 and the Qualitative analysis of proton inelastic scattering for diquark searching

Vladimir V. Bytev,a1 Stepan. S. Shimanskiya

a Joint Institute for Nuclear Research, 141980 Dubna (Moscow Region), Russia

Abstract In this paper we discuss exclusive reactions which analysis can be used to receive direct indication of diquark existence. We make estimations of diquark scattering process measurement in inelastic proton-proton collisions. It was shown that putting special restrictions over kinematics and particles in ﬁnal state of process it will be possible to enhance potential diquark contribution to scattering up to 104. We put qualitative characteristics of process with diquark and ways to distinguish it from quark scattering in model-independent way. arXiv:2109.12025v1 [hep-ph] 24 Sep 2021

1E-mail: [email protected] IPPP/21/31, ZU-TH 45/21

Next-to-leading order QCD corrections to diphoton-plus-jet production through gluon fusion at the LHC

Simon Badgera, Thomas Gehrmannb, Matteo Marcolib, Ryan Moodiec

aDipartimento di Fisica and Arnold-Regge Center, Universitàdi Torino, and INFN, Sezione di Torino, Via P. Giuria 1, 10124 Torino, Italy bPhysik-Institut, UniversitätZürich, Wintherturerstrasse 190, CH-8057 Zürich, Switzerland cInstitute for Particle Physics Phenomenology, Department of Physics, Durham University, South Road, Durham, DH1 3LE, UK

Abstract We compute the next-to-leading order (NLO) QCD corrections to the gluon-fusion subprocess of diphoton-plus-jet production at the LHC. We compute fully diﬀerential distributions by combining two-loop virtual corrections with one-loop real radiation using antenna subtraction to cancel infrared divergences. We observe signiﬁcant corrections at NLO which demonstrate the importance of combining these corrections with the quark-induced diphoton-plus-jet channel at next-to-next- to-leading order (NNLO). Keywords: NLO corrections

1. Introduction Recent breakthroughs in two-loop amplitude technology are opening up a new range of precision two-to-three scattering problems. Diphoton-plus-jet production has been one of the ﬁrst predictions to appear at NNLO in QCD [1]. This progress is extremely timely given the continually improving experimental measurements of diphoton signatures [2]. Predictions for pure diphoton production have been known to NNLO accuracy for some time [3–7] and the steps towards N3LO are being taken with the completion of the three-loop amplitudes [8]. Diphoton-plus-jet signa-

arXiv:2109.12003v1 [hep-ph] 24 Sep 2021 tures are of particular importance at the LHC since they form the largest background to Higgs production at high transverse momenta. The extra jet is necessary to ensure a non-zero transverse momenta in the diphoton system. The recently computed NNLO corrections [1] of diphoton-plus-jet production display a good perturbative convergence, except in regions where the loop-mediated gluon-fusion process (which contributes to the cross section only from NNLO onwards) is numerically sizable compared to other contributions. In order to capture the full eﬀects of the QCD corrections, it is important to include loop-induced gluon-fusion channels from at least one order higher in the perturbative series. These corrections are the subject of this article. Preprint submitted to Physics Letters B September 27, 2021 An alternative scheme for eﬀective range corrections in pionless EFT

M. Ebert,1, ∗ H.-W. Hammer,1, 2, † and A. Rusetsky3, 4, ‡ 1Technische UniversitätDarmstadt, Department of Physics, 64289 Darmstadt, Germany 2ExtreMe Matter Institute EMMI and Helmholtz Forschungsakademie Hessen fürFAIR (HFHF), GSI Helmholtzzentrum fürSchwerionenforschung GmbH, 64291 Darmstadt, Germany 3HISKP and BCTP, Rheinische Friedrich-Wilhelms UniversitätBonn, 53115 Bonn, Germany 4Tbilisi State University, 0186 Tbilisi, Georgia Abstract We discuss an alternative scheme for including effective range corrections in pionless effective field theory. The standard approach treats range terms as perturbative insertions in the T -matrix. In a finite volume this scheme can lead to singular behavior close to the unperturbed energies. We consider an alternative scheme that resums the effective range but expands the spurious pole of the T -matrix created by this resummation. We test this alternative expansion for several model potentials and observe good convergence. arXiv:2109.11982v1 [hep-ph] 24 Sep 2021

∗ [email protected] † [email protected] ‡ [email protected]

1 SciPost Physics Submission

MCNET-21-13

Accelerating Monte Carlo event generation – rejection sampling using neural network event-weight estimates

K. Danziger1, T. Janßen2, S. Schumann2, F. Siegert1

1 Institut fürKern- und Teilchenphysik, TU Dresden, Dresden, Germany 2 Institut fürTheoretische Physik, Georg-August-UniversitätGöttingen,Göttingen, Germany

September 27, 2021

Abstract

The generation of unit-weight events for complex scattering processes presents a severe challenge to modern Monte Carlo event generators. Even when using sophisticated phase-space sampling techniques adapted to the underlying transition matrix elements, the efficiency for generating unit-weight events from weighted samples can become a limiting factor in practical applications. Here we present a novel two-staged unweighting procedure that makes use of a neural-network surrogate for the full event weight. The algorithm can significantly accelerate the unweighting process, while it still guarantees unbiased sampling from the correct target distribution. We apply, validate and benchmark the new approach in high-multiplicity LHC production processes, including Z/W +4 jets and tt¯+3 jets, where we find speed-up factors up to ten. arXiv:2109.11964v1 [hep-ph] 24 Sep 2021

1 5.1 Introduction

A theory vade mecum for PSI experiments

G. Colangelo1, F. Hagelstein2, A. Signer2,3? and P.Stoffer4

1 Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, University of Bern, Switzerland 2 Paul Scherrer Institut, 5232 Villigen PSI, Switzerland 3 University of Zurich, Physik-Institut, 8057 Zurich, Switzerland 4 University of Vienna, Faculty of Physics, Boltzmanngasse 5, 1090 Vienna, Austria

? [email protected]

Review of Particle Physics at PSI doi:10.21468/SciPostPhysProc.5

Abstract

This article gives a compact introduction and overview of the theory underlying the experiments described in the rest of this review.

Copyright G. Colangelo et al. Received 17-05-2021 This work is licensed under the Creative Commons Accepted 21-07-2021 Check for Attribution 4.0 International License. Published 06-09-2021 updates Published by the SciPost Foundation. doi:10.21468/SciPostPhysProc.5.005

5.1 Introduction The purpose of this article is to give a broad overview of the theory background to the experiments that have been and are carried out at the Paul Scherrer Institute. Space limitations make it

arXiv:2109.11948v1 [hep-ph] 24 Sep 2021 impossible to go into depth or provide a self-contained theoretical summary. Much more modestly, we aim to put the experiments into context and provide key references for further reading. The experiments we refer to are listed in Table 5.1 and they will be described in greater detail in separate sections/articles of the Review of Particle Physics at PSI [1–23]. These experiments either lead to precise determinations of physical parameters required as input for other experiments (e.g., muon life time, pion mass), or search for physics beyond the Standard Model (BSM). The BSM searches proceed along different frontiers. One way to search for new physics is to consider physical observables whose Standard Model (SM) contributions either vanish or are too small to be experimentally accessible. In other words, they are identical to zero for practical purposes. Examples are charged lepton-ﬂavor violating (cLFV) muon decays or a permanent neutron electric dipole moment (EDM). To put constraints on the branching ratios of BSM decays, one has to

1 Prepared for submission to JHEP

Conversions of propagation eigenstates of supernova neutrinos by atomic electrons

Motohiko Kusakabea,1 aSchool of Physics, and International Research Center for Big-Bang Cosmology and Element Gen- esis, Beihang University 37, Xueyuan Rd., Haidian-qu, Beijing 100083 China

E-mail: [email protected]

Abstract: Electron number densities in stars and the Earth are inhomogeneous because of atomic electrons. The large inhomogeneities on atomic-scale tend to form at tops of respective layers of stars, and 1s electrons of O locally produce weak potentials higher than that of the high MSW resonance. Then, supernova neutrinos experience vast numbers of non-adiabatic transitions. This inhomogeneous electron potential generates finite amplitudes of all three propagation eigenstates, and wave packets effectively separate. Then, spectral differences between three flavors significantly diminishes after propagation. arXiv:2109.11942v1 [hep-ph] 24 Sep 2021

1Corresponding author. arXiv:2109.11932v1 [hep-ph] 24 Sep 2021 cin([].Frteitrcino on rjciewt ag nucle large a with projectile point a colors of of number interaction large e the a the For of in approximation summarized [2]). and diagr ( a (’reggeons’) the [1] action gluons of e.g. summation reggeized (see on of approach based interaction JIMWLK approach BFKL-Bartels so-called su the the the of been of and has context QCD the the in in behavior both high-energy the years many For Introduction 1 nlrpdte n omc oedffiut silsrtdb ut e at is few quite collisions by illustrated AB c 10]. as boundary difficult, 9, for with more [8, much equation one so the analogous f and rather generalized rapidities but the final been rapidity equation in has evolution equation BK for and this the equation 7] functions 6, Unlike Green [5, BFKL nuclei. papers the our which In with 4], vertex. [3, constructed diagrams equation (BK) pomeron Balitski-Kovchegov the equation, closed culyteiflec fteqarpepmrninteraction O- pomeron for quadruple estimates Numerical the of small. 25. influence to the up actually rapidities and 0.2 order ihsml rigfntosalw ofidtemnmmo the of minimum the Th find cut. to allows infrared an functions introduce trying to simple compels with and region infrared neato srltvl ml.As trto ftepart the of iteration Also small. relatively is interaction h ubro oosaecntutd h etxcnit of consists vertex The const coupling constructed. the rapidity are in in derivative colors orders of first) number and the (zeroth leading the of ouino h atri osbeol o aiiissmall rapidities for only possible is constant motion latter coupling of the equations and of action solution the of part corresponding The h orpmrnvre ssuidi h etraieQCD. perturbative the in studied is vertex four-pomeron The ..BanDp fHg nryphysics, Energy High of Dep. Braun M.A. an-eesugSaeUiest,Russia University, State Saint-Petersburg α s orpmrnvertex Four-pomeron fteodro rae hnuiy hnteqarpepomer quadruple the when unity, than greater or order the of , ∂ y n h te ihteBK neato ewe pomerons. between interaction BFKL the with other the and etme 7 2021 27, September Abstract N 1 c nbt prahsoearvsa simple a at arrives one approaches both in with rta n ut large quite and 2 than er n n udmnn in subdominant and ant r on.Teiterative The found. are un u ob rather be to out turns w em,oewt a with one terms, two aitoa approach variational e olsosso that show collisions O ∂ t oiaigterms dominating Its rcliin ftoheavy two of collisions or y cinat action m osrce o the for constructed ams jc fitniestudy intensive of bject niin tiiiland initial at onditions susal nthe in unstable is drfrne therein) references nd etv C reggeon QCD ffective culysm h fan the sums actually epsa t solution its at tempts h rpepomeron triple the α sadi the in and us s omr an more no fthe of on Bottom-flavored inclusive emissions in the variable-flavor number scheme: a high-energy analysis

Francesco Giovanni Celiberto ,1, 2, 3, ∗ Michael Fucilla ,4, 5, † Dmitry Yu. Ivanov ,6, ‡ Mohammed M.A. Mohammed ,4, 5, § and Alessandro Papa 4, 5, ¶ 1European Centre for Theoretical Studies in Nuclear Physics and Related Areas (ECT*), I-38123 Villazzano, Trento, Italy 2Fondazione Bruno Kessler (FBK), I-38123 Povo, Trento, Italy 3INFN-TIFPA Trento Institute of Fundamental Physics and Applications, I-38123 Povo, Trento, Italy 4Dipartimento di Fisica, Universit`adella Calabria, I-87036 Arcavacata di Rende, Cosenza, Italy 5Istituto Nazionale di Fisica Nucleare, Gruppo collegato di Cosenza I-87036 Arcavacata di Rende, Cosenza, Italy 6Sobolev Institute of Mathematics, 630090 Novosibirsk, Russia

Abstract

We propose the inclusive semi-hard production, in proton-proton collisions, of two bottom- flavored hadrons, as well as of a single bottom-flavored hadron accompanied by a light jet, as novel channels for the manifestation of stabilization effects of the high-energy resummation under next-to-leading-order corrections. Our formalism relies on a hybrid high-energy and collinear factorization, where the BFKL resummation of leading and next-to-leading energy logarithms is used together with collinear factorization. We present results for cross sections and azimuthal correlations differential in rapidity, which are widely recognized as standard observables where to hunt for distinctive signals of the BFKL dynamics. We propose the study of double differential distributions in the transverse momenta of final-state particles as a common basis to investigate the interplay of different kinds of resummation mechanisms.

arXiv:2109.11875v1 [hep-ph] 24 Sep 2021 Keywords: QCD phenomenology, high-energy resummation, heavy ﬂavor, bottom production

∗Electronic address: [email protected] †Electronic address: [email protected] ‡Electronic address: [email protected] §Electronic address: [email protected] ¶Electronic address: alessandro.papa@ﬁs.unical.it Prepared for submission to JHEP

Boosted decision trees in the era of new physics: a smuon analysis case study

Alan S. Cornella, Wesley Doorsamyb, Benjamin Fuksc, Gerhard Harmsena, and Lara Masona,d aDepartment of Physics, University of Johannesburg, PO Box 524, Auckland Park 2006, South Africa bInstitute for Intelligent Systems, University of Johannesburg, PO Box 524, Auckland Park 2006, South Africa cLaboratoire de Physique Théoriqueet Hautes Energies (LPTHE), UMR 7589, Sorbonne Université et CNRS, 4 place Jussieu, 75252 Paris Cedex 05, France dUniversitéde Lyon, F-69622 Lyon, France: UniversitéLyon 1, Villeurbanne CNRS/IN2P3, UMR5822, Institut de Physique des 2 Infinis de Lyon E-mail: [email protected], [email protected], [email protected], [email protected], [email protected]

Abstract: Machine learning algorithms are growing increasingly popular in particle physics analyses, where they are used for their ability to solve difficult classification and regres- sion problems. While the tools are very powerful, they may often be under- or mis-utilised. In the following, we investigate the use of gradient boosting techniques as applicable to a generic particle physics problem. We use as an example a Beyond the Standard Model smuon collider analysis which applies to both current and future hadron colliders, and we compare our results to a traditional cut-and-count approach. In particular, we interrogate the use of metrics in imbalanced datasets which are characteristic of high energy physics problems, offering an alternative to the widely used area under the curve (auc) metric through a novel use of the F-score metric. We present an in-depth comparison of feature selection and investigation using a principal component analysis, Shapley values, and feature permutation methods in a way which we hope will be widely applicable to future

arXiv:2109.11815v1 [hep-ph] 24 Sep 2021 particle physics analyses. Moreover, we show that a machine learning model can extend the 95 % confidence level exclusions obtained in a traditional cut-and-count analysis, while potentially bypassing the need for complicated feature selections. Finally, we discuss the possibility of constructing a general machine learning model which is applicable to probe a two-dimensional mass plane. arXiv:2109.11785v1 [hep-ph] 24 Sep 2021 ln ihtecrepnigrnraie prtr ndimensional 7]. in [6, operators schemes renormalized (MS) corresponding the with along n rc at,mtvtdb h atta hs w at on do parts two these dec that energy on fact four-term Focusing the a obtain by renormalization. to motivated under step parts, origin key the trace A on and anomaly. insight trace fundamental the “some by miss decompositions energy iei n oeta nris he-emdcmoiinof decomposition three-term A se by energies. refined be qu potential a can and into contribution mass kinetic quark hadron the the of particles, of massive aspects decomposition being two natural a these fact Keeping in decompos obtains two equilibrium. above mechanical the of that concluded constraint which [5] Ref. in presented eew aeafehlo ttevra hoe ntecneto qua the of virial context the the system, 14–16]. closed in [5, a theorem Refs. for in virial that, forward the show put at we equilibrium r particular, look result In fresh a quarter, a 13]. remaining take the we provides part Here trace the and mass hsclmaigi aiu otxsi pedxA n ieabifac brief a give presen and We A, VII. Appendix Section in in contexts various summarized in then preser meaning are of physical interpre findings physical importance providing Our the when VI careful theory. tra particularly Section the by be in played should de show role one energy the We on in four-term comment energy”. operators the also We EMT anomalous in the 18]. terms 17, of mass [6–10, renormalization recently different the the in V of decomposition Section structure mass in t nucleon review of the We interpretation of physical 10]. problem the the on for III consequences Section in elaborating After otiuin ae ntecomponent the on based contributions ut ieetpcue ftencenms,tigrn eae wit debates triggering mass, nucleon deco the meaning. both physical of operators, pictures renormalized different underlying quite the of form precise sn h rc fteeeg-oetmtno ET operator (EMT) tensor energy-momentum the of trace the using h u fiscnttetmse.Udrtnighwhdo mas hadron how Understanding masses. constituen constituent elementary its these of of sum states the bound exclusively finds one hl geigwt h ahmtclapcso es 57,i sclaime is it [5–7], Refs. of aspects mathematical the with agreeing While h rsn oki raie sflos nScinI esuyi de in study we II Section In follows. as organized is work present The nodrt lrf h iuto,agnrlPoincaré-covariant a general a situation, the clarify to order In ogao h ulo ashsbe eopsdi frame-indepe a in decomposed been has mass nucleon the ago, Long gluon and quarks prevents confinement of phenomenon the QCD In 1 PT NS cl oyehiu,Isiu oyehiu d Polytechnique Institut Polytechnique, Ecole CNRS, CPHT, htteqatmaoaoseeg sntagniecontribut genuine a not th is Examining symmetry. energy translation mass. of anomalous nucleon consequence the quantum to the contribution an that “quantum new th of a as concept as well the works as minim study components, in also its We tensor Poincaré of symmetry. energy-momentum interpretation mechan physical and QCD the the decomposition the virial mass of the nucleon renormalization particular the in for discuss implications We decomposition. mass erve n xmn ndti eetdvlpet regardi developments recent detail in examine and review We 2 Cédric Lorcé, ulo asdcmoiinadmcaia equilibrium mechanical and decomposition mass nucleon eateto hsc,SR,Tml nvriy Philadelp University, Temple SERC, Physics, of Department 3 4 iatmnod iia nvri` el td iPva I Pavia, di Studi Università degli Fisica, di Dipartimento siuoNzoaed iiaNcer,Szoed ai,I- Pavia, di Sezione Nucleare, Fisica di Nazionale Istituto 1 nrymmnu esri QCD: in tensor Energy-momentum nra Metz, Andreas T 00 T ti on httetaeespr rvdstreqatr fth of quarters three provides part traceless the that found is it , 00 nters rm a enpooe nRf.[,4.Laigaieth aside Leaving 4]. [3, Refs. in proposed been has frame rest the in Dtd etme 7 2021) 27, September (Dated: .INTRODUCTION I. 2 abr Pasquini, Barbara T ai,Rued aly 12 aasa,France Palaiseau, 91128 Saclay, de Route Paris, e opsto,wihhsbe ne oediscussion some under been has which composition, 00 fre oa vra hoe”i es 3 ,8–10]. 4, [3, Refs. in theorem” “virial a as to eferred s arnmse a hrfr agl ie from differ largely therefore can masses Hadron ts. pstosaemteaial orc u provide but correct mathematically are mpositions tosatal i nomto bu aswt the with mass about information mix actually itions eETcmoet,w nlz nScinI the IV Section in analyze we components, EMT he dshm-needn nlsshsrcnl been recently has analysis scheme-independent nd aint prtr pern nalattice-regulated a in appearing operators to tation aaigters nry(..qakms)fo the from mass) quark (i.e. energy rest the parating tmxudrLrnztasomtosadhence and transformations Lorentz under mix ot e rs steeoeo totimportance. utmost of therefore is arise ses r otiuinadaguncnrbto.Quarks contribution. gluon a and contribution ark moiini osprt h M notraceless into EMT the separate to is omposition rmcicdswt h osrito mechanical of constraint the with coincides orem oepae ytetaeaoayand anomaly trace the by played role e g i h arncpyiscmuiyaottheir about community physics hadronic the hin fti omi hthsbe icse recently, discussed been has what is form this of h arncbudsaepyissprtd one separated, physics state bound hadronic the h ih fteagmnspeetdi es [8– Refs. in presented arguments the of light the euaiain(R nminimal-subtraction-type in (DR) regularization µν eaoayadterltdcneto “quantum of concept related the and anomaly ce rmi unu edter n its and theory field quantum in orem dn a noaqakadaguncontribution gluon a and quark a into way ndent igtasainsmer,adw ru that argue we and symmetry, translation ving ute eal ntevra hoe n its and theorem virial the on details further t oapa ntepyia pcrm Instead, spectrum. physical the in appear to s cleulbim ercnie the reconsider We equilibrium. ical ftencenms” aeyterl played role the namely mass”, nucleon the of ,4 3, T tmfil hoy(F) e loRf.[11– Refs. also see (QFT); theory field ntum mlu nry rpsdi some in proposed energy” omalous on fteD prahi pedxB. Appendix in approach DR the of count µν aiu ruet,w conclude we arguments, various e alteQTvrino h iiltheorem. virial the of version QFT the tail lsbrcintp cee and schemes al-subtraction-type o otems u ue sa as rule, sum mass the to ion gteqeto ftenucleon the of question the ng n ioeRodini Simone and Ri Stp cee n h operator the and schemes MS-type in DR n[]ta h w-emadthree-term and two-term the that [8] in d 1 ] ae,adcmoiinit four into decomposition a Later, 2]. [1, 210Pva Italy Pavia, -27100 i,P 92,USA 19122, PA hia, 70 ai,Italy Pavia, 27100 ,4 3, nucleon e e arXiv:2109.11754v1 [hep-ph] 24 Sep 2021 hroe eas h aetp of type same the because thermore, h M cl o ihrodrtrsdcesswt the with decreases of terms precision higher-order ex- the for are scale then, PMC And scales the PMC form. perturbative determined in the pressed orders, different at tec re r nrdcdi re oasr different absorb to order in scales of introduced PMC categories are distinct order which each in at approach, multi-scale a as terms. renormalon convergent divergent of highly elimination are the to series due perturbative and resultant independent, scheme the invariant, scale are Γ( predictions has width dealing decay for It well the works approach with ambiguity. PMC the scale that found renormalization been such the inate Γ( to correction pQCD precise achieving of pC)cretoso ea it Γ( width decay of corrections (pQCD) insae nteltrtr,tepicpeo maximum of principle the literature, [ the (PMC) renormaliza- conformality In of choice the scale. to tion independent observ- be physical should by a able required invariance, As group renormalization predictions. sys- the theoretical important an for with as error treated scale tematic usually renormalization is range guessed arbitrary the an uncertainty using the them, by Among caused possible. as much as tainty re (N next-to-next-to-next-to-next-to-leading order to up calculated losi nipratcanlfrsuyn h Higgs the studying for channel [ important two phenomenology new into an the decays searching boson is Higgs for gluons The helpful SM. also the Standard beyond is the physics it of test and precision (SM), in Model role important an plays ∗ § ‡ † h ig-oo decay Higgs-boson The [email protected] [email protected] [email protected] [email protected] ti oe httePCwsoiial introduced originally was PMC the that noted is It o h ups,w edt euetepC uncer- pQCD the reduce to need we purpose, the For nqatmcrmdnmc QD,teHgsboson Higgs the (QCD), chromodynamics quantum In 4 O [ LO) trigpitadesrstesaeivrac tec orde each at invariance scale the ensures and point starting oeacrt QDpeito ntedcywdh ..Γ( i.e. width, decay hi the be on can prediction terms pQCD higher-order accurate unknown more to due perturba dependence in scale ambiguity scale renormalization conventional re C orcin yuigtenwysgetdinfinite- suggested newly (PMC the conformality using intrinsic by corrections QCD order hr h rosaesurdaeae faltemninder mentioned the all of averages squared are errors the where hnhiKyLbrtr o atcePyisadCosmology and Physics Particle for Laboratory Key Shanghai { 1 2 β nti ae,w nlz h oa ea it fteHgsde Higgs the of width decay total the analyze we paper, this In 2 eateto hsc,CogigUiest,Cogig40 Chongqing University, Chongqing Physics, of Department hnqn e aoaoyfrSrnl ope hsc,Cho Physics, Coupled Strongly for Laboratory Key Chongqing i – } 1 11 trsit corresponding into -terms .A rsn,teprubtv QCD perturbative the present, At ]. .W r hsfcn h opportunity the facing thus are We ]. h-inGao Chu-Tian hnhiJa ogUiest,Saga 020 epesR People’s 200240, Shanghai University, Tong Jiao Shanghai 12 3 – NA,KyLbrtr o atceAtohsc n Cosmol and Astrophysics Particle for Laboratory Key INPAC, H 15 → a ensgetdt elim- to suggested been has ] gg [ ) H 16 1 , – ∗ → H { 18 β ∞ nteitiscconformality intrinsic the on igGn Wu Xing-Gang → i ,eg h PMC the e.g. ], gg ,wihue h rnil fmxmmcnomlt PC as (PMC) conformality maximum of principle the uses which ), } trsemerge -terms gg H α sn h nnt-re cl-etn prahbased approach scale-setting infinite-order the using Dtd etme 7 2021) 27, September (Dated: s aebeen have ) → n fur- And . gg ). 1 , 2 , † uDn Huang Xu-Dong em ilaettemgiueof magnitude the affect will terms etraiecoefficients perturbative etn prahcle stePMC the as called scale- approach infinite-order setting an consideration, into property (iCF) a ihtehl fRE h MMshm sguede- gauge scales is PMC mMOM-scheme the The [ RGE. fixing pendent of of help ambiguities the the with avoid momentum to physical [ that the (mMOM-scheme) under scheme ones subtraction space the into Refs.[ them form in argued been become could dependence scale scale [ residual large PMC very such the of weak, series is perturbative the of convergence fewer since orders perturbative of increment nteltrtr [ literature the in esteeetv clsa ahodrb eurn l the all requiring but by same, order each the at scale-dependent are scales effective coefficients the conformal sets the PMC, the efrhripoe.I h rsn ae,w ilre- will we paper, Γ( present width the decay indicates In the This analyze could improved. predictions PMC further eliminated. previous be the are of scales precision the PMC na- that the perturbative original of the resid- to eachture the due at and dependence form, scales scale perturbative PMC ual in fixed no-longer newly are the orders way, this Via actly. tional hr em constant Fermi where u oukonprubtv em [ dependence terms scale perturbative residual unknown has to PMC due the Thus known. are orcin yuigtePMC the using by corrections w losa the at gluons two s eety yfrhrtkn h nrni conformality intrinsic the taking further by Recently, rcial,tedcywdho h ig easinto decays Higgs the of width decay the Practically, = Γ( ieQDi lmntd n h residual the and eliminated, is QCD tive H α Sshm a era rmRefs.[ from read be can MS-scheme s .B sn h PMC the using By r. → / o sources. ror re cl-etn prahbsdon based approach scale-setting order H 29 4 colo hsc n Astronomy, and Physics of School , 20 π hysprse.W hsoti a obtain thus We suppressed. ghly gg ,adfrdfiiees eaotteLandau the adopt we definiteness, for and ], 31 epesRpbi fChina of Republic People’s 1331, → gig413,PR hn and China P.R. 401331, ngqing .Ti sbcueteukonhigher-order unknown the because is This ]. and = ) a channel cay gg 1 { , ) ‡ 21 β M µ | α 36 PMC i pbi fChina of epublic n u Zeng Jun and r s 6 .TePMC The ]. } H 3 √ odrlvlcnb xrse as expressed be can level -order trsa ahodrt aihex- vanish to order each at -terms sternraiainsae The scale. renormalization the is G 16 2 ∞ g (MOE), ogy π F G C – 336 = F 18 " H k k X ∈ H ,i sipratt trans- to important is it ], =0 ∞ 4 [0 1 = → , ∞ 4] → . C prah the approach, 42 gg ( ∞ k approach. . M +7 − 16638 3 ( gg µ § α pto up 6 H prahsat from starts approach ∞ . . 19 r 01 92 s pt N to up ) ne h conven- the under ) ) a tec order. each at a enproposed been has ,adi h pQCD the if and ], KeV, s k × +2 α the s 6 2 ( 10 - µ – 11 r 22 − { )) 4 5 β .A has As ]. OQCD LO – # i 27 GeV } , -terms such ] − (1) 2 , arXiv:2109.11722v1 [hep-ph] 24 Sep 2021 ryfrhdoi litcflw[–0 nintermediate in [8–10] [5– flow meson elliptic to hadronic prop- baryon for (NCQ) of erty scaling ratio quark phenom- number-of-constituent enhanced surprising and the 7] several as found such 2000s Collider ena the Ion Heavy in heavy-ion Relativistic (RHIC) that at un- recall the experiments We on collision hadronization. inspiration of new provide derstanding often reactions ergy pro- popular hadron reactions. of the energy describe kinds high to in three used duction often are are and which 4] models [2] [3, fragmentation recombination cluster quark present. [1], at fragmentation models non- phenomenological String by a process described is is cre- (QCD) Hadronization and gluons chromodynamics and reactions. quantum quarks energy perturbative state high final in from ated hadrons of tion atnsse ssiluulydsrbdb fragmentation by described usually mechanism. still is system parton h arnzto fpro e ihhigh hand, with other jet the parton created On of matter hadronization quark collisions. the bulk heavy-ion of relativistic quark hadronization in the of [11– mechanism for study coalescence 16] the parton prompt or observations (re-)combination These range. r sal eaddt ecoeyrltdt h for- the of to Observation in related (QGP). phenomena plasma closely these quark-gluon be phenomena of to These mation regarded usually an [21–24]. collisions heavy-ion are meson relativistic in observed to been enhance have baryon 20], [19, of long- collectivity or ratio 18], pro- ridge [17, as hadron correlation such in Collider events range phenomena multiplicity Hadron new high in Large of duction series at a available found (LHC) energies at lisions ∗ shaofl@mail.sdu.edu.cn xeietldt fhdo rdcini ihen- high in production hadron of data Experimental forma- the of process the to refers Hadronization ntels eae xeiet of experiments decade, last the In inl fqakcmiaina arnzto in hadronization at combination quark of Signals 1 colo hsclSineadItlietEgneig Ji Engineering, Intelligent and Science Physical of School cln rpryfor property scaling that ursfrom quarks oa obte nesadtepoet fsalpro syste parton small of property the understand better to as so ehns.Tetidi hteprmna aafor data experimental that is third The mechanism. olso energies. collision √ C acltos hs eut niaeasmlrt betw similarity a indicate results These calculations. QCD sn h pcrmo pdw ursfrom quarks up/down of spectrum the using pcrmrto.W ugs ytmtcmaueet in measurements systematic suggest We ratios. spectrum arn tmdrpdt in mid-rapidity at hadrons 2 colo hsc n hsclEgneig uuNra Uni Normal Qufu Engineering, Physical and Physics of School s efidsgaso ur obnto thdoiainfo th from hadronization at combination quark of signals find We 200 = I. p T INTRODUCTION pcr of spectra pp e n hta H nris epredict We energies. LHC at that and GeV φ and aaadta fu/onqak from quarks up/down of that and data p Λ olsosteeoein- therefore collisions A , p Ξ T u Song, Jun − pcr of spectra and , pp pp olsosat collisions K ∗ p n Acol- pA and 0 1 T Ω a esl-ossetydsrbduigtesetu fst of spectrum the using described self-consistently be can a-ogLi, Hai-hong n small and − and √ φ s p T 200 = d d n htfor that and p aaadta fcamqak rmperturbative from quarks charm of that and data 1 obnto f(ni)urswt h qa velocity. equal the with (anti-)quarks of combination at ac fqakcmiainin combination perfor- the in quark study of to findings tuned interesting mance is our with it of [31–37], PHTHIA energies view LHC as In such generators parameters. calcula- event with of compared tions usually were [38–44]. data experiments RHIC Experimental of years early in collaboration e ie upiigindication. surprising a gives per e oa ofidtesmlrt rdffrnei hadron in difference of or study The similarity the in find production to as so GeV imo osiun ursadatqak a systemat- can antiquarks describe and ically quarks mecha- constituent (EVC) of combination nism equal-velocity [31–37] an studies recent that Our found studies creation [25–30]. theoretical aspects possible intensive different from the attracts This i.e., mini-QGP. question, of interesting an voke h e etr fhdo rdcinin indicate production studies hadron of our feature successful, new par- usually the small and to system applied that ton often viewpoint is traditional mechanism the fragmentation with Compared hadrons. itiuino arn( hadron a of distribution pos- the in to related mini-QGP be of may formation This sible energies. LHC at lisions f(anti-)quarks of H energies. LHC hi qa eoiycmiaincncretyconstruct correctly can combination that velocity so equal (anti-)quarks their constituent are (anti-)quarks Here, n eglnShao Feng-lan and e.Tefis stecnttetqaknumber quark constituent the is first The GeV. p norECmdl[2 3,ahdo sfre ythe by formed is hadron a 33], [32, model EVC our In h rdcino dnie arn in hadrons identified of production The p √ f T T f M s B p pcr fsnl-hr arn n their and hadrons single-charm of spectra pcrmof spectrum igUiest,Sadn 715 China 273155, Shandong University, ning i i e arnpouto in production hadron een = pp II. aai h qa-eoiycombination equal-velocity the in data ( ( p p p T 0 e a ytmtclymaue ySTAR by measured systematically was GeV 200 rae in created m T T olsosat collisions pcr of spectra = ) = ) xeietldt of data experimental e UR UBRSAIGOF SCALING NUMBER QUARK pp est,Sadn 715 China 273165, Shandong versity, HADRONIC p pp p κ κ T T olsosat collisions B M pcr fietfidhdosi hspa- this in hadrons identified of spectra olsosa w olso nryscales. energy collision two at collisions i D i pcr flgtflvradsingle-charm and light-flavor of spectra f f q ∗ q 2, 1 + 1 p √ ( pp ∗ ( x and r lowl described well also are s x 1 1 olsosa different at collisions 200 = p dN/dp p T T ρ ) P ) 0 f T h eodis second The . f q pp q ¯ 2 e nfuture in GeV SPECTRA 2 T p pp ( ( √ pp olsosat collisions stepouto those of product the is ) x T x 2 2 pcr of spectra olsosat collisions s p p and T T 200 = ) ) range f . q p pp 3 pp olsosat collisions A ( x pp olsosat collisions and 3 GeV p T collisions √ ) p s , col- A = 200 (1) (2) p T Preprint no. NJU-INP 051/21 Revealing pion and kaon structure via generalised parton distributions

K. Raya,1, 2, 3 Z.-F. Cui,4, 5 L. Chang,1, ∗ J. M. Morgado,6 C. D. Roberts,4, 5, † and J. Rodr´ıguez-Quintero6, ‡ 1School of Physics, Nankai University, Tianjin 300071, China 2Departamento de F´ısica Teórica y del Cosmos, Universidad de Granada, E-18071, Granada, Spain 3Instituto de Ciencias Nucleares, Universidad Nacional Autónomade México, Apartado Postal 70-543, CDMX 04510, México 4School of Physics, Nanjing University, Nanjing, Jiangsu 210093, China 5Institute for Nonperturbative Physics, Nanjing University, Nanjing, Jiangsu 210093, China 6Department of Integrated Sciences and Center for Advanced Studies in Physics, Mathematics and Computation; University of Huelva, E-21071 Huelva; Spain. (Dated: 2021 September 23) Clear windows onto emergent hadron mass (EHM) and modulations thereof by Higgs boson interactions are provided by observable measures of pion and kaon structure, many of which are accessible via generalised parton distributions (GPDs). Beginning with algebraic GPD Ansätze, constrained entirely by hadron-scale π and K valence-parton distribution functions (DFs), in whose forms both EHM and Higgs boson influences are manifest, numerous illustrations are provided. They include the properties of electromagnetic form factors, impact parameter space GPDs, gravitational form factors and associated pressure profiles, and the character and consequences of all-orders evolution. The analyses predict that mass-squared gravitational form factors are stiffer than electromagnetic form factors; reveal that K pressure profiles are tighter than π profiles, with both mesons sus- taining near-core pressures at magnitudes similar to that expected at the core of neutron stars; deliver parameter-free predictions for π and K valence, glue, and sea GPDs at the resolving scale ζ = 2 GeV; and predict that at this scale the fraction of meson mass-squared carried by glue and sea combined matches that lodged with the valence degrees-of-freedom, with a similar statement holding for mass-squared radii. Keywords: continuum Schwinger function methods; emergence of hadron mass; Nambu-Goldstone modes – pions and kaons; nonperturbative quantum field theory; parton distributions; strong interactions

I. INTRODUCTION is ill-founded; and this realisation eliminates any model- independent link between vector-meson photoproduction Uncovering the source of the vast bulk of visible mass and the in-proton scale anomaly [13, 14]. in the Universe is one of the highest priorities in mod- In contrast, increasingly many robust connections are ern physics. It is the quest to explain the emergence being drawn between EHM and the properties of pions of hadron mass (EHM) [1–8]. Two questions define the and kaons, revealed in both phenomenological sketches heart of this problem [9]: how does a proton, built from and theoretical predictions [5–8, 15–25]. Importantly, three valence quarks whose masses are comparable to the links can be verified empirically [1–4]. Such is the that of an electron, acquire a mass mp ' 1 GeV; and context for the analysis of π- and K-meson generalised how does the pion, built from similar valence degrees-of- parton distributions (GPDs) [26–28] described herein, freedom, remain nearly massless? Within the Standard which explores their capacity to reveal novel expressions Model, the answers are expected to lie within quantum of EHM in measures of pion and kaon structure and is chromodynamics (QCD) and relate to the physical ex- also both a modernisation and extension of the studies pressions of the scale anomaly in this theory. in Refs. [29, 30]. It has been argued that low-energy elastic scattering of Our report is arranged as follows. SectionII recapitu- arXiv:2109.11686v1 [hep-ph] 23 Sep 2021 the heavy-mesons V = J/ψ, Υ from the proton can pro- lates known properties of meson light-front wave func- vide insights into connections between the proton mass tions (LFWFs) and explains features that will subse- and QCD’s scale anomaly [10]. Direct access to such re- quently be exploited. Of particular importance are the actions is difficult. So, proposals for measurements of connections drawn between LFWFs, distribution ampli- the scale anomaly have been based on the supposition tudes (DAs), and distribution functions (DFs). This is that the photoproduction of such mesons from the pro- followed in Sec.III by elucidation of our approach to the ton, e+p → e0 +V +p, provides a route to V +p → V +p construction of realistic Ansätze for such LFWFs. [11, 12]. However, it is now known that this assumption SectionIV exploits the overlap representation [31] to define GPDs, H(x, ξ, t), on the DGLAP domain (|x| ≥ ξ) [32–35]; uses it to develop algebraic forms for π and K GPDs; and explicates some of their consequences, with ∗ [email protected] results that include predictions for meson electromag- † [email protected] netic form factors. ‡ [email protected] Impact parameter space (IPS) GPDs [36] for the π Constraining Time Dependent Dark Matter Signals from the Sun

Mohammadreza Zakeria,b, Yu-Feng Zhoua,c,d [email protected], [email protected] a CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China b Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506, USA c School of Fundamental Physics and Mathematical Sciences, Hangzhou Institute for Advanced Study, UCAS, Hangzhou 310024, China d International Centre for Theoretical Physics Asia-Paciﬁc, Beijing/Hangzhou, China

Abstract Dark matter (DM) particles captured by the Sun can produce high energy electrons outside the Sun through annihilating into meta-stable mediators. The corresponding cosmic-ray electron signals observed by the space-based experiments will be time dependent due to the orbital motion of the space-based detectors. The shape of this time dependence is predictable given the orbital information of the detectors. Since the high-energy CR electron (with energy E ą 100 GeV) ﬂuxes are expected to be constant in time, non-observation of such time variation can be used to place upper limits on the DM annihilation cross section. We analyze the time dependence of dark matter cosmic-ray signals in three space-based experiments: AMS-02, DAMPE and CALET. Under the assumption that no time dependent signal is observed, we derive the 95% C.L. exclusion limits on the signal strength from the current data. We map our limits onto the parameter space of the dark photon model and ﬁnd that the constraints are comparable with that derived from the supernova SN1987A.

1 Introduction

The origin, acceleration and propagation of cosmic rays (CRs) in the Milky Way have been the subject of cosmic ray measurements by detectors such as AMS, DAMPE, CALET, and Fermi-LAT. Supernova remnants are considered [1] as the main sources for Galactic CRs. It’s also broadly accepted [2] that the CRs are accelerated by the diffusive Fermi acceleration mechanism [3] as the supernova shock wave moves through the interstellar medium1. One of the main sources of errors in analyses of CRs comes from the uncertainties in the propagation models. Therefore, constraining the propagation parameters using our observations will reduce this uncertainty. For example, it was arXiv:2109.11662v1 [hep-ph] 23 Sep 2021 shown [6] that the boron to carbon ratio (B{C) data combined with the proton flux measurements from AMS-02 can be used to determine the propagation parameters. There are also ambiguities in modeling the acceleration and propagation effects in the heliosphere especially for low energy CRs. Studying the solar modulation effects on the CRs can advance our understanding of the cosmic ray problems. As an example, the Voyager-1 spacecraft measurements of low energy CRs and their time dependence [7,8] are very significant in analyzing the local interstellar spectra and solar modulation effects on CRs. We can also investigate the propagation of different particle species in the Milky Way and the charge-sign dependency of the solar modulation effects by studying the differences of modulation effects among electrons, positrons, and nuclei [9]. On the high-energy (HE) side of the spectrum, since the Sun is not expected to be a source of HE CRs, probing

1For a critical review of this standard paradigm for the origin of the Galactic CRs check [4], and for an example of the calculation of the spectra for HE protons and nuclei see [5].

1 Multicomponent scalar dark matter at high-intensity proton beam experiments

Amalia Betancur,1, ∗ Andr´esCastillo,2, 3, 4, † Guillermo Palacio,1, ‡ and Juan Suarez1, § 1Grupo F´ısica Te´orica y Aplicada, Universidad EIA, A.A. 7516, Envigado, Colombia

2Grupo de F´ısica Experimental, Big Data y Anal´ıtica, Universidad Sergio Arboleda, Bogot´aDC, Colombia

3Instituto de Astrof´ısica de Canarias, C/ V´ıa L´actea, s/n E38205 La Laguna (Tenerife), Spain

4Universidad de La Laguna, Departamento de Astrof´ısica, La Laguna, Tenerife, Spain La Laguna (Tenerife), Spain

We study a scalar dark matter (DM) model with two DM species coupled to the Standard Model (SM) particles via a sub-GeV dark photon. In this model, we find that DM conversion occurs through the dark photon and it plays a fundamental role in setting the observed relic abundance. Furthermore, the two DM candidates can be produced at fixed-target experiments a la Beam- Dump. Detailed predictions for signal and backgrounds are obtained with the help of MadDump and NuWro Montecarlo generators. We explore the potential reach on the sensitivity of DUNE near detector and SHiP experiment, and we find that portions of the parameter space will be within reach of the two experiments.

I. INTRODUCTION

Despite the increasing evidence of the existence of dark matter (DM), its direct observation continues to be elusive. Many searches for DM have been carried out in a variety of experiments, in indirect detection looking for signals of DM annihilation in outer space, underground facilities looking for DM interaction with a detector, and even at the Large Hadron Collider, all yielding null results so far. The DM program has mostly focused on the search for a DM particle at the electroweak scale, because it has been considered that such mass scale is natural for a Weakly Interactive Massive Particle (WIMP), one of the preferred DM candidates. Nevertheless, the null results so far, have prompted the community to consider diﬀerent candidates and/or mass ranges. For instance, DM candidate could still be a WIMP with thermal decoupling but with a mass that lies well below the electroweak scale and is therefore not within reach of the usual direct detection experiments sensitive to weak scale WIMPS. It is then natural that such DM candidate has escaped observation. For this DM to communicate with the SM a new portal is needed, which could be scalar, vector or fermion [1]. One of the most widely studied models is the vector portal, where the mediator is a new gauge boson, similar to the photon, which communicates to the SM through a kinetic mixing [2–14]. Such new gauge boson, the dark photon, has preaviously been proposed to explain the observation of gamma-ray lines in the center of the Milky Way [2] and the discrepancy between the theoretical prediction and the measured value of the muon anomalous magnetic moment ∆aµ [15, 16], as is included in the historical account of [17]. In this work we will consider a dark photon that arises from arXiv:2109.11586v1 [hep-ph] 23 Sep 2021 a U(1)D symmetry, which is broken, hence yielding a massive vector boson. This dark photon will not be considered for the aforementioned reported anomalies, but as a portal to communicate the DM sector with the SM sector. In regards to the DM particles that communicate with the SM through the dark photon, the preferred candidates have been singlets, either fermion or scalar particles that are only charged under the U(1)D symmetry. Nevertheless, it is important to consider models that go beyond the simplest scenarios and that have more potential for rich phenomenology. For instance, the total relic abundance could be saturated by several DM candidates, which is natural considering that 5% of the matter-energy content of the Universe is made of a myriad of particles. This idea, dubbed multicomponent DM has gained more attention because it introduces viable candidates while also allowing for a rich phenomenology that could be within reach of actual and near-future experiments. Such models allow

∗ [email protected] † [email protected] ‡ [email protected] § [email protected] Prepared for submission to JHEP

MIT–CTP 5327

Probing a dilute short lived Quark Gluon Plasma medium with jets

Varun Vaidya Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, U.S.A. E-mail: [email protected]

Abstract: I look at propagation of jets that act as a probe of the medium in the phenomenolog- ically relevant case of a short lived dilute Quark Gluon Plasma(QGP) created during heavy ion collisions. Working in the regime where the lifetime of the medium is of the order or smaller than the formation time of the energetic jet, I derive a factorization formula with a manifest separation of scales for an illustrative jet substructure observable using the Open Quantum system Effective Field Theory(EFT) approach developed in [1], in terms of a medium structure function and a medium induced jet function . The medium structure function or the "PDF of the medium" that describes the observable independent physics of the QGP remains identical to that in a dilute long lived medium considered in [2], while the medium induced jet function is modified and in- corporates quantum interference between the hard interaction that creates the jet and subsequent medium interactions. The medium jet function continues to enjoy a BFKL rapidity evolution just as for a long lived medium, albeit with a modified renormalization scale that now depends on the

arXiv:2109.11568v1 [hep-ph] 23 Sep 2021 medium lifetime. A ﬁrst order dark SU(2)D phase transition with vector dark matter in the light of NANOGrav 12.5 yr data

Debasish Borah,1, ∗ Arnab Dasgupta,2, 3, † and Sin Kyu Kang4, ‡ 1Department of Physics, Indian Institute of Technology Guwahati, Assam 781039, India 2Pittsburgh Particle Physics, Astrophysics, and Cosmology Center, Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15206, USA 3Institute of Convergence Fundamental Studies , Seoul-Tech, Seoul 139-743, Korea 4School of Liberal Arts, Seoul-Tech, Seoul 139-743, Korea Abstract We study a dark SU(2)D gauge extension of the standard model (SM) with the possibility of a strong ﬁrst order phase transition (FOPT) taking place below the electroweak scale in the light of NANOGrav 12.5 yr data. As pointed out recently by the NANOGrav collaboration, gravitational waves (GW) from such a FOPT with appropriate strength and nucleation temperature can explain their 12.5 yr data. We impose a classical conformal invariance on the scalar potential of SU(2)D sector involving only a complex scalar doublet with negligible couplings with the SM Higgs. While a FOPT at sub-GeV temperatures can give rise to stochastic GW around nano-Hz frequencies being in agreement with NANOGrav ﬁndings, the SU(2)D vector bosons which acquire masses as a result of the FOPT in dark sector, can also serve as dark matter (DM) in the universe. The relic abundance of such vector DM can be generated in a non-thermal manner from the SM bath via scalar portal mixing. We also discuss future sensitivity of gravitational wave experiments to the model parameter space.

Introduction: The NANOGrav collaboration has re- try. More specifically, in our recent work [10], a complex cently searched for a gravitational wave (GW) signal scalar singlet charged under such Abelian gauge sym- produced from a first order phase transition (FOPT) in metry while simultaneously obeying classical conformal 45 pulsars data set spanning over 12.5 yr [1]. Accord- invariance led to a strong FOPT at bubble nucleation ing to their analysis, the 12.5 yr data can be explained temperature much below electroweak scale. For earlier in terms of a FOPT occurring at a temperature below works on FOPT within such Abelian gauge extended sce- the electroweak (EW) scale. It should however be noted narios, please refer to [13–19] and references therein. In that similar signals can also be generated by supermas- this work, we consider the possibility of realising such a sive black hole binary (SMBHB) mergers. In 2020, the low scale FOPT within a dark non-Abelian gauge model, NANOGrav collaboration also reported their evidence for particularly focusing on a SU(2)D gauge extension. Sev- an isotropic stochastic GW background with a power law eral earlier works on such dark FOPT with non-Abelian spectrum having frequencies around the nano-Hz regime gauge symmetries have already studied the consequences [2]. Although the statistical significance for a stochastic for GW at eLISA type experiments [20–22]. On the other GW brackground discovery claim was low, it still led to hand, in this work we focus on the possibility of having a several interesting new physics explanations like cosmic FOPT within such non-Abelian dark sectors at sub-EW string origins [3–5], FOPT origins [5–10] etc. While more scale such that the resulting stochastic GWs can have data are required to settle these issues, the pulsar timing frequencies typically within the NANOGRav or other arrays (PTAs) like NANOGrav sensitive to GW of ex- present as well as future PTA type experiments. tremely low frequencies do offer a complementary probe of GW background to future space-based interferometers While such dark phase transition of strongly first order like eLISA [11, 12]. Future data have the potential to and resulting GWs have been discussed earlier as well, we probe many of the above-mentioned new physics expla- study this possibility within a dark non-Abelian gauge nations for such low frequency GW background. sector for the first time after NANOGrav collaboration arXiv:2109.11558v1 [hep-ph] 23 Sep 2021 Motivated by the results from NANOGrav collabora- analysed their 12.5 year data in the context of GWs from tion explained in terms of a FOPT characterized by the the FOPT at a sub-EW scale [1]. Another motivation preferred ranges for strength (α∗) as well as temperature for such dark gauge symmetries like U(1)D, as discussed (T∗) of the phase transition as shown in [1], we study a in earlier work [10], is that it can also be motivated simple model to achieve such a strong FOPT below elec- from DM point of view where a singlet field charged troweak scale. The idea of such strong FOPT at low scale under this gauge symmetry plays the role of DM. On has been explored within the context of Abelian gauge the other hand, a dark non-Abelian gauge symmetry like extended models with or without additional dark mat- SU(2)D, as we adopt in this work, naturally contains ter (DM) fields charged under the same gauge symme- a DM candidate in the form of the massive vector 1 boson . We consider such a minimal scenario of SU(2)D

∗Electronic address: [email protected] † 1 Electronic address: [email protected] Note that U(1)D vector boson can also be made a viable DM ‡ Electronic address: [email protected] candidate by tuning the kinetic mixing U(1)D and U(1)Y of

1 Dark Conﬁnement and Chiral Phase Transitions: Gravitational Waves vs Matter Representations

1, 2, 3, 4, 3, 5, 6, Manuel Reichert, ∗ Francesco Sannino, † Zhi-Wei Wang, ‡ and Chen Zhang § 1Department of Physics and Astronomy, University of Sussex, Brighton, BN1 9QH, U.K. 2Scuola Superiore Meridionale, Largo S. Marcellino, 10, 80138 Napoli NA, Italy 3CP3-Origins, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark 4Dipartimento di Fisica “E. Pancini”, Universitàdi Napoli Federico II — INFN sezione di Napoli, Complesso Universitario di Monte S. Angelo Edificio 6, via Cintia, 80126 Napoli, Italy 5Department of Astronomy and Theoretical Physics, Lund University, 22100 Lund, Sweden 6INFN Sezione di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italy We study the gravitational-wave signal stemming from strongly coupled models featuring both, dark chiral and confinement phase transitions. We therefore identify strongly coupled theories that can feature a first-order phase transition. Employing the Polyakov-Nambu-Jona-Lasinio model, we focus our attention on SU(3) Yang-Mills theories featuring fermions in fundamental, adjoint, and two-index symmetric representations. We discover that for the gravitational-wave signals analysis, there are significant differences between the various representations. Interestingly we also observe that the two-index symmetric representation leads to the strongest first-order phase transition and therefore to a higher chance of being detected by the Big Bang Observer experiment. Our study of the confinement and chiral phase transitions is further applicable to extensions of the Standard Model featuring composite dynamics.

I. INTRODUCTION 1. We have systematically investigated the chiral and confinement phase transitions including their inter- Future gravitational-wave (GW) observatories provide play for SU(3) Yang-Mills theory with fermions new opportunities to investigate the existence of dark in fundamental, adjoint, and two-index symmet- sectors that are currently inaccessible. Because the ex- ric representations via the Polyakov-Nambu-Jona- periments will be sensitive to strong first-order phase Lasinio (PNJL) model. transitions in the early universe, it is paramount to un- 2. Representations matter: fermions in the two-index derstand the landscape of theories that can lead to GW symmetric representation increase the strength of signals. Within the landscape of theories, asymptotically the first-order confinement phase transition while free gauge-fermion systems are privileged, being already the fermions in the adjoint representation decrease chosen by nature to constitute the backbone of the Stan- it. Notably, for the adjoint and two-index sym- dard Model (SM). It is therefore reasonable to expect metric representation cases with one Dirac flavour, them also to appear in the dark sector [1–11]. the chiral phase transition is a second-order phase In [12], we embarked on a systematic investigation of transition. the composite landscape by providing the first compre- hensive study of the dark confinement phase transition 3. For all representations (chiral and confinement stemming from pure gluonic theories. There, we first ex- phase transitions), the inverse duration of the phase tended the state-of-the-art knowledge of the confinement transition is large, β˜ ∼ O(104). We discuss the phase transition to arbitrary number of colours by com- large value of β˜ in the thin-wall approximation bining effective approaches with lattice results and then where it is given by a competition of the surface determined their GW imprints. tension with the latent heat. This sheds light arXiv:2109.11552v1 [hep-ph] 23 Sep 2021 Using our work [12] on pure gauge dynamics as a step- on generic features of non-abelian gauge-fermion ping stone, we now investigate the dark dynamics stem- systems and helps finding models with stronger ming from the confinement and chiral phase transitions gravitational-wave signals. arising when adding fermions in different matter repre- 4. The confinement phase transition with fermions sentations to gauge theories. For previous works on chiral in the two-index symmetric representation has the and confinement phase transitions see e.g. [13–20]. We best chance of being detected by the Big Bang Ob- highlight below the main findings of our work: server (BBO) with a potential signal-to-noise ratio SNR ∼ O(10). 5. We further provide an outlook of which strongly ∗ Email: [email protected]; ORCID: 0000-0003-0736-5726 coupled models can potentially lead to a first-order † Email: [email protected]; ORCID: 0000-0003-2361-5326 chiral phase transition. The methods used in this ‡ Email: [email protected]; ORCID: 0000-0002-5602-6897 work can be readily applied to some of these mod- § Email: chen.zhang@fi.infn.it; ORCID: 0000-0003-2649-8508 els as well as other dark and bright extensions of Coexisting normal and intruder configurations in 32Mg

N. Kitamuraa,1,∗, K. Wimmerb,c,d,e, A. Povesf, N. Shimizua, J. A. Tostevinh, V. M. Badere,g, C. Bancroftd, D. Barofskyd, T. Baughere,g, D. Bazine, J. S. Berrymane, V. Bildsteini, A. Gadee,g, N. Imaia, T. Kr¨ollj, C. Langere, J. Lloydd, E. Lunderberge,g, F. Nowackik, G. Perdikakisd,e, F. Recchiae, T. Redpathd, S. Saenzd, D. Smalleye, S. R. Stroberge,g, Y. Utsunol,a, D. Weisshaare, A. Westerbergd

aCenter for Nuclear Study, University of Tokyo, Wako, Saitama 351-0198, Japan bInstituto de Estructura de la Materia, CSIC, 28006 Madrid, Spain cDepartment of Physics, University of Tokyo, Bunkyo, Tokyo 113-0033, Japan dDepartment of Physics, Central Michigan University, Mt. Pleasant, Michigan 48859, USA eNational Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, Michigan 48824, USA fDepartamento de F´ısica Teórica, Universidad Autónomade Madrid, 28049 Madrid, Spain gDepartment of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA hDepartment of Physics, University of Surrey, Guildford, Surrey GU2 7XH, United Kingdom iDepartment of Physics, University of Guelph, Guelph, Ontario N1G 2W1, Canada jInstitut fürKernphysik, Technische UniversitätDarmstadt, 64289 Darmstadt, Germany kInstitut Pluridisciplinaire Hubert Curien, 67037 Strasbourg, France lAdvanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan

Abstract Situated in the so-called “island of inversion,” the nucleus 32Mg is considered as an archetypal example of the disappearance of magicity at N = 20. We report on high statistics in-beam spectroscopy of 32Mg with a unique approach, in that two direct reaction probes with different sensitivities to the underlying nuclear structure are employed at the same time. More specifically, states in 32Mg were populated by knockout reactions starting from 33Mg and 34Si, lying inside and outside the island of inversion, respectively. The momentum distributions of the reaction residues and the cross sections leading to the individual final states were confronted with eikonal-based reaction calculations, yielding a significantly updated level scheme for 32Mg and spin-parity assignments. By fully exploiting observables obtained in this measurement, a variety of structures coexisting in 32Mg was unraveled. Comparisons with theoretical predictions based on shell-model overlaps allowed for clear discrimination between different structural models, revealing that the complete theoretical description of this key nucleus is yet to be achieved. Keywords: In-beam γ-ray spectroscopy, Island of inversion, Direct reactions, Radioactive beams, Shell model

While shell structure in atomic nuclei is the foundation of inversion, which has received significant attention over of nuclear physics, the breakdown of the traditional magic the years. numbers observed in unstable nuclei has revealed modifica- The nucleus 32Mg has been studied by various ap- tions in the shell structure and the important roles played proaches, e.g., β decay of 32Na, Coulomb excitation, frag- by underlying nuclear interactions [1–5]. In the “island of mentation reactions, and inclusive one-neutron knockout inversion” in the nuclear chart, ground states of neutron- from 33Mg [13–18]. Most importantly, following the ob- rich sd-shell nuclei exhibit strong admixtures of intruder servation of the low excitation energy of the first 2+ state configurations involving the fp-shell orbitals, resulting in back in 1979 [13], the large quadrupole collectivity pro- arXiv:2109.12090v1 [nucl-ex] 24 Sep 2021 the disappearance of the canonical N = 20 magic num- vided evidence for the strong deformation of the ground ber [6–8]. This concept was extended to other regions in state and the quenched N = 20 gap [16]. In 2010, the the nuclear chart. To date, islands of inversion have been first-excited 0+ state was identified in the t(30Mg,p) reac- proposed at N = 8, 20, 28, 40, and 50 [1, 9–12], and tion [19], and later the existence of this state was confirmed universal understanding of these islands and their driving in an in-beam γ-ray spectroscopy experiment [20]. From a + mechanism has been one of the main topics in modern nu- simplified point of view, the 02 was interpreted as a coun- clear physics. This Letter focuses on the nucleus 32Mg terpart of the near-spherical ground states of 30Mg [21] (N = 20 and Z = 12), at the heart of the original island and 34Si [22], and it was taken as a manifestation of shape coexistence, which is proposed to be the universal phe- ∗Corresponding author nomenon in the islands of inversion [23]. However, it was Email address: [email protected] (N. Kitamura) later shown that its structure is much more complex. Key 1Present address: Department of Physics and Astronomy, Uni- properties of the structural change between the 0+ states versity of Tennessee, Knoxville, Tennessee 37996, USA

Preprint submitted to Physics Letters B September 27, 2021