Issue 
A&A
Volume 601, May 2017



Article Number  A87  
Number of page(s)  7  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201630229  
Published online  10 May 2017 
The jetdisk symbiosis without maximal jets: 1D hydrodynamical jets revisited
^{1} Anton Pannekoek Institute for Astronomy, University of Amsterdam, PO Box 94249, 1090 GE, Amsterdam, The Netherlands
email:
p.k.crumley@uva.nl
^{2} Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA
^{3} Mathematical Sciences and STAG Research Centre, University of Southampton, SO17 1BJ, UK
Received: 11 December 2016
Accepted: 2 March 2017
In this work we discuss the recent criticism by Zdziarski (2016, A&A, 586, A18) of the maximal jet model derived in Falcke & Biermann (1995, A&A, 293, 665). We agree with Zdziarski that in general a jet’s internal energy is not bounded by its restmass energy density. We describe the effects of the mistake on conclusions that have been made using the maximal jet model and show when a maximal jet is an appropriate assumption. The maximal jet model was used to derive a 1D hydrodynamical model of jets in agnjet, a model that does multiwavelength fitting of quiescent/hard state Xray binaries and lowluminosity active galactic nuclei. We correct algebraic mistakes made in the derivation of the 1D Euler equation and relax the maximal jet assumption. We show that the corrections cause minor differences as long as the jet has a small opening angle and a small terminal Lorentz factor. We find that the major conclusion from the maximal jet model, the jetdisk symbiosis, can be generally applied to astrophysical jets. We also show that isothermal jets are required to match the flat radio spectra seen in lowluminosity Xray binaries and active galactic nuclei, in agreement with other works.
Key words: black hole physics / galaxies: jets / hydrodynamics / stars: black holes / galaxies: active
© ESO, 2017
1. Introduction
In this work we reexamine the maximal jet model of accreting blackholes with disks and jets derived in Falcke & Biermann (1995). The idea behind a maximal jet is that there is a strict upper limit on the power carried by a jet in terms of its mass flux. As first pointed out by Zdziarski (2016), the maximal jet is the result of an erroneous conclusion that the internal energy density of a gas must be less than or equal to the restmass energy density of the gas. The maximal jet model was used to link the power carried by the jet to the power in the accretion disk, arriving at the highly influential and widely used “jetdisk symbiosis”. The main result of the jetdisk symbiosis states that the total power in the jet, L_{j}, can be related to the mass accretion rate of the disk, Ṁ_{disk}, through the jet’s Lorentz factor γ_{j} and an efficiency, η< 1. .
Without appealing to a maximal jet, we argue in this paper that the jetdisk symbiosis is reasonable for astrophysical accreting blackholes in general. In fact, it is reasonable to estimate the power of the jet to be , as long as one takes η ≲ a few, as opposed to being strictly less than one. It is intuitive why this conclusion should hold in a BlandfordPayne type jet where the disk itself is powering the jet (Blandford & Payne 1982), but the conclusion should also hold for jets that are powered by the black hole’s rotational energy via the BlandfordZnajek mechanism (Blandford & Znajek 1977). In the BlandfordZnajek mechanism, the jet power is proportional to the poloidal magnetic flux at the black hole, and the magnetic flux that can be carried to the black hole is ultimately limited by the massaccretion rate (Narayan et al. 2003; Tchekhovskoy et al. 2011). This explains why the jetdisk symbiosis has been such a successful concept.
However, the maximal jet conclusion results solely from an algebraic mistake in Falcke & Biermann (1995) and cannot be applied broadly to accreting blackholes with jets. We argue in this work that if the jet is efficiently accelerated and has a small terminal Lorentz factor, the initial enthalpy should be roughly equal to the restmass energy, in agreement with a maximal jet. A jet with a large terminal Lorentz factor will start with a large enthalpy, but if the jet is efficiently accelerated, the internal energy of the jet will be approximately equal to the restmass energy density after the jet has reached its final Lorentz factor (Vlahakis & Königl 2003; McKinney 2006). There are some cases when this approximation breaks down. For a steadystate, axisymmetric, magneticallyaccelerated outflow to be efficiently accelerated, it must stay causally connected in the transverse direction. However, it is unclear if this requirement is an actual impediment to the magnetic acceleration of astrophysical jets, or an artifact of the symmetries imposed. (For a concise review of magnetic acceleration of jets, see Komissarov2011.) A radiation or thermal pressuredriven jet will be efficiently accelerated even if it is conical and freestreaming. We assume the jet opening angle is small enough to ensure it remains in causal contact with the external medium throughout this paper.
In addition, we correct algebraic errors made in the derivation of the jet’s Lorentz factor as a function of distance in Falcke (1996). This Lorentz factor profile is used to calculate the dynamics of a jet in agnjet, an outflow dominated model of lowluminosity accreting black holes that has been applied to several different lowluminosity active galactic nuclei and Xray binaries (Markoff et al. 2005, 2008, 2015; Maitra et al. 2009; van Oers et al. 2010; Plotkin et al. 2015; Connors et al. 2017). The physics behind agnjet is presented in detail in Markoff et al. (2005) & Maitra et al. (2009). We characterize how the changes to the Lorentz factor profile effect the resulting spectral energy distribution calculated by agnjet.We find that the aforementioned algebraic mistakes have a negligible effect on the radiation from the outflow as long as the jet is roughly isothermal, has a small Lorentz factor, and is launched with an aspect ratio of order unity.
Our paper is organized as follows: in Sect. 2 we outline the mistake made when the maximal jet was derived. Then, we argue why the main conclusion of the maximal jet model, the jetdisk symbiosis, still holds. We also describe the systems for which the maximal jet model can be applied. In Sect. 3, we reexamine 1D pressuredriven jets, relaxing the maximal jet requirement. In Sect. 4, we make the dynamics of agnjet selfconsistent, and find the combined effects of all the changes to the dynamics on the calculated spectral energy distribution are small. We end by summarizing and discussing our results.
2. Bernoulli’s equation and maximal jets
The total power of an axisymmetric, conical jet at height z from the launching point with an opening angle θ is equal to the jet’s Lorentz factor γ_{j} times the enthalpy flux, (1)where ω is the enthalpy. In ideal magnetohydrodynamics (MHD), a baryonic jet with a comoving number density of protons n, has an enthalpy given by (2)U_{j}, P_{j} are the total energy density and pressure of the jet, which can be broken down into a gas component (U_{th}, P_{th}) and a magnetic component (U_{B} = P_{B} = B^{2}/ 8π)^{1}. The jet’s gas pressure can be related to the internal energy of the gas via the adiabatic index, Γ, P_{th} = (Γ − 1)U_{th}. We define the magnetization parameter, σ, as B^{2}/ (4πnm_{p}c^{2}).^{2} The enthalpy then simplifies to (3)In Falcke & Biermann (1995), the authors assume that the magnetic fields are isotropically turbulent, and therefore can be treated as an ideal gas with adiabatic index Γ. The authors then write the enthalpy in terms of the total jet internal energy density U_{j} as (4)The authors then rewrite the enthalpy in terms of the sound speed, β_{s}, adiabatic index, and density. The sound speed is (5)(see e.g., Königl 1980). From Eqs. (4) and (5), one can derive the well known result that the maximal sound speed is . Substituting Eq. (5) into Eq. (4) and solving for ω yields (6)Compare the above equation to the equivalent equation in Falcke & Biermann (1995). They give a formula that is the approximation when β_{s} is small, (7)The mistake is that Falcke & Biermann (1995) then use Eq. (7) to argue ω must be less than 2nm_{p}c^{2} because β_{s} must be less than . Clearly this is wrong, because when using the correct formula, Eq. (6), ω diverges to infinity as . This mistake was first pointed out by Zdziarski (2016).
Since the total jet power is equal to the Lorentz factor times the enthalpy flux, the mistake leads to a maximal jet power , where Ṁ_{j} is the mass flux through the jet. Protons are not created at the jet base, so Falcke & Biermann (1995) argue Ṁ_{j} ≤ Ṁ_{disk}, where Ṁ_{disk} is the mass accretion rate of the disk. Therefore, Falcke & Biermann (1995) conclude , η ≤ 1, i.e., the jetdisk symbiosis. The mass flux through the jet could be larger than Ṁ_{disk} if the jet has significant baryon loading from winds. If the disk drives a strong wind, the disk wind may be able to increase the baryon loading, but it is generally thought that in Xray binaries a strong disk wind is associated with the soft states, i.e., when the jet is not present (Neilsen & Lee 2009). New evidence shows simultaneous winds and jets in the high state, but there is no evidence that shows a strong wind during the low hardstate (MuñozDarias et al. 2016, 2017; Kalemci et al. 2016). Stellar winds could also play a role in increasing Ṁ_{j} in AGN or highmass Xray binaries (Komissarov 1994). At any rate, the baryon loading would only significantly increase the power carried by the jet when the power of the intercepted baryons exceeds that of the jet.
It is a reasonable approximation that the power of the jet does not significantly exceed even if the power in the jet is supplied by the spin of the black hole via the BlandfordZnajek mechanism. In the BlandfordZnajek mechanism, the power extracted is proportional to the magnetic flux carried to the black hole (Blandford & Znajek 1977). The magnetic flux carried to the blackhole is proportional to Ṁ_{BH}, and GRMHD simulations find that the resultant jet power exceeds by a factor of at most ~ 3 (Tchekhovskoy et al. 2011). However, there is no general requirement that the resultant jet’s power is dominated by the jet’s kinetic energy throughout the jet. For instance, in a Poyntingdominated jet, the power is and σ can be ≫ 1. Clearly, the enthalpy flux can be much larger than the mass flux, but in the next subsection, we will argue from completely different grounds than Falcke & Biermann (1995) that while a jet can start with any initial enthalpy, the jet will reach an ω/nm_{p}c^{2} ~ 2 at the modified magnetosonic fast point if it is efficiently accelerated. Furthermore, in efficiently accelerated jets that have mildly relativistic terminal Lorentz factors, i.e., γ_{j}β_{j} ~ 1, the initial enthalpy will not exceed the rest mass energy density by a significant amount. The reason why this is true can be seen most easily via Bernoulli’s equation.
2.1. Bernoulli’s equation
In a steadystate, conservative jet, the total power carried by the jet is constant along the jet. Dividing the power by another conserved quantity, the particle number flux through a crosssection of the jet, one rederives the relativistic Bernoulli’s equation (see e.g., Königl 1980): (8)Bernoulli’s equation is simply a statement that the energy of the fluid per particle must not change as the particles travel along a streamline, as long as particles are not created or destroyed. In this paper we are only considering selfsimilar jets with no gradient of the pressure in the toroidal direction of the jet. If a jet is launched with an initial Lorentz factor γ_{0} with an initial enthalpy per particle of ω_{0}/n_{0}, then because ω/n ≥ m_{p}c^{2}, it is clear from Eqs. (8) and (3) there is a maximal Lorentz factor: (9)If a jet achieves γ_{max}, it means it has converted 100% of its Poynting and thermal energy into bulk kinetic energy.Doing so in a steadystate, Poyntingdominated jet requires that the jet be causally connected in the transverse direction^{3}. If the jet goes out of causal contact in the transverse direction, the magnetic pressure may not be able to accelerate the jet further. If the flow is completely spherical with σ_{0} ≫ 1, the outflow reaches a terminal Lorentz factor (Goldreich & Julian 1970). If instead the jet has an opening angle θ_{j}, the terminal Lorentz factor is (Kumar & Zhang 2015).When the jet is thermally or radiatively driven, Eq. (9) should hold regardless of the jet’s geometry. Vlahakis & Königl (2003), Komissarov et al. (2007) show that in some Poynting dominated jets, electromagnetic fields will naturally selfcollimate the jet and ensure that the jet remains causally connected up to the modified magnetosonic fast point. In doing so the jet reaches rough equipartition between kinetic flux and Poynting flux, and the jet reaches a final Lorentz factor equal to .^{4} Therefore, using Bernoulli’s equation, the jet should have ω/nm_{p}c^{2} ~ 2 at the modified magnetosonic fast point, right after the jet has finished accelerating. Furthermore, ω_{0}/n_{0}m_{p}c^{2} ≈ 2 is likely to be a good assumption for mildlyrelativistic outflows with γ_{j}β_{j} of order unity, because if ω_{0} ≫ n_{0}m_{p}c^{2} the jets would reach γ_{j} that are too large. Mildly relativistic jets are expected to occur in quiescent/hard state blackhole Xray binaries^{5} (Gallo et al. 2003; MillerJones et al. 2007). A mildly relativistic jet may be launched by the supermassive blackhole at the center of our galaxy, Sgr A* (Falcke et al. 2009; Brinkerink et al. 2015).
Therefore, ω_{0}/n_{0} ~ 2m_{p}c^{2} should be a fine assumption as long as we restrict ourselves to mildly relativistic outflows with small enough opening angles. In jets with large Lorentz factors, γ_{j}β_{j} ≫ 5, like blazars, BL Lacs, relativistic tidal disruption events, and gammaray bursts, ω_{0}/n_{0} ≫ 2m_{p}c^{2}, but if the jet is accelerated efficiently, at the modified magnetosonic fast point ω/nm_{p}c^{2} ~ 2.^{6} The previous jet model of agnjet, and the one described in this paper, are not capable of reproducing the jets in these objects.
In the process of examining the maximal jet model, we have discovered additional minor algebraic mistakes in the Lorentz factor profile derived in Falcke (1996) and used in agnjet. In the rest of the paper, we show that these mistakes are minor and do not affect the conclusions drawn from fitting the model to low Lorentz factor sources. We then show that isothermal or nearly isothermal jets are required to have a flat radio spectrum.
3. Pressuredriven conical jets
In this section we reproduce with minor corrections the derivation of the onedimensional propagation of a quasiisothermal hydrodynamic jet from Falcke (1996). We find that the corrected quasiisothermal acceleration profile agrees within ~20% with the result in Falcke (1996), and the corrected Lorentz factor profile is much closer to the profile in a perfectly isothermal flow (see Fig. 1).
Fig. 1 Different Lorentz factor profiles: The profile derived in Falcke (1996) and used in agnjet cf. Markoff et al. (2005, yellow solid line), the profile derived in this work after correcting algebraic errors (Eq. ( 18), red solid line), and the profiles found by solving the 1D Euler equation of a conical jet assuming the gas is isothermal (green dashed line) and adiabatic (lilac solid line). If the jet follows Bernoulli’s equation, γ_{max}β_{max} ≈ 2.39 when γ_{0}β_{0} ≈ 0.485 and Γ = 4 / 3precisely the terminal value in the adiabatic jet Lorentz factor profile. 

