Large-scale Jet Research Articles

Star and jet multiplicity in the high-mass star forming region IRAS 05137+3919

Context. We present a study of the complex high-mass star forming region IRAS 05137+3919 (also known as Mol8), where multiple jets and a rich stellar cluster have been described in previous works. Aims. Our goal is to determine the number of jets and shed light on their origin, and thus determine the nature of the young stars powering these jets. We also wish to analyse the stellar clusters by resolving the brightest group of stars. Methods. The star forming region was observed in various tracers and the results were complemented with ancillary archival data. The new data represent a substantial improvement over previous studies both in resolution and frequency coverage. In particular, adaptive optics provides us with an angular resolution of 80 mas in the near IR, while new mid- and far-IR data allow us to sample the peak of the spectral energy distribution and thus reliably estimate the bolometric luminosity. Results. Thanks to the near-IR continuum and millimetre line data we can determine the structure and velocity field of the bipolar jets and outflows in this star forming region. We also find that the stars are grouped into three clusters and the jets originate in the richest of these, whose luminosity is ~ 2.4 × 104 L ⊙ . Interestingly, our high-resolution near-IR images allow us to resolve one of the two brightest stars (A and B) of the cluster into a double source (A1+A2). Conclusions. We confirm that there are two jets and establish that they are powered by B-type stars belonging to cluster C1. On this basis and on morphological and kinematical arguments, we conclude that the less extended jet is almost perpendicular to the line of sight and that it originates in the brightest star of the cluster, while the more extended one appears to be associated with the more extincted, double source A1+A2. We propose that this is not a binary system, but a small bipolar reflection nebula at the root of the large-scale jet, outlining a still undetected circumstellar disk. The gas kinematics on a scale of ~0.2 pc seems to support our hypothesis, because it appears to trace rotation about the axis of the associated jet.

Real-time evolution of a large-scale relativistic jet

Context. Astrophysical jets are ubiquitous in the Universe on all scales, but their large-scale dynamics and evolution in time are hard to observe since they usually develop at a very slow pace. Aims. We aim to obtain the first observational proof of the expected large-scale evolution and interaction with the environment in an astrophysical jet. Only jets from microquasars offer a chance to witness the real-time, full-jet evolution within a human lifetime, since they combine a 8 short ′ , few parsec length with relativistic velocities. Methods. The methodology of this work is based on a systematic recalibraton of interferometric radio observations of microquasars available in public archives. In particular, radio observations of the microquasar GRS 1758−258 over less than two decades have provided the most striking results. Results. Significant morphological variations in the extended jet structure of GRS 1758−258 are reported here that were previously missed. Its northern radio lobe underwent a major morphological variation that rendered the hotspot undetectable in 2001 and reappeared again in the following years. The reported changes confirm the Galactic nature of the source. We tentatively interpret them in terms of the growth of instabilities in the jet flow. There is also evidence of surrounding cocoon. These results can provide a testbed for models accounting for the evolution of jets and their interaction with the environment.

THE COMPLEX NORTH TRANSITION REGION OF CENTAURUS A: RADIO STRUCTURE

We present deep radio images of the inner 50 kpc of Centaurus A, taken with the Karl G. Jansky Very Large Array (VLA) at 90cm. We focus on the Transition Regions between the inner galaxy - including the active nucleus, inner radio lobes, and star-forming disk - and the outer radio lobes. We detect previously unknown extended emission around the Inner Lobes, including radio emission from the star-forming disk. We find that the radio-loud part of the North Transition Region, known as the North Middle Lobe, is significantly overpressured relative to the surrounding ISM. We see no evidence for a collimated flow from the Active Galactic Nucleus (AGN) through this region. Our images show that the structure identified by Morganti et al. (1999) as a possible large-scale jet appears to be part of a narrow ridge of emission within the broader, diffuse, radio-loud region. This knotty radio ridge is coincident with other striking phenomena: compact X-ray knots, ionized gas filaments, and streams of young stars. Several short-lived phenomena in the North Transition Region, as well as the frequent re-energization required by the Outer Lobes, suggest that energy must be flowing through both Transition Regioins at the present epoch. We suggest that the energy flow is in the form of a galactic wind.

