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 detailed study of the lowest-lying 1/2(-) and 3/2(-) Lambda Q resonances both in the heavy 2 2 quark (bottom and charm) and the strange sectors. We have paid special attention to the interplay between the constituent quark-model and chiral baryon-meson degrees of freedom, which are coupled using a unitarized scheme consistent with leading-order heavy quark symmetries. We show that the Lambda(b)(5912) [J(P) = 1/2(-)], Lambda(b)(5920) [J(P) = 3/2(-)] and the Lambda(c)(2625) [J(P) = 3/2-], and the Lambda(1520) [J(P) = 3/2(-)] admitting larger breaking corrections, are heavyquark spin-flavor siblings. They can be seen as dressed quark-model states with Sigma Q(()*()) pi molecular components of the order of 30%. The J(P)=1(-) Lambda(2595) has, however, a higher molecular 2 probability of at least 50%, and even values greater than 70% can be easily accommodated. This is because it is located almost on top of the threshold of the Sigma(c)pi pair, which largely influences its properties. Although the light degrees of freedom in this resonance would be coupled to spin-parity 1(-) as in the Lambda(b)(5912), Lambda(b)(5920) and Lambda(c)(2625), the Lambda(c)(2595) should not be considered as a heavy-quark spin-flavor partner of the former ones. We also show that the Lambda(1405) chiral two-pole pattern does not have analogs in the 1 – charmed and bottomed sectors, because the 2 N D-(*()) and N (B) over bar (()*()) channels do not play for heavy quarks the decisive role that the N (K) over bar does in the strange sector, and the notable influence of the bare quark-model states for the charm and bottom resonances. Finally, we predict the existence of two Lambda(b)(6070) and two Lambda(c)(2765) heavy-quark spin and flavor sibling odd parity states.

Abstract: We perform a study of the (B*+B0), (BB+)-B-*0 correlation functions using an extension of the local hidden gauge approach which provides the interaction from the exchange of light vector mesons and gives rise to a bound state of these components in I = 0 with a binding energy of about 21 MeV. After that, we face the inverse problem of determining the low energy observables, scattering length and effective range for each channel, the possible existence of a bound state, and, if found, the couplings of such a state to each (B*+B0), (BB+)-B-*0 component as well as the molecular probabilities of each of the channels. We use the bootstrap method to determine these magnitudes and find that, with errors in the correlation function typical of present experiments, we can determine all these magnitudes with acceptable precision. In addition, the size of the source function of the experiment from where the correlation functions are measured can be also determined with a high precision.