Open with DEXTER 
The onedimensional propagation of a supersonic jet in the z direction follows the Euler equation given in Pomraning (1973), Falcke et al. (2009): (10)Following Falcke & Biermann (1995), we assume that the jet is welldescribed by a fluid with adiabatic index 4/3 and use the enthalpy from Eq. (4) (11)so Eq. (10) becomes: (12)Particle number conservation along the jet forces the number density to be (13)where θ is the opening angle of the jet (i.e., the crosssectional radius of the jet divided by the height of the jet).
The jet is launched at an initial height of z_{0}, and is assumed to be traveling at the sound speed. The sound speed in the jet’s rest frame, β_{s}, is (14)Instead of requiring U_{j} = nm_{p}c^{2} throughout the jet, we introduce a new parameter ζ, which is the ratio between the initial internal energy of the jet and the initial restmass energy density, i.e., U_{0} = ζn_{0}m_{p}c^{2}. ζ = 1 corresponds to the maximal jet conditions in Falcke (1996). The initial sound speed at the base of the jet is (15)For ζ = 1, Γ = 4 / 3 the sound speed is ≈ 0.43. Since the jet is assumed to be launched traveling at the sound speed, the initial velocity of the jet is (16)For Γ = 4 / 3, ζ = 1, γ_{0}β_{0} ≈ 0.485.
The jet Lorentz factor profile from Falcke (1996) can be derived by treating the jet as a conical jet (θ = θ_{0}) and using particle number conservation to get a z dependence on the density (see Eq. (13)). To fix the z dependence on the internal energy, we need a prescription of how the temperature changes with a change in density. If the jet is isothermal, T_{j} is constant and U_{e} ∝ n. If the jet is adiabatic, U_{e} ~ nkT_{j} ∝ n^{Γ}, and therefore T_{j} ∝ n^{Γ − 1} ∝ (γ_{j}β_{j})^{1 − Γ}z^{2 − 2Γ}. Falcke (1996) assumes that the gas is only able to do PdV work in the zdirection, hence the only adiabatic losses are due to the jet’s acceleration. We call this assumption quasiisothermal. If the jet is quasiisothermal, the temperature T_{j}, is proportional to (γ_{j}β_{j})^{1 − Γ}. It is difficult to understand how exactly the gas is prevented from doing PdV work in the lateral direction. It is far more realistic to assume there is continuous particle acceleration to counteract the adiabatic losses due to the expansion, as Blandford & Königl (1979) use to explain their isothermal jet model. For a heating mechanism to recover the quasiisothermal temperature dependence, it means it must be capable of compensating the large adiabatic losses due to the lateral expansion, but not the comparatively small adiabatic losses from the acceleration. Why the heating mechanism would do this is unclear. We retain T_{j} ∝ (γ_{j}β_{j})^{1 − Γ} here for historical reasons, and we note that when using the correct Euler equation, the difference between the quasiisothermal case and the isothermal case is negligible for the small Lorentz factors achieved in jets with U_{j,0} ~ n_{0}m_{p}c^{2} and Γ = 4 / 3, as assumed in this work.
In a quasiisothermal jet, U_{j} is (17)Substituting Eqs. (17) and (13) into Eq. (10), and assuming the jet is launched with an initial γ_{0}β_{0} equal to the sound speed (Eq. (16)), the 1D Euler equation that results is The above equation should reduce to the jet Lorentz factor profile used in Falcke (1996), Markoff et al. (2005) when ζ = 1. However, it differs from Eq. (2) in Falcke (1996): (20)(21)The difference between our equation and the equation in Falcke (1996) can be accounted for as follows: the − Γγ_{j}β_{j} term in Eq. (18) results from a neglected term, the difference in the exponent in ξ results from an arithmetic error, and finally the difference in the inside of the parenthesis of ξ terms is from setting instead of using the proper value given in Eq. (16). The difference between the solutions of Eqs. (18) and (20) are small and shown in Fig. 1. In Fig. 1, we also include solutions to the 1D Euler equations when the jet is isothermal (T_{j} = const., i.e., Eq. 20 with ξ = 1) and adiabatic (T_{j} ∝ (γ_{j}β_{j})^{1 − Γ}z^{2 − 2Γ}, see Eq. (25)).
The above quasiisothermal and isothermal solutions do not conserve energy, nor do they follow the Bernoulli equation. The violation of Bernoulli’s equation is clear by looking at maximal Lorentz factor on any pressure driven jet that conserves energy, with U_{0} = ζn_{0}m_{p}c^{2}, Eq. (9) becomes (22)or γ_{max} = 7γ_{0}/ 3 for a relativistic gas starting with equal parts internal energy density and rest mass energy density, U_{0} = n_{0}m_{p}c^{2}. All solutions except for the adiabatic solution eventually reach a γ_{j} that exceeds γ_{max}, which is easily seen in Fig. 1.
The total amount of heating needed to explain the solution in agnjet is equivalent to the increase in the jet power. Using Eq. (1) for the jet’s total power, the increase of power in a quasiisothermal jet is (23)The required heating to power a quasiisothermal jet is shown in Fig. 2. The heating must come from some internal process in the jet, but the heating mechanism is not capable of being captured by our timeindependent, laminar flow treatment in this work. One possibility is that heating originates from internal shocks (Malzac 2013). Internal shocks would do more than just heat the gas, they would also change the momentum and hence the dynamics. Magnetic reconnection can convert magnetic energy into thermal energy to keep the jet’s electrons isothermal, but reconnection would not increase the total power carried by the jet.
Fig. 2 Total power carried by a quasiisothermal jet as a function of distance from the black hole z normalized to its initial power. The power of the jet increases as it propagates (see Eq. (23)). One possibility to explain how the jet is heated is via some internal process not captured by the timeindependent Euler equation, such as internal shocks. 

