| BZs >gnuplot.ps Acta Polytechnica Vol. 51 No. 2/2011 On the Intrinsic Simplicity of Spectral Variability of GRBs A. Chernenko Abstract In this paper we present a Multi-Scale Correlation Analysis (MSCA) of the light curves of gamma-ray bursts recorded in different energy ranges. This analysis allows us to identify time intervals where emission variability can be reduced to a single physical parameter and can therefore be robustly attributed to a single physical emitter. The properties of these intervals can then be investigated separately, and the spectral properties of individual emitters can be analysed. The signatures of hidden dynamical relations between individual emitters are also discussed. Keywords: Gamma-ray bursts, variability, multi-scale, correlations, random matrices, eigenvalues. 1 Introduction: dimension of spectral variability Relativistically expanding sources produce radiation that manifests a rich spectral evolution and diverse light curve structures. Even if individual emitting regions have a similar and simple geometry and rest frame spectra, spectral and temporal evolution within any period of an observer’s time will be very diverse, since many individual emitters contribute to this period at different viewing angles and Lorentz factors, at different stages of development. In reality, the physical properties and even the physical nature of individual emitters within a single source can also be quite diverse. In this situation we believe that the very first stage of spectroscopy should be model independent, and should give a clear understanding of the follow- ing: 1. The dimension of variability for the en- tire period of the observations — how many spectral parameters are in fact independently variable during the period of the observations. 2. Segmentation into low-dimension peri- ods — can the entire period of the observations be split into sub-periods for each of which the internally measured dimension of variability is considerably lower than the dimension for the entire period 3. Hidden dynamics — if there is evidence that the emission properties at any moment of the observations are dependent on its properties at earlier moment(s) of time. If such signatures are found, can the dimension of variability be fur- ther decreased by taking into account the hidden dynamics. In this paper I will focus on the first two items on this list. 2 Measuring the dimension of variability by correlating light curves Correlation analysis of emission light curves recorded in different energy ranges (wave lengths) is a nat- ural, and model-independent way of estimating the dimension of spectral variability. The rank RCM of the corresponding correlation matrix CM is exactly its numerical measure. Figure (1) presents the light curves of GRB911118. These light curves measured in p = 8 energy ranges by BATSE/COMPTON Large Area Detectors with time resolution of 0.016 s pro- duced a data matrix with dimension p×n = 8×1 000. Figure (2) presents the corresponding correlation matrix. The spectrum of its sorted eigenvalues EV (i) (i = 1, . . . , 8) is also presented. It is clear that there are 2 non-zero eigenvalues, and therefore the rank of the correlation matrix RCM = 2. In this case, we would say that spectral variability is two- dimensional. This original result was first obtained by Cher- nenko and Mitrofanov in [1] and was interpreted as evidence for two non-correlated emission components related to physically distinct emission sources. Later similar results were obtained in a number of papers, e.g. [2]. There is also an alternative interpretation involving a single emission component with two in- dependently variable parameters [3]. The next logical step in the investigation, as was noted in the introduction, is segmentation, i.e. to look at the GRB through time windows of different widths and positions in order to find sub-intervals where the spectral variability is just one-dimensional, i.e. the covariance matrix has just a single non-zero eigenvalue. We call this approach Multi-Scale Cor- relation Analysis, and will present its results later 61 Acta Polytechnica Vol. 51 No. 2/2011 Fig. 1: Gamma-ray light curves of gamma-ray burstGB911118 detected by theBATSE/COMPTONexperiment are pre- sented asmeasured in 8 broad energy channels between 37 keVand approx. 2MeVwith time resolution of 0.016 seconds. The amplitudes of the light curves are arbitrarily scaled to emphasize their shapes Fig. 2: Correlation matrix calculated for the 8 light curves of GB011118 presented in Figure (1) is graphically shown in the left panel. In the right panel, the spectrum of its sorted eigenvalues EV (i) is presented. The eigenvalues are normalized by the largest EV (1). Only the two largest eigenvalues are evidently non-zero, so the rank of the correlation matrix RCM =2 in this paper. But before that two important state- ments should be put forward: • As our goal is to study sub-intervals of the GRB time history on multiple time scales down to the shortest ones where covariance can be calculated (i.e. n → 3), we need to deal with cases when n < p and the covariance matrix is p − n + 1 times singular. • In the real world, times series are “contami- nated” by Poisson noise. Therefore, to deter- mine the number of eigenvalues which are sig- nificantly larger than zero we need to know how to propagate Poisson noise from raw time series to noise in eigenvalues. 3 Statistical properties of eigenvalues of correlation matrices for Poisson noise time series of arbitrary duration Let us consider a data matrix with dimensions n × p, where p is number of time series and n is number of measurements (duration) for each of the time series. Let us assume that each of the time series is a station- ary Poisson process with parameter λ. Then, with- out loss of generalization we can transform the Pois- son process by means of the Anscombe formula [4] to 62 Acta Polytechnica Vol. 51 No. 2/2011 a nearly Gaussian process with standard deviation σ = 1: A : x ⇒ 2 √ x + 3 8 (1) The correlation matrices for Gaussian time series with σ = 1 are known in random matrices theory as Wishart matrices. Analytical results on their eigen- values are basically limited to the quasi-continuous density distribution of the eigenvalues in the asymp- totic case when p ≤ n → ∞, which is given by the Marcenko-Pastur density function [5]. Analytical re- sults on the distribution of individual eigenvalues are basically limited to the density distribution of the largest eigenvalue, which is described by the Tracy- Wisdom law [7]. Although this law was found to be reasonably accurate for n and p as small as 5 [6], for our investigation this is not sufficient since i) we need to know the statistical properties of smaller eigenval- ues in order to determine the dimension of spectral variability and ii) the Tracy-Wisdom law holds for n ≥ p, which poses an unnecessary restriction on our analysis. To sum up, in the case of n � ∞ and espe- cially when n <∼ p � ∞, when the correlation ma- trix becomes singular, existing analytical results are insufficient. Therefore, we attempted a number of Monte-Carlo simulations to derive the sample distri- bution for eigenvalues for data matrices with different n : [4, 8, . . . , 1 024] and p : [4, 8, 12, 16, 32, 64]. For each realization of data matrix D we com- puted, by means of SVD, the spectrum of eigenvalues. For each pair of n and p, the number of such random realizations was 107. In Figure (3) we present the probability distributions for non-singular eigenvalues for some values of n and p. The goal of this Monte-Carlo experiment was to estimate the one-sided confidence limits EV pi for sorted eigenvalues EVi, i : [1, n] for a typical set of α- values: α : [0.1, 10−2, 10−3, 10−4, 10−5]. In Figure (4) these confidence limits are presented as functions of n. 4 The quest for one-dimensional spectral variability Again, mathematically speaking, one-dimensional spectral variability during some time interval [t1, t2] takes place if a correlation matrix built from multi- channel time histories over this time interval has only one non-zero eigenvalue. In a real case when the Pois- son background is present, there will be no zero eigen- values. Replacing the “one non-zero eigenvalue” by “one significant eigenvalue” would not be correct be- cause if the S/N ratio is low, then even if the first few eigenvalues are quite comparable, only the first one could eventually be significant. Fig. 3: Probability distributions for individual sorted non-singular eigenvalues of sample correlation matrices for Poisson noise time series are presented for some values of time series number p and lengths n. The leftmost panels correspond to cases n < p when the correlation matrices are p − n +1 times singular and only n − 1 largest eigenvalues are non-zero. The rightmost panels correspond to asymptotic cases when P � n → ∞. In this case, the eigenvalues are densely and symmetrically grouped around one, as expected from the Marcenko-Pastur density function [5] 63 Acta Polytechnica Vol. 51 No. 2/2011 Fig. 4: Upper confidence limits are presented for the individual sorted non-singular eigenvalues of sample correlation matrices as derived from the corresponding probability distribution presented in Figure (3). The solid lines correspond to p = 0.5 representing therefore the most probable values of the eigenvalues, while the two dotted lines in each plot correspond to p =10−3 and p =10−5 Thus, for the purposes of this analysis the follow- ing definition of one-dimensionality will be used: 1. the first two eigenvalues are significant with a predefined significance level of α and 2. their ratio EV (1/2) = EV (1)/EV (2) is greater than a predefined level of �. None of the BATSE Gamma-ray bursts that we have studied manifested one-dimensional spectral variability during its entire duration [3]. Therefore, following the idea proposed in the Introduction, we attempted segmentation — an analysis of the sub- periods to find out whether one-dimensional variabil- ity exists on shorter times scales somewhere during a gamma-ray burst. 4.1 Segmentation: looking for intervals of one-dimensional spectral variability Gamma-ray bursts may last from a fraction of a sec- ond to hundreds of seconds, manifesting complicated temporal structures on different time scales. To iden- tify periods of one-dimensional spectral variability in such a complex phenomenon we have, in principle, to use a truly redundant multi-scale approach. In the case of multi-scale correlation analysis this means that the correlation matrices should be sampled for all possible time windows [ti, tj ], where i, j : [1..N ], N being the total number of time intervals for a given GRB. More precisely, the procedure of Multi-Scale Cor- relation Analysis was as follows: • For a moment t and time scale τ , the correlation matrix CM (t, τ ) size n× n is calculated for light curves in n energy channels for the time interval [t − τ /2, t + τ /2] • First two largest eigenvalues EV (1, t, τ ) and EV (2, t, τ ) of CM (t, τ ) are calculated • If both EV (1, t, τ ) and EV (2, t, τ ) are significant then the ratio EV (1/2, t, τ ) = EV (1)/EV (2) is calculated as a local measure of the variability dimension at moment t. If EV (1/2) � 1 we conclude that the dimension of variability at mo- ment t and time scale τ is equal to 1. For the time scale of τ = 1 s the result of this procedure is illustrated in Figure (5). This procedure is then performed for a number of time scales, and for each moment t the scale τmax, which maximizes EV (1/2, t, τ ) is determined: EV (1/2, t)max = max τ (EV (1/2, t, τ )) (2) Peaks of EV (1/2, t)max mark time intervals [t − τmax/2, t + τmax/2] where one-dimensional variabil- ity is most prominent. In Figure (6) the time histo- ries of EV (1/2, t)max and τmax(t) are presented to- gether with the original multi-channel light curves 64 Acta Polytechnica Vol. 51 No. 2/2011 Fig. 5: Segmentation of the GB911118 light curves into sub-intervals with one-dimensional variability is demonstrated using Multi-Scale Correlation Analysis on the single time scale τ = 1 s. The upper panel presents the time histories of the two largest eigenvalues EV (1, t) and EV (2, t). The ratio EV (1/2, t) for the same time scale τ =1 s is presented in the middle panel. The peaks of this ratio mark the center points of sub-intervals with one-dimensional variability. The bottom panel shows three most pronounced sub-intervals with one-dimensional variability in the original multi-channel light curves of GB911118 by means of dashed lines 65 Acta Polytechnica Vol. 51 No. 2/2011 Fig. 6: Segmentation of the GB911118 light curves into sub-intervals with one-dimensional variability on multiple time scales using Multi-Scale Correlation Analysis. Peaks of the EV (1/2, t)max curve that is presented in the upper panel mark the middle points of these intervals. The corresponding time scale τmax, which is presented in the middle panels, gives the width of these intervals. This is illustrated in the bottom panel, which shows the original multi-channel light curves of GB911118 by means of dashed lines for 3 arbitrarily selected peaks of EV (1/2, t)max of GRB911118. The time history of τmax(t) is step- wise because for computational reasons the calcula- tions were performed for the limited number of log- arithmically spaced values: τ = 0.064 · 2n (s) where n = 0, . . . , 7. 4.2 The nature of the second and higher dimensions and hidden dynamics of spectral variability Figure (5) shows that there are sub-intervals of GB911118 (e.g. around t = 3 s and t = 7 s) where the spectral variability is not one-dimensional at any time scale. There are two alternatives: i) emis- sion during such time intervals is of a different, two-dimensional, nature or ii) emission with one- dimensional spectral variability is persistent along the entire GRB but there is an additional component also with one-dimensional variability which emerges in various parts of the light curve and then becomes negligible. In principle, a possible physical connection be- tween emitters within one source should manifest it- self via a correlation between their light curves. How- 66 Acta Polytechnica Vol. 51 No. 2/2011 ever, in the case of relativistic expansion these light curves would likely be shifted in time and energy in the observer’s frame. This could add extra dimen- sions to otherwise low-dimensional spectral variabil- ity. To deal with this problem, we should allow lags between the individual light curves and look for the CM with the lowest rank RCM over the continuum of the lags. Such an analysis was done for the entire du- ration of GB911118 [8], and we found a combination of lags that makes RCM ∼ 1 for the entire burst. 5 Conclusions We can see that the spectral variability in GRBs is at most two-dimensional, and it may be possible to reduce it even to one-dimension by introducing time lags between the light curves in different energy chan- nels. This brings us to the hypothesis of a single emission source within a GRB with a single domi- nating variable physical parameter and some kind of “echoes” shifted in time and energy for geometrical reasons. A natural physical framework for modeling this behavior would be an optically thick expanding photosphere. However, this picture is in contradic- tion with models that involve more than one emission mechanism, e.g. Ryde et al [9]. References [1] Chernenko, A., Mitrofanov, I.: MNRAS, 274, 361–368 (1995). [2] Bagoly, Z., et. al.: A & A, 493, Issue 1, 51–54 (2009). [3] Chernenko, A.: in the Proceedings of the IAU 8th Asian-Pacific Regional Meeting, edited by S. Ikeuchi, J. Hearnshaw, and T. Hanawa, the As- tronomical Society of Japan, 2002, pp. 321–322. [4] Anscombe, F. J.: Biometrika, 35, 246–254 (1948). [5] Marcenko, V. A., Pastur, L. A.: Math. USSR-Sb. 1 507–1 536 (1967). [6] Iain, M., Johnstone: Ann. Statist. Volume 29, Number 2 (2001), 295–327. [7] Tracy, C. A., Widom, H.: Comm. Math. Phys. 177 727–754. [8] Chernenko, A.: in the Proceedings of XXIèmes Rencontres de Blois Windows on the Universe, June 21–26, 2009, Blois, France, (2009). [9] Ryde, F., et al.: ApJLett, 625, L95–L98 (2005). Anton Chernenko E-mail: anton@cgrsmx.iki.rssi.ru Space Research Institute 117 997, 84/32 Profsoyuznaya Str, Moscow, Russia 67