The jet-ISM interaction in the outer filament of Centaurus A

The interaction between the radio plasma ejected by the active nucleus of a galaxy and the surrounding medium is a key process that can have a strong impact on the interstellar medium of the galaxy and hence on galaxy evolution. The closest laboratory where we can observe and investigate this phenomenon is the radio galaxy Centaurus A. About 15 kpc north-east of this galaxy, a particularly complex region is found: the so-called Outer Filament where jet-cloud interactions have been proposed to occur. We investigate the presence of signatures of jet-ISM interaction by a detailed study of the kinematics of the ionized gas, expanding on previous results obtained from the HI. We observed two regions of the outer filament with VLT/VIMOS in the IFU observing mode. Emission from Hbeta and [OIII]4959,5007\AA\ is detected in both pointings. We found two distinct kinematical components of ionized gas that well match the kinematics of the nearby HI cloud. One component follows the regular kinematics of the rotating gas while the second shows similar velocities to those of the nearby HI component thought to be disturbed by an interaction with the radio jet. We suggest that the ionized and atomic gas are part of the same dynamical gas structure originating as result of the merger that shaped Centaurus A and which is regularly rotating around Centaurus A as proposed by other authors. The gas (ionized and HI) with anomalous velocities is tracing the interaction of the Large-Scale radio Jet with the ISM, suggesting that, although poorly collimated as structure, the jet is still active. However, we can exclude that a strong shock is driving the ionization of the gas. It is likely that a combination of jet entrainment and photoionization by the UV continuum from the central engine is needed in order to explain both the ionization and the kinematics of the gas in the Outer Filament.

MODELING MHD ACCRETION-EJECTION: EPISODIC EJECTIONS OF JETS TRIGGERED BY A MEAN-FIELD DISK DYNAMO

We present MHD simulations exploring the launching, acceleration and collimation of jets and disk winds. The evolution of the disk structure is consistently taken into account. Extending our earlier studies, we now consider the self-generation of the magnetic field by an $\alpha^2\Omega$ mean-field dynamo. The disk magnetization remains on a rather low level, that helps to evolve the simulations for $T > 10,000$ dynamical time steps on a domain extending 1500 inner disk radii. We find a magnetic field of the inner disk similar to the commonly found open field structure, favoring magneto-centrifugal launching. The outer disk field is highly inclined and predominantly radial. Here, differential rotation induces a strong toroidal component that plays a key role in outflow launching. These outflows from the outer disk are slower, denser, and less collimated. If the dynamo action is not quenched, magnetic flux is continuously generated, diffuses outward the disk, and fills the entire disk. We have invented a toy model triggering a time-dependent mean-field dynamo. The duty cycles of this dynamo lead to episodic ejections on similar timescales. When the dynamo is suppressed as the magnetization falls below a critical value, the generation of the outflows and also accretion is inhibited. The general result is that we can steer episodic ejection and {\em large-scale jet knots} by a {\em disk-intrinsic dynamo} that is time-dependent and regenerates the jet-launching magnetic field.

Large-scale jets from active galactic nuclei as a source of intracluster medium heating: cavities and shocks

The evolution of powerful extragalactic jets is not only interesting by itself, but also for its impact on the evolution of the host galaxy and its surroundings. We have performed long-term axisymmetric numerical simulations of relativistic jets with different powers to study their evolution through an environment with a pressure and density gradient. Our results show key differences in the evolution of jets with different powers in terms of the spatial and temporal scales of energy deposition. According to our results, the observed morphology in X-ray cavities requires that an important fraction of the jet's energetic budget is in the form of internal energy. Thus, light, lepton-dominated jets are favoured. In all cases, heating is mainly produced by shocks. Cavity overpressure is sustained by an important population of thermal particles. Our simulations reproduce the cool-core structure in projected, luminosity-weighted temperature. We have performed an additional simulation of a slow, massive jet and discuss the differences with its relativistic counterparts. Important qualitative and quantitative differences are found between the non-relativistic and the relativistic jets. Our conclusions point towards a dual-mode of AGN kinetic feedback, depending on the jet power.