Open with DEXTER 
The jet will conserve energy if instead of assuming T_{j} ∝ (γ_{j}β_{j})^{1 − Γ}, we use an adiabatic jet where U_{j} ∝ n^{Γ}, or T_{j} ∝ (γ_{j}β_{j})^{1 − Γ}z^{2 − 2Γ}. In an adiabatic conical jet the z dependence of the internal energy is: (24)and the full 1D Euler equation is (25)(26)The solution to the adiabatic 1D Euler equation when ζ = 1 is shown in Fig. 1. Unlike the quasiisothermal case, the jet reaches the maximal Lorentz factor equal to that predicted by the Bernoulli equation γ_{0}(1 + Γ).
4. Radiation, collimation, and a more selfconsistent agnjet
The agnjet model was developed in Markoff et al. (2005), Maitra et al. (2009) as way of fitting multiwavelength spectra of blackhole jets. The model is described in full detail in the aforementioned papers, but we will give a brief description here. In agnjet the gas is assumed to be moving at a velocity equal to the sound speed though a nozzle with constant radius r_{0} that ends at z_{0}. At z_{0}, the jet is allowed to expand at constant velocity equal to the sound speed^{7} in the crosssection radial direction, and weakly accelerated in the z direction by the jet’s pressure, i.e., it follows the 1D Euler equation in the z direction. Electrons initially have a MaxwellJuttner thermal distribution in the nozzle and the jet, and a fraction of the electrons are accelerated into a nonthermal population at a height z_{acc}. The electrons radiate via synchrotron and inverse Compton processes.
Fig. 3 How a different initial aspect ratio of the jet, r_{0}/z_{0}, changes the dynamics of a selfcollimated quasiisothermal jet compared to the conical jet considered in Fig. 1. As long as r_{0}/z_{0} ~ 1, the difference between a conical and selfcollimated jet is much less than a factor of 2. 

