Relativistic viscosity has the following applications: change in the mean particle momentum and spreading around the mean for the cosmic rays (CR) ; acceleration of CR in shear flows, such as active galactic nuclei (AGN), gamma ray burst (GRB) and jets  ; interaction of a neutral particle with the microwave background radiation (CMB)  ; and generation of the CMB in the expanding universe . The Lagrangian and a Hamiltonian for a relativistic particle moving in a dissipative medium characterized by a force that depends on the square of the velocity of the particle have been derived . This paper is a highly idealized attempt to model SN light curves by assuming that the resistivity of the ambient interstellar medium is quadratic to the velocity of the SN envelope. Section 2 derives a relativistic equation of motion in the presence of viscosity proportional to the square of the velocity. Section 3 applies the relativistic results to the motion of SN 1993J, to the light curve of GRB 130427A and GRB 060729. The CMB is not related with the model that is presented here.
2. The Equation of Motion
We assume a one dimensional motion with a resistive force of Newtonian type, , where B is a constant, is the considered mass and is the velocity. In the following, we will consider only positive and decreasing velocities. Newton’s second law in special relativity is:
where F is the force, p is the relativistic momentum, m is the relativistic mass, is the rest mass, c is the velocity of light, is the velocity and is the Lorentz factor; see equation (7.16) in . The first order differential equation in the velocity that governs the motion is
The solution to this first order differential equation for in an implicit form is
where is the velocity at . An explicit solution for the velocity can be obtained by considering the physical solution of the previous algebraic equation of fourth degree in v
Before we continue, we will introduce the following simplification.
Conjecture 1. In the presence of more than one solution for the temporal evolution of the velocity/space, we select the physical one that has a positive decreasing/increasing behavior. According to this statement, the physical solution is
We already know that at . Once we know that at , it is possible to derive the unknown parameter B from the previous formula
At this moment, we are unable to obtain an analytical solution for . Therefore, the trajectory is obtained by a numerical integration of Equation (8).
An approximation for the trajectory is represented by a Taylor series expansion about of order 2
which means the following approximate velocity as function of time
We now analyze the mildly relativistic case for which the relativistic momentum in the case of a unit mass is
which means the following first order differential equation for the velocity
The velocity in the mildly relativistic case is
The indefinite integral of the velocity is
which leads to the following trajectory in the mildly relativistic case
The unknown parameter B in the mildly relativistic case is
The Relativistic Luminosity
The relativistic transfer of energy through a surface, A, is
wherep is the pressure. For sake of simplicity, we take p = 0, and is the Lorentz factor, see eqn. A31 in  or eqn. (43.44) in . In the case of a spherical cold expansion
We now assume the following power-law behavior for the density in the advancing layer of radius r
which has the following temporal scaling
where is given by the Taylor series represented by Equation (10). The mechanical relativistic luminosity is
Once we insert formulae (10) and (11) in this equation, we obtain an approximate expression for the mechanical relativistic energy, which has a complicated expression that is not here reported. We now parameterize the presence of the absorption introducing a the optical thickness . The observed luminosity is assumed to be
where , which is here assumed to be constant in the interval of time considered here, allows the match between theory and observations. The optical thickness takes the value in the case of an optically thin medium or can be function of time to simulate the complex behavior of the observed luminosity; more details can be found in . The observed absolute magnitude is
where is a constant of match between theory and observations.
3. Astrophysical Applications
The chosen astrophysical units are pc for length and yr for time. With these units, the initial velocity is .
As a first example, we apply the above results to the deduction of the parameter B to SN 1993J for which observational times and velocities are available  , see Table 1. A test for the quality of the fits is represented by the merit function
where , and are the theoretical radius, the observed radius and the observed uncertainty, respectively. Once B is derived on an observational basis, Figure 1 reports the analytical velocity and Figure 2 reports the numerical trajectory, with data as given in Table 1.
A comparison between the numerical solution with the Runge Kutta method and a Taylor series is reported in Figure 3.
The second example is applied to the light curve (LC) of GRB 130427A, which was the most luminous gamma-ray burst in the last 30 years, see Figure 1 in . Figure 4 reports the X-flux as a function of time and the relative theoretical
Table 1. Numerical values of the parameters for the fit and the theoretical model applied to SN 1993J.
Figure 1. Analytical velocity versus time (full line) for SN 1993J as given by formula (8).
Figure 2. Numerical radius (full line) and astronomical data of SN 1993J with vertical error bars.
data with data as in Table 2.
The third example is dedicated to the LC for GRB 060729 as observed by the Ultraviolet and Optical Telescope (UVOT) in the time interval [10-2-26] days, see Figure 1 in . Figure 5 presents the LC of UVOT (U) apparent magnitude for GRB 060729 with data as in Table 3, and Figure 6 presents the temporal
Figure 3. Numerical radius (blue line) and Taylor solution (red line) as represented by Equation (10).
Figure 4. Flux in the X-ray as function of time in seconds for GRB 130427A (empty stars) and theoretical curve as given by Equation (25) (full line) when with data as in Table 2.
Table 2. Numerical values of the parameters for the theoretical model applied to GRB 130427A.
Table 3. Numerical values of the parameters for the theoretical model applied to GRB 060729.
Figure 5. The LC of UVOT (U) + HST (F330W) for GRB 060729 (empty stars) and theoretical curve with radius as given by the relativistic numerical model with data as in Table 3. The theoretical magnitude is given by Equation (26) (full line).
Figure 6. The time dependence of (empty stars) for GRB 060729 and a logarithmic polynomial approximation of degree 5 (full line). Parameters as in Table 3.
behavior of the optical depth.
Physics. We have derived velocity as a function of time for a relativistic neutral particle that moves in a dissipative medium in the presence of friction, which depends on the square velocity of the particle. The trajectory (i.e., space as a function of time) can be deduced by numerical integration of the velocity or by a Taylor series of the differential equation of the second order.
Astrophysics. An application of the obtained results to SN 1993J allows us to derive the constant B and draw a comparison between the observed and theoretical trajectory, see Table 1 and Figure 2. The LCs of GRB 130427A and GRB 060729 were simulated, see Figure 4 and Figure 5.
 Rieger, F. and Duffy, P. (2007) Cosmic-Ray Acceleration and Viscosity. In: Aschenbach, B., Burwitz, V., Hasinger, G. and Leibundgut, B., Eds., Relativistic Astrophysics Legacy and Cosmology—Einstein’s, Springer, Berlin, 477-479.
 Webb, G.M., Barghouty, A.F., Hu, Q. and le Roux, J.A. (2018) Particle Acceleration Due to Cosmic-Ray Viscosity and Fluid Shear in Astrophysical Jets. The Astrophysical Journal, 855, 31.
 Dedkov, G. and Kyasov, A. (2010) Tangential Force and Heating Rate of a Neutral Relativistic Particle Mediated by Equilibrium Background Radiation. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 268, 599.
 Kolekar, S., Shankaranarayanan, S. and Chitre, S.M. (2020) Generation of CMB and Cosmological Constant via Bulk Viscosity. General Relativity and Gravitation, 52, 98.
 López, G.V., Montes, G.C. and Zanudo, J.G.T. (2015) One Dimensional Relativistic Free Particle in a Quadratic Dissipative Medium. Journal of Modern Physics, 6, 121-125.
 Zaninetti, L. (2021) Energy Conservation in the Thin Layer Approximation: IV. The Light Curve for Supernovae. International Journal of Astronomy and Astrophysics, 11, 37.
 Marcaide, J.M., Mart-Vidal, I., Alberdi, A. and Pérez-Torres, M.A. (2009) A Decade of SN 1993J: Discovery of Radio Wavelength Effects in the Expansion Rate. A&A, 505, 927.
 Mart-Vidal, I., Marcaide, J.M., Alberdi, A., Guirado, J.C., Pérez-Torres, M.A. and Ros, E. (2011) Radio Emission of SN1993J: The Complete Picture. II. Simultaneous Fit of Expansion and Radio Light Curves A & A, 526, A143.
 De Pasquale, M., Page, M., Kann, D., Oates, S., Schulze, S., Zhang, B., Cano, Z., Gendre, B., Malesani, D., Rossi, A., Gehrels, N., Troja, E., Piro, L., Boër, M. and Stratta, G. (2017) Challenging the Forward Shock Model with the 80 Ms Follow up of the X-Ray Afterglow of Gamma-Ray Burst 130427A. Galaxies, 5, 6.