Abstract: Optimisation problems are ubiquitous in particle and astrophysics, and involve locating the optimum of a complicated function of many parameters that may be computationally expensive to evaluate. We describe a number of global optimisation algorithms that are not yet widely used in particle astrophysics, benchmark them against random sampling and existing techniques, and perform a detailed comparison of their performance on a range of test functions. These include four analytic test functions of varying dimensionality, and a realistic example derived from a recent global fit of weak-scale supersymmetry. Although the best algorithm to use depends on the function being investigated, we are able to present general conclusions about the relative merits of random sampling, Differential Evolution, Particle Swarm Optimisation, the Covariance Matrix Adaptation Evolution Strategy, Bayesian Optimisation, Grey Wolf Optimisation, and the PyGMO Artificial Bee Colony, Gaussian Particle Filter and Adaptive Memory Programming for Global Optimisation algorithms.

Abstract: The Xi(++)(cc) -> Xi('+)(c)pi(+) decay is observed using proton-proton collisions collected by the LHCb experiment at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 5.4 fb(-1). The Xi(++)(cc) -> Xi('+)(c)pi(+) decay is reconstructed partially, where the photon from the Xi('+)(c) -> Xi(+)(c)gamma decay is not reconstructed and the pK(-)pi(+) final state of the Sc+ baryon is employed. The Xi(++)(cc) -> Xi('+)(c)pi(+) branching fraction relative to that of the Xi(++)(cc) -> Xi('+)(c)pi(+) decay is measured to be 1.41 +/- 0.17 +/- 0.10, where the first uncertainty is statistical and the second systematic.

Abstract: We present a novel benchmark application of a quantum algorithm to Feynman loop integrals. The two on-shell states of a Feynman propagator are identified with the two states of a qubit and a quantum algorithm is used to unfold the causal singular configurations of multiloop Feynman diagrams. To identify such configurations, we exploit Grover's algorithm for querying multiple solutions over unstructured datasets, which presents a quadratic speed-up over classical algorithms when the number of solutions is much smaller than the number of possible configurations. A suitable modification is introduced to deal with topologies in which the number of causal states to be identified is nearly half of the total number of states. The output of the quantum algorithm in IBM Quantum and QUTE Testbed simulators is used to bootstrap the causal representation in the loop-tree duality of representative multiloop topologies. The algorithm may also find application and interest in graph theory to solve problems involving directed acyclic graphs.

Abstract: We generalize the Chern-Simons modified gravity to the metric-affine case and impose projective invariance by supplementing the Pontryagin density with homothetic curvature terms which do not spoil topologicity. The latter is then broken by promoting the coupling of the Chern-Simons term to a (pseudo)-scalar field. The solutions for torsion and nonmetricity are derived perturbatively, showing that they can be iteratively obtained from the background fields. This allows us to describe the dynamics for the metric and the scalar field perturbations in a self-consistent way, and we apply the formalism to the study of quasi normal modes in a Schwarzschild black hole background. Unlike in the metric formulation of this theory, we show that the scalar field is endowed with dynamics even in the absence of its kinetic term in the action. Finally, using numerical methods we compute the quasinormal frequencies and characterize the late-time power law tails for scalar and metric perturbations, comparing the results with the outcomes of the purely metric approach.

Abstract: Scatterings of galactic dark matter (DM) particles with the constituents of celestial bodies could result in their accumulation within these objects. Nevertheless, the finite temperature of the medium sets a minimum mass, the evaporation mass, that DM particles must have in order to remain trapped. DM particles below this mass are very likely to scatter to speeds higher than the escape velocity, so they would be kicked out of the capturing object and escape. Here, we compute the DM evaporation mass for all spherical celestial bodies in hydrostatic equilibrium, spanning the mass range [10(-)(10) – 10(2)] M-circle dot, for constant scattering cross sections and s-wave annihilations. We illustrate the critical importance of the exponential tail of the evaporation rate, which has not always been appreciated in recent literature, and obtain a robust result: for the geometric value of the scattering cross section and for interactions with nucleons, at the local galactic position, the DM evaporation mass for all spherical celestial bodies in hydrostatic equilibrium is approximately given by E-c/T-chi similar to 30, where E-c is the escape energy of DM particles at the core of the object and T-chi is their temperature. In that case, the minimum value of the DM evaporation mass is obtained for super-Jupiters and brown dwarfs, m(ev)(ap) similar or equal to 0.7 GeV. For other values of the scattering cross section, the DM evaporation mass only varies by a factor smaller than three within the range 10(-41) cm(2) <= sigma(p) <= 10(-31) cm(2), where sigma(p) is the spin-independent DM-nucleon scattering cross section. Its dependence on parameters such as the galactic DM density and velocity, or the scattering and annihilation cross sections is only logarithmic, and details on the density and temperature profiles of celestial bodies have also a small impact.