Phase-space Integrals Research Articles

Parallel superposition of small-amplitude oscillatory shear flow upon steady shear flow from rotarance theory

The power of a macromolecular theory for the transport properties of a polymeric liquid increases with the number of analytical expressions for its most important material functions. In this work, we add another of these canonical function to our recent series of material function derivations for rotarance theory. By rotarance theory, we mean the explanation of the elasticity of polymeric liquids by use of (i) the diffusion equation to get the orientation distribution in Euler coordinates, and (ii) the integration in phase space using this distribution to get the target material function. In this paper, we target parallel superposition of oscillatory shear flow upon steady shear flow. We arrive at analytical expressions for both parts of the complex viscosity in parallel superposition. We find that these explain the classic experimental observations in parallel superposition: (a) the maximum in the real part of the complex viscosity, and (b) the negative values of minus its imaginary part, and (c) the independence of the steady mean shear stress from the superposed oscillation.

Two-loop master integrals for leading-color pp→tt¯H\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ pp\ o t\\overline{t}H $$\\end{document} amplitudes with a light-quark loop

We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in tt¯H\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ t\\overline{t}H $$\\end{document} production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations in ϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Three-point energy correlators in hadronic Higgs boson decays

We present the analytic calculation of the leading-order three-point energy correlator (EEEC) in hadronic Higgs decays, including both the gluon-initiated channel H→gg+X and quark-initiated channel H→qq¯+X. The phase-space integration is evaluated directly using Mandelstam variables sij=(pi+pj)2, and the appearing square roots can be rationalized by either conformal ratios or celestial coordinate variables. Throughout the calculation, we observe the same transcendental function space as in N=4 super–Yang-Mills theory and e+e−→hadrons. Different infrared limits are also explored using the full analytic result, offering the fixed-order data for EEEC factorization and resummation. Given its nontrivial shape dependence, the EEEC presents an excellent opportunity to explore the dynamics of gluon jets originating from the H→gg decay channel at future lepton colliders. Published by the American Physical Society 2024

Next-to-Next-to-Leading Order QCD Corrections to Semi-Inclusive Deep-Inelastic Scattering.

We present the first results for the next-to-next-to leading order (NNLO) corrections to the semi-inclusive deep-inelastic scattering process in perturbative quantum chromodynamics. We consider the quark initiated flavor nonsinglet process and obtain the complete contributions analytically at leading color. All relevant virtual and real emission Feynman diagrams have been computed using integration-by-parts reduction to master integrals and two approaches for their subsequent evaluation (parametric phase-space integration and method of differential equations). The numerical analysis demonstrates the significance of the NNLO corrections and their great impact on the reduction of the residual scale dependence.

One-loop corrections to the double-real emission contribution to the zero-jettiness soft function at N3LO in QCD

We present an analytic calculation of the one-loop correction to the double-real emission contribution to the zero-jettiness soft function at N3LO in QCD, accounting for both gluon-gluon and quark-antiquark soft final-state partons. We explain all the relevant steps of the computation including the reduction of phase-space integrals to master integrals in the presence of Heaviside functions, and the methods we employed to compute them.

Noether invariance theory for the equilibrium force structure of soft matter

We give details and derivations for the Noether invariance theory that characterizes the spatial equilibrium structure of inhomogeneous classical many-body systems, as recently proposed and investigated for bulk systems (Sammüller et al 2023 Phys. Rev. Lett. 130 268203). Thereby an intrinsic thermal symmetry against a local shifting transformation on phase space is exploited on the basis of the Noether theorem for invariant variations. We consider the consequences of the shifting that emerge at second order in the displacement field that parameterizes the transformation. In a natural way the standard two-body density distribution is generated. Its second spatial derivative is thereby balanced by two further and different two-body correlation functions, which respectively introduce thermally averaged force correlations and force gradients in a systematic and microscopically sharp way into the framework. Separate exact self and distinct sum rules express this balance. We exemplify the validity of the theory on the basis of computer simulations for the Lennard–Jones gas, liquid, and crystal, the Weeks–Chandler–Andersen fluid, monatomic Molinero–Moore water at ambient conditions, a three-body-interacting colloidal gel former, the Yukawa and soft-sphere dipolar fluids, and for isotropic and nematic phases of Gay–Berne particles. We describe explicitly the derivation of the sum rules based on Noether’s theorem and also give more elementary proofs based on partial phase space integration following Yvon’s theorem.

Analytic auxiliary mass flow to compute master integrals in singular kinematics

