Abstract: We present a search for magnetic monopoles and high-electric-charge objects using LHC Run 2 root s = 13TeV proton-proton collisions recorded by the ATLAS detector. A total integrated luminosity of 138 fb(-1) was collected by a specialized trigger. No highly ionizing particle candidate was observed. Considering the Drell-Yan and photon-fusion pair production mechanisms as benchmark models, cross-section upper limits are presented for spin-0 and spin-1/2 magnetic monopoles of magnetic charge 1g(D) and 2g(D) and for high-electric-charge objects of electric charge 20 <= vertical bar z vertical bar <= 100, for masses between 200 GeV and 4000 GeV. The search improves by approximately a factor of three the previous cross-section limits on the Drell-Yan production of magnetic monopoles and high-electric charge objects. Also, the first ATLAS limits on the photon-fusion pair production mechanism of magnetic monopoles and high-electric-charge objects are obtained.

Abstract: We present a novel nonperturbative approach for calculating the form factors of the quark-gluon vertex in terms of an unknown three-point function, in the Landau gauge. The key ingredient of this method is the exact all-order relation connecting the conventional quark-gluon vertex with the corresponding vertex of the background field method, which is Abelian-like. When this latter relation is combined with the standard gauge technique, supplemented by a crucial set of transverse Ward identities, it allows the approximate determination of the nonperturbative behavior of all 12 form factors comprising the quark-gluon vertex, for arbitrary values of the momenta. The actual implementation of this procedure is carried out in the Landau gauge, in order to make contact with the results of lattice simulations performed in this particular gauge. The most demanding technical aspect involves the approximate calculation of the components of the aforementioned (fully dressed) three-point function, using lattice data as input for the gluon propagators appearing in its diagrammatic expansion. The numerical evaluation of the relevant form factors in three special kinematical configurations (soft-gluon and quark symmetric limit, zero quark momentum) is carried out in detail, finding qualitative agreement with the available lattice data. Most notably, a concrete mechanism is proposed for explaining the puzzling divergence of one of these form factors observed in lattice simulations.

Abstract: We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that forces the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting function. The method relies on a set of hyperparameters that control how gradients are treated in the physical loss. The performance of our methodology is demonstrated by solving Riemann problems in special relativistic hydrodynamics, extending earlier studies with PINNs in the context of the classical Euler equations. The solutions obtained with the GA-PINN model correctly describe the propagation speeds of discontinuities and sharply capture the associated jumps. We use the relative l(2) error to compare our results with the exact solution of special relativistic Riemann problems, used as the reference ''ground truth'', and with the corresponding error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by the GA-PINN model is comparable to that obtained with a shock-capturing scheme, achieving a performance superior to that of the baseline PINN algorithm in general. An additional benefit worth stressing is that our PINN-based approach sidesteps the costly recovery of the primitive variables from the state vector of conserved variables, a well-known drawback of grid-based solutions of the relativistic hydrodynamics equations. Due to its inherent generality and its ability to handle steep gradients, the GA-PINN methodology discussed in this paper could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions.

Abstract: We present a novel method for computing the nonperturbative kinetic term of the gluon propagator from an ordinary differential equation, whose derivation hinges on the central hypothesis that the regular part of the three-gluon vertex and the aforementioned kinetic term are related by a partial Slavnov-Taylor identity. The main ingredients entering in the solution are projection of the three-gluon vertex and a particular derivative of the ghost-gluon kernel, whose approximate form is derived from a Schwinger-Dyson equation. Crucially, the requirement of a pole-free answer determines the initial condition, whose value is calculated from an integral containing the same ingredients as the solution itself. This feature fixes uniquely, at least in principle, the form of the kinetic term, once the ingredients have been accurately evaluated. In practice, however, due to substantial uncertainties in the computation of the necessary inputs, certain crucial components need be adjusted by hand, in order to obtain self-consistent results. Furthermore, if the gluon propagator has been independently accessed from the lattice, the solution for the kinetic term facilitates the extraction of the momentum-dependent effective gluon mass. The practical implementation of this method is carried out in detail, and the required approximations and theoretical assumptions are duly highlighted.

Abstract: We present a novel analysis of the gluon gap equation, where its full nonlinear structure is duly taken into account. In particular, while in previous treatments the linearization of this homogeneous integral equation introduced an indeterminacy in the scale of the corresponding mass, the current approach determines it uniquely, once the value of the gauge coupling at a given renormalization point is used as input. A crucial ingredient for this construction is the “kinetic term” of the gluon propagator, whose form is not obtained from the complicated equation governing its evolution, but is rather approximated by suitable initial Ansatze, which are subsequently improved by means of a systematic iterative procedure. The multiplicative renormalization of the central equation is carried out following an approximate method, which is extensively employed in the studies of the standard quark gap equation. This approach amounts to the effective substitution of the vertex renormalization constants by kinematically simplified form factors of the three- and four-gluon vertices. The resulting numerical interplay, exemplified by the infrared suppression of the three-gluon vertex and the mild enhancement of the four-gluon vertex, is instrumental for obtaining positive-definite and monotonically decreasing running gluon masses. The resulting gluon propagators, put together from the gluon masses and kinetic terms obtained with this method, match rather accurately the data obtained from large-volume lattice simulations.