NuSTARREVEALS THE COMPTONIZING CORONA OF THE BROAD-LINE RADIO GALAXY 3C 382

Broad-line radio galaxies (BLRGs) are active galactic nuclei that produce powerful, large-scale radio jets, but appear as Seyfert 1 galaxies in their optical spectra. In the X-ray band, BLRGs also appear like Seyfert galaxies, but with flatter spectra and weaker reflection features. One explanation for these properties is that the X-ray continuum is diluted by emission from the jet. Here, we present two NuSTAR observations of the BLRG 3C 382 that show clear evidence that the continuum of this source is dominated by thermal Comptonization, as in Seyfert 1 galaxies. The two observations were separated by over a year and found 3C 382 in different states separated by a factor of 1.7 in flux. The lower flux spectrum has a photon-index of $\Gamma=1.68^{+0.03}_{-0.02}$, while the photon-index of the higher flux spectrum is $\Gamma=1.78^{+0.02}_{-0.03}$. Thermal and anisotropic Comptonization models provide an excellent fit to both spectra and show that the coronal plasma cooled from $kT_e=330\pm 30$ keV in the low flux data to $231^{+50}_{-88}$ keV in the high flux observation. This cooling behavior is typical of Comptonizing corona in Seyfert galaxies and is distinct from the variations observed in jet-dominated sources. In the high flux observation, simultaneous Swift data are leveraged to obtain a broadband spectral energy distribution and indicates that the corona intercepts $\sim 10$% of the optical and ultraviolet emitting accretion disk. 3C 382 exhibits very weak reflection features, with no detectable relativistic Fe K$\alpha$ line, that may be best explained by an outflowing corona combined with an ionized inner accretion disk.

The North Middle Lobe of Centaurus A

AbstractWe present new, deep VLA 327 MHz, GALEX Far-UV, and Hα images of the inner ~50 kpc of Centaurus A. We find the structure identified by Morganti et al. 1999 as a possible “large scale jet” is part of a knotty, linear feature within a broader region of diffuse radio emission. The linear feature is coincident with a narrow ribbon of Far-UV and Hα emission that extends 6-35 kpc from the galaxy core, as well with a similar ridge of soft X-ray emission. The Far-UV image also shows that a strong starburst is occurring in the central dusty disk, with a star-formation rate of ~ 2M⊙ yr−1. We suggest that the various peculiar phenomena seen to the NE of the galaxy can be explained by a wind from the starburst disk, enhanced by energy input from the AGN.

Observational aspects of AGN jets at high energy

AbstractFor the last two decades, significant and dramatic progress has been made in understanding astrophysical jet sources, particularly in the X-ray and gamma-ray energy bands. For example, the Chandra X-ray observatory reveals a number of AGN jets extending from kpc to Mpc scales. More recently, the Fermi Gamma-ray Space Telescopes launched in 2008 started monitoring the gamma-ray sky with excellent sensitivity of about ten times greater than that of EGRET onboard CGRO, and has detected more than 2,000 sources (mostly AGNs) as of 2014. Moreover, Fermi-LAT has discovered gamma-ray emissions not only from blazars but from a dozen radio galaxies not previously known to emit gamma-rays. Closer to home, the Fermi-bubbles were discovered to extend 50 degrees above and below the Galactic center. These large scale diffuse gamma-ray structures are similar in structure to AGN lobes such as those seen in Cen A and provide evidence for past activity in our Galactic center. In this review, I will first summarize recent highlights of large scale jets in radio galaxies, specifically resolved by the Chandra X-ray observatory. Next I will move on to the gamma-ray sky to present some highlights from Fermi-LAT observations of “misaligned” blazars, namely radio galaxies. I will discuss a unification scheme connecting blazars and misaligned radio galaxies. In the last part, I will also briefly comment on recent multiband observations of the Fermi-bubble and possible impacts on the AGN jet physics in the near future.

Possible proton synchrotron origin of X-ray and gamma-ray emission in large-scale jet of 3C 273

