The equation of motion plays a relevant role in our understanding of the physics of the Herbig-Haro objects (HH) after  . A common example is to evaluate the velocity of the jet in HH34 as 300 km/s, see  , 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 . It is, therefore, possible to speak of proper motions of young stellar outflows, see    .
The first set of theoretical efforts exclude the magnetic field:  have modeled the slowing down of the HH 34 superjet as a result of the jet’s interaction with the surrounding environment,  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,  found some constraints on the physical and chemical parameters of the clump ahead of HHs and  reviewed some important understanding of outflows from young stars.
The second set of theoretical efforts include the magnetic field:  analysed the HH 1-2 region in the L1641 molecular cloud and found a straight magnetic field of about 130 micro-Gauss;  analysed HH 211 and found field lines of the magnetic field with different orientations;  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.
The velocity evolution of the HH34 jet has recently been analysed in [SII] 2, (672 nm), frames and Table 1 in  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
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
where is the position at . The fitted trajectory, distance versus time, is
and the fitted velocity as function of time is
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
from which we can conclude that the velocity decreases with increasing distance, see Figure 1.
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
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
The velocity as function of time is
where is the initial velocity. The distance at time t is
The time as function of distance is obtained by the inversion of this equation
The velocity as a function of space is
The numerical value of B is
3.2. Newton’s Behaviour
The equation of motion is
The velocity as function of time is
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
The time as function of distance is obtained by the inversion of the above equation
The velocity as function of the distance is
The numerical value of A is
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
where r is the radius of the jet. Section A at position is
where is the opening angle and is the initial position on the x-axis. At position x, we have
The conservation of energy flux states that
where is the velocity at position x and is the velocity at position , see Formula A28 in . More details can be found in  . The density is assumed to decrease as a power law
where is the density at and a positive parameter. The differential equation that models the energy flux is
The velocity as a function of the position, x,
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.
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,
and the velocity as function of time
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
and the astrophysical version is
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
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
By inserting in this equation the velocity to first order as given by Equation (25), the radiative losses, , are
The sum of the radiative losses between and x is given by the following integral, L,
The conservation of the flux of energy in the presence of the back-reaction due to the radiative losses is
The real solution of the cubic equation for the velocity to the second order, , is
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,
The solution for x of the above minimum allows us to derive the jet’s length, ,
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  ), we should include the precession of the source and motion of the host star, following a scheme outlined in . The various coordinate systems are, , ,. The vector representing the motion of the jet is represented by the following matrix:
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:
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
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
On assuming, for the sake of simplicity, that and, the translation matrix becomes
The final matrix A representing the “motion law” can be found by composing the four matrices already described;
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 . 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.
Figure 7. Continuous trajectory of the superjet connected with HH34: the three Eulerian angles characterising the point of view are, and. The precession is characterised by the angle and by the angular velocity. The star has velocity, the considered time is 29,000 yr and the other parameters are as in Figure 3.
Figure 8. Continuous trajectory of the superjet connected with HH34: the three Eulerian angles characterising the point of view are, and. The other parameters as in Figure 7.
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.
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 . 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  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,  or  , is
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
where is the optical depth at frequency:
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:
where K is a constant. The intensity is now
which in the case of constant density, C, is
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
where Q, the radiative losses, is given by Equation (33). The explicit form of this equation is
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 . The length of this locus in the optically thin layer approximation is
The number density is constant in the circle of radius a and therefore the intensity of the radiation is
Figure 9. Observational points of intensity at 4.5 μm (empty stars) and theoretical curve as given by Equation (53) (full line). when, , and other parameters as in Figure 3.
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 10. 1D cut of the intensity, I, when.
Figure 11. 2D map of the intensity of a jet which has length 0.1 pc and radius of 0.01 pc.
The orientation of the object is characterised by the Euler angles and therefore by a total rotation matrix, E, see . The matrix point is represented by the following matrix, B,
• The intensity 2D map is obtained by summing the points of the rotated images.
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:
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 14. 3D view of the trefoil when the three Eulerian angles which characterise the point of view are, and.
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.
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.
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.
Credit for Figure 13 is given to ESO.
 Raga, A.C., Reipurth, B. and Noriega-Crespo, A. (2019) The HH34 Jet/Counterjet System at 1.5 and 4.5 μm. Revista Mexicana de Astronomia y Astrofisica, 55, 117-123.
 Raga, A.C., Noriega-Crespo, A., Rodrguez-González, A., et al. (2012) The Kinematics of HH 34 from HST Images with a Nine-Year Time Baseline. The Astrophysical Journal, 748, 103.
 Raga, A.C., Noriega-Crespo, A., Carey, S.J. and Arce, H.G. (2013) Proper Motions of Young Stellar Outflows in the Mid-Infrared with Spitzer (IRAC). I. The NGC 1333 Region. Astronomical Journal, 145, 28.
 Noriega-Crespo, A., Raga, A.C., Moro-Martn, A., Flagey, N. and Carey, S.J. (2014) Proper Motions of Young Stellar Outflows in the Mid-Infrared with Spitzer II HH 377/Cep E. New Journal of Physics, 16, Article ID: 105008.
 Guzmán, A.E., Garay, G., Rodrguez, L.F., Contreras, Y., Dougados, C. and Cabrit, S. (2016) A Protostellar Jet Emanating from a Hypercompact H II Region. The Astrophysical Journal, 826, Article No. 208.
 Raga, A.C., Reipurth, B., Esquivel, A., Castellanos-Ramrez, A., Velázquez, P.F., Hernández-Martnez, L., Rodrguez-González, A., Rechy-Garca, J.S., Estrella-Trujillo, D. and Bally, J. (2017) Proper Motions of the HH 1 Jet. Revista Mexicana de Astronomia y Astrofisica, 53, 485-495.
 López-Martn, L., Raga, A.C., López, J.A. and Meaburn, J. (2001) Theory and Observations of a Jet in the σ Orionis Region: HH 444. Revista Mexicana de Astronomia y Astrofisica Conference Series, 10, 61-64.
 Viti, S., Girart, J.M., Garrod, R., Williams, D.A. and Estalella, R. (2003) The Molecular Condensations Ahead of Herbig-Haro Objects. II. A Theoretical Investigation of the HH 2 Condensation. Astronomy & Astrophysics, 399, 187-195.
 Raga, A.C., Reipurth, B., Cantó, J., Sierra-Flores, M.M. and Guzmán, M.V. (2011) An Overview of the Observational and Theoretical Studies of HH 1 and 2. Revista Mexicana de Astronomia y Astrofisica, 47, 425-437.
 Kwon, J., Choi, M., Pak, S., Kandori, R., Tamura, M., Nagata, T. and Sato, S. (2010) Magnetic Field Structure of the HH 1-2 Region: Near-Infrared Polarimetry of Point-Like Sources. The Astrophysical Journal, 708, 758-769.
 Lee, C.F., Rao, R., Ching, T.C., Lai, S.P., Hirano, N., Ho, P.T.P. and Hwang, H.C. (2014) Magnetic Field Structure in the Flattened Envelope and Jet in the Young Protostellar System HH 211. The Astrophysical Journal, 797, L9.
 Lee, C.F., Hwang, H.C. and Li, Z.Y. (2016) Angular Momentum Loss in the Envelope-Disk Transition Region of the HH 111 Protostellar System: Evidence for Magnetic Braking? The Astrophysical Journal, 826, Article No. 213.
 Zaninetti, L. (2016) Classical and Relativistic Flux of Energy Conservation in Astrophysical Jets. Journal of High Energy Physics, Gravitation and Cosmology, 2, 41-56.
 Zaninetti, L. (2009) Scaling for the Intensity of Radiation in Spherical and Aspherical Planetary Nebulae. Monthly Notices of the Royal Astronomical Society, 395, 667-691.