arxiv: v1 [astro-ph.co] 27 Jul 2011

Transcription

1 Understanding Giant Radio Galaxy J : Large-Scale Morphology, Environment, and Energetics J. Machalski 1, M. Jamrozy 1, L. Stawarz 2,1, and D. Kozie l-wierzbowska 1 arxiv: v1 [astro-ph.co] 27 Jul Astronomical Observatory, Jagellonian University, ul. Orla 171, Kraków, Poland 2 Institute of Space and Astronautical Science, JAXA, Yoshinodai, Sagamihara, Kanagawa, , Japan machalsk@oa.uj.edu.pl, jamrozy@oa.uj.edu.pl, stawarz@astro.isas.jaxa.jp, eddie@oa.uj.edu.pl ABSTRACT In this paper we consider the possibility that the structure of the largest radio galaxy J is formed by a restarted rather than a primary jet activity. This hypothesis is motivated by the unusual morphological properties of the source, suggesting almost ballistic propagation of powerful jets in a particularly low-density environment. New radio observations of J confirm its morphology consisting of only two narrow lobes; no trace of any outer low-density cavity due to the previous jet activity is therefore detected. Different model fits performed using the newly accessed radio data imply relatively young age of the source, its exceptionally high expansion velocity, large jet kinetic power, and confirm particularly low-density environment. We find that it is possible to choose a realistic set of the model parameters for which the hypothetical outer lobes of J are old enough so that their expected radio surface brightness is substantially below the rms noise level of the available radio maps. On the other hand, the extremely low density of the gas surrounding the J lobes is consistent with the mean density of the baryonic matter in the Universe. This suggests that the source may be instead located in a real void of the galaxy and matter distribution. In both cases the giant radio lobes of J are expected to modify substantially the surrounding matter by driving strong shocks and heating the gas located at the outskirts of the filamentary galactic distribution. Finally, we also find that the energetic requirements for the source are severe in terms of the total jet power and the total energy deposited by the outflows far away from the central engine. Subject headings: radiation mechanisms: non-thermal galaxies: active galaxies: jets galaxies: individual (J )

2 2 1. Introduction Giant radio galaxies (hereafter GRGs) are luminous but low-surface brightness sources dominated at radio wavelengths by the emission of extended lobes with linear sizes 1Mpc. Morphologically, most of GRGs resemble powerful Fanaroff-Riley type II objects (FR IIs; Fanaroff & Riley 1974), and hence dynamical models developed for classical doubles are typically used to infer their intrinsic parameters based on the multifrequency radio data obtained using a variety of radio instruments (e.g., Konar et al. 2008; Jamrozy et al. 2008; Machalski et al. 2009). The key difficulty in modeling and understanding of such sources is however the fact that the main characteristics of the ambient medium surrounding giant lobes are hardly known. In fact, there are almost no observational constraints on the structure, temperature, or density of the intergalactic gas at Mpc distances from the centers of poor clusters or groups in which powerful radio galaxies, including GRGs, are expected to be located (Jamrozy et al. 2004; Saripalli et al. 2005). But one may also turn the argument around, and argue that GRGs can indeed be used as unique probes of the intergalactic medium (IGM) and of its cosmological evolution at the outskirts, or even outside of groups and clusters of galaxies. This idea was discussed and pushed forward by, e.g., Subrahmanyan et al. (2008) and Safouris et al. (2009). Clearly, using GRGs as cosmological probes requires a consensus regarding the nature of those objects. Are therefore GRGs just evolved FR II galaxies, for which large linear sizes are strictly due to their advanced ages of the order of 100Myr? Or are they instead intrinsically different from classical doubles(as suggested by, e.g., Kaiser & Alexander 1999)? A detailed analysis of various samples of GRGs favors the former possibility, suggesting however in addition that densities of the medium surrounding the discussed objects are lower than average for the whole FRII population (Mack et al. 1998; Schoenmakers et al. 2000b; Machalski & Jamrozy 2006). It is interesting to note in this context an analogy with the sources located at the other extremum of the size distribution of radio galaxies, namely Compact Symmetric Objects (CSOs), which spectroscopically are often classified as GigaHertz-Peaked Spectrum (GPS) sources. For these, the long-debated question was if their linear sizes kpc are due to young ages ( 0.1Myr) of otherwise regular FRIIlike radio structures, or instead are due to exceptionally dense ambient medium frustrating evolved jets and confining their lobes within host galaxies (see O Dea 1998, for a review). Or do they constitute yet another intrinsically distinct population of extragalactic radio sources? The merging agreement is that compact radio galaxies are indeed newly born precursors of classical doubles (and also of low-power FR I sources; see the discussion in Stawarz et al. 2008, and references therein), even though some of them may reside in particularly overdense environment (see, e.g., García-Burillo et al. 2007).