The computation of master integrals from their differential equations requires boundary values to be supplied by an independent method. These boundary values are often desired at singular kinematical points. We demonstrate how the auxiliary mass flow technique can be extended to compute the expansion coefficients of master integrals in a singular limit in an analytical manner, thereby providing these boundary conditions. To illustrate the application of the method, we re-compute the phase space integrals relevant to initial-final antenna functions at NNLO, now including higher-order terms in their ϵ-expansion in view of their application in third-order QCD corrections.

On the Partition Temperature of Massless Particles in High-Energy Collisions

Although partition temperature derived using the Darwin–Fowler method is exact for simple scenarios, the derivation for complex systems might reside in specific approximations whose viability is not ensured if the thermodynamic limit is not attained. This work elaborates on a related problem relevant to relativistic high-energy collisions. On the one hand, it is simple enough that closed-form expressions can be obtained precisely for the one-particle distribution function. On the other hand, the resulting expression is not an exponential form, and therefore, it is not straightforward that the notion of partition function could be implied. Specifically, we derive the one-particle distribution function for massless particles where the phase-space integration is performed exactly for the underlying canonical ensemble consisting of a given number of particles. We discuss the viability of the partition temperature in this case. Possible implications of the obtained results regarding the observed Tsallis distribution in transverse momentum spectra in high-energy collisions are also addressed.

Efficient phase-space generation for hadron collider event simulation

We present a simple yet efficient algorithm for phase-space integration at hadron colliders. Individual mappings consist of a single t-channel combined with any number of s-channel decays, and are constructed using diagrammatic information. The factorial growth in the number of channels is tamed by providing an option to limit the number of s-channel topologies. We provide a publicly available, parallelized code in C++ and test its performance in typical LHC scenarios.

MadNIS - Neural multi-channel importance sampling

Theory predictions for the LHC require precise numerical phase-space integration and generation of unweighted events. We combine machine-learned multi-channel weights with a normalizing flow for importance sampling, to improve classical methods for numerical integration. We develop an efficient bi-directional setup based on an invertible network, combining online and buffered training for potentially expensive integrands. We illustrate our method for the Drell-Yan process with an additional narrow resonance.

Steady elongational flow from rotarance theory

The physics of macromolecular orientation has been used to explain the elasticity of polymeric liquids. Specifically, by first sculpting a rigid bead–rod likeness of the macromolecule, we can then derive its hydrodynamic resistance to orientation. The solution for the orientation distribution function has then been used, by integration in phase space, to get rheological material functions in both (i) small- and (ii) large-amplitude oscillatory shear flow, including its limiting case, and (iii) steady shear flow. However, rheological material functions in steady homogeneous extension from rigid bead–rod theory remain elusive. In this paper, we derive the orientation distribution function, and the rheological material functions, for suspensions of general rigid bead–rod structures. We focus on the time-steady viscosities in extension, and we first do so for general extensional kinematics. We then obtain the viscosities in steady extension for (iv) uniaxial extension, (v) planar extension, and (vi) biaxial extension. We close with a worked example, in which we use our new result for the steady uniaxial extensional viscosity to build a bridge between the macromolecular theory and the Oldroyd framework for rheological constitutive models. We, thus, arrive at a constitutive equation whose parameters are deducible from the moments of inertia of the macromolecule, and thus, deducible from macromolecular architecture alone. Our model is accurate up to third order for time-independent flows and is accurate to second order for time-dependent ones.

Two-loop infrared renormalization with on-shell methods

Within the framework proposed by Caron-Huot and Wilhelm, we give a recipe for computing infrared anomalous dimensions purely on-shell, efficiently up to two loops in any massless theory. After introducing the general formalism and reviewing the one-loop recipe, we extract a practical formula that relates two-loop infrared anomalous dimensions to certain two- and three-particle phase space integrals with tree-level form factors of conserved operators. We finally provide several examples of the use of the two-loop formula and comment on some of its formal aspects, especially the cancellation of ‘one-loop squared’ spurious terms. The present version of the paper is augmented with a detailed treatment of the structure of infrared divergences in massless theories of scalars and fermions up to two loops. In the calculation we encounter divergent phase space integrals and show in detail how these cancel among each other as required by the finiteness of the anomalous dimension. As a non-trivial check of the method, we also perform the computation with a standard diagrammatic approach, finding perfect agreement.

PSGen, a generator of phase space parameterizations for the multichannel Monte Carlo integration