Open with DEXTER 
The assumption of a lateral expansion at constant speed while accelerating vertically results in a jet that is selfcollimating. The effect of the selfcollimation on the dynamics was not previously considered in agnjet. We show here that the effect is small as long as the initial crosssectional radius of the jet is roughly equal to the launch height.
The crosssectional radius of the jet, r, is assumed to follow (27)and conservation of number density of particles is (28)For the quasiisothermal case, i.e., what is currently used in agnjet, the internal energy profile is (29)(ζ = 1 in agnjet). With the crosssectional radius assumed in Eq. (27), the 1D Euler equation becomes: (30)and as before, (31)There is now another free parameter in the velocity profile: the initial aspect ratio of the jet: r_{0}/z_{0}. We show the effects of this new parameter and compare our results to a quasiisothermal conical jet in Fig. 3.
To divvy up the total internal energy into an electron energy density and magnetic energy density, agnjet uses a free parameter k, defined as the ratio of the magnetic energy density to the electron energy density. The dependence of the magnetic field on the height is (32)U_{j} has a height dependence that is different if the jet is assumed to be isothermal, quasiisothermal, or adiabatic. In the isothermal jet, B ∝ r^{1}(γ_{j}β_{j})^{− 1 / 2}, which is slightly slower than the expected dependence if there is flux conservation of a toroidal magnetic field (∝ r^{1}(γ_{j}β_{j})^{1}). In the adiabatic or quasiisothermal case, the magnetic field decreases faster than in an isothermal jet. The electron’s characteristic Lorentz factor γ_{e} has the same distance and bulk Lorentz factor dependence as the temperature. γ_{e,0} is a free parameter. The density profile of the electrons is determined by number conservation, but the initial number of electrons and positrons is fixed by requiring that U_{j,0} = n_{p}m_{p}c^{2}, (33)The above relationship will break down if k is too large, and it results in a nonphysical, electrostatically charged jet with n_{e}<n_{p}. k must satisfy the following inequality (34)If agnjet requires a k large enough to violate the above inequality to fit the spectrum of an object, the k is inconsistent with the model of agnjet. In this case, it likely means that the jet is Poynting dominated, i.e., U_{B}>U_{p}, a scenario which is not considered in agnjet. We note that versions of agnjet prior to 2014 did not use Eq. (33) to set the number of electrons and positrons, and instead simply set n_{e} = n_{p}. In the earlier version, agnjet was only selfconsistent if 1 + k was equal to the right hand side of Eq. (34). If we force n_{e} = n_{p}, as in the earlier version of agnjet, when 1 + k< 110 [ T_{e,0}/ (10^{11} K) ] ^{1}, the difference from a selfconsistent solution can be estimated by solving the Euler equation in Eq. (30) with a ζ = (1 + k) [ T_{e,0}/ (10^{11} K) ] / 110. The differences are likely to be minor even if ζ ≪ 1 because the jet is not accelerated very much even when ζ = 1. If in the previous version of agnjet a solution was found with 1 + k ≫ 110 [ T_{e,0}/ (10^{11} K) ] ^{1}, the jet is Poynting dominated, and its dynamics are not correctly calculated by agnjet.
We plot an example multiwavelength spectralenergy distribution using agnjet in Fig. 4. In the figure it is clear that the quasiisothermal, isothermal, and previous agnjet with minor errors all give roughly the same result: a nearly flat, selfabsorbed synchrotron spectrum at frequencies below the thermal bump, and a Compton hump in the Xrays. In Fig. 4, we assume the jet to be selfcollimating when calculating the dependence of the internal energy of the jet with height, but we use the velocity profiles from Fig. 1. Unlike the models with a constant or nearlyconstant temperature throughout the jet, the adiabatic model shows a steep rise from the selfabsorbed radio emission to the thermal synchrotron bump. We find that to have a flat radio spectra the jet must be kept at a nearly constant temperature, a similar conclusion as found by Blandford & Königl (1979) when fitting the radio cores of AGNs or as found by Mościbrodzka & Falcke (2013) when fitting the radio spectrum of Sgr A*, or as found by Pe’er & Casella (2009) when fitting the radio emission of Xray binaries.
5. Summary and discussion
In this work, we have reanalyzed the hydrodynamical jets derived in Falcke (1996). When deriving the Lorentz factor profile, Falcke (1996) used the maximal jet model from Falcke & Biermann (1995), that contains an algebraic error. The maximal jet assumption has two main conclusions – the jet’s power is dominated by its kinetic power, and the jet power must be less than the mass accretion rate of the disk onto the black hole times the bulk Lorentz factor of the jet. We argued in this letter that the second conclusion is likely true in general for astrophysical jets. Even in jets that extract their energy from the blackhole spin, the energy in the jet does not greatly exceed because the amount of magnetic flux that can be carried to the black hole depends on the mass accretion rate. The maximal efficiency of jet production found in GRMHD simulations is ≲ 300% (Tchekhovskoy et al. 2011; Nemmen & Tchekhovskoy 2015). Furthermore, while it is true that the jet’s power need not be dominated by its kinetic power, we argued that this is a good approximation in a jet with a sufficiently small opening angle and a small terminal bulk Lorentz factor. Xray binaries and lowluminosity AGN likely host such jets, so models based on these assumptions are wellfounded for such objects. In addition, we corrected minor algebraic mistakes made in the derivation of the Lorentz factor profile from Falcke (1996). The effects of correcting these errors are to make the quasiisothermal jet behave more similarly to a jet where the temperature is kept completely constant.
Fig. 4 Different spectral energy distribution (SED) calculated for conical jets in agnjet. The colors correspond to the same Lorentz factor profiles in Fig. 1. The models where the temperature is constant or nearly constant all show similar SED, while the adiabatic model has a steep rise below the thermal synchrotron bump. 