3 3 During the recent years it became clear that understanding the evolution of extragalactic radio sources including the youngest and the oldest examples of this class of astrophysical objects involves however not only accessing the properties of the ambient medium into which they evolve, or the main jet parameters (like jet total kinetic power, lifetime, etc.), but also characterizing a complex duty cycle of the jet activity. In particular, population studies of radio galaxies and quasars indicated that the jet activity must be intermittent, or at least highly modulated, in addition on different timescales (e.g., Reynolds & Begelman 1997; Alexander 2000; Schoenmakers et al. 2001; Czerny et al. 2009). The direct observational evidences for such an intermittency consist of different manifestations of the interrupted and/or restarted jets in radio-loud Active Galactic Nuclei (AGN), including interestingly some GRGs with compact inner structures of the CSO/GPS type (e.g., Schoenmakers et al. 1999; Konar et al. 2004, 2008; Saikia et al. 2007). This is an important point indeed, since it indicates that the appearance of large-scale jets and lobes is shaped not exclusively by the jet lifetime and the jet kinetic power relative to the density of the ambient medium, as believed previously, but is a complex convolution of these factors with the jet duty cycle. Note in particular that the restarted jets do not propagate within the undisturbed interstellar/intergalactic medium, but instead within the environment modified by the passage of the outflow during the previous epoch of the jet activity, and this may affect substantially the observed properties of the newly formed lobes. Double-double radio galaxies (hereafter referred to as DDRGs) are the most striking and clear examples of a highly intermittent activity of the jet engine with no realignment of the outflows direction in the two subsequent activity epochs (Lara et al. 1999; Schoenmakers et al. 2000a). These are the sources in which an extended radio structure comprises at least two pairs of lobes with a common core, both strongly co-aligned to within a few degrees to each other. At present, about 15 of such systems have been identified from either radio or X-ray observations (for a review see Saikia & Jamrozy 2009). The outer and also the inner lobes of DDRGs are uniquely characterized by an edge-brightened morphology of the FRII type, what led at the first place to the interpretation of the inner structures in DDRGs as regular lobes indeed, inflated by a plasma left behind a head (a termination shock) of a new-born jet which propagates through an outer cocoon formed during the previous epoch of the jet activity (Kaiser et al. 2000). The alternative interpretation for the inner structures of DDRGs would involve ballistic and naked (lobeless) jets observed in a few performed 2D and 3D hydrodynamical simulations of restarting outflows in AGN (Clarke & Burns 1991; Clarke 1997). Bulk of the observations and analysis carried till now for several best known DDRGs are in general consistent with the former scenario (see, e.g., Saripalli et al. 2002, 2003; Konar et al. 2006; Safouris et al. 2008; Machalski et al. 2010), even though some potential problems of the model has been recently discussed in

4 4 Brocksopp et al. (2007, 2011). The problems arise when the exact evolution of the inner lobes of DDRGs is analyzed. These inner lobes appear rather slim on the high-dynamics radio maps, displaying in particular axial ratios in the range of 9 16, to be compared with the range of 2 5 established for the outer lobes of same systems, and also for the lobes in regular FRII radio galaxies. This indicates that the restarted jets in DDRGs propagate within the environment much rarified with respect to the original IGM in which the outer lobes are evolving, as in fact expected in the framework of the jet intermittency scenario, but yet much denser than the one provided simply by the injection of a jet plasma to the outer lobes in the previous epoch of the jet activity (see, e.g., the detailed analysis by Safouris et al. 2008). Thus required contamination of the remnant cocoons by the additional matter is one of the central questions regarding the nature of DDRGs (see the discussion in Kaiser et al. 2000). It is not clear, however, if any entrainment process operates on sufficiently short timescales to enrich substantially the outer cocoons with the required amount of the cold matter from the IGM during the lifetime of a radio source. For this reason, Brocksopp et al. (2007, 2011) argued that radio emission observed from the inner structures of DDRGs is not due to the lobes of newly born jets, but predominantly due to bow shocks induced within the extended lobes by only briefly interrupted relativistic outflows. In this scenario no contamination of the extended/outer cocoon by the additional matter would be needed. In this paper we consider and analyze a possibility that the structure of the largest radio galaxyknownuptodate J withthetotalprojected linear sizeofd = 4.69Mpc, discovered by Machalski et al. (2008) is formed by a restarted, rather than a primary jet activity. The motivation for the performed analysis and new radio observations is due to the fact that the axial ratio for the lobes in this outstanding example of a giant, R T 12, is of the same order or even larger than axial ratios for the inner lobes in all well studied DDRGs, and therefore very different from the typical axial ratios characterizing normal size FR II sources or other known GRGs. Since in the framework of the standard dynamical model for classical doubles the axial ratio of the lobes depends predominantly on the jet kinetic power relative to the density of the surrounding, J must be a powerful source evolving in a particularly low-density environment, as indeed concluded in Machalski et al. (2008). The question is if the low-density of the matter surrounding giant lobes of J is due to a unique location of J in a void region, or due to some previous jet activity epoch which took place before the formation of the currently observed lobes, and which rarified the IGM thereby (in analogy with DDRGs). In the former case, J may be used as a particularly interesting (due to its enormous size) probe of a warm-hot phase of the intergalactic medium (hereafter WHIM) within a filamentary galaxy distribution far away from any group or cluster of galaxies. In the latter case, J would be the most