PSGen is a new general purpose Fortran program which has been written to facilitate the Monte Carlo phase space integration of the S matrix element of any 2→n scattering process, with n=2,...,9, provided by the user. The program is written in Fortran 90/95. It uses a new very fast algorithm that automatically generates calls to Fortran subroutines containing different phase space parameterizations of the considered class of processes. The parameterizations take into account mappings of poles due to the Feynman propagators of unstable heavy particles decaying into 2 or 3 on shell final state particles according to predefined patterns, possible single or double t-channel poles and peaks due to one on shell photon or gluon radiation. The individual subroutines are organized in a single multichannel kinematics subroutine which can be easily called while computing the phase space integral of the S matrix element as a function of generated particle four momenta, in either the leading or higher orders of the perturbation series. The particle four momenta can be used in a quadruple precision version, if necessary. Program summaryProgram Title:PSGenCPC Library link to program files:https://doi.org/10.17632/y7nprysy8g.1Developer's repository link:https://www.kk.us.edu.pl/psgen.htmlLicensing provisions: GPLv3Programming language: Fortran 90/95Nature of problem: Automatic generation of Fortran routines which allow for the multichannel Monte Carlo integration of user's provided S matrix element as a function of generated particle four momenta of any 2→n, with n=2,...,9, scattering process. The S matrix can be defined in either the leading or higher orders.Solution method: A new very fast algorithm written in Fortran 90/95 generates calls to different handwritten subroutines containing phase space parameterizations which, among others, take into account poles due to the Feynman propagators of unstable heavy particles decaying into 2 or 3 on shell final state particles according to predefined patterns, possible single or double t-channel poles and peaks due to one on shell photon or gluon radiation. The generated routines are moved to a directory, where they are combined in a single properly normalized multichannel phase space parameterization. The corresponding kinematics subroutine is self-consistent and can be readily combined with any program which calculates the S matrix element in either the leading or higher orders of the perturbation series. The particle four momenta can be generated in a quadruple precision version, if necessary. The automatically generated subroutine containing calls to different kinematics subroutines can be easily supplemented with calls to user's own made kinematics subroutines.

Complete two-loop QCD amplitudes for tW production at hadron colliders

We have calculated the complete two-loop QCD amplitudes for hadronic tW production by combining analytical and numerical techniques. The amplitudes have been first reduced to master integrals of eight planar and seven non-planar families, which can contain at most four massive propagators. Then a rational transformation of the master integrals is found to obtain a good basis so that the dimensional parameter decouples from the kinematic variables in the denominators of reduction coefficients. The master integrals are computed by solving their differential equations numerically. We find that the finite part of the two-loop squared amplitude is stable in the bulk of the phase space. After phase space integration and convolution with the parton distributions, it increases the LO cross section at the 13 TeV LHC by about 3%.

Decoherence effects in reactor and Gallium neutrino oscillation experiments: a QFT approach

We adopt the quantum field theoretical method to calculate the amplitude and event rate for a neutrino oscillation experiment, considering neutrino production, propagation and detection as a single process. This method allows to take into account decoherence effects in the transition amplitude induced by the quantum mechanical uncertainties of all particles involved in the process. We extend the method to include coherence loss due to interactions with the environment, similar to collisional line broadening. In addition to generic decoherence induced at the amplitude level, the formalism allows to include, in a straightforward way, additional damping effects related to phase-space integrals over momenta of unobserved particles as well as other classical averaging effects. We apply this method to neutrino oscillation searches at reactor and Gallium experiments and confirm that quantum decoherence is many orders of magnitudes smaller than classical averaging effects and therefore unobservable. The method used here can be applied with minimal modifications also to other types of oscillation experiments, e.g., accelerator based beam experiments.

Constructing the general partial wave and renormalization in effective field theory

We construct the general partial wave amplitude basis for the $N\ensuremath{\rightarrow}M$ scattering, which consists of Poincar\'e Clebsch-Gordan coefficients, in a Lorentz invariant form given in terms of the spinor-helicity variables. The inner product of the Clebsch-Gordan coefficients is defined, which converts on-shell phase space integration into an algebraic problem. We also develop the technique of partial wave expansions of arbitrary amplitudes, including those with infrared divergence. These are applied to the computation of anomalous dimension matrix for general effective operators, where unitarity cuts for the loop amplitudes, with an arbitrary number of external particles, are obtained via partial wave expansion.

Matsubara dynamics approximation for generalized multi-time correlation functions.