Abstract The large-scale jet of quasar 3C 273 has been observed in radio to gamma-ray frequencies. Earlier the X-ray emission from knot A of this jet has been explained with inverse Compton scattering of the cosmic microwave background radiations by the shock accelerated relativistic electrons in the jet. More recently it has been shown that this mechanism overproduces the gamma-ray flux at GeV energy and violates the observational results from Fermi LAT. We have considered the synchrotron emission from a broken power-law spectrum of accelerated protons in the jet to explain the observed X-ray to gamma-ray flux from knot A. The two scenarios discussed in our work are (i) magnetic field is high, synchrotron energy loss time of the protons is shorter than their escape time from the knot region and the age of the jet and (ii) their escape time is shorter than their synchrotron energy loss time and the age of the jet. These scenarios can explain the observed photon spectrum well for moderate values of Doppler factor. The required jet luminosity is high ∼1046 erg s−1 in the first scenario and moderate ∼1045 erg s−1 in the second, which makes the second scenario more favourable.

On the spine-layer scenario for the very high-energy emission of NGC 1275

We discuss the $\gamma$-ray emission of the radiogalaxy NGC 1275 (the central galaxy of the Perseus Cluster), detected by Fermi-LAT and MAGIC, in the framework of the "spine-layer" scenario, in which the jet is assumed to be characterized by a velocity structure, with a fast spine surrounded by a slower layer. The existence of such a structure in the parsec scale jet of NGC 1275 has been recently proved through VLBI observations. We discuss the constraints that the observed spectral energy distribution imposes to the parameters and we present three alternative models, corresponding to three different choices of the angles between the jet and the line of sight ($\theta_{\rm v}=6^{\circ}, 18^{\circ}$ and 25$^{\circ}$). While for the the case with $\theta_{\rm v}=6^{\circ}$ we obtain an excellent fit, we consider this solution unlikely, since such small angles seems to be excluded by radio observations of the large-scale jet. For $\theta_{\rm v}=25^{\circ}$ the required large intrinsic luminosity of the soft (IR--optical) component of the spine determines a large optical depth for $\gamma$-rays through the pair production scattering $\gamma \gamma\rightarrow e^+ e^-$, implying a narrow cut--off at $\sim50$~GeV. We conclude that intermediate angles are required. In this case the low frequency and the high--energy emissions are produced by two separate regions and, in principle, a full variety of correlations is expected. The correlation observed between the optical and the $\gamma$-ray flux, close to linearity, is likely linked to variations of the emissivity of the spine.

A close-pair binary in a distant triple supermassive black hole system.

Galaxies are believed to evolve through merging, which should lead to some hosting multiple supermassive black holes. There are four known triple black hole systems, with the closest black hole pair being 2.4 kiloparsecs apart (the third component in this system is at 3 kiloparsecs), which is far from the gravitational sphere of influence (about 100 parsecs for a black hole with mass one billion times that of the Sun). Previous searches for compact black hole systems concluded that they were rare, with the tightest binary system having a separation of 7parsecs (ref. 10). Here we report observations of a triple black hole system at redshift z = 0.39, with the closest pair separated by about 140 parsecs and significantly more distant from Earth than any other known binary of comparable orbital separation. The effect of the tight pair is to introduce a rotationally symmetric helical modulation on the structure of the large-scale radio jets, which provides a useful way to search for other tight pairs without needing extremely high resolution observations. As we found this tight pair after searching only six galaxies, we conclude that tight pairs are more common than hitherto believed, which is an important observational constraint for low-frequency gravitational wave experiments.

Three-dimensional general relativistic radiation magnetohydrodynamical simulation of super-Eddington accretion, using a new code harmrad with M1 closure

Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant application using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20\%$, the large-scale electromagnetic jet efficiency is of order $10\%$, and the total radiative efficiency that reaches large distances remains low at only order $1\%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs.\ disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $z\sim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.

EVIDENCE OF PARSEC-SCALE JETS IN LOW-LUMINOSITY ACTIVE GALACTIC NUCLEI