5 5 extreme example of a restarting radio galaxy, providing interesting constraints on the jet duty cycle and cosmological evolution of radio-loud AGN. New observations of J , conducted with the Giant Metrewave Radio Telescope (GMRT) and Very Large Array(VLA), confirm its morphology consisting of only two narrow lobes discerned in the VLA 1400MHz sky surveys FIRST (Becker et al. 1995) and NVSS (Condon et al. 1998), together with the radio core coincident with a parent galaxy at the redshift of z = The whole radio structure is highly collinear and symmetric: the arm ratiois1.08only, andthemisalignment angleis1.3. No traceofanyouter cocoonis however detected in our new observations. We modeled J using the DYNAGE algorithm of Machalski et al. (2007), along with the three bona-fide DDRGs treated as comparison and control sources. The resulting radio luminosity linear size (P ν D) evolutionary diagrams are used to determine conditions allowing the hypothetical outer lobes in the discussed object to dim in radio fluxes to a level precluding the detection at the given sensitivity limits. The archival data regarding the exemplary DDRGs used in the modeling, as well as the new GMRT and VLA observations of J resulting in high-resolution low- and high-frequency maps of the source, are presented in 2 of the paper. The applied model and the obtained model results are described in 3. Further discussion regarding the jet energetics and the properties of the environment of J are given in 4. Final conclusions are summarized in 5. Distances, linear sizes, volumes, and luminosities of the analyzed sources are calculated for the modern ΛCDM cosmology (Ω m = 0.27, Ω Λ = 0.73, and H 0 = 71kms 1 Mpc 1 ). 2. The Data In order to test the hypothesis that the observed radio structure of J results from a secondary jet activity in the parent AGN, we gathered new GMRT and VLA data. These allow us to model the source using the DYNAGE algorithm (see 3.1), along with the three carefully selected DDRGs, treated as comparison sources, for which the outer and inner lobes are observed and well mapped at a number of radio frequencies. The results of the newly conducted observations of J as well as the details regarding the selected DDRGs are given below. Main properties of the studied sources are summarized in Table 1.

6 GMRT and VLA Observations of J The analysis of J presented in this paper is based on the radio observations conducted with the GMRT and VLA, as well as on the archival VLA data from the NRAO 1400MHzFIRSTandNVSSsurveys (Becker et al.1995;condon et al.1998). Theobserving log for the new GMRT and VLA observations is given in Table2, including the name of the telescope and the array configuration, the frequency of observations and the primary beamwidth, as well as the dates of the exposures. The MHz VLA image of the entire source is presented in Figure1. The enlarged GMRT and VLA images of the northern and southern lobes at four different observing frequencies are shown in Figure 2. The images at 1400 MHz are reproduced from the FIRST survey, and the instrumental characteristics for all the radio maps included in this paper are collected in Table3. The obtained radio fluxes of J and the resulting luminosities are listed in Table 4. The low-frequency GMRT observations at MHz and MHz were made in the standard manner, with each observation of the target source interspersed with observations of calibrator sources. The phase calibrators PKS B (3C 298) and PKS B were observed at MHz and MHz, respectively, after each of several 20 minute-long exposures of the target centered on the core position. 3C286 was used as the flux density and bandpass calibrator based on the scale of Baars et al. (1977). The total observing times on the target sources were about 200 min (at 328.8MHz) and 180 min (at 617.3MHz). The data were edited and reduced using the NRAO AIPS package, and then self-calibrated to produce the best possible images. The high-frequency MHz observations conducted with the VLA in its D configuration were made in two runs aimed at a different lobe of J each (see Table2). The interferometric phases were calibrated every 20 min with the phase calibrator PKS B The source 3C286 was used as the primary flux density calibrator. 41 and 93 minute-long exposures were conducted for the fields centered on the northern and the southern lobes, respectively. This allowed us to reach an rms noise values of about and 0.036mJybeam 1 for the two exposures respectively. The data were edited, reduced and self-calibrated with the NRAO AIPS package. Images produced from the resulting datasets were corrected for the primary beam pattern Comparison Sample of DDRGs The three DDRGs selected for this project are J , J , and J The archival data for these used in our modeling are summarized in Table5.

7 7 The source J ( B ) is one of the best studied DDRG(Schoenmakers et al. 2000a; Kaiser et al. 2000; Konar et al. 2006). Multifrequency radio maps of this object presentedinkonar et al.(2006)showbroad,diffuseouterlobesextendedoutofthecollinearpair of narrow inner lobes. The dynamical modeling of the source discussed in Machalski et al. (2009) and in Brocksopp et al. (2011) resulted in significantly different values for the fitted model parameters. This is most likely due to different model setups and different assumptions involved in the modeling. Hence J deserves further consideration. The source J ( PKS B ) was analyzed by Saripalli et al. (2003), Safouris et al. (2008), and Machalski et al. (2010). The outer lobes of this DDRG extend back to the parent galaxy forming a common cocoon with no significant gaps in radio brightness. The inner narrow lobes are therefore fully embedded within the outer lobes, all the way from the radio core to their termination points. The ages and other characteristics of both pairs of lobes, as well as the properties of their surrounding environments, have been previously determined using the DYNAGE algorithm by Machalski et al. (2010). The DDRG J was discovered by Saikia et al. (2006). In spite of the lack of spectroscopic redshift determination for the parent galaxy, the source is taken into consideration in our work because its inner lobes are brighter at radio wavelengths than the outer lobes. This is quite unusual characteristic in the whole sample of all known DDRGs. The other morphological properties of the source are however typical for its class. An uncertainty in the photometric redshift estimate for J affects the determination of its physical size and luminosity, but should not influence much the main conclusions. 3. Modeling The standard and widely considered scenario for the time evolution of FR II type radio sources is based on the pioneering work by Scheuer (1974) and Blandford & Rees (1974), modified subsequently by a number of authors (e.g., Begelman & Cioffi 1989; Falle 1991; Kaiser & Alexander 1997; Kaiser et al. 1997; Blundell et al. 1999; Manolakou & Kirk 2002). This scenario enables for a self-consistent evaluation of the source (lobes ) length and radio luminosity for a given jet power, jet lifetime, and density of the ambient medium, accounting for the cumulative effects of energy losses of relativistic particles injected by the jet into the evolving cocoon due to the synchrotron radiation, inverse Compton scattering of CMB photons, and adiabatic expansion of the lobes. Our analysis presented below is based on a particular dynamical description by Kaiser & Alexander (1997) with the radiative processes properly taken into account following Kaiser, Dennett-Thorpe & Alexander (1997; hereafter KDA model), as summarized in Kaiser & Best (2007).