We introduce a semi-classical approximation for calculating generalized multi-time correlation functions based on Matsubara dynamics, a classical dynamics approach that conserves the quantum Boltzmann distribution. This method is exact for the zero time and harmonic limits and reduces to classical dynamics when only one Matsubara mode is considered (i.e., the centroid). Generalized multi-time correlation functions can be expressed as canonical phase-space integrals, involving classically evolved observables coupled through Poisson brackets in a smooth Matsubara space. Numerical tests on a simple potential show that the Matsubara approximation exhibits better agreement with exact results than classical dynamics, providing a bridge between the purely quantum and classical descriptions of multi-time correlation functions. Despite the phase problem that prevents practical applications of Matsubara dynamics, the reported work provides a benchmark theory for the future development of quantum-Boltzmann-preserving semi-classical approximations for studies of chemical dynamics in condensed phase systems.

GNN-assisted phase space integration with application to atomistics

Overcoming the time scale limitations of atomistics can be achieved by switching from the state-space representation of Molecular Dynamics (MD) to a statistical-mechanics-based representation in phase space, where approximations such as maximum-entropy or Gaussian phase packets (GPP) evolve the atomistic ensemble in a time-coarsened fashion. In practice, this requires the computation of expensive high-dimensional integrals over all of phase space of an atomistic ensemble. This, in turn, is commonly accomplished efficiently by low-order numerical quadrature. We show that numerical quadrature in this context, unfortunately, comes with a set of inherent problems, which corrupt the accuracy of simulations—especially when dealing with crystal lattices with imperfections. As a remedy, we demonstrate that Graph Neural Networks, trained on Monte-Carlo data, can serve as a replacement for commonly used numerical quadrature rules, overcoming their deficiencies and significantly improving the accuracy. This is showcased by three benchmarks: the thermal expansion of copper, the martensitic phase transition of iron, and the energy of grain boundaries. We illustrate the benefits of the proposed technique over classically used third- and fifth-order Gaussian quadrature, highlight the impact on time-coarsened atomistic predictions, and discuss the computational efficiency. The latter is of general importance when performing frequent evaluation of phase space or other high-dimensional integrals, which is why the proposed framework promises applications beyond the scope of atomistics.

Compression and information entropy of binary strings from the collision history of three hard balls

We investigate how to measure and define the entropy of a simple chaotic system, three hard spheres on a ring. A novel approach is presented, which does not assume the ergodic hypothesis. It consists of transforming the particles’ collision history into a sequence of binary digits. We then investigate three approaches which should demonstrate the non-randomness of these collision-generated strings compared with random number generator created strings: Shannon entropy, diehard randomness tests and compression percentage. We show that the Shannon information entropy is unable to distinguish random from deterministic strings. The Diehard test performs better, but for certain mass-ratios the collision-generated strings are misidentified as random with high confidence. The zlib and bz2 compression algorithms are efficient at detecting non-randomness and low information content, with compression efficiencies that tend to 100% in the limit of infinite strings. Thus ‘compression algorithm entropy’ is non-extensive for this chaotic system, in marked contrast to the extensive entropy determined from phase-space integrals by assuming ergodicity.

Quark and gluon two-loop beam functions for leading-jet pT and slicing at NNLO

We compute the complete set of two-loop beam functions for the transverse momentum distribution of the leading jet produced in association with an arbitrary colour-singlet system. Our results constitute the last missing ingredient for the calculation of the jet-vetoed cross section at small veto scales at the next-to-next-to-leading order, as well as an important ingredient for its resummation to next-to-next-to-next-to-leading logarithmic order. Our calculation is performed in the soft-collinear effective theory framework with a suitable regularisation of the rapidity divergences occurring in the phase-space integrals. We discuss the occurrence of soft-collinear mixing terms that might violate the factorisation theorem, and demonstrate that they are naturally absorbed into the beam functions at two loops in the exponential rapidity regularisation scheme when performing a multipole expansion of the measurement function. As in our recent computation of the two-loop soft function, we present the results as a Laurent expansion in the jet radius R. We provide analytic expressions for all flavour channels in x space with the exception of a set of R-independent non-logarithmic terms that are given as numerical grids. We also perform a fully numerical calculation with exact R dependence, and find that it agrees with our analytic expansion at the permyriad level or better. Our calculation allows us to define a next-to-next-to-leading order slicing method using the leading-jet pT as a slicing variable. As a check of our results, we carry out a calculation of the Higgs and Z boson total production cross sections at the next-to-next-to-leading order in QCD.