Relativistic Motion with Viscosity. I Newton’s Law of Resistance ()
1. Introduction
Relativistic viscosity has the following applications: change in the mean particle momentum and spreading around the mean for the cosmic rays (CR) [1]; acceleration of CR in shear flows, such as active galactic nuclei (AGN), gamma ray burst (GRB) and jets [2] [3]; interaction of a neutral particle with the microwave background radiation (CMB) [4] [5]; and generation of the CMB in the expanding universe [6]. 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 [7]. 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:
(1)
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 [8]. The first order differential equation in the velocity that governs the motion is
(2)
The solution to this first order differential equation for
in an implicit form is
(3)
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
(4)
with
(5)
and
(6)
where
(7)
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
(8)
We already know that
at
. Once we know that
at
, it is possible to derive the unknown parameter B from the previous formula
(9)
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
(10)
which means the following approximate velocity as function of time
(11)
We now analyze the mildly relativistic case for which the relativistic momentum in the case of a unit mass is
(12)
which means the following first order differential equation for the velocity
(13)
The velocity in the mildly relativistic case is
(14)
where
(15)
The indefinite integral of the velocity is
(16)
where
(17)
which leads to the following trajectory in the mildly relativistic case
(18)
The unknown parameter B in the mildly relativistic case is
(19)
The Relativistic Luminosity
The relativistic transfer of energy through a surface, A, is
(20)
wherep is the pressure. For sake of simplicity, we take p = 0, and
is the Lorentz factor, see eqn. A31 in [9] or eqn. (43.44) in [10]. In the case of a spherical cold expansion
(21)
We now assume the following power-law behavior for the density in the advancing layer of radius r
(22)
which has the following temporal scaling
(23)
where
is given by the Taylor series represented by Equation (10). The mechanical relativistic luminosity is
(24)
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
(25)
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 [11]. The observed absolute magnitude is
(26)
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 [12] [13], 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 [14]. 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 [15]. 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.
4. Conclusions
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.