The nuclear radio emission of low-luminosity active galactic nuclei (LLAGN) is often associated with unresolved cores. In this paper we show that most LLAGN present extended jet radio emission when observed with sufficient angular resolution and sensitivity. They are thus able to power, at least, parsec-scale radio jets. To increase the detection rate of jets in LLAGN, we analyze subarcsecond resolution data of three low-ionization nuclear emission regions. This yields the detection of extended jet-like radio structures in NGC\,1097 and NGC\,2911 and the first resolved parsec-scale jet of NGC\,4594 (Sombrero). The three sources belong to a sample of nearby LLAGN for which high-spatial-resolution spectral energy distribution of their core emission is available. This allows us to study their accretion rate and jet power ($Q_\mathrm{jet}$) without drawing on (most) of the ad hoc assumptions usually considered in large statistical surveys. We find that those LLAGN with large-scale radio jets ($>$100 pc) have $Q_\mathrm{jet} > 10^{42}$ erg s$^{-1}$, while the lowest $Q_\mathrm{jet}$ correspond to those LLAGN with parsec-scale ($\leq$100 pc) jets. The $Q_\mathrm{jet}$ is at least as large as the radiated bolometric luminosity for all LLAGN, in agreement with previous statistical studies. Our detection of parsec-scale jets in individual objects further shows that the kinematic jet contribution is equally important in large- or parsec-scale objects. We also find that the Eddington-scaled accretion rate is still highly sub-Eddington ($<10^{-4}$) when adding the $Q_\mathrm{jet}$ to the total emitted luminosity (radiated plus kinetic). This indicates that LLAGN are not only inefficient radiators but that they also accrete inefficiently or are very efficient advectors.

Spatiotemporal behaviors of absolute density of atomic oxygen in a planar type of Ar/O2 non-equilibrium atmospheric-pressure plasma jet

A method combining of two-photon absorption laser-induced fluorescence and vacuum ultra-violet absorption spectroscopy is constructed and used to get the fine spatial distribution of O generated by a planar type of Ar/O2 non-equilibrium atmospheric-pressure plasma jet. The O density in the quasi-uniform region, which occupies about 80% of the total plasma width, is as high as 1015 cm−3. The lifetime of O is estimated to be about 360 ± 60 µs. The loss of O atoms is due to the three-body reaction of O + O2 + Ar → O3 + Ar confirmed by a simple calculation with a few rate equations. The results are very useful for the simulation of Ar/O2 plasmas, the design of large-scale planar plasma jets and the development of their potential applications.

Axis-switching of a micro-jet

In this study, it is shown that free microjets can undergo complex transitions similar to large-scale free jets despite relatively low Reynolds numbers. Using an original experimental method allowing for the 3D reconstruction of the instantaneous spatial organization of the microjet, the axis-switching of a micro-jet is observed for the first time. This is the first experimental evidence of such complex phenomena for free micro-jets. Combining these experimental results with Direct Numerical Simulations it is shown that the mechanism responsible for the axis-switching is the deformation of a micro-vortex ring due to induction by the corner vortices, as it occurs in large scale non-circular jets.

Lagrangian reconstructions of temperature and velocity in a model of surface ocean turbulence

The characterization of submesoscale dynamics is crucial to apprehend their impact on the global ocean properties. Direct measurements of fine structures over the world oceans, nevertheless, are at present severely limited by the spatial resolution of available satellite products. In this work we numerically investigate the possibility to reconstruct tracer fields, like surface temperature, at small scales, from low-resolution data using a Lagrangian technique based on the properties of chaotic advection. The capabilities of the method are explored in the context of a forced Surface Quasi Geostrophic (SQG) turbulent flow representing a large-scale meandering jet and smaller-scale eddies. Both qualitative and quantitative comparisons are performed between the original (high-resolution) fields and their reconstructions that use only low-resolution data. Good agreement is found for filamentary structures, even in the presence of a large-scale forcing on the tracer dynamics. The statistics of tracer gradients, which are relevant for assessing the possibility to detect fronts, are found to be accurately reproduced. Exploiting SQG theory, the reconstruction technique is also extended to obtain the velocity field in three dimensions when temperature is the tracer. The results indicate that relevant features of dynamical quantities at small scales may be adequately deduced from only low-resolution temperature data. However, the ability to reconstruct the flow is critically limited by the energetic level of submesoscales. Indeed, only structures generated by non-local mesoscale features can be well retrieved, while those associated to the local dynamics of submesoscale eddies cannot be recovered.