8 DYNAGE Algorithm Similarly to our previous work(machalski et al. 2008, 2009, 2010), the analysis presented in this paper is performed using the DYNAGE algorithm of Machalski et al. (2007). This numerical code allows us to solve the inverse problem, i.e. to determine for a given set of observables four main free parameters of the KDA model, namely (i) the jet power Q j, (ii) the density of the ambient gaseous medium near the radio core ρ 0, (iii) the age of the lobes t, and (iv) the effective injection spectral index α inj. The latter parameter approximates the slope of the emission continuum of the electron energy distribution freshly injected at the jet head jet termination shock (hotspot) into the expanding lobes. Note that the term effective stands here for the average over a broad energy range and over the lifetime of the source, since the energy spectrum of particles accelerated at the hotspot may change with time, and may be far from a single power-law approximation (see the discussion in Machalski et al. 2007). A relatively robust determination of these four free parameters is enabled by fitting the applied dynamical/radiative model to the main observables, namely (i) the projected linear length of the cocoon D, (ii) the total volume of the cocoon V, and (iii) the lobes radio power at different frequencies P νi (where i = 1,2,3,...). As in the KDA model, we assume an inclination θ = 90 and a cylindrical geometry for each source. This, for the total length D and width w, gives the total volume V = πd 3 /(4AR 2 ), where AR 2R T = D/w is the axial ratio of the whole radio structure (in the case of roughly symmetrical lobes). The other model parameters (such as adiabatic indices for the plasma inside and outside of the cocoon, ˆγ c and ˆγ IGM, respectively) and model assumptions are discussed in Machalski et al. (2007). The original KDA model assumes a constant and uninterrupted jet flow delivering ultrarelativistic gas to the expanding lobes throughout the whole lifetime of a source. As such, it is well suited for describing relatively young or normal-size FR II radio galaxies, but may not be directly applicable for modeling the evolved objects, i.e. remnant/outer lobes of GRGs and DDRGs, whose observed properties indicate in many cases a faded jet activity at the time of observations. Therefore, in order to analyze such systems in the framework of the discussed scenario, it seems necessary to introduce one more free parameter to the model, namely the jet lifetime t j, which may be substantially different than the age of the source t. Note in this context that when the jet activity is terminated, the adiabatic expansion of the cocoon inflated by the jets is expected to slow down, or even to stop completely (see Brocksopp et al. 2007, 2011). Hence two different phases of the dynamical evolution of such lobes should in fact be considered, the first one corresponding to the ongoing feeding of the lobes by the active jets (source ages t t j ), and the second one corresponding to the case when the jet flow is ceased (t > t j ). As discussed in Brocksopp et al. (2011), a few/several sound-crossing times through the lobes are needed to establish the new expansion regime

9 9 during the second phase of the evolution, and this timescale may be comparable to the age of the source in the case of large-scale outer lobes of GRGs and DDRGs. Following the above argument, in our fitting procedure we neglect the effect of a slower expansion phase when evaluating the dynamical parameters of a source. This allows us in particular to estimate the length of a jet as and the longitudinal expansion velocity of a cocoon as ( ) 1/(5 β) Qj l j (t) = c 1 t 3/(5 β), (1) ρ 0 a β 0 v h (t) = dl j dt = 3c ( ) 1/(5 β) 1 Qj t (β 2)/(5 β). (2) 5 β ρ 0 a β 0 In the above c 1 is a numerical constant depending on a value of R T, a 0 is the radius of the central plateau in the density distribution of the surrounding gas, and β is the exponent in the anticipated ambient density profile ρ(r) = ρ 0 (r/a 0 ) β. The total projected length of the entire structure is D = 2l j. The ratio of the pressure within the lobes heads to that within the rest of the cocoon, P hc p h /p c, is taken as P hc = ( β)R β T following Kaiser (2000). We note that only if t j is significantly shorter than t, thus evaluated model parameters would correspond to the lower limit on the ratio Q j /ρ 0. But since in our discussion below ( 4) we emphasize large jet kinetic powers and low ambient medium densities required to fit the observed spectrum of J (and of the other giant sources), equations 1 2 constitute a very conservative choice indeed. Besides, in the presented model fitting the jet lifetimes are never less than 70% of the source ages. Importantly, the evaluation of the lobes luminosities is done taking into account only those relativistic particles which are injected into the lobes during the first phase of the evolution, corresponding to the ongoing jet activity. The integration over the source spectrum is performed following the description provided in Kaiser & Cotter (2002) in the context of relict radio sources. In the KDA model the minimum-energy condition is applied, being determined by accounting for the energies stored in the lobes magnetic field and in relativistic radiating particles only. However, the extended lobes of luminous radio galaxies may not fulfill exactly the minimum-energy condition, and may contain especially at large stages of their evolution a significant fraction of non-radiating (possibly mildly-relativistic) particles, as indicated by the recent X-ray observations of a number of such systems as well as by some theoretical considerations. Yet the departures from the energy equipartition between magnetic field and radiating electrons within the FR II lobes are never that large(kataoka & Stawarz 2005;

10 10 Croston et al. 2005), as demonstrated also for several particular cases of GRGs (Konar et al. 2009; Isobe et al. 2009, 2011a,b). Hence, in our analysis we apply the standard minimumenergy condition, but allow for a non-negligible contribution of non-radiating particles (ultrarelativistic or only mildly-relativistic protons) to the total pressure inside the cocoons. The ratio of the energy density stored in the lobes magnetic field B to that stored in all the lobes particles is then R eq = (1+s inj )/[4(k +1)], where s inj = 2α inj +1 is the energy slope of the electrons (with minimum and maximum Lorentz factors γ min and γ min, respectively) freshly accelerated at the jet termination shock, and k is the energy/pressure fraction of the non-radiating particles Results The observational data summarized in 2 are fitted using the DYNAGE algorithm described in 3.1. The results of the modeling are given in Table6. In the fitting procedure we fix the following model parameters: adiabatic indices for the plasma within and outside of the cocoons ˆγ c = ˆγ IGM = 5/3, minimum electron Lorentz factor γ min = 1 and the maximum electron Lorentz factor γ max = 10 7, as well as the inclination of the lobes θ = 90 (see Machalski et al. 2007, for more discussion regarding these model assumption). In the anticipated density profile characterizing the ambient medium into which the lobes evolve, ρ(r) = ρ 0 (r/a 0 ) β, we fix the characteristic radius of the density plateau a 0 = 10kpc, and in the case of the outer cocoons propagating through an undisturbed IGM the power-law index β = 1.5. Such a density profile is expected to be a good approximation of the hot gaseous medium in groups of galaxies hosting luminous radio sources (Kaiser et al. 1997). When modeling the inner lobes, which propagate within a uniform plasma of the outer cocoons, we set instead β = 0 or 0.1 (as determined by a quality of the fits). Note that the assumption regarding a uniform distribution of the matter within the lobes is justified by the fact that the sound speed thereby has to be in general high when compared to the dynamical timescale involved (see Kaiser et al. 2000). Finally, we assumed the fraction of the non-radiating particles k = 0 for the inner lobes, and k = 10 for the outer lobes, in order to account for the expected entrainment processes enriching the evolved cocoons with the matter. The model fits return the following free parameters: the total jet kinetic power Q j, the density of the central plateau in the distribution of the ambient medium ρ 0, the injection spectral index α inj, the age of the lobes t, and the jet lifetime t j. These allow us to evaluate further some other parameters of the analyzed sources, including the advance velocity of the jets v h (t)/c (see equation2), the average expansion velocity of the source D/(2ct), the