Open with DEXTER 
The Lorentz factor profile from Falcke (1996) was used in the agnjet model described in Markoff et al. (2005) & Maitra et al. (2009). agnjet has been used to fit the multiwavelength spectra of outflow dominated Xray binaries and nearby lowluminosity active galactic nuclei. In agnjet, the jet is assumed to be selfcollimating, but the collimation was not selfconsistently applied to the jet’s dynamics. We show that the effects of the selfcollimation are negligible for quasiisothermal jets as long as the aspect ratio of the jet at the launching point is of order unity. We examined the effects of the new Lorentz factor profile on the spectral energy distribution calculated by agnjet. We find there is very little difference between the assumed quasiisothermal jet and a completely isothermal jet, however, we find a large difference in the spectral energy distribution between an adiabatic jet and a jet where the temperature is kept roughly constant. We find that isothermal jets are required to match the flat radio spectra seen in hard/quiescent state XRB and lowluminosity AGN. We do not selfconsistently account for how the jet is kept hot. If the jet is heated internally, the heating mechanism would change the Lorentz factor profile from the one calculated in this work. If the gas is shock heated, the shocks will change the jet’s momentum, and if the jet is Poynting dominated, the dissipation of the magnetic field to heat the particles will change the magnetic pressure gradient.
Finally, we note that in addition to this work, there is an ongoing effort to derive a MHDconsistent jet model that will be capable of calculating the jet properties selfconsistently, e.g., Polko et al. (2010, 2013, 2014), Ceccobello et al. (in prep.). Such models will be able to address many of the shortcomings of the models described in this work.
The requirement that Poyntingdominated jets remain causally connected in the transverse direction can be relaxed if the outflow is impulsive (Granot et al. 2011).
Although see Heinz & Merloni (2004), MillerJones et al. (2006) regarding the difficulties in placing a strong upper limit on the Lorentz factor of Xray binaries.
Acknowledgments
P.C. would like to thank Heino Falcke, Peter Biermann, Sera Markoff, Pawan Kumar, Rodolfo Barniol Duran, and Thomas D. Russell for valuable discussions. P.C. acknowledges financial support from the WARP program of The Netherlands Organisation for Scientific Research (NWO). C.C. and R.C. acknowledge NWO/Nova. Y.C. is suported by the European Union’s Horizon 2020 research and innovation programme under the Marie SklodowskaCurie Global Fellowship grant agreement No. 703916.
References
 Blandford, R. D., & Königl, A. 1979, ApJ, 232, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433 [NASA ADS] [CrossRef] [Google Scholar]
 Brinkerink, C. D., Falcke, H., Law, C. J., et al. 2015, A&A, 576, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Connors, R. M. T., Markoff, S., Nowak, M. A., et al. 2017, MNRAS, 466, 4121 [NASA ADS] [Google Scholar]
 Falcke, H. 1996, ApJ, 464, L67 [NASA ADS] [CrossRef] [Google Scholar]
 Falcke, H., & Biermann, P. L. 1995, A&A, 293, 665 [NASA ADS] [Google Scholar]
 Falcke, H., Markoff, S., & Bower, G. C. 2009, A&A, 496, 77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gallo, E., Fender, R. P., & Pooley, G. G. 2003, MNRAS, 344, 60 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., & Julian, W. H. 1970, ApJ, 160, 971 [NASA ADS] [CrossRef] [Google Scholar]
 Granot, J., Komissarov, S. S., & Spitkovsky, A. 2011, MNRAS, 411, 1323 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Heinz, S., & Merloni, A. 2004, MNRAS, 355, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Kalemci, E., Begelman, M. C., Maccarone, T. J., et al. 2016, MNRAS, 463, 615 [NASA ADS] [CrossRef] [Google Scholar]
 Komissarov, S. S. 1994, MNRAS, 269, 394 [NASA ADS] [CrossRef] [Google Scholar]
 Komissarov, S. S. 2011, Mem. Soc. Astron. It., 82, 95 [NASA ADS] [Google Scholar]
 Komissarov, S. S., Barkov, M. V., Vlahakis, N., & Königl, A. 2007, MNRAS, 380, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Königl, A. 1980, Phys. Fluids, 23, 1083 [NASA ADS] [CrossRef] [Google Scholar]
 Kumar, P., & Zhang, B. 2015, Phys. Rep., 561, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Maitra, D., Markoff, S., Brocksopp, C., et al. 2009, MNRAS, 398, 1638 [NASA ADS] [CrossRef] [Google Scholar]
 Malzac, J. 2013, MNRAS, 429, L20 [NASA ADS] [CrossRef] [Google Scholar]
 Markoff, S., Nowak, M. A., & Wilms, J. 2005, ApJ, 635, 1203 [NASA ADS] [CrossRef] [Google Scholar]
 Markoff, S., Nowak, M., Young, A., et al. 2008, ApJ, 681, 905 [NASA ADS] [CrossRef] [Google Scholar]
 Markoff, S., Nowak, M. A., Gallo, E., et al. 2015, ApJ, 812, L25 [NASA ADS] [CrossRef] [Google Scholar]
 McKinney, J. C. 2006, MNRAS, 368, 1561 [NASA ADS] [CrossRef] [Google Scholar]
 MillerJones, J. C. A., Fender, R. P., & Nakar, E. 2006, MNRAS, 367, 1432 [NASA ADS] [CrossRef] [Google Scholar]
 MillerJones, J. C. A., Rupen, M. P., Fender, R. P., et al. 2007, MNRAS, 375, 1087 [NASA ADS] [CrossRef] [Google Scholar]
 Mościbrodzka, M., & Falcke, H. 2013, A&A, 559, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 MuñozDarias, T., Casares, J., Mata Sánchez, D., et al. 2016, Nature, 534, 75 [NASA ADS] [CrossRef] [Google Scholar]
 MuñozDarias, T., Casares, J., Mata Sánchez, D., et al. 2017, MNRAS, 465, L124 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2003, PASJ, 55, L69 [NASA ADS] [Google Scholar]
 Neilsen, J., & Lee, J. C. 2009, Nature, 458, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Nemmen, R. S., & Tchekhovskoy, A. 2015, MNRAS, 449, 316 [NASA ADS] [CrossRef] [Google Scholar]
 Pe’er, A., & Casella, P. 2009, ApJ, 699, 1919 [NASA ADS] [CrossRef] [Google Scholar]
 Plotkin, R. M., Gallo, E., Markoff, S., et al. 2015, MNRAS, 446, 4098 [NASA ADS] [CrossRef] [Google Scholar]
 Polko, P., Meier, D. L., & Markoff, S. 2010, ApJ, 723, 1343 [NASA ADS] [CrossRef] [Google Scholar]
 Polko, P., Meier, D. L., & Markoff, S. 2013, MNRAS, 428, 587 [NASA ADS] [CrossRef] [Google Scholar]
 Polko, P., Meier, D. L., & Markoff, S. 2014, MNRAS, 438, 959 [NASA ADS] [CrossRef] [Google Scholar]
 Pomraning, G. C. 1973, The equations of radiation hydrodynamics (Oxford: Pergamon Press) [Google Scholar]
 Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS, 418, L79 [NASA ADS] [CrossRef] [Google Scholar]
 van Oers, P., Markoff, S., Rahoui, F., et al. 2010, MNRAS, 409, 763 [NASA ADS] [CrossRef] [Google Scholar]
 Vlahakis, N., & Königl, A. 2003, ApJ, 596, 1080 [NASA ADS] [CrossRef] [Google Scholar]
 Zdziarski, A. A. 2016, A&A, 586, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Figures
