Abstract: We present a search for nine lepton-number-violating and three lepton-flavor-violating neutral charm decays of the type D-0 -> h'(-) h(-) l'(+) l(+) and D-0 -> h'(-) h(+) l'(+/-) l(-/+), where h and h' represent a K or pi meson and l and l' an electron or muon. The analysis is based on 468 fb(-1) of e(+) e(-) annihilation data collected at or close to the Upsilon(4S) resonance with the BABAR detector at the SLAC National Accelerator Laboratory. No significant signal is observed for any of the twelve modes, and we establish 90% confidence level upper limits on the branching fractions in the range (1.0-30.6) x 10(-7). The limits are between 1 and 3 orders of magnitude more stringent than previous measurements.

Abstract: We present a precise lattice computation of pseudoscalar and vector heavy-light meson masses for heavy-quark masses ranging from the physical charm mass up to similar or equal to 4 times the physical b-quark mass. We employ the gauge configurations generated by the European Twisted Mass Collaboration (ETMC) with N-f = 2 + 1 + 1 dynamical quarks at three values of the lattice spacing (a similar or equal to 0.062; 0.082; 0.089 fm) with pion masses in the range M-pi similar or equal to 210-450 MeV. The heavy-quark mass is simulated directly on the lattice up to similar or equal to 3 times the physical charm mass. The interpolation to the physical b-quark mass is performed using the ETMC ratio method, based on ratios of the meson masses computed at nearby heavy-quark masses, and adopting the kinetic mass scheme. The extrapolation to the physical pion mass and to the continuum limit yields m(b)(kin) (1 GeV) = 4.61(20) GeV, which corresponds to (m) over bar (b) ((m) over bar (b)) 4.26(18) GeV in the (MS) over bar scheme. The lattice data are analyzed in terms of the heavy-quark expansion (HQE) and the matrix elements of dimension-four and dimension-five operators are extracted with a good precision, namely,(Lambda) over bar = 0.552(26) GeV, mu(2)(pi) = 0.321(32) GeV2, and mu(2)(G)(m(b)) = 0.253(25) GeV2. The data also allow for a rough estimate of the dimension-six operator matrix elements. As the HQE parameters play a crucial role in the inclusive determination of the Cabibbo-Kobayashi-Maskawa matrix elements V-ub and V-cb, their precise determination on the lattice may eventually validate and improve the analyses based on fits to the semileptonic moments.

Abstract: We present a novel theory of gravity by considering an extension of symmetric teleparallel gravity. This is done by introducing, in the framework of the metric-affine formalism, a new class of theories where the nonmetricity Q is nonminimally coupled to the matter Lagrangian. More specifically, we consider a Lagrangian of the form L similar to f(1)(Q) + f(2)(Q)L-M, where f(1) and f(2) are generic functions of Q, and L-M is the matter Lagrangian. This nonminimal coupling entails the nonconservation of the energy-momentum tensor, and consequently the appearance of an extra force. The formulation of the gravity sector in terms of the Q instead of the curvature may result in subtle improvements of the theory. In the context of nonminimal matter couplings, we are therefore motivated to explore whether the new geometrical formulation in terms of the Q, when implemented also in the matter sector, would allow more universally consistent and viable realizations of the nonminimal coupling. Furthermore, we consider several cosmological applications by presenting the evolution equations and imposing specific functional forms of the functions f(1)(Q) and f(2)(Q), such as power-law and exponential dependencies of the nonminimal couplings. Cosmological solutions are considered in two general classes of models, and found to feature accelerating expansion at late times.

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.