FERMI RULES OUT THE INVERSE COMPTON/CMB MODEL FOR THE LARGE-SCALE JET X-RAY EMISSION OF 3C 273

The X-ray emission mechanism in large-scale jets of powerful radio quasars has been a source of debate in recent years, with two competing interpretations: either the X-rays are of synchrotron origin, arising from a different electron energy distribution than that producing the radio to optical synchrotron component, or they are due to inverse Compton scattering of cosmic microwave background photons (IC/CMB) by relativistic electrons in a powerful relativistic jet with bulk Lorentz factor Gamma ~ 10-20. These two models imply radically different conditions in the large-scale jet in terms of jet speed, kinetic power, and maximum energy of the particle acceleration mechanism, with important implications for the impact of the jet on the large-scale environment. A large part of the X-ray origin debate has centered on the well-studied source 3C 273. Here we present new observations from Fermi which put an upper limit on the gamma-ray flux from the large-scale jet of 3C 273 that violates at a confidence greater that 99.9% the flux expected from the IC/CMB X-ray model found by extrapolation of the UV to X-ray spectrum of knot A, thus ruling out the IC/CMB interpretation entirely for this source when combined with previous work. Further, this upper limit from Fermi puts a limit on the Doppler beaming factor of at least delta < 9, assuming equipartition fields, and possibly as low as delta < 5, assuming no major deceleration of the jet from knots A through D1.

Aeroacoustics of volcanic jets: An overview

Explosive volcanic eruptions can inject large volumes of ash into heavily traveled air corridors; they pose a significant societal and economic hazard. They also generate large amplitude atmospheric infrasound waves (~0.01-20 Hz), which can be recorded at thousands of kilometers from the eruption and can provide detailed information on the timing, duration, and relative vigor of the volcanic explosions. In order to provide more detail about the eruption process based on acoustic signals, a quantitative model for the acoustic source process within the volcanic eruption column is needed. Volcanic eruption columns are modeled by a momentum-driven jet flow, transitioning with altitude into a thermally buoyant plume. Infrasound recordings from such activity resemble the large-scale turbulence similarity spectrum, indicating that large-scale volcanic jet flows generate an infrasonic form of jet noise. However, volcanic jet noise deviates from pure-air laboratory jet noise because of complexities such as multiphase flow (especially loading with ash particles); nozzle/crater geometry and roughness; buoyancy effects; and high temperature and density effects. We propose a new framework for understanding acoustic sources at volcanoes based on aeroacoustics research, which is being developed through multi-disciplinary integration of field, numerical, and laboratory studies.

Hints of an axion-like particle mixing in the GeV gamma-ray blazar data?

Axion-Like Particles (ALPs), if exist in nature, are expected to mix with photons in the presence of an external magnetic field. The energy range of photons which undergo strong mixing with ALPs depends on the ALP mass, on its coupling with photons as well as on the external magnetic field and particle density configurations. Recent observations of blazars by the Fermi Gamma-Ray Space Telescope in the 0.1–300 GeV energy range show a break in their spectra in the 1–10 GeV range. We have modeled this spectral feature for the flat-spectrum radio quasar 3C454.3 during its November 2010 outburst, assuming that a significant fraction of the gamma rays convert to ALPs in the large scale jet of this blazar. Using theoretically motivated models for the magnetic field and particle density configurations in the kiloparsec scale jet, outside the broad-line region, we find an ALP mass ma ∼ (1−3)⋅10−7 eV and coupling gaγ ∼ (1−3)⋅10−10 GeV−1 after performing an illustrative statistical analysis of spectral data in four different epochs of emission. The precise values of ma and gaγ depend weakly on the assumed particle density configuration and are consistent with the current experimental bounds on these quantities. We apply this method and ALP parameters found from fitting 3C454.3 data to another flat-spectrum radio quasar PKS1222+216 (4C+21.35) data up to 400 GeV, as a consistency check, and found good fit. We find that the ALP-photon mixing effect on the GeV spectra may not be washed out for any reasonable estimate of the magnetic field in the intergalactic media.