Fig. 1 Different Lorentz factor profiles: The profile derived in Falcke (1996) and used in agnjet cf. Markoff et al. (2005, yellow solid line), the profile derived in this work after correcting algebraic errors (Eq. ( 18), red solid line), and the profiles found by solving the 1D Euler equation of a conical jet assuming the gas is isothermal (green dashed line) and adiabatic (lilac solid line). If the jet follows Bernoulli’s equation, γ_{max}β_{max} ≈ 2.39 when γ_{0}β_{0} ≈ 0.485 and Γ = 4 / 3precisely the terminal value in the adiabatic jet Lorentz factor profile. 

Open with DEXTER  
In the text 
Fig. 2 Total power carried by a quasiisothermal jet as a function of distance from the black hole z normalized to its initial power. The power of the jet increases as it propagates (see Eq. (23)). One possibility to explain how the jet is heated is via some internal process not captured by the timeindependent Euler equation, such as internal shocks. 

Open with DEXTER  
In the text 
Fig. 3 How a different initial aspect ratio of the jet, r_{0}/z_{0}, changes the dynamics of a selfcollimated quasiisothermal jet compared to the conical jet considered in Fig. 1. As long as r_{0}/z_{0} ~ 1, the difference between a conical and selfcollimated jet is much less than a factor of 2. 

Open with DEXTER  
In the text 
Fig. 4 Different spectral energy distribution (SED) calculated for conical jets in agnjet. The colors correspond to the same Lorentz factor profiles in Fig. 1. The models where the temperature is constant or nearly constant all show similar SED, while the adiabatic model has a steep rise below the thermal synchrotron bump. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.