Issue |
A&A
Volume 599, March 2017
|
|
---|---|---|
Article Number | A55 | |
Number of page(s) | 7 | |
Section | Astrophysical processes | |
DOI | https://doi.org/10.1051/0004-6361/201628810 | |
Published online | 28 February 2017 |
Time lag in transient cosmic accreting sources
1 Physics Department, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185 Rome, Italy
2 National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Kashirskoe Shosse 31, 115409 Moscow, Russia
3 Space Research Institute of Russian Academy of Sciences, Profsoyuznaya 84/32, 117997 Moscow, Russia
e-mail: gkogan@iki.rssi.ru
4 INAF–Istituto di Astrofisica e Planetologia Spaziali, via del Fosso del Cavaliere, 100, 00133 Roma, Italy
e-mail: franco.giovannelli@iaps.inaf.it
Received: 28 April 2016
Accepted: 18 November 2016
Context. We develop models for time lag between the maxima of the source brightness in different wavelengths during a transient flash of luminosity that is connected with a short-period increase of the mass flux onto the central compact object.
Aims. We derive a simple formula for finding the time delay among events in different wavelengths which is valid in general for all disk-accreting cosmic sources. We quantitatively also discuss a model for time-lag formation in active galactic nuclei (AGNs).
Methods. In close binaries with accretion disks, the time lag is connected with effects of viscosity that define a radial motion of matter in the accretion disk. In AGN flashes, the infalling matter has a low angular momentum, and the time lag is defined by the free-fall time to the gravitating center.
Results. We show the validity of these models by means of several examples of galactic and extragalactic accreting sources.
Key words: accretion, accretion disks / galaxies: nuclei
© ESO, 2017
1. Introduction
The optical behavior of the Be star in the high-mass X-ray transient A0535+26/HDE245770, as discussed by Giovannelli & Sabau-Graziati (2011), shows that the luminosity at periastron is typically enhanced by ~ 0.02 to a few tenths magnitude, and the X-ray outburst occurs eight days after the periastron. Indeed, an increase of the mass flux occurs at periastron. This flush reaches the external part of the temporary accretion disk around the neutron star and moves to the hot central parts of the accretion disk and the neutron star surface. The time necessary for this way is dependent on the turbulent viscosity in the accretion disk, as discussed by Giovannelli et al. (2013; GBK13). Giovannelli et al. (2015) discussed the behaviour of A0535+26/HDE245770 during 2014 by using the GBK13 ephemeris – JDopt−outb = JD0(2 444 944) ±n(111.0 ± 0.4) days, in order to describe the passages of the neutron star at the periastron at the 106th and 108th cycles after the optical event occurred on December 5, 1981 (811205-E), which suggested that the incoming X-ray outburst occurred on December 13, 1981 (Nagase et al. 1982). The results of Giovannelli et al. (2015) and the derived relationship ΔVmag vs. Ix can predict not only the arrival time of X-ray outbursts, but also its intensity Ix.
GBK13 have constructed a quantitative model for explaining the time delay between optical and X-ray outbursts. The mechanism proposed by GBK13 to explain the X-ray – optical delay in A 0535+26/HDE245770 is based on an enhanced mass flux propagation through the viscous accretion disk. The observed time delay is related to the motion of a high-mass flux region from the outer boundary of the neutron star Roche lobe to the Alfvén surface that is due to the action of the α-viscosity. This mechanism, known as the UV – optical delay (the delay of the extreme ultraviolet (EUV) flash with respect to the optical flash), was observed and modeled for cataclysmic variables (e.g., Smak 1984, 1998; Lasota 2001).
Lags between the optical and X-ray maxima in flashes had also been observed in active galactic nuclei (AGNs). Nandra et al. (1998) found a delay of ~4 days between UV and X-ray emissions in NGC 7469; Maoz et al. (2000) found a delay of ~100 days between optical and X-ray emissions in the Seyfert galaxy NGC 3516; Marshall et al. (2008) found a delay of ~15 days between optical and X-ray emissions in Mkr 509, and Doroshenko et al. (2009) found a delay of ≈10 days between R, I and X-ray luminosities in the Seyfert galaxy 3C 120; Shemmer et al. (2003) found a delay of 2.4 ± 1.0 days between optical and X-ray emissions in NGC 4051. The time lags observed in AGNs are much shorter than we would expect from the scaling in formulae (7) and (19). In the disk accretion model we expect time lags orders of magnitude larger than in galactic X-ray transients, and judging from observations, these lags are of the same order. We suggest that this is connected with a much faster contraction of accreting matter, which occurs in the case of quasi-spherical accretion of matter with low angular momentum.
In this paper we derive a simple formula for finding the time delay among events in different wavelengths, which is valid in general for all disk-accreting cosmic sources, and quantitatively discuss a model for time-lag formation in AGNs.
2. Equations of the accretion disk structure
When the characteristic time of variability of the mass flux along the accretion disk is longer than the relaxation time of the local disk equilibrium, it is possible to use the approximation of local equilibrium Shakura (1972), see also Bisnovatyi-Kogan (2011), to calculate the transient disk structure. The equilibrium along a radius of the accretion disk around a star with a mass M is determined by the Keplerian rotational velocity ΩK(1)Writing the equation of the vertical equilibrium in approximate algebraic form, we obtain
(2)where
is the speed proportional to the sound velocity, P and ρ are the (gas + radiation) pressure and density at the symmetry plane of the accretion disk, and h is the semi-thickness of the accretion disk. The specific angular momentum l of the matter in the accretion disk is connected to the rotation velocity as
(3)The mass flux through the disk at radius r is connected to the radial velocity vr as
(4)We use an α approximation for the turbulent viscosity (Shakura 1972) when the (rφ) component of the stress tensor trφ is written as
(5)where the phenomenological non-dimensional parameter α ≤ 1. The condition of stationarity of the angular momentum, in which the outward viscous radial flux of the angular momentum is balanced by the angular momentum of the inward flux of the mass, is written as (see, e.g., Bisnovatyi-Kogan 2011)
(6)The main input into the time lag comes from the outer regions of the disk with l ≫ lin. Then we have from Eqs. (4) and (6) the expression for the radial velocity in the form
(7)We also define the surface density Σ, and write Eq. (6) in light of Eq. (3), using condition l ≫ lin, in the form
(8)The equation of the local thermal balance in the accretion disk, when the heat produced by viscosity Q+ is entirely emitted through the sites of the optically thick accretion disk with a total flux Q−, at l ≫ lin is written as (see, e.g., Bisnovatyi-Kogan 2011)
(9)Here T is the temperature in the symmetry plane of the accretion disk, a is the constant of the radiation energy density, c is the speed of light, and ϰ is the Thompson (scattering) opacity of the matter.
3. Calculation of the time lag in the disk accretion model for galactic X-ray sources
Based on Eqs. (8) and (2), Eq. (9) may be written as (10)From the second relation in Eq. (8), using Eq. (2), we obtain the relation
(11)The time lag in transient X-ray sources is typically formed in the regions of the accretion disk that are dominated by gas pressure with P = Pg, and the scattering opacity ϰ= 0.2(1 + X), where X is the hydrogen mass fraction. The gas constant ℛ for the hydrogen-helium fully ionized plasma and the gas pressure are written as (Landau & Lifshitz 1980)
(12)Here k is the Boltzmann constant and mp is the proton mass. By multiplying Eqs. (10) and (11), we obtain, based on Eq. (12), the relation
(13)Our main goal is to express the local parameters of the accretion disk through the observed parameters: mass of the compact star (or a black hole) M, mass flux into the star Ṁ measured through its bolometric luminosity, and the effective temperature Teff of the object, evaluated from its spectra. The radius where radiation is the most effective during the transient accretion may be found from the thermal equilibrium condition by accepting that the radiation is emitted from the sides of the disk with an effective temperature Teff
(14)giving
(15)Using Eqs. (1) and (15), Eq. (13) is written in the form
(16)The radial velocity from Eq. (7) may be written as
(17)Based on Eqs. (15) and (16), we obtain
(18)For a given Ṁ in the flash, the value of Teff changes in time, while the wave of a large Ṁ moves to the star along the accretion disk (GBK13; Giovannelli et al. 2015; Bisnovatyi-Kogan et al. 2014). The equation for the radius of the hot layer flash-wave moving to the star is written as
(19)We can reduce it to the variable Teff instead of r, using Eqs. (15) and (18), when it is written in the form
(20)Writing Eq. (20) as
(21)we obtain its solution in the form
(22)For the initial condition Teff(0) = T0, we find the time τ, when the effective temperature Teff = T1, by the relation
(23)If T0 corresponds to the maximum of the optics and T1 corresponds to the X-ray maximum, then the value of τ represents the time lag between these two maxima, which for the transient X-ray source.
A0535+26/HDE245770 is close to eight days (Giovannelli & Sabau-Graziati 2011). For T1 ≫ T0, using Eq. (21), we have (24)By substituting numbers into Eq. (24) and introducing dimensionless values
(25)we obtain the value of the time lag between optical and X-ray maxima in the transient Be/X-ray source, resulting from a rapid increase in the mass flux of matter with X = 0.7 through the accretion disk, using Eq. (12), in the form
(26)This formula is useful to form an approximate idea about the time lag between optical and X-ray maximum emissions in a cosmic disk-accreting source, regardless of whether it is a white dwarf, a neutron star, or a black hole. The time lag depends on the mass as (~ m2/3), and only very weakly on the mass flux, ṁ1/15. The dependence on viscosity α and on T4 is stronger.
For A0535+26/HDE245770, using the viscosity α = 0.1, as derived by the GBK13 model, m = 1.4 and ṁ = 1, we obtain from Eq. (26) days. Taking into account the experimental time lag τexp ≃ 8 days, the derived temperature is T4 = 2.78. Owing to the very weak dependence of τ on the accretion rate, the time lag in this source remains practically the same for flashes of different intensity, like the flash observed by de Martino et al. (1989), with a mass accretion rate Ṁ ≃ 7.7 × 10-7M⊙ yr-1, ṁ = 77. When we use the neutron star mass M = 1.5 M⊙ (Giovannelli et al. 2007), the viscosity coefficient α = 0.15, and T4 ≃ 2.8 (De Loore et al. 1984), the time delay computed with Eq. (26), τ = 8.06 days, coincides with the experimental time delay (X-ray – optical) of ~ 8 days.
When we assume the mass of a white dwarf M=0.97 M⊙ (Giovannelli et al. 1983), a mass accretion rate Ṁ=4×1017 g s-1≈6.3×10-9M⊙ yr-1, ṁ=0.63 (Giovannelli & Sabau-Graziati 2012), a viscosity α = 0.2 (e.g. Smak 1998), and T4 = 4 (Hameury et al. 1999), the time delay in the cataclysmic variable SS Cygni computed with Eq. (26) is τ ≃ 1.8 days. The experimental time delay (UV – optical) is 0.9–1.4 days (Wheatley et al. 2003), and for α = 0.3, we obtain τ = 1.35 days, which is well inside the error box.
When we assume the mass of a neutron star M = 1.4 M⊙ (Waterhouse et al. 2016), a mass accretion rate Ṁ = 4 × 1017 g s-1 ≈ 6.3 × 10-9M⊙ yr-1, ṁ = 0.63 (Yamaoka et al. 2011; Meshcheryakov et al. 2013), a viscosity α = 0.2, and T4 ≃ 2.8, the time delay in the X-ray source Aql X-1 computed with Eq. (26) is τ ≃ 4.4 days. The observed time delay (X-ray – optical) of ~ 3 days (Shahbaz et al. 1998) is better reproduced for α = 0.3, when τ = 3.2 days follows from Eq. (26).
When we assume the mass of a black hole M ≃ 7 M⊙ (Orosz & Bailyn 1997), and the approximate relationship M/Ṁ ~ 107 yr (Karakula et al. 1984) in the black hole X-ray transient GRO J1655-40, we have Ṁ ~ 7 × 10-7M⊙ yr-1. For T4 = 3, we obtain from Eq. (26) a time delay τ ≃ 15.6 days at α = 0.2, and τ ≃ 11.3 days at α = 0.3. The optical precursor was observed ~ 6 days before the X-ray flash (Orosz et al. 1997), therefore Eq. (26) overestimates the observed value. This time delay was observed in the X-ray flash of the LMXB X-ray nova with a black hole, which is typically connected with instability that develops in the accretion disk after a sufficiently high mass is accumulated in the quiescent state. When an instability develops, it leads to non-stationary behavior with higher viscosity. This explains why Eq. (26), which has been obtained for a stationary accretion disk, overestimates the value of the time lag. However, when we use a higher viscosity value, for example, α = 0.6, the time lag computed with Eq. (26) is τ ≃ 6.5 days, in agreement with the observed time delay.
4. Time lags in AGN observations
Supermassive black holes in AGNs radiate as a result of accretion, and time lags similar to galactic X-ray sources are observed during their sporadic variability, see Marshall et al. (2008), Nandra et al. (1998), and other references below. The region of the optical emission in AGNs is situated somewhat closer to the black hole than in galactic sources, and the choice between radiative and gas pressure is not evident. To make this choice, we compared local solutions for the accretion disk structure and determined the boundary between the radiation-supported and gas-supported regions. The solution of the accretion disk structure in a local approximation gives the expressions for the pressure (in CGS units) in the equatorial plane for radiation-dominated (Pr), and gas-dominated (Pg) regions as (Bisnovatyi-Kogan 2011)
(27)At the boundary between two regions these pressures are equal, which determines the value of xb as
(28)Here and below, we use the notations
(29)The disk effective temperature at the boundary x = xb follows from Eqs. (14), (28), and (29) as
(30)For the characteristic temperature of the optical emission T4 = 2.5, we obtain from Eq. (30) the relation at which this temperature is situated at the border between gas- and radiation-dominated regions as
(31)for α = 0.2. The accretion rate at the critical Eddington luminosity corresponds to ṁ ≈ 16 for a Schwarzschild BH and to ṁ ≈ 2.5 for the limiting Kerr metric. In most AGNs, the optical emission therefore comes from the radiation-dominated regions, in contrast to the accretion into a binary galactic BH.
The formula for the time lag for the optics coming from a radiation-dominated disk region may be obtained in the same way as for the gas-dominated regions. Now (32)By multiplying Eqs. (10) and (11), we obtain the relation
(33)By dividing Eqs. (10) and (11), we obtain
(34)Using Eq. (32), we obtain
(35)When we use Eq. (35) in Eq. (33), we obtain
(36)When we use Eqs. (1), (15), and (36) in Eq. (7), we obtain the radial velocity in the accetion disk as a function of the effective temperature in the form
(37)Similar to the gas-dominated case above, we obtain the equation for Teff as
(38)Writing this equation as
(39)we obtain the expression for the time lag τr between the effective temperature Teff = T0, and Teff = T1 ≫ T0, in the radiation-dominated disk region as
(40)Using notations from Eq. (25), we obtain the expression for the time lag τr in the form
(41)For AGNs we have used the value of T4 = 2.5 when no direct data were available. For the α-viscosity we have used a value α = 0.2 when no other indications were available.
The application of Eq. (41) to some AGNs gives results that are much longer than those observed. This means that the mechanism responsible for the time lag in AGNs is different from that of galactic transient accreting sources. In particular, we have the estimations of the time lag of the objects listed below.
In Mrk 509 the mass of the supermassive black hole (SMBH) is estimated as M ≃ 1.4 × 108M⊙ (Peterson et al. 2004), and the mass flux Ṁ ~ 14 M⊙ yr-1 (Karakula et al. 1984). For the viscosity coefficient α = 0.2 and T4 = 2.5, the time delay computed with the formula (41) is τ ≃ 2.2 × 104 days, while the observational delay is τobs ≃ 15 days (Marshall et al. 2008). Mehdipour et al. (2011) observed no significant time lag between optical and X-ray variability. In our model of the time lag in AGNs, which is presented in the next section, the time lag with optics preceding X-rays is connected with the tidal disruption of stars approaching an SMBH. The type of the flare and a real time lag strongly depend on the properties of the disrupted star and may be very different. Therefore the absence of any time lag found by Mehdipour et al. (2011) may indicate a different (probably more compact) disrupted star, or any other physical mechanism that might have caused the flare.
In NGC 7469 the mass of SMBH is estimated as M ≃ 107M⊙ (Peterson et al. 2004; Onken et al. 2004), and the mass flux Ṁ ~ 1 M⊙ yr-1 (Karakula et al. 1984). For the viscosity coefficient α = 0.2, and T4 = 2.5 the time delay computed with Eq. (41) is τ ≃ 3.4 × 104 days, while the observational time delay (UV – optical) is τobs ≃ 4 days (Nandra et al. 1998).
In 3C 120 the mass of SMBH is estimated as M ≃ 5.5 × 107M⊙ (Onken et al. 2004; Vestergaard & Peterson 2006), and the mass flux Ṁ ~ 5.5 M⊙ yr-1 (Karakula et al. 1984). For the viscosity coefficient α = 0.2 and T4 = 2.5, the time delay computed with Eq. (41) is τ ≃ 2.5 × 104 days, while the observational time delays in this object are much shorter.
Different time lags (X-ray–R,I,B) in 3C 120 have been studied by Doroshenko et al. (2009). The results of this paper are rather controversial, and cannot be explained by a unique model. The authors report that the Rc- and Ic-band flux variations lag significantly behind the B-band fluxes by 3.9 and 6.2 days, respectively, when the position of the CCF centroid is considered at 0.8rmax. In addition, the X-ray variability on a long timescale also lags behind the B-band variations by 5.3 days. However, the confidence level of this estimate is only 87%. A more detailed analysis of the correlation between the X-ray and optical emissions revealed a fairly complex picture: the degree of correlation between the optical and X-ray flux variations is different at different times. They also claim in the discussion that according to their data, the X-ray emission in 3C 120 lags behind rather than leads the optical emission when they take the strong CCF asymmetry into account. This fact is also difficult to reconcile with the simple reprocessing model. We explain this particular type of time lag by the model described in the next section, where we use for our estimations τlag ≃ 10 days. Other types of the time lag mentioned in this paper could be explained by a reprocessing model or any other model. Doroshenko et al. (2009) conclude: “according to our data, the X-ray events may sometimes lead the optical ones (episode 3), sometimes occur almost simultaneously (episode 1) or may lag behind”. This clearly shows that the results of this paper cannot be interpreted in only one way.
In NGC 3516 the mass of SMBH is estimated as M ≃ 3.17 × 107M⊙ (Denney et al. 2010), and the mass flux Ṁ ~ 3.17 M⊙ yr-1 (Karakula et al. 1984). For the viscosity coefficient α = 0.2 and T4 = 2.5, the time delay computed with Eq. (41) is τ ≃ 2.8 × 104 days, while the experimental time delay (X-ray–optical) is τobs ≃ 100 days (Maoz et al. 2000).
In NGC 4051 the mass of SMBH is estimated as M ≃ 5 × 105M⊙ (Shemmer et al. 2003), and the mass flux Ṁ ~ 0.05 M⊙ yr-1 (Karakula et al. 1984). For the viscosity coefficient α = 0.2 and T4 = 2.5, the time delay computed with Eq. (41) is τ ≃ 5.6 × 104 days, while the experimental time delay (X-ray–optical) is τobs = 2.4 ± 1.0 days (Shemmer et al. 2003). In later observations, Alston et al. (2013) did not find any correlation between optical and X-ray light curves. This behavior is similar to the observational properties of MRK 509 that we discussed above, and may be explained in the same way.
The tidal disruption of a star by an SMBH was discovered by ASAS-SN with robotic telescopes in the center of the galaxy PGC 043234 and was followed by optical and X-ray instruments, including Swift and XMM (Miller et al. 2015). The discovery magnitude indicated a substantial flux increase over archival optical images of this galaxy. Archival X-ray studies rule out the possibility that PGC 043234 has a standard active nucleus that could produce bright flaring, because it was not detected in the ROSAT All-sky survey (Voges et al. 1999). As was noted, the multiwavelength light curves in Fig. 1 of Miller et al. (2015) clearly indicate a tidal disruption event. The authors noted that the X-ray points in this figure carry relatively large errors, but it is distinctly visible in this figure that the maximum of the X-ray flux occurs about five days later than the maximum in the optical band. This means that the lime lag in the outburst in the nucleus of the galaxy PGC 043234 has an evident similarity with the outbursts in the five nuclei described above. The timescale of the lag in this nucleus is also short and cannot be explained by the processes within a radiation-dominated region in the model of a non-stationary accretion disk.
We note note that a brief qualitative explanation of these lags in AGNs, very similar to the GBK13 quantitative model, was suggested by Marshall et al. (2008). The conclusions of this model evidently contradict observational data. Therefore we consider another model for the formation of the time delay in optical – X-ray flashes from AGN.
5. Model of a time lag formation in AGNs
The accretion onto an SMBH in AGN takes place from surrounding gas, presumably formed by stellar winds of surrounding bulge stars. The angular momentum of stars in a quasi-spherical bulge is low, therefore formation of an accretion disk may not occur, and accretion could take place in the form of a spherical flow. The process of the tidal disruption of the star that approaches an SMBH was investigated numerically in different approximations, see, for instance, Carter & Luminet (1983), Frolov et al. (1994), Diener et al. (1995), Marck et al. (1996), Ivanov et al. (2003), and Ivanov & Chernyakova (2006).
The flashes in AGNs, which are presumably connected to tidal disruptions of the surrounding stars in a close encounter with an SMBH, are accompanied by rapid ejection of matter with the formation of a jet flowing outside, and another rapid jet directed toward the SMBH. A large part of the inner jet moves to the SMBH with a high velocity that is on the order of the free-fall velocity. The tidal disruption leads to the optical flash, and the X-ray flash starts when the matter of the inner jet is sufficiently heated to radiate in the X-ray region. It is not clear from the numerical simulations how the angular momentum is distributed over the matter falling onto the SMBH, nor which part of the matter falls with a small impact parameter onto the SMBH. We suggest that in the observed events this part is sufficiently large to create a strong initial flash.
The time delay between the optical and X-ray flashes follows from observations, listed above; the radius at which the optical flash occurs is calculated for the motion with free-fall velocity vff as (42)Taking τff equal to observational time delay τobs, we obtain a radius of the optical flash ropt as
(43)where τobs is expressed in days and the SMBH mass m is given in solar masses.
A tidal disruption of a star occurs when the tidal force from the SMBH Ft becomes comparable to the gravitational force of the star Fs at the radius of a star Rs with a mass Ms. The radius of the tidal disruption rt when these two forces become equal is written as, see, for example, Bisnovatyi-Kogan et al. (1982),
Properties of stars tidally disrupted by SMBH in AGN’s.
We see that the flashes in AGN that originate in tidal disruption may occur when a giant star with a radius of between a few tens and a few hundreds of solar radii enters the tidal radius. The observational time lag between the optical and X-ray flashes is quite consistent with the model we considered here, in which the optical flash occurs at the radius of the tidal disruption and the X-ray flash occurs when the matter that is accreted with the free-fall velocity becomes hot enough as a result of adiabatic heating.
We note that another possibility to interpret the short time delays in AGNs is based on the irradiation model, see, for example, Ulrich et al. (1997). This model could explain recent extensive observations of NGC 5548 in X-rays (SWIFT), UV, and optical light (HST) in the “reverberation mapping” campaign, see Edelson et al. (2015) and Fausnaugh et al. (2016). In this object the UV and optical light curves are lagging the X-rays in short (one- to two-day) intervals, so that the lag time increases for longer wavelengths. The irradiation model qualitatively explains this behavior quite well, but the estimated size of the disk is a factor of three larger than the prediction from standard thin-disk theory (Fausnaugh et al. 2016). A sample of 21 active galactic nuclei was analyzed with data from the Swift satellite to study the variability properties of the population in the X-ray, UV, and optical band (Buisson et al. 2016). A correlated variability between the emission in X-rays and UV is significantly detected for 9 of the 21 sources and is consistent with the UV lagging the X-rays. This behavior would be seen if the correlated UV variations were produced by the reprocessing of X-ray emission. In all cases the observed UV lags are somewhat longer than expected for a standard thin disk. This may be connected to the incompleteness of the local accretion disk theory, see Novikov & Thorne (1973).
6. Conclusions
We developed models of time lags between optical and X-ray flashes for close-binary galactic sources with accretion disks and for an AGN with an SMBH that is embedded in a quasi-spherical bulge.
The time lag in disk-accreting galactic close-binary sources is based on a sudden increase in the accretion flow that starts at the disk periphery and is related to the optical maximum. The massive accretion layer propagates to the central compact source as a result of the turbulent viscosity. The X-ray flash occurs when this massive layer reaches the inner hot regions of the accretion disk and falls into the central compact object. The matter in the accretion disk moves inside with a speed that is determined by the turbulent viscosity. We described this model quantitatively and derived an analytic formula that determines the value of the time lag. This formula gives results that agree well with observational values.
The flashes in an AGN are considered in the model when a disruption of a star that is in the evolution phase of a giant enters the radius of strong tidal forces. The matter with a low angular momentum that is released by the star falls into the SMBH in the form of a quasi-spherical flow with a velocity that is close to the free-fall velocity. An X-ray flash occurss when the falling matter reaches the hot inner regions. The time lag observed in these sources is identified with the time of the matter falling from the tidal radius onto the central region. The values of the tidal radius that we calculated in this model were compared with the theoretical radii of a tidal disruption that depends on the masses of the SMBH and of the star, and on the radius of the star. Knowing the SMBH masses from observations, and making a reasonable suggestion for the stellar mass that is on the order of one solar mass, we obtained that the radii of the disrupted star are between a few tens and a few hundreds of R⊙. These radii are characteristic of stars of moderate mass on the giant phase of evolution, see, for instance, Bisnovatyi-Kogan (2011).
The matter with larger angular momentum that appeared in the disruption of the star is expected to form an accretion disk through which the matter will move to the center as a result of turbulent viscosity, similarly to flashes in close galactic binaries. This motion is much slower than free-fall velocity and may last for many years. After such a flash in AGNs, we therefore expect a long-duration irregular variability in the whole electromagnetic spectrum. As may be seen in Fig. 1 of Miller et al. (2015), the X-ray and optical luminosity enter the non-steady plateau with a nonzero signal, which may be connected with a transition to the disk-accretion stage.
The variability properties observed in many AGNs, where optical and UV emission lags the X-ray light curve, may be explained by the model in which an X-ray flash in the center of AGN is followed by reradiation of the surrounding accretion disk.
Acknowledgments
The work of GSBK was partially supported by the Russian Foundation for Basic Research Grant No. 14-02-00728 and No. OFI-M 14-29-06045, the Russian Federation President Grant for Support of Leading Scientific Schools, Grant No. NSh-261.2014.2, and by Basic Research Program P-7 of the Presidium of the Russian Academy of Sciences. We are indebted to NASAs ADS for its magnificent literature and bibliography serving.
References
- Alston, W. N., Vaughan, S., & Uttley, P. 2013, MNRAS, 429, 75 [NASA ADS] [CrossRef] [Google Scholar]
- Bisnovatyi-Kogan, G. S. 2011, Stellar Physics 2: Stellar Evolution and Stability. (Berlin-Heidelberg: Springer-Verlag) [Google Scholar]
- Bisnovatyi-Kogan, G. S., Churaev, R. S., & Kolosov, B. I. 1982, A&A, 113, 179 [NASA ADS] [Google Scholar]
- Bisnovatyi-Kogan, G. S., Klepnev, A. S., & Giovannelli, F. 2014, Int. J. Mod. Phys. Conf. Ser., 28, 1460186 [CrossRef] [Google Scholar]
- Buisson, D. J. K., Lohfink, A. M., Alston, W. N., & Fabian, A. C. 2016, MNRAS, 464, 3194 [Google Scholar]
- Carter, B., & Luminet, J.-P. 1983, A&A, 121, 97 [NASA ADS] [Google Scholar]
- De Loore, C., Giovannelli, F., van Dessel, E. L., et al. 1984, A&A, 141, 279 [NASA ADS] [Google Scholar]
- de Martino, D., Waters, L. B. F. M., Giovannelli, F., & Persi, P. 1989, in Two Topics in X-ray Astronomy, Vol. 1: X-ray Binaries, Vol. 2: AGN and the X-ray Background, eds. J. Hunt, & B. Battrick, ESA SP, 296, 519 [Google Scholar]
- Denney, K. D., Peterson, B. M., Pogge, R. W., et al. 2010, ApJ, 721, 715 [NASA ADS] [CrossRef] [Google Scholar]
- Diener, P., Kosovichev, A. G., Kotok, E. V., Novikov, I. D., & Pethick, C. J. 1995, MNRAS, 275, 498 [NASA ADS] [Google Scholar]
- Doroshenko, V. T., Sergeev, S. G., Efimov, Y. S., Klimanov, S. A., & Nazarov, S. V. 2009, Astron. Lett., 35, 361 [NASA ADS] [CrossRef] [Google Scholar]
- Edelson, R., Gelbord, J. M., Horne, K., et al. 2015, ApJ, 806, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Fausnaugh, M. M., Denney, K. D., Barth, A. J., et al. 2016, ApJ, 821, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Frolov, V. P., Khokhlov, A. M., Novikov, I. D., & Pethick, C. J. 1994, ApJ, 432, 680 [NASA ADS] [CrossRef] [Google Scholar]
- Giovannelli, F., & Sabau-Graziati, L. 2011, Acta Polytechnica, 51, 21 [Google Scholar]
- Giovannelli, F., & Sabau-Graziati, L. 2012, Acta Polytechnica, 52, 11 [NASA ADS] [Google Scholar]
- Giovannelli, F., Gaudenzi, S., Rossi, C., & Piccioni, A. 1983, Acta Astron., 33, 319 [NASA ADS] [Google Scholar]
- Giovannelli, F., Bernabei, S., Rossi, C., & Sabau-Graziati, L. 2007, A&A, 475, 651 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Giovannelli, F., Bisnovatyi-Kogan, G. S., & Klepnev, A. S. 2013, A&A, 560, A1 (GBK13) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Giovannelli, F., Bisnovatyi-Kogan, G. S., Bruni, I., et al. 2015, Acta Astron., 65, 107 [NASA ADS] [Google Scholar]
- Hameury, J.-M., Lasota, J.-P., & Dubus, G. 1999, MNRAS, 303, 39 [NASA ADS] [CrossRef] [Google Scholar]
- Ivanov, P. B., & Chernyakova, M. A. 2006, A&A, 448, 843 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ivanov, P. B., Chernyakova, M. A., & Novikov, I. D. 2003, MNRAS, 338, 147 [NASA ADS] [CrossRef] [Google Scholar]
- Karakula, S., Tkaczyk, W., & Giovannelli, F. 1984, Adv. Space Res., 3, 335 [NASA ADS] [CrossRef] [Google Scholar]
- Landau, L. D., & Lifshitz, E. M. 1980, Statistical physics. Pt.1 (Oxford: Pergamon Press) [Google Scholar]
- Lasota, J.-P. 2001, New Astron. Rev., 45, 449 [NASA ADS] [CrossRef] [Google Scholar]
- Maoz, D., Edelson, R., & Nandra, K. 2000, AJ, 119, 119 [NASA ADS] [CrossRef] [Google Scholar]
- Marck, J. A., Lioure, A., & Bonazzola, S. 1996, A&A, 306, 666 [NASA ADS] [Google Scholar]
- Marshall, K., Ryle, W. T., & Miller, H. R. 2008, ApJ, 677, 880r [NASA ADS] [CrossRef] [Google Scholar]
- Mehdipour, M., Branduardi-Raymont, G., Kaastra, J. S., et al. 2011, A&A, 534, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meshcheryakov, A., Khamitov, I., Eselevich, M., et al. 2013, ATel, 5114, 1 [NASA ADS] [Google Scholar]
- Miller, J. M., Kaastra, J. S., Miller, M. C., et al. 2015, Nature, 526, 542 [Google Scholar]
- Nagase, F., Hayakawa, S., Kunieda, H., et al. 1982, ApJ, 263, 814 [NASA ADS] [CrossRef] [Google Scholar]
- Nandra, K., Clavel, J., Edelson, R. A., et al. 1998, ApJ, 505, 594 [NASA ADS] [CrossRef] [Google Scholar]
- Novikov, I. D., & Thorne, K. S. 1973, in Black Holes, eds. C. Dewitt, & B. S. Dewitt (Les Astres Occlus), 343 [Google Scholar]
- Onken, C. A., Ferrarese, L., Merritt, D., et al. 2004, ApJ, 615, 645 [NASA ADS] [CrossRef] [Google Scholar]
- Orosz, J. A., & Bailyn, C. D. 1997, ApJ, 477, 876 [NASA ADS] [CrossRef] [Google Scholar]
- Orosz, J. A., Remillard, R. A., Bailyn, C. D., & McClintock, J. E. 1997, ApJ, 478, L83 [NASA ADS] [CrossRef] [Google Scholar]
- Peterson, B. M., Ferrarese, L., Gilbert, K. M., et al. 2004, ApJ, 613, 682 [NASA ADS] [CrossRef] [Google Scholar]
- Shahbaz, T., Bandyopadhyay, R. M., Charles, P. A., et al. 1998, MNRAS, 300, 1035 [NASA ADS] [CrossRef] [Google Scholar]
- Shakura, N. I. 1972, AZh, 49, 921 [Google Scholar]
- Shemmer, O., Uttley, P., Netzer, H., & McHardy, I. M. 2003, MNRAS, 343, 1341 [NASA ADS] [CrossRef] [Google Scholar]
- Smak, J. 1984, Acta Astron., 34, 161 [NASA ADS] [Google Scholar]
- Smak, J. I. 1998, Acta Astron., 48, 677 [Google Scholar]
- Ulrich, M.-H., Maraschi, L., & Urry, C. M. 1997, ARA&A, 35, 445 [NASA ADS] [CrossRef] [Google Scholar]
- Vestergaard, M., & Peterson, B. M. 2006, ApJ, 641, 689 [NASA ADS] [CrossRef] [Google Scholar]
- Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389 [NASA ADS] [Google Scholar]
- Waterhouse, A. C., Degenaar, N., Wijnands, R., et al. 2016, MNRAS, 456, 4001 [NASA ADS] [CrossRef] [Google Scholar]
- Wheatley, P. J., Mauche, C. W., & Mattei, J. A. 2003, MNRAS, 345, 49 [NASA ADS] [CrossRef] [Google Scholar]
- Yamaoka, K., Krimm, H. A., Tomsick, J. A., et al. 2011, ATel, 3686, 1 [NASA ADS] [Google Scholar]
All Tables
Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.