1. Introduction
The equation of motion plays a relevant role in our understanding of the physics of the Herbig-Haro objects (HH) after [1] [2]. A common example is to evaluate the velocity of the jet in HH34 as 300 km/s, see [3] , without paying attention to its spatial or temporal evolution. A precise evaluation of the evolution of the jet’s velocity with time in HH34 has been done, for example, by [4]. It is, therefore, possible to speak of proper motions of young stellar outflows, see [5] [6] [7] [8].
The first set of theoretical efforts exclude the magnetic field: [9] have modeled the slowing down of the HH 34 superjet as a result of the jet’s interaction with the surrounding environment, [10] have shown that a velocity profile in the jet beam is required to explain the observed acceleration in the position-velocity diagram of the HH jet, [11] found some constraints on the physical and chemical parameters of the clump ahead of HHs and [12] reviewed some important understanding of outflows from young stars.
The second set of theoretical efforts include the magnetic field: [13] analysed the HH 1-2 region in the L1641 molecular cloud and found a straight magnetic field of about 130 micro-Gauss; [14] analysed HH 211 and found field lines of the magnetic field with different orientations; [15] analysed HH 111 and found evidence for magnetic braking.
These theoretical efforts to understand HH objects leave a series of questions unanswered or partially answered, as follows:
• Is it possible to find a law of motion for turbulent jets in the presence of a medium with a density that decreases as a power law?
• Is it possible to introduce the back reaction into the equation of motion for turbulent jets to model the radiative losses?
• Can we model the bending of the super-jet connected with HH34?
• Can we explain the bow shock visible in HH34 with the theory of the image?
To answer these questions, this paper reviews in Section 2 the velocity observations of HH34 at a 9 yr time interval, Section 3 analyses two simple models as given by the Stoke’s and Newton’s laws of resistance, Section 4 applies the conservation of the energy flux in a turbulent jet to find an equation of motion, Section 5 models the extended region of HH34, the so called “superjet”, and Section 6 reports some analytical and numerical algorithms that allow us to build the image of HH34.
2. Preliminaries
The velocity evolution of the HH34 jet has recently been analysed in [SII] 2, (672 nm), frames and Table 1 in [4] reports the Cartesian coordinates, the velocities, and the dynamical time for 18 knots in 9 years of observations. To start with time, t, equal to zero, we fitted the velocity versus distance with the following power law
(1)
where v and x are the velocity and the length of the jet,
is the velocity at
and
with its relative error is a parameter to be found with a fitting procedure. The integration of this equation gives the time as a function of the position, x, as given by the fit
(2)
where
is the position at
. The fitted trajectory, distance versus time, is
(3)
and the fitted velocity as function of time is
(4)
Table 1. Numerical values for the physical parameters of HH34.
The adopted physical units are pc for length and year for time, and the useful conversion for the velocity is
.
The fit of Equation (1) when x is expressed in pc gives
(5)
from which we can conclude that the velocity decreases with increasing distance, see Figure 1.
The time is derived from Equation (2) and Table 1 reports the basic parameters of HH34. This time is more continuous in respect to the dynamical time reported in column 6 of Table 1 in [4].
3. Two Simple Models
When a jet moves through the interstellar medium (ISM), a retarding drag force
, is applied. If v is the instantaneous velocity, then the simplest model assumes
(6)
where n is an integer. Here, the case of
and
is considered. In classical mechanics,
is referred to as Stoke’s law of resistance and
is referred to as Newton’s law of resistance.
Figure 1. Observational points of velocity in pc/yr versus distance in pc (empty circles) and best fit as given by Equation (1) (full line).
3.1. Stoke’s Behaviour
The equation of motion is given by
(7)
The velocity as function of time is
(8)
where
is the initial velocity. The distance at time t is
(9)
The time as function of distance is obtained by the inversion of this equation
(10)
The velocity as a function of space is
(11)
The numerical value of B is
(12)
where
is the velocity at point
; the data of Table 1 gives B = 0.0009549 (Figure 2).
3.2. Newton’s Behaviour
The equation of motion is
(13)
The velocity as function of time is
(14)
Figure 2. Observational points of velocity in pc/yr versus distance in pc (empty circles), best fit as given by Equation (1) (full line), Stokes behaviour as given by Equation (11) (dashed line) and Newton behaviour as given by Equation (17) (dot-dash-dot-dash).
where
is the initial velocity. The distance at time t is
(15)
The time as function of distance is obtained by the inversion of the above equation
(16)
The velocity as function of the distance is
(17)
The numerical value of A is
(18)
where
is the velocity at point
; the data of Table 1 gives A = 5.68381834.
4. Energy Flux Conservation
The conservation of the energy flux in a turbulent jet requires a perpendicular section to the motion along the Cartesian x-axis, A
(19)
where r is the radius of the jet. Section A at position
is
(20)
where
is the opening angle and
is the initial position on the x-axis. At position x, we have
(21)
The conservation of energy flux states that
(22)
where
is the velocity at position x and
is the velocity at position
, see Formula A28 in [16]. More details can be found in [17] [18]. The density is assumed to decrease as a power law
(23)
where
is the density at
and
a positive parameter. The differential equation that models the energy flux is
(24)
The velocity as a function of the position, x,
(25)
Figure 3 reports the velocity as a function of the distance and the observed points.
We now have four models for the velocity as a function of time and Table 2 reports the merit function
, which is evaluated as
Figure 3. Observational points of velocity in pc/yr versus distance in pc (empty stars). The theoretical fit is given by Equation (25) (full line) with parameters
,
and
.
Table 2. The values of the
for four models of velocity of HH34.
(26)
where
represents the observed value at position i and
the theoretical value at position i. A careful analysis of Table 2 allows us to conclude that the turbulent jet performs better in respect to the Stokes’s and Newton’s behavior.
The trajectory, i.e. the distance as function of the time,
(27)
and the velocity as function of time
(28)
Figure 4 reports the trajectory as a function of time and of the observed points.
The rate of mass flow at the point x,
, is
(29)
and the astrophysical version is
(30)
where
is the opening angle in rad, x and
are expressed in pc, n is the number density of protons at
expressed in particles cm-3,
is the solar mass and
is the initial velocity at point
expressed in km/s. This rate of mass flow as function of the distance x increases when
, is constant when
, and decreases when
.
Figure 4. Observational points of distance in pc versus time in years (empty stars). The theoretical curve is given by Equation (27) (full line) with the same parameters as in Figure 3.
The Back Reaction
Let us suppose that the radiative losses are proportional to the flux of energy
(31)
where
is a constant that is thought to be
. By inserting in the above equation the considered area,
, and the power law density here adopted the radiative losses are
(32)
By inserting in this equation the velocity to first order as given by Equation (25), the radiative losses,
, are
(33)
The sum of the radiative losses between
and x is given by the following integral, L,
(34)
The conservation of the flux of energy in the presence of the back-reaction due to the radiative losses is
(35)
The real solution of the cubic equation for the velocity to the second order, , is
(36)
Figure 5 reports the effect of introducing the losses on the velocity as function of the distance for a given value of, i.e. the velocity decreases more quickly.
The presence of the back-reaction allows us to evaluate the jet’s length, which can be derived from the minimum in the corrected velocity to second order as a function of x,
(37)
which is
(38)
The solution for x of the above minimum allows us to derive the jet’s length, ,
(39)
Figure 6 reports an example of the jet’s length as a function of the parameter.
Figure 5. Velocity to the second order as function of the distance, see Equation (36), when (full line) and (dashed line), other parameters as in Figure 3.
Figure 6. Length of the jet as a function of, and the other parameters are as in Figure 3.
5. The Extended Region
To deal with the complex shape of the continuation of HH34 (e.g. see the new region HH173 discovered by [19] ), we should include the precession of the source and motion of the host star, following a scheme outlined in [20]. The various coordinate systems are, , ,. The vector representing the motion of the jet is represented by the following matrix:
(40)
where the jet motion L(t) is considered along x axis.
The jet axis, x, is inclined at an angle relative to an axis, and therefore the matrix, which represents a rotation through z axis, is given by:
(41)
The jet is undergoing precession around the axis and is the angular velocity of precession expressed in radians per unit time. The transformation from the coordinates fixed in the frame of the precessing jet to the nonprecessing coordinate is represented by the matrix
(42)
The last translation represents the change of the framework from (), which is co-moving with the host star, to a system (), in comparison to which the host star is in a uniform motion. The relative motion of the origin of the coordinate system is defined by the Cartesian components of the star velocity, and the required matrix transformation representing this translation is
(43)
On assuming, for the sake of simplicity, that and, the translation matrix becomes
(44)
The final matrix A representing the “motion law” can be found by composing the four matrices already described;
(45)
The three components of the previous matrix A represent the jet’s motion along the Cartesian coordinates as given by an observer who sees the star moving in a uniform motion. The point of view of the observer can be modeled by introducing the matrix E, which represents the three Eulerian angles, see [21]. A typical trajectory is reported in Figure 7 and a particularised point of view of the same trajectory is reported in Figure 8 in which a loop is visible.
6. Image Theory
This section summarises the continuum observations of HH34, reviews the transfer equation with particular attention to the case of an optically thin layer, analyses a simple analytical model for theoretical intensity, reports the numerical algorithm that allows us to build a complex image and introduces the theoretical concept of emission from the knots.
6.1. Observations
The system of the jet and counter jet of HH34 has been analysed at 1.5 μm and 4.5 μm, see Figure 3 in [3]. The intensity is almost constant, for the first 12" of the jet and for the first 20" of the counter jet. At larger distances, the intensity drops monotonically. At a distance of 414 pc as given by [4] the conversion between physical and angular distance is. For example, at 1.5 μm, the emission is mainly due to the [Fe II]1.64 μm line.
6.2. The Transfer Equation
For the transfer equation in the presence of emission only see, for example, [22] or [23] , is
(46)
where is the specific intensity, s is the line of sight, is the emission coefficient, is a mass absorption coefficient, is the density of mass at position s, and the index denotes the frequency of emission. The solution to Equation (46) is
(47)
where is the optical depth at frequency:
(48)
We now continue to analyse a case of an optically thin layer in which is very small (or is very small) and where the density is replaced by the concentration of the emitting particles:
(49)
where K is a constant. The intensity is now
(50)
which in the case of constant density, C, is
(51)
The increase in brightness is proportional to the concentration of particles integrated along the line of sight.
6.3. Theoretical Intensity
The flux of observed radiation along the centre of the jet, is assumed to scale as
(52)
where Q, the radiative losses, is given by Equation (33). The explicit form of this equation is
(53)
This relation connects the observed intensity of radiation with the rate of energy transfer per unit area. A typical example of the jet of HH34 at 4.5 μm is reported in Figure 9.
6.4. Emission from a Cylinder
A thermal model for the image is characterised by a constant temperature and density in the internal region of the cylinder. Therefore, we assume that the number density C is constant in a cylinder of radius a and then falls to 0, see the simplified transfer Equation (51). The line of sight when the observer is situated at the infinity of the x-axis and the cylinder’s axis is in the perpendicular position is the locus parallel to the x-axis, which crosses the position y in a Cartesian x-y plane and terminates at the external circle of radius a. A similar treatment for the sphere is given in [24]. The length of this locus in the optically thin layer approximation is
(54)
The number density is constant in the circle of radius a and therefore the intensity of the radiation is
(55)
A typical example of this cut is reported in Figure 10 and the intensity of all the cylinder is reported in Figure 11.
6.5. Numerical Image
The numerical algorithm that allows us to build a complex image in the optically thin layer approximation is now outlined.
• An empty, value = 0, memory grid which contains 4003 pixels is considered.
• The points which fill the jet in a uniform way to simulate the constant density in the emitting particles are inserted, value = 1, in
• Each point of has spatial coordinates which can be represented by the following matrix, A,
Figure 11. 2D map of the intensity of a jet which has length 0.1 pc and radius of 0.01 pc.
(56)
The orientation of the object is characterised by the Euler angles and therefore by a total rotation matrix, E, see [21]. The matrix point is represented by the following matrix, B,
(57)
• The intensity 2D map is obtained by summing the points of the rotated images.
A typical result of the simulation is reported in Figure 12, which should be compared with the observed image as given by Figure 13.
Figure 12. 2D intensity map of HH34, parameters as in Figure 8.
Figure 13. Three-color composite image of the young object HH34.
6.6. The Mathematical Knots
The trefoil knot is defined by the following parametric equations:
(58)
(59)
(60)
with. The visual image depends on the Euler angles, see Figure 14.
The image in the optically thin layer approximation can be obtained by the numerical method developed in Section 6.5 and is reported in Figure 15.
This 2D map in the theoretical intensity of emission shows an enhancement where two mathematical knots apparently intersect.
Figure 15. Image of the trefoil with parameters as in Figure 14, the side of the box in pc is 1 and the radius of the tube in pc is 0.006.
7. Conclusions
Laws of motion:
We analysed two simple models for the law of motion in HH objects as given by the Stoke’s and Newton’s behaviour, see Section 3. A third law of motion is used for turbulent jets in the presence of a medium whose density decreases with a power law, as given by Equation (23). The model that is adopted for the turbulent jets conserves the flux of energy. For example, Equation (25) reports the velocity as function of the position. The analysis for observed theoretical velocity as function of time/space, see Table 2, assigns the smaller value to the turbulent jet.
Back reaction:
The insertion of the back reaction in the equation of motion allows us to introduce a finite rather than infinite jet’s length, see Equation (39).
The extended region:
The extended region of HH34 is modeled by combining the decreasing jet’s velocity with the constant velocity and precession of the central object, see the final matrix (45).
The theory of the image:
We have analysed the case of an optically thin layer approximation to provide an explanation for the so called “bow shock” that is visible in HH34. This effect can be reproduced when two emitting regions apparently intersect on the plane of the sky, see the numerical simulation as given by Figure 12. This curious effect of enhancement in the intensity of emission can easily be reproduced when the image theory is applied to the mathematical knots, see the example of the trefoil in Figure 15.
Acknowledgements
Credit for Figure 13 is given to ESO.