11 11 density of the surrounding medium around the lobes heads ρ(d), the cocoon pressure p c, and the magnetic field intensity within the lobes B. The radio spectra for the control sample of DDRGs fitted by the applied model are presented in Figure3. The quality of the fits is quantified by the reduced χ 2 values given in Table6, and also in Figure The Case of J The high-resolution GMRT image of J at 617 MHz confirms very narrow lobes of the source, as already suggested by the 1400MHz map available from the FIRST survey. The axial ratio of the lobes, R T 12, is of the same order or even bit higher than that characterizing the inner lobes in basically all the well-studied DDRGs, and much higher than that established for normal-size FRIIs, outer lobes of DDRGs, or GRGs. Also the age of the observed radio structure of J found from the preliminary model fitting presented by Machalski et al. (2008), namely 50Myr, as well as the average expansion speed evaluated by the authors, 0.1c, are totallydiscordant with the corresponding values determined for the other GRGs. Therefore, the idea that the observed structure of J might be due to a restarted rather than a primary jet activity seems well justified. Unfortunately, all the observations made till now do not show a trace of any relict (older) cocoonembedding the observed lobes. In the first step of the modeling, we repeated the calculations of Machalski et al. (2008) using the new radio data presented in 2. The observed lobes of J are thus modeled as structures formed by a primary jet activity, with the standard ambient medium density profile appropriate for group of galaxies assumed (a 0 = 10kpc and β = 1.5), along with k = 1 (model Primary in Table6). It is worth noting however that the assumption regarding the ambient medium density may not be strictly valid in the particular case of the analyzed object, since the optical galaxy counts around J do not reveal any excess over the background level (see in this context 4.1). Besides, even if some gaseous halo exists around the host galaxy of the source, similar to the virialized hot IGM of elliptical galaxies and groups of galaxies, its extension up to >Mpc distances from the central AGN in a form of a steep power-law density distribution seems questionable. Thus, the parameters returned by the Primary model fit have to be taken with caution. In addition to this, we also model the lobes of J as structures inflated by restarting jets, in analogy to the inner lobes of DDRGs, setting in this case β = 0.1 and k = 0 (model Inner in Table6). For illustration, the Inner model fit to the comoving spectrum of the entire structure of J is shown in Figure 4. Both Primary and Inner models return comparable physical parameters of the stud-

