We present a 3D special-relativistic radiation hydrodynamics code. It uses the radiative inversion scheme with the M1-closure relation for the radiation equations, which allows the treatment of a wide range of optical depth, temperature and opacity. The radiation field is treated in the grey-body approximation. We present the standard 1D and 2D tests that include both optically thin and thick scenarios, as well as hydrodynamical and radiation pressure dominated configurations. As an application in 3D, we show the evolution of a jet driven by radiation-hydrodynamics with a helical perturbation. The code is expected to allow the exploration of scenarios in high-energy astrophysics where the radiation is important, like sources of GRBs.

We present the three-dimensional numerical simulation of jets modeled with Relativistic Radiation Hydrodynamics (RRH), that evolve across two environments: i) a stratified surrounding medium and ii) a 16TI progenitor model. We consider opacities consistent with various processes of interaction between the fluid and radiation, specifically, free-free, bound-free, bound-bound and electron scattering. We explore various initial conditions, with different radiation energy densities of the beam in hydrodynamical and radiation pressure dominated scenarios, considering only ultra-relativistic jets. In order to investigate the impact of the radiation field on the evolution of the jets, we compare our results with purely hydrodynamical jets. Comparing among jets driven by RRH, we find that radiation pressure dominated jets propagate slightly faster than gas pressure dominated ones. Finally, we construct the luminosity Light Curves (LCs) associated with all these cases. The construction of LCs uses the fluxes of the radiation field which is fully coupled to the hydrodynamics equations during the evolution. The main properties of the jets propagating on the stratified surrounding medium are that the LCs show the same order of magnitude as the gamma-ray luminosity of typical Long Gamma-Ray Bursts $10^{50}-10^{54}$erg/s and the difference between the radiation and gas temperatures is of nearly one order of magnitude. The properties of jets breaking out from the progenitor star model are that the LCs are of the order of magnitude of low-luminosity GRBs $10^{46}-10^{49}$ erg/s, and in this scenario the difference between the gas and radiation temperature is of four orders of magnitude, which is a case far from thermal equilibrium.

We present a method to classify initial conditions of a long gamma ray bursts model sourced by a single relativistic shock. It is based on the use of artificial neural networks (ANNs) that are trained with light curves (LC) generated with radiation relativistic hydrodynamics simulations. The model we use consists in a single shock with a highly relativistic injected beam into a stratified surrounding medium with profile $1/r^2$. In the process we only consider the bremsstrahlung radiation and Thomson scattering process. The initial conditions we use to train the ANN are three: the rest mass density, Lorentz factor and radiation energy density of the beam that produces the relativistic shock, together with the LC generated during the process. The classification selects the location of a box in the 3d parameter space that better fits a given LC, and in order to decrease the uncertainty of the parameters this box is refined and the classification selects a new box of smaller size.

We explore the possibility that a single relativistic shock, where the gas dynamics is coupled with radiation, can fit the light curves of long GRBs. For this we numerically solve the one dimensional relativistic radiation hydrodynamics equations with a single initial shock. We calculate light curves due to the evolution of this shock in terms of the velocity of the shock, the opacity of the gas, mass density and density of radiated energy. We explore how the variation of each of these parameters provides different features in the light curves. As examples we include the fitting of two long GRBs.

We present the dynamics of rotating Bose Condensate galactic dark matter halos, made of an ultralight spinless boson. We restrict to the case of adding axisymmetric rigid rotation to initially spherically symmetric structures and show there are three regimes: i) small angular momentum, that basically retains the drawbacks of spherically symmetric halos related to compactness and failure at explaining galactic RCs, ii) an intermediate range of values of angular momentum that allow the existence of long-lived structures with acceptable RC profiles, and iii) high angular momentum, in which the structure is dispersed away by rotation. We also present in detail the new code used to solve the Gross-Pitaevskii Poisson system of equations in three dimensions.

In this paper we show that spherically symmetric BEC dark matter halos, with the $\sin r/r$ density profile, that accurately fit galactic rotation curves and represent a potential solution to the cusp-core problem are unstable. We do this by introducing back the density profiles into the fully time-dependent Gross-Pitaevskii-Poisson system of equations. Using numerical methods to track the evolution of the system, we found that these galactic halos lose mass at an approximate rate of half of its mass in a time scale of dozens of Myr. We consider this time scale is enough as to consider these halos are unstable and unlikely to be formed. We provide some arguments to show that this behavior is general and discuss some other drawbacks of the model that restrict its viability.