12 12 iedobject inparticularq j ergs 1, α inj 0.55, t 35Myr, v h (t)/c D/(2ct) 0.2, p c ( ) dyncm 2, and B 0.6µG except of the central density of the ambient medium ρ 0. We note here relatively young age of the source implied by the model fits (in between of the ages evaluated for the inner and outer lobes of the analyzed DDRGs), its relatively high expansion velocity (ten times higher that the corresponding values found for the outer lobes in the control sample of DDRGs), and the magnetic field intensity exactly within the range established for other GRGs through the X-ray observations (Konar et al. 2009; Isobe et al. 2009, 2011a,b). The goodness of the fits in both Primary and Inner cases is practically the same; the test χ 2 red cannot be therefore used to distinguish between the models. Importantly, the only difference between the two models considered, i.e. the emerging value of the ρ 0 parameter, is just a direct consequence of the different ambient density profiles assumed (β = 1.5 versus β = 0.1); indeed, the density of the surrounding medium around the lobes heads, which is of a main importance in determining the other model parameter for the evolved/giant system like J , is basically identical in both scenarios, and is in addition very low, namely ρ(d) gcm 3. Hence the emerging conclusion is that regardless on the unknown density profile of the ambient medium surrounding the giant lobes of J and regardless on the unknown contribution from the non-radiating particles to the internal pressure of its cocoon (at least within the explored range) the lobes of J must evolve in a very rarified surrounding. We return to this point in the next section 4. If the low density medium surrounding J is due to the previous jet activity, i.e., if this source is indeed of a double-double type, the question is where are its outer lobes. As already emphasized above, such structures are not detected in the newly performed radio observations, but this may be only because of their advanced ages precluding any detection due to substantial cooling of the radiating electrons. In order to explore this idea in more detail, we evaluate the expected properties of the hypothetical outer lobes, within the framework of the general model for the evolution of luminous radio galaxies anticipated in this paper, by assuming the jet power in the previous activity epoch three times larger than the power emerging from the Inner modeling of J , namely Q j ergs 1. This is in analogy to the control sample of DDRGs, for which the previous jet activity epochs seem to be universally characterized by larger kinetic powers of the outflows when compared to the ongoing jet activity epochs. We also fix the ambient medium density profile β = 1.5, and set the central density as ρ 0 = gcm 3, which is the average determined for 30 GRGs studied by Machalski et al. (2009) and Machalski (2011). Finally, we set the aspect ratio for the hypothetical relict lobes of J as half of the value characterizing the observed radio structure, R T = 6. The expected time evolution of the length and luminosity of the outer lobes is then calculated in a number

13 13 of models differentiated by values of t j and t. Searching for a model providing relict lobes length not shorter than that of the currently observed structure, as well as radio fluxes below the current sensitivity limits, we find a possible solution for D = 4.70Mpc, t j = 106Myr and t out = 178Myr (see model Outer in Table6). Figure5 presents the expected observed spectrum of the hypothetical outer lobes corresponding to the model solution Evolutionary Tracks In the framework of the considered model and for a given set of the observables, the whole time evolution can be easily calculated for each analyzed source. Figure 6 shows such an evolution, depicted in the luminosity size evolutionary tracks, P ν D, for the outer and inner lobes of the three exemplary DDRGs modeled in this paper. Here the star symbols indicate the actual ages of the lobes emerging from the modeling, while the consecutive age markers are connected with dotted lines. The diagrams indicate that very soon after a termination of the jet flow the outer lobes decrease quickly in radio luminosity, and slow down their expansion. The inner lobes may then overtake the outer cocoon both in the luminosity and in size. Strictly speaking, at the point when the new-born structures grow beyond the older ones the model applied here to describe the evolution of the inner lobes is not valid anymore; hence the segments of the evolutionary curves shown in Figure 6 which corresponds to such a case do not describe any realistic situation. However, our calculations indicate also that this is expected to happen only when the inner lobes reach an age of about or over 100Myr. The analogous evolutionary tracks for J are shown in Figure 7. The tracks corresponding to the Primary and Inner models for the observed radio structure are denoted by solid and dashed curves, respectively. The evolutionary track for the hypothetical outer lobes in J , corresponding to the Outer model parameters, is given by short-dashed curve. The actual ages of the observed and of the hypothetical/outer radio structures emerging from the model fits are marked by the filled and open star symbols, respectively. As can be inferred from the presented diagrams 6 7, for the ages of the outer lobes t < t j corresponding to D(t) a 0 the lobes luminosities decrease slowly with time. This is because of the ongoing synchrotron and adiabatic cooling of ultrarelativistic electrons evolving within slowly expanding cocoons. We note that in the framework of the applied model the adiabatic losses do not depend on the ageof a source, due to the assumed constant ratio between the lobes pressures P hc p h /p c throughout the lifetime of a system. Also, inverse-compton energy losses of the radiating particles are mostly negligible: for all the

14 14 outer lobes modeled in this paper the energy densities stored in the lobes magnetic field exceed the energy densities of the CMB photon field at the source redshifts, u B B 2 /8π > u CMB (1 + z) 4 ergcm 3, up to late phases of the evolution t 100Myr. On the other hand, once the jet flow decays and the source grow old (ages t > t j ), combined synchrotron and inverse-compton energy losses quickly cool most of the particles radiating within the relict lobes, causing a rapid decrease in the source luminosity. The inner lobes evolve in a different manner. In the uniform, low-density but highpressure environment their adiabatic (sideway) expansion is hampered, and as a result the lobes luminosity increases with time. However, the simultaneously decreasing energy density of the lobes magnetic field falls at some point below the energy density of the CMB photon field, causing the onset of the luminosity decline due to the dominant inverse-compton energy losses of the radiating electrons. It is worth noting that such an effect, predicted already by Baldwin (1982), strongly depends on the density of the source environment. This is illustrated in Figure8, which shows the P ν D evolutionary tracks of a fiducial source (as considered by Kaiser et al. 1997) evolving within a uniform ambient medium (β=0) characterized by three different values of the density: ρ 0 = gcm 3 (solid line, reproducingfigure4inkaiser et al.1997), gcm 3 (dashedline), and gcm 3 (dotted line). As illustrated, the effect is stronger and stronger for lower and lower density of the ambient medium, but in all the cases the inverse-compton dominance and the related onset of the luminosity decline starts around the source ages of about t 100Myr, when the source sizes are D 0.1 1Mpc. The time evolution of length and luminosity of the outer and inner lobes of the exemplary DDRG J is shown in Figure9. 4. Discussion The main observational argument for a restarted jet activity in J is the very high axial ratio of the observed radio structure. This finding suggests almost ballistic propagation of powerful jets in a particularly low-density ambient medium. If this low density of the source environment is due to the previous jet activity indeed, or rather due to a unique location of J in a void region of the galaxy distribution, is a central question we aim to address in the presented analysis. Extremely large linear size of the studied object, D 4.7 Mpc, makes the investigation particularly relevant. Note that radio cocoons exceeding a few Mpc in size are quite unexpected at the redshift of 0.3, since all the other extended FRII systems at z 0.3, including outer lobes of DDRGs, are substantially shorter than this. Regardless on the details of the jet duty cycle and of the origin of the lobes in the system, however, the energetic requirements for the outflows in J are

15 15 severe, as discussed below Galaxy Counts As an independent check of the IGM density in the vicinity of J we have performed simple galaxy counts around the position of the source, and compared the resulting values with the mean background number density of galaxies determined by Tyson (1988), log N bg = (0.39 R 4.8) deg 2 mag 1, (3) for the isophotal red magnitudes R. We note that the isophotal limits are sufficiently low for most of the galaxies so that R can be considered as total magnitudes to the precision of the faint-galaxy photometry. However, in order to perform reliable counts of galaxies surrounding the analyzed target, it is necessary to choose proper selection criteria in order to remove from the list stellar objects, keeping in mind that the host galaxy of J itself is very faint: R = 19.65mag and B = 21.82mag (Machalski et al. 2008). Luminous extragalactic radio sources are hosted by elliptical galaxies located in the centers of poor clusters and groups. A standard measure of cluster richness can therefore be used to check whether J shares this characteristic or nor. Here we follow a particularly simple and frequently applied procedure of counting the objects within 3 Mpc radius around a given R magnitude galaxy, which are brighter than R +3 (e.g., Hill & Lilly 1991). For this purpose we use the Minnesota Automated Plate Scanner (MAPS) Catalog of the first-epoch Palomar Observatory Sky Survey (POSS-I), which was derived from digitized scans of glass copies of the blue and red plates of the original POSS-I. The limiting magnitudes of this database preclude however detections of the systems as faint as R = = 22.65mag. For this reason, we define N as a number of galaxies with red POSS-I magnitudes E ± 1.5, which are located within 3 Mpc radius around J (corresponding to the angular separation of 11.1 ). Two different selection criteria have been considered: (i) 0.5 galnod 1.0 at both E and O plates indexing non-stellar objects, and (ii) colors O E > 1.4mag providing rejection of bluer objects, i.e., by assumption, quasars. As a result, 68 and 72 objects fulfilling the above criteria have been found, respectively. The emerging numbers are, quite surprisingly, fully compatible with each other. Diving the corresponding mean value 70 by a circular sky area of 3Mpc radius at the position of the target source gives us N 650 galaxies per deg 2. The relevant number of background galaxies implied by equation3 is N bg Hence, the resulting ratio N/N bg 0.25 indicates that J is indeed located in distinctly poor region of the galaxy distribution, as further confirmed by investigating the number density of galaxies within a few larger concentric rings, up to 16Mpc from the target; this number density occurs to be almost constant.

16 Relict Lobes? The P ν D diagrams shown in Figure6 illustrate that after a termination of the jet activity, luminosities of the outer lobes in DDRGs decrease quickly with time, falling eventually below the inner lobes luminosities. This is further depicted in Figure 9 for the particular case of DDRG J Therefore, if the hypothetical outer lobes of J are old enough, one may expect their radio surface brightness to fall substantially below the rms noise level of the available radio maps. Figure7 shows that the 1400MHz luminosity of these lobes, calculated for the Outer model parameters (Table 6), is indeed lower than that of the currently observed structure in the system. On the other hand, one should expect at the same time a substantial steepening of the radio continuum for the dying cocoon, and hence possibly non-negligible fluxes at low observing frequencies. Figure 5 indicates in particular that the total flux density of the presumed outer lobes in J might be comparable to the observed flux density at frequencies 300MHz. Around 74MHz, the frequency of the VLSS survey (Cohen et al. 2007), the estimated flux is about 1.4Jy. Even though relatively high, this flux is still not enough for a robust detection of the relicts in the archival dataset. We note that the VLSS survey is characterized by the average sensitivity of σ 0.1Jy per beam of 80 80, and hence that the detection limit of about 5σ for an unresolved source reads as 0.5Jybeam 1. Assuming therefore that each of the two putative outer lobes of J has an area of about 2 3, i.e. about 3 4 beams, a 74MHz flux of about Jy would be required for the detection in the VLSS. This is 2 3 times the 0.7Jy flux expected for each of the lobes in the framework of the Outer model. Future high-sensitivity VLA or Low Frequency Array (LOFAR) observations at frequencies below 100 MHz could be however more promising for investigating a presence of the hypothetical outer cocoon in J Interestingly, this structure could be also possibly probed at X-ray frequencies, since inverse-compton scattering of the CMB photons to the X-ray range involves electron energies lower than those involved in producing observable radio synchrotron emission of the lobes. Searching for analogous relict ( ghost ) structures in other evolved radio galaxies and radio-loud quasars using current X-ray instruments has been in fact already proposed and considered in the literature (e.g., Fabian et al. 2009; Blundell & Fabian 2011; Mocz et al. 2011) Age and Expansion Velocity In the previous section 3.3 we have noted a relatively young age of the observed radio structure of J emerging from the model fits, t 35Myr, which is in between of the ages evaluated for the inner and outer lobes of the analyzed DDRGs. Correspondingly,

17 17 a relatively high expansion velocity v h (t)/c D/(2ct) 0.2 has been found for the source, about ten times higher that the analogous values estimated for the outer lobes in the control sample of DDRGs, and a factor of at least a few higher than the expansion velocity of the inner lobes in such systems. Are however these somehow extreme values realistic at all, or are they just artifacts of the incorrect scenario applied? In all the kinematic models for FR II-type sources, characteristic asymmetries in lengths and luminosities between the two opposite lobes are expected due to the light time travel effects (e.g., Longair & Riley 1979; Best et al. 1995; Scheuer 1995; Ryś 2000). The size ratio between the longer (jet side) and shorter (counter-jet side) lobes is in particular q = [1+(v h /c)cosθ]/[1 (v h /c)cosθ], giving the advance speed of the jets ( ) q 1 ( c ) v h =. (4) q +1 cosθ In the case of J we have estimated q 1.078, leading to v h c/cosθ and thus implying θ < 88. If we further assume lobes inclination in the source θ 80, the expansion velocity reads as v h 0.216c, and hence the source age as t D/(2v h ) 35.4Myr. Thus, the kinematic estimates of the age and of the advance velocity of the observed radio structure in J are fully consistent with the values resulting from the dynamical modeling presented in 3.3. We note that the jet inclination θ 80 is in fact more probable than θ 90 adopted in the fitting procedure also for the other DDRGs analyzed. However, the resulting model parameters are not sensitive to such relatively minor changes in the assumed value of θ, since a physical size of the lobes increases only by 2% when the source inclination is decreased from 90 to 80. Finally, we comment on the fact that even higher advance velocity of the jet head has been estimated for the radio-loud quasar 3C 273 by Stawarz (2004), namely v h 0.85c, under the hypothesis that the observed large-scale structure in this source is also a manifestation of the restarted jet activity. The model proposed by Stawarz (2004), which was inspired similarly by somewhat unusual morphological properties of 3C 273 namely by the extremely narrow and only one-sided large-scale lobe in the system is in fact analogous to the one discussed here for J Density of the Ambient Medium An upper limit for the thermal gas density within cocoons of FRII sources may be estimated from the internal depolarization (lack of any signatures for such) of the observed lobes radio emission. In particular, Garrington & Conway (1991) found that the product of

18 18 the cold gas number density and the magnetic field intensity within the studied lobes is on average n th B < µgcm 3. This, for the typical range 3µG< B < 30µG found by means of multiwavelength analysis of the non-thermal lobes emission (including the radio and X-ray data for normal-size FRII systems; e.g., Kataoka & Stawarz 2005; Croston et al. 2005, and references therein), gives very roughly n th < (1 10) 10 4 cm 3, or equivalently ρ th m p n th < (0.3 3) gcm 3. Below we assume that same limit holds also for the cocoons of GRGs, and the outer lobes of DDRGs. Ontheother hand, theupper limit ontheplasma density within theextended lobesmay be provided by analyzing the total amount of ultrarelativistic (or at least mildly relativistic) plasma deposited by the jets into the expanding cocoons. Neglecting the work done by such expanding cocoons on the ambient medium, the total energy density deposited by a pair of jets with kinetic energy Q j over the source lifetime t into the cavity with volume V = πd 3 /(16R 2 T ) is simply u rel 2Q j t/v. This energy density is shared between the lobes magnetic field, electron-positron pairs, and also possibly protons, u tot = u B +u e + u p. For the purpose of rough estimates, below we assume an energy equipartition between different plasma species within the evolved lobes, as well as mixed(hadronic-leptonic) jet composition, leading to u rel 3u p 3 γ p m p c 2 n rel, where γ p is a mean Lorentz factor of the protons within the lobes, and n rel is their number density. Noting that in the case of an electronproton jet propagating with relativistic bulk velocities, for which a significant/dominant fraction of the kinetic energy is carried by the particles (as in fact suggested by the most recent analysis for luminous quasar sources; e.g., Sikora & Madejski 2000; Kataoka et al. 2008), protons injected into the lobes through the jet termination shock are expected to be characterized by mean energies γ p Γ j, where Γ j is the bulk Loretz factor of the outflow (see, e.g., the discussion in Kino & Takahara 2004; Stawarz et al. 2007). Hence the minimum number density of the lobes reads as n rel 2Q j t 3V Γ j m p c 2. (5) For the parameters of the hypothetical relict lobes of J implied by the Outer model (Table6), namely Q j = ergs 1, t = 178Myr, D = 4.7Mpc and R T = 6 (giving the total volume of the cocoon V cm 3 ), assuming in addition conservatively large-scale jet bulk Lorentz factor Γ j 3 (see Wardle & Aaron 1997), we obtain n rel cm 3, or equivalently ρ rel m p n rel gcm 3. This value is in a very good agreement with the analogous one estimated by Safouris et al. (2008) for the outer lobes of DDRG PKS B (J in this paper). To sum up, the dynamical modeling presented in 3 implies that the density of the ambient medium around the observed lobes of J is ρ J1420 inn gcm 3.