arxiv: v2 [astro-ph.he] 18 Jul 2018

Transcription

1 Draft version July 19, 2018 Typeset using LATEX default style in AASTeX61 SMOKE AND MIRRORS: SIGNAL-TO-NOISE AND TIME-REVERSED STRUCTURES IN GAMMA-RAY BURST PULSE LIGHT CURVES Jon Hakkila, 1, 2 Stephen Lesage, 2 Stanley McAfee, 2 Eric Hofesmann, 2 Corinne Maly Taylor, 2, 3 and Thomas Cannon 2 arxiv: v2 [astro-ph.he] 18 Jul The Graduate School, University of Charleston, SC at the College of Charleston, 66 George St., Charleston, SC , USA 2 Department of Physics and Astronomy, College of Charleston, 66 George St. Charleston, SC , USA 3 Department of Physics and Astronomy, Clemson University, Clemson, SC 29631, USA (Received July 19, 2018; Revised; Accepted) Submitted to ApJ ABSTRACT We demonstrate that the smoke of limited instrumental sensitivity smears out structure in gamma-ray burst (GRB) pulse light curves, giving each a triple-peaked appearance at moderate signal-to-noise and a simple monotonic appearance at low signal-to-noise. We minimize this effect by studying six very bright GRB pulses (signal-to-noise generally > 100), discovering surprisingly that each exhibits complex time-reversible wavelike residual structures. These mirrored wavelike structures can have large amplitudes, occur on short timescales, begin/end long before/after the onset of the monotonic pulse component, and have pulse spectra that generally evolve hard to soft, re-hardening at the time of each structural peak. Among other insights, these observations help explain the existence of negative pulse spectral lags, and allow us to conclude that GRB pulses are less common, more complex, and have longer durations than previously thought. Because structured emission mechanisms that can operate forwards and backwards in time seem unlikely, we look to kinematic behaviors to explain the time-reversed light curve structures. We conclude that each GRB pulse involves a single impactor interacting with an independent medium. Either the material is distributed in a bilaterally symmetric fashion, the impactor is structured in a bilaterally symmetric fashion, or the impactor s motion is reversed such that it returns along its original path of motion. The wavelike structure of the time-reversible component suggests that radiation is being both produced and absorbed/deflected dramatically, repeatedly, and abruptly from the monotonic component. Keywords: gamma-ray burst: general, stars: jets, methods: data analysis, methods: statistical Corresponding author: Jon Hakkila hakkilaj@cofc.edu

2 2 Hakkila et al. 1. INTRODUCTION It has become apparent that pulses are the basic radiative units of gamma-ray burst (GRB) prompt emission (e.g., see Hakkila & Preece (2011); Hakkila et al. (2015, 2018) and references contained therein), containing key information about the physical mechanisms of GRBs and the environments in which they occur. However, GRB pulses do not conform to the standard definition of single vibrations or short bursts of sound, electric current, light, or other waves (Oxford English Dictionary). Instead, a variety of observed properties characterize GRB pulses, including durations ranging from milliseconds to hundreds of seconds (Norris et al. 1996), longer decay times than rise times (pulse asymmetry) (Norris et al. 1996), longer durations at lower energies than at higher energies (Richardson et al. 1996), durations that anti-correlate with their peak fluxes (Hakkila & Preece 2011), fluences that correlate with their durations (Hakkila & Preece 2011), emissions that start nearly simultaneously in different energy bands (Hakkila & Nemiroff 2009), hard-to-soft spectral evolution (Norris et al. 2005; Hakkila et al. 2015), more pronounced hard-to-soft evolution for hard asymmetric pulses than for soft symmetric ones (Hakkila et al. 2015, 2018), similar correlated behaviors regardless of the burst class to which they belong (long, intermediate, or short; Hakkila & Preece (2011); Hakkila et al. (2015, 2018)), and more symmetric, longer duration light curves at lower energies than at higher energies (Hakkila et al. 2015). Additionally, the interpulse separations in short GRBs increase with the durations of both the initial and subsequent pulses (Hakkila et al. 2018). This effect has not been apparent in long GRB pulses (Ramirez-Ruiz & Fenimore 2000). GRB pulses are non-monotonic in nature (Hakkila & Preece 2014; Hakkila et al. 2015, 2018). They undergo brighter intensity fluctuations than are expected from only Poisson noise, and these fluctuations are in phase with the underlying monotonic pulse structure. The fluctuations are thus not indicative of noise, nor of isolated pulses, but rather of pulse components. This non-monotonicity is critically important in characterizing GRB pulses, because otherwise every observed non-stochastic fluctuation in a pulse light curve can be interpreted as being a separate pulse, leading to improper conclusions about the number of pulses, pulse characteristics, and inter-pulse separations. The most basic non-monotonic pulse fluctuations are smoothly-varying, triple-peaked structures that can be viewed by fitting a GRB pulse light curve with a monotonic function, subtracting the fit from the data, and examining the remaining pulse residuals. Although pulses generally exhibit hard-to-soft evolution, those with noticeable triple-peaked structures tend to undergo spectral re-hardening just before or at the time of a peak, and this can give some pulses (particularly those that are softer and more symmetric) the appearance of having intensity-tracking spectral evolution. Although some pulses exhibit smooth monotonic or triple-peaked structures, others exhibit more complex, chaotic structures that are often particularly noticeable at higher energies. Pulses observed in noisy data show little structure and are best described as monotonic. We use the word structure to describe variations exceeding those commonly expected from a monotonic GRB pulse light curve. This structure definition refers to any intensity variation that results in an inadequate monotonic fit to an otherwise distinct emission episode. As such, it is instrument-dependent, because bright GRB pulses often exhibit structure whereas faint pulses rarely do. We have recently developed a classification scheme for Short GRB pulses (Hakkila et al. 2018) which demonstrates that the amount of structure in GRB pulses increases as the signal-to-noise (S/N) increases, and as the temporal bin size decreases. In other words, structure seems to be an inherent feature of GRB pulses that is not always apparent in data because it can get washed out by instrumental effects. Structure is less pronounced in faint GRB pulses than in bright ones (Hakkila & Preece 2014; Hakkila et al. 2015, 2018), and S/N plays an important role in this. Signal-to-noise can decrease for a variety of instrumental reasons, including inefficient photon detection, small detector surface area, decreased temporal bin size, decreased spectral range, increased spectral resolution, and detection at lower (noisier) energies. It seems intuitively obvious that structure and noise should become indistinguishable from one another when they have comparable amplitudes. However, structure and an underlying smooth pulse profile are both present in light curves of GRB pulses. As signal-to-noise decreases, does structure disappear before or after the smoothly-varying remainder of the pulse disappears? How does this affect the properties of pulses measured at low signal-to-noise relative to those at high signal-to-noise? Is it possible that other pulse properties might accompany these changes at low signal-to-noise (e.g., duration, asymmetry, and spectral evolution), and what might such changes, if observed, tell us about the influence of instrumental biases on GRB pulse measurements? We wish to explore these questions further using pulses from bright Long GRBs observed by BATSE (the Burst And Transient Source Experiment; Horack (1991)) on NASA s Compton Gamma Ray Observatory (CGRO). BATSE had one of the largest surface areas of any GRB experiment and thus was able to observe many GRBs at large S/N ratios.

3 Smoke and Mirrors 3 2. EMPIRICAL PULSE AND RESIDUAL MODEL 2.1. The Pulse Model In order to explore structure as a function of S/N, we must first define and characterize structureless pulses. We do this by using the monotonic asymmetric intensity function of Norris et al. (2005): I(t) = Aλe [ τ1/(t ts) (t ts)/τ2], (1) where t is time since trigger, A is the pulse amplitude, t s is the pulse start time, τ 1 is the pulse rise parameter, τ 2 is the pulse decay parameter, and the normalization constant λ is λ = exp [2(τ 1 /τ 2 ) 1/2 ]. The intensity is the count rate obtained from the total counts observed in BATSE energy channels 1 (20 50 kev), 2 ( kev), 3 ( kev), and 4 (300 kev 1 MeV) at 64 ms temporal resolution. Poisson statistics and a two-parameter background counts model of the form B = B 0 + BS t are assumed (where B is the background counts in each bin and B 0 and BS are constants denoting the mean background (counts) and the rate of change of the mean background (counts/s)). This model can be used to produce observable pulse parameters such as pulse peak time (τ peak = t s + τ 1 τ 2 ), pulse duration (w) and pulse asymmetry (κ). As a result of the smooth rapid rise and more gradual fall of the pulse model, w and κ are measured relative to some fraction of the peak intensity. Following (Hakkila & Preece 2011), we measure w and κ at I meas /I peak = e 3 (corresponding to 0.05I peak ), or where µ = τ 1 /τ 2, and w = τ 2 [9 + 12µ] 1/2, (2) κ [1 + 4µ/3] 1/2. (3) Asymmetries range from symmetric (κ = 0 when τ 1 τ 2 and µ ) to asymmetric (κ = 1 when τ 1 τ 2 and µ 0) with longer decay than rise times The Residual Model Residuals to the Norris et al. (2005) model are generated by obtaining the best monotonic pulse fit for the observed pulse light curve and subtracting this fit from the data. Small yet distinct deviations in the residuals are found to be systematically in phase with the monotonic light curve (Hakkila & Preece 2014), and these deviations are needed to accurately describe GRB pulse shapes. Although the deviations are closely aligned with the pulse duration, they are not always contained within it. Thus we have defined a larger fiducial time interval w f as w f = t end t start = 4.4τ 2 ( 1 + µ/2 + 1) + τ 1 τ 2, (4) with the fiducial end time t end given by t end = w 2 (1 + κ) + t s + τ peak (5) and the fiducial start time t start given by t start = t s 0.1[ w 2 (1 + κ) τ peak]. (6) The fiducial time interval is a unitless interval introduced to demonstrate the generic nature of GRB pulse residual structures, regardless of the durations and asymmetries of the underlying pulses containing them. Residuals are obtained using this time interval so that normalized pulse residual structures can be directly compared with one another, and so that goodness-of-fit measurements used in making these comparisons are not based on the amount of background sampled. The residuals exhibit variations that can be fitted with an empirical function (Hakkila & Preece 2014) in the fiducial time interval f where 0 t 1: aj 0 ( Ω f [t 0;f t 0.005]) if t < t 0;f res(t ) = a if t 0;f t t 0;f aj 0 ( s f Ω f [t t 0;f 0.005]) if t > t 0;f (7)

4 4 Hakkila et al. Here, J 0 (x) is an integer Bessel function of the first kind, t 0;f is the unitless time of the residual peak, a f is the amplitude of the residual peak, Ω f is the unitless Bessel function s angular frequency that defines the timescales of the residual wave (a large Ω f corresponds to a rapid rise and fall), and s f is a dimensionless scaling factor that relates the fraction of time which the function pre-t 0;f has been compressed relative to its time-inverted form post-t 0;f. The time during which the pulse intensity is a maximum is required to be a plateau instead of a peak because otherwise the width of the Bessel function s first lobe would not be broad enough to capture the extent of the central peak. Therefore, a duration of w plateau 0.010w fid is hard-coded into the fitting routine to create a small plateau structure in the function. Without evidence in the pulse shape that the Bessel function continued beyond the third zero (following the second half-wave), Hakkila & Preece (2014) truncated the function at the third zero J 0 (x = ±8.654), giving the fitting function a triple-peaked appearance. Because the central plateau most commonly aligns with the peak of the monotonic pulse, the three residual function peaks were named the precursor peak, the central peak, and the decay peak, respectively (Hakkila & Preece 2014). In some pulses the central peak of the residual function does not closely align with the peak of the monotonic pulse (see Hakkila et al. (2018) and Section 4 of this manuscript). The reason for this misalignment is still not well-understood but it appears to signify that the residual function and the monotonic function are separate pulse components. After fitting, the unitless fiducial values are converted back to the observed time interval using: and a = a f (8) s = s f, (9) t 0 = t 0;f (t end t start ) + t start, (10) Ω = Ω f /(t end t start ), (11) where t start and t end are the real time values corresponding to the start and end of the fiducial time interval. The Bessel function frequency Ω is now given in units of s 1 and the residual peak time t 0 is now given in units of s. The flux contribution that the residual function makes to the combined pulse fit, as well as its relative temporal alignment, is capable of significantly altering a pulse shape from that of monotonicity. Although the precursor peak and the decay peak of the residual function typically fall during the rise and decay portions of the monotonic fitting function, in some cases the precursor peak can occur before the formal rise of the monotonic pulse. When the precursor peak is harder than the central peak, the hard-to-soft pulse evolution is pronounced. When the precursor peak is softer than the central peak, the pulse begins evolving hard-to-soft but then re-hardens to a maximum value around the time of the central peak, and the pulse hardness appears to track with intensity. Additionally, the central residual peak and the peak of the monotonic pulse do not always align with one another, producing an extra smooth pulse structure in some cases. The 64 ms S/N is a measure of the peak brightness of each GRB, relative to its background count measured with 64 ms temporal resolution. The S/N is S/N = (P 64 B)/ P 64, (12) where P 64 is the 64 ms peak counts and B is the mean background count. This signal-to-noise ratio is primarily appropriate when analyzing GRB pulses fit on the 64 ms time scale. Using the approach of Hakkila et al. (2018), we classify GRB pulses based on the amount of structure they exhibit, as defined relative to the χ 2 ν p values of the pulse fits. This method first fits each pulse with the monotonic Norris et al. (2005) fitting function and analyzes the goodness-of-fit with a χ 2 test. We define χ 2 over the fiducial timescale and ν from the Norris et al. (2005) model as the number of temporal bins minus the number of pulse- and background-fit parameters two for the background and four for a single pulse. The p value associated with χ 2 ν is the probability that a value of χ 2 or less having ν degrees of freedom will be measured in a random distribution. We consider good fits (indicating relatively smooth light curves) to be those having standard best-fit p values of p best Next, this monotonic pulse fit is subtracted from the light curve leaving a residual light curve. The residual light curve is then fitted with the Hakkila & Preece (2014) residual fitting model. Finally the residual fit is added to the monotonic pulse fit to produce a total fit, and this is tested to see whether or not adding the residual fitting function significantly improved the overall pulse fit. A χ 2 test is used to indicate whether or not the residual function should be included

5 Smoke and Mirrors 5 in the fit: χ 2 is the difference in χ 2 obtained from the Norris et al. (2005) model minus that obtained from the Norris et al. (2005) model combined with the Hakkila & Preece (2014) residual model. The difference in the number of degrees of freedom between these fits is four for a single pulse. If adding the residual fit significantly improved the pulse fit, then the combined models act as the pulse fit. But if the addition of the residual fit does not improve the overall pulse fit, then only the original monotonic pulse fit is used to represent the final pulse fit. We require a χ 2 p value of p 10 3 for the model to be considered significantly improved and therefore use the combined pulse models to represent the final pulse fit. We use this approach to classify GRB pulses into one of four groups. Simple pulses are those best characterized by Equation 1 alone (p best ). Blended pulses are fits improved using Equation 7 (p 10 3 and p best ). Structured pulses have many characteristics that can be explained by the pulse and residual models, but also have significant deviations from these structures (p 10 3 and 10 5 p best < ). Complex pulses have complicated light curves (p best < 10 5 ). 3. ANALYSIS We hypothesize that bright, complex GRB emission episodes are really structured pulses that can be reduced to triple-peaked and monotonic shapes with the loss of S/N. We test this hypothesis by analyzing four of BATSE s brightest GRBs. Specifically, we study BATSE triggers 143 (GRB ), 249 (GRB ), 5614 (GRB ), and 7301 (GRB B) (Paciesas et al. 2000). All four bursts have been observed in BATSE energy channels 1 4 (20 kev to 1 MeV) with 64 ms temporal resolution. We initially limited our sample to five pulses with S/N > 100. However, because BATSE 7301 contains two bright pulses, we also chose to analyze the second pulse of BATSE 143 (S/N = 55.9) for completeness. As we will show in Section 3.1.1, any pulse with S/N > 25 can be considered to be bright for the purposes of this analysis. BATSE trigger 143 contains two distinct emission episodes peaking at 1.7 seconds and 48 seconds after the trigger. The brighter first event (143p1) has a S/N = 131.5, a duration of w = 9.9 seconds, an asymmetry of κ = 0.75, and has considerable temporal structure, while the smoother and fainter second event (143p2) has S/N = 55.9, w = 11.2 seconds, κ = 0.59, and shows clear evidence of the triple-peaked structure. BATSE trigger 249 is a complex, time-symmetric emission episode with S/N = peaking 21.9 seconds after the trigger and containing at least 5 separate emission peaks. The entire episode lasts some 38 seconds, although the fitted pulse lasts only w = 20.8 seconds and has an asymmetry of κ = BATSE trigger 5614 is a bright (S/N = 177.2) single-peaked asymmetric emission episode (κ = 0.69) peaking 8.4 seconds after the trigger with a duration of w = 6.7 seconds. Although the burst has the appearance of a fairly smooth single pulse, the main part of the pulse rise does not occur until roughly 7 seconds after the trigger. Of the six pulses in the sample, this one most clearly exhibits the triple-peaked pulse shape described in previous papers. BATSE trigger 7301 is a complex burst with what appear to be four separate complex emission episodes. For reasons we will discuss later, we identify only the two brightest episodes as pulses 7301p1 and 7302p2 (these peak 15.8 seconds and seconds after the trigger, respectively) while excluding the fainter emission episodes starting at the trigger (t = 0 seconds) and at around t = 32 seconds. Pulse 7301p1 is a bright, complex episode with S/N = 100.9, w = 10.7 seconds and κ = Pulse 7301p2 is a bright, complex episode with S/N = 144.3, w = 29.3 seconds and κ = 0.52; the pulse contains three bright peaks each lasting a few seconds which occur during fainter emission lasting roughly 50 seconds. We initially assume that each of the bright, complex emission episodes in BATSE triggers 143, 249, 5614, and 7301 is composed of a single structured pulse. This is completely antithetical to the standard pulse-fitting approach in which each pulse would be treated as a light-curve variation that could not be attributed to noise. Using such an approach, 143p1, 249, 7301p1, and 7301p2 would be best fit by dozens of pulses. Even the apparently smoothly-varying 143p2 and 5614 would consist of at least half a dozen pulses because the burst is so structured at this S/N level. Fits to the six pulses using the monotonic model described by Equation 1 and the residual model described by Equation 7 are shown in Figures 1 6. The left panels of Figures 1 6 show the residual structure while the right panels show both the monotonic fit (dotted line) and the combined fit (solid line). Each pulse is classified as complex using the formal criteria described in Section 2.2, indicating that the combination of the monotonic pulse model and the residual model do not adequately describe the pulse light curves. However, in each case it can be seen that the models contribute to only an approximate first-order fit by creating envelopes that contain the most noticeable structures, and that the residual fits are only rough approximations to the actual residual structure.

6 6 Hakkila et al. Figure 1. BATSE trigger 143p1 observed at S/N = The left panel shows the Hakkila & Preece (2014) residual structure, while the right panel shows both the Norris et al. (2005) pulse fit (dotted line) and the combined Norris et al. (2005) pulse plus Hakkila & Preece (2014) residual fit (solid line). In this and all subsequent plots, residual count rate and count rate refer to the number of counts per 64 ms bin. Figure 2. BATSE trigger 143p2 observed at S/N = The left panel shows the Hakkila & Preece (2014) residual structure, while the right panel shows both the Norris et al. (2005) pulse fit (dotted line) and the combined Norris et al. (2005) pulse plus Hakkila & Preece (2014) residual fit (solid line). Figure 3. BATSE trigger 249 observed at S/N = The left panel shows the Hakkila & Preece (2014) residual structure, while the right panel shows both the Norris et al. (2005) pulse fit (dotted line) and the combined Norris et al. (2005) pulse plus Hakkila & Preece (2014) residual fit (solid line).

7 Smoke and Mirrors 7 Figure 4. BATSE trigger 5614 observed at S/N = The left panel shows the Hakkila & Preece (2014) residual structure, while the right panel shows both the Norris et al. (2005) pulse fit (dotted line) and the combined Norris et al. (2005) pulse plus Hakkila & Preece (2014) residual fit (solid line). Figure 5. BATSE trigger 7301p1 observed at S/N = The left panel shows the Hakkila & Preece (2014) residual structure, while the right panel shows both the Norris et al. (2005) pulse fit (dotted line) and the combined Norris et al. (2005) pulse plus Hakkila & Preece (2014) residual fit (solid line). Figure 6. BATSE trigger 7301p2 observed at S/N = The left panel shows the Hakkila & Preece (2014) residual structure, while the right panel shows both the Norris et al. (2005) pulse fit (dotted line) and the combined Norris et al. (2005) pulse plus Hakkila & Preece (2014) residual fit (solid line).

8 8 Hakkila et al Smoke: Pulse Properties as Functions of Signal-to-Noise We estimate the properties of these six BATSE pulses at lower S/N by reducing their assumed intrinsic counts distributions and adding stochastic instrumental noise. The pulse plus residual fits described thus far have been based on the summed counts in BATSE s four energy channels Since some pertinent pulse properties are energy-dependent (e.g., spectral lags), we use a process that produces reductions in each individual energy channel rather than only in the summed four-channel light curve. For each of the six BATSE pulses listed, we estimate their properties at lower S/N values by removing the background from each energy channel and assuming that the remaining photon flux distribution is a fair representation of the intrinsic flux distribution. When summed together, these channel-dependent counts allow us to reproduce the intrinsic summed four-channel counts observed by the instrument. We then define roughly ten equal logarithmic four-channel count intervals between the observed peak count rate and a measure close to background (ostensibly near S/N 3). For each energy channel we iteratively decrease the peak counts proportionally through these logarithmic decrements, add energy-dependent Poisson noise to the results, and recombine the energy channel-dependent light curves to obtain a reduced-intensity four-channel light curve. This process assumes that the reduced pulse would have been detected with a similar background rate if it had been detected at a lower S/N, an assumption that is only partially valid because the background rate varied slowly and regularly throughout the Compton Gamma-Ray Observatory s mission in response to the satellite s eccentric orbit. The noisification process also assumes that there is no systematic coupling between detector energy channels, and that the detector response does not depend significantly on count rate. There are potential problems with both of these assumptions: the non-diagonal energy response of BATSE s Large Area Detectors (LADs) created a coupling between energy channels, and low S/N observations were affected by spacecraft orientation and relative direction between the Earth and a GRB (e.g., see Pendleton et al. (1996)). However, we believe that both assumptions are valid as first-order approximations Pulse Structure As S/N decreases, pulse structure becomes less noticeable. At high S/N, pulse structure is easily distinguishable and separable from Poisson noise. At low S/N, structure becomes indistinguishable from noise, while the underlying smooth monotonic pulse remains identifiable and measurable. This becomes an issue when determining the number of pulses potentially present in a GRB light curve: at low S/N the monotonic pulse dominates over structure, and the light curve appears to be that of a single pulse. At moderate and high S/N these structures are clearly not the result of noise, and can even appear to be so pronounced that they can be mistaken for separate pulses. However, confusion in the number of pulses can be reduced by recognizing that the structures are linked in phase with the main body of emission, either sitting atop the monotonic pulse (e.g., 143p1; see Figure 1) or a connected to it as a continuum of episodes like cutout paper doll chains (e.g., 249; see Figure 3). The importance of this relationship will be discussed further in Section 4. Figure 7 demonstrates two stages in the intensity reduction of BATSE trigger 249. If this pulse had been detected at S/N = 5.2 (left panel), then it would be classified as blended because both monotonic and residual pulse models are needed to completely account for the pulse light curve variations. If the pulse had been detected at S/N = 4.3 (right panel), then a fit to the residual function would not be statistically significant (as measured by p ) and the pulse would be classified as simple. Figure 8 summarizes the effect of S/N on pulse classification for the pulses in this study. All pulses are classified as complex at high S/N, and all pulses are classified as simple at low S/N. The S/N values at which these transitions occur (complex to structured, structured to blended, and blended to simple) vary from pulse to pulse depending on the intrinsic pulse structure, the background at the time the burst was detected, and the way in which the random background fluctuations have interacted with the reduced pulse. It is useful to compare the S/N distributions of pulses populating each class to the overall pulse S/N distribution. To do this we must estimate the numbers of pulses in each class because each pulse must be fitted before it can be classified, and this is a time-consuming procedure. Even estimating the overall pulse S/N distribution must be performed with care because measuring the mean background for each burst first requires us to separate the background from the pulse (and from multiple pulses when the burst contains more than one); this again requires us to fit each pulse before we can characterize its S/N.

9 Smoke and Mirrors 9 Figure 7. Two stages in the intensity reduction of 249. The left panel shows what the pulse can look like at S/N = 5.2, where the pulse and residual fits combine to place the pulse in the blended class. The right panel shows what the pulse can look like at S/N = 4.3, where the residuals do not statistically improve the fit and the pulse is classified as a simple pulse. Although the addition of the residual function does not significantly improve the fit at S/N = 4.3, the combined pulse and residual fits are shown for visual comparison. 4 Complex (4) pulse class 3 2 Structured (3) Blended (2) 1 Simple (1) signal-to-noise Figure 8. GRB pulse classification as S/N is reduced. Indicated are the general paths taken during S/N reduction by each of the five pulses shown: 5614 (red dotted line), 143p1 (blue dashed-dotted line), 143p2 (gray dashed line), 249 (orange dasheddouble-dotted line), and 7301p1 (green solid line). The reduced classification of 7301p2 is not plotted, but is similar to that of 249. Simple pulses (class 1) are best characterized by monotonic pulses, blended pulses (class 2) are best characterized by a monotonic pulse overlaid by a smooth triple-peaked residual function, structured pulses (class 3) are improved but not adequately characterized by a pulse plus the residual function, and complex pulses (class 4) are too structured to be described by these simple pulse models. In order to carry out our estimate we first note that there is a strong correlation between S/N values and published BATSE 1024 ms C max /C min values (the count rate at which a GRB triggers relative to the minimum count rate for a trigger to occur): S/N = 4.606(C max /C min ) 0.669, (13) with a correlation coefficient of R 2 = From this correlation we can estimate the S/N distribution of the brightest pulses in each BATSE burst using C max /C min (again, we cannot identify the S/N distribution of all pulses without going through the time-consuming process of fitting all pulses). Second, we obtain the distribution of BATSE GRB pulses fitted in previous papers (e.g., Hakkila & Preece (2014); Hakkila et al. (2015)). We make sure to exclude GRBs with T 90 durations less than 2 seconds (T 90 is the duration spanned by 90% of the flux) because 64 ms binning distorts any structure present; Hakkila et al. (2018)). We overlay the estimated pulse distribution with the observed distributions of simple, blended, structured, and complex pulses in Figure 9. It is apparent that the current fitted

10 10 Hakkila et al. Figure 9. The estimated S/N distribution of the brightest BATSE pulse in each burst (solid black line). The plot is overlaid with the measured S/N values of classified BATSE pulses (primarily from Hakkila & Preece (2014); Hakkila et al. (2015)) Simple pulses (blue dotted line) are found at the lowest S/N ( log(s/n) = 0.83 or S/N = 6.7), blended pulses (dashed green line) have ( log(s/n) = 1.03 or S/N = 10.7), structured pulses (dotted-dashed orange line) have ( log(s/n) = 1.15 or S/N = 14.1), and complex pulses (solid red line) have ( log(s/n) = 1.51 or S/N = 32.7). The large S/N of the pulses in this sample (S/N > 100) can be contrasted with the rest of the distribution. Note also that all pulses with S/N 25 have been classified as complex. pulse distribution favors pulses exhibiting the triple-peaked structure (e.g., blended and structured) because of our bias towards looking for evidence of it. Based on the classification criteria of these pulses and on the distribution of BATSE pulse S/N values, we estimate that at least 50% of pulses from long BATSE GRBs are found at S/N 8. Most of these are likely to be classified as either simple or blended, although pulses containing more structure are also found at these S/N values. This result explains both why Hakkila & Preece (2014) were able to establish a simple residual model using low-s/n pulses and why the residual model did not suffice when applied to bright pulses. It is entirely likely that all GRB pulses exhibit complex structures, but only a small fraction of pulses exhibit these structures without having them smoothed by instrumental effects related to S/N. However, caution must be exerted in generalizing and interpreting why GRB pulse structure exists at low S/N. Not all faint pulses are simple; some have been classified as blended, structured, and even complex. Smeared-out, complex, intrinsic pulse structures can appear indistinguishable from overlapping pulses at low S/N Pulse Duration and Asymmetry Although observing a GRB pulse at reduced S/N smears out its pulse structure, the pulse s duration (w) and asymmetry (κ) extracted from the monotonic model shown in Equation 1 do not appear to be systematically affected. Figure 10 demonstrates this for the bright pulses in this sample as S/N is reduced. Measured w and κ values are, however, affected statistically by large uncertainties at low S/N (S/N < 6). We also note that the pulses with the largest amounts of structure preceding and following the main monotonic component (249, 7301p1, and 7301p2) have the largest duration and asymmetry uncertainties at all S/N values. Pulses in which structure overlays the monotonic component (143p1, 143p2, and 5614) have smaller uncertainties that actually serve to guide the monotonic pulse fit and decrease uncertainty in measuring the pulse fit parameters (τ 1, τ 2, t s, and A) Pulse Spectral Lags Spectral lags indicate the time difference between the optimal alignment of pulse s high-energy and low-energy fluxes. Lags are measured using the cross-correlation function (CCF). The CCF is either applied to the data (data lag); e.g., Cheng et al. (1995); Norris et al. (1996); Band (1997)) or to the fitted pulses (fit lag; Norris et al. (2005); Hakkila & Preece (2011)). Spectral lags correlate with GRB luminosity (Norris et al. 1996) and provide a simple measurement of spectral properties. Although there appears to be a general relationship between both types of lag and pulse evolution (e.g., Kocevski et al. (2003); Schaefer (2004); Ryde (2005); Hakkila & Preece (2011); Hakkila et al. (2015, 2018)), some pulses have negative data lags (Ukwatta et al. 2010; Roychoudhury et al. 2014; Chakrabarti et al. 2018) and data lags are often inconsistent with fit lags (e.g., Ukwatta et al. (2010)).

11 Smoke and Mirrors 11 Figure 10. Fitted durations (left panel) and asymmetries (right panel) of pulses 143p1 (plusses), 143p2 (crosses), 249 (triangles), 5614 (diamonds), 7301p1 (asterisks), and 7301p2 (squares) as their S/N values are reduced. Fitted pulse durations and asymmetries do not change systematically with decreased S/N, with the exception of 7301p2 (due to insufficient contributions to the fit from the loss of emission in the pulse s broad wings ). Some measurements appear to undergo discontinuous shifts as S/N drops; this happens when pulse structure abruptly disappears into the noise and the fit parameters have to fit a more simply-structured pulse. Uncertainties increase for all pulse measurements when S/N < 6. We predict that data lags obtained from structured pulse data will be different than fit lags of the same pulses. Short-duration pulses have shorter lags (of both types) than long-duration pulses (e.g., Norris et al. (2005); Peng et al. (2007); Hakkila et al. (2008)). This is not terribly surprising since a pulse s lag must be shorter than the pulse s duration since the flux distributions generally overlap (Hakkila & Preece 2011). Similarly, structured pulses must have short data lags, since structure often occurs on short timescales, and since structure is generally bright enough to contribute significantly to the CCF peak. Fit lags measured by applying the CCF to multi-channel Norris et al. (2005) pulse fits, however, are not constrained to be short, because the fits effectively smear out any energy-dependent structures by redistributing them. Since the broad CCF distributions made from these smeared light curves smoothly span the entire pulse duration, the extracted fit lags can exceed the timescales of the structures. Additionally, pronounced structure in one energy channel but not another can shift the CCF, and therefore the lag obtained from it. We calculate the data lags from the CCF for several energy channels for each pulse, and compare these to the lags of the smooth fitted monotonic pulse components (fit lags), the results of which can be found in Table 1. We have chosen to demonstrate lags between channels 3 and 1, channels 3 and 2, and channels 4 and 1 because these are commonly cited energy channels. As discussed above, we find that the greatest similarities between data lags and Norris et al. (2005) fit lags (obtained by cross-correlating the energy-dependent pulse fits) occur in the smoothest, least-structured pulses (143p2 and 5614), while the greatest differences are found in the most-structured pulses (143p1, 249, 7301p1, and 7301p2). The data lags and Norris et al. (2005) fit lags of 143p1 are relatively similar, even though this pulse shows pronounced structure. The reason for this appears to be that the rapidly-varying structure is primarily embedded within the pulse decay and does not markedly change the fit parameters of any of the energy-dependent pulse fits. We can also see if the spectral lags of bright pulses change as they are observed closer to the instrumental noise limit. Pulse data lag measurements across a variety of S/N values are shown in Figures Figure 11 shows the data lags of pulses 143p1 (left panel) and 143p2 (right panel), Figure 12 shows the data lags of pulses 249 (left panel) and 5614 (right panel), and Figure 13 shows the data lags of pulses 7301p1 (left panel) and 7301p2 (right panel). The data lags of most pulses in our sample do not change systematically when viewed at lower S/N, although the uncertainties in measuring these values do increase and some lags become unmeasurable due to the loss of signal in a particular energy channel (usually channels 1 or 4). Possible exceptions are the channel 3 vs. 1 and channel 4 vs. 1 lags of pulses 249 and 5614; these exhibit slight systematic increases not measured in their channel 3 vs. 2 lags. These increases are caused by the loss of pulse structure in noisy channel 1 but not in the other channels; without this structure the lag measures a cross-correlation between structured channel 3 or 4 emission and a now smoother channel 1 light curve. The resulting data lags are larger than the data lags of the original unreduced pulses (which are sensitive to structure) and smaller than the fit lags of the unreduced pulses (which are insensitive to structure). With these caveats, the loss of pulse structure at low S/N does not appear to introduce systematic changes in spectral lag measurements.

12 12 Hakkila et al. Table 1. Pulse Spectral Lags Pulse 31 data lag 31 fit lag 32 data lag 32 fit lag 41 data lag 41 fit lag 143p ± ± ± ± ± ± p ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± p ± ± ± ± ± ± p ± ± ± ± ± ± 0.05 Note Data lag is the time delay between emission in a high-energy channel and a low-energy channel, measured using the CCF from the time at which the maximum value of the cross-correlation function occurs. Fit lag refers to a lag obtained by cross-correlating the Norris et al. (2005) pulse fits in each energy channel, and is thus a smoothed version of the standard CCF lag p p2 Lag (seconds) 0 Lag (seconds) ms signal-to-noise 64 ms signal-to-noise Figure 11. Data lags of pulses 143p1 (left panel) and 143p2 (right panel). Black dotted lines indicate lags between BATSE channels 3 and 2, blue solid lines indicate lags between BATSE channels 3 and 1, and orange dashed lines indicate lags between BATSE channels 4 and 1. Lags in both 143p1 and 143p2 decrease systematically with decreasing S/N Lag (seconds) 0 Lag (seconds) ms signal-to-noise 64 ms signal-to-noise Figure 12. Data lags of pulses 249 (left panel) and 5614 (right panel). Black dotted lines indicate lags between BATSE channels 3 and 2, blue solid lines indicate lags between BATSE channels 3 and 1, and orange dashed lines indicate lags between BATSE channels 4 and 1. Lags in 249 remain constant with decreasing S/N while those of 5614 increase systematically with decreasing S/N.

13 Smoke and Mirrors p p2 Lag (seconds) 0 Lag (seconds) ms signal-to-noise ms signal-to-noise Figure 13. Data lags of pulses 7301p1 (left panel) and 7301p2 (right panel). Black dotted lines indicate lags between BATSE channels 3 and 2, blue solid lines indicate lags between BATSE channels 3 and 1, and orange dashed lines indicate lags between BATSE channels 4 and 1. Lags in 7301p1 do not undergo systematic changes, while lags in 7301p2 increase systematically by large amounts with decreasing S/N. Pulse lag and duration are both luminosity indicators (Hakkila et al. 2008), yet neither seems to be systematically affected by decreased S/N, although lag results can be marginally affected when one lag channel is significantly noisier than the other. Thus, measurements of these characteristics at low S/N are generally reliable, with uncertainties mostly produced by small number photon counting statistics Mirrors: Temporally Reflected Residuals Because instrumental noise has been shown to smear out GRB pulse structure, we are forced to confront the fact that our best chance to learn about GRB pulse structure is to study the spectro-temporal evolutionary properties of bright, complex GRBs. The study of complex GRB pulse characteristics has generally been avoided primarily because authors have disagreed on the number of pulses represented in bursts containing structure if the structure cannot be attributed to noise. However, as a result of our analysis in Section 3.1.1, we are able to approach each bright episode under the assumption that it represents a single structured pulse. The six pulses in our sample suggest that bright pulse residuals exhibit a wavelike function similar to that described by the residual model, but with a larger number of complete waves than the three used to develop Equation 7. The complex residual structures of these pulses can be seen in the left panels of Figures 1 6. Hakkila & Preece (2014) assumed that the wave used in the empirical fitting function did not extend outside the Bessel function s third zeros at J 0 (x = ±8.654) because the low- to moderate-s/n pulse residuals showed no evidence of this. That assumption is not true at the higher S/N values of pulses in this sample, where 143p1 appears to exhibit nine sets of up-and-down residual fluctuations, 143p2 has five, 249 has at least seven, 5614 has six, 7301p1 has roughly six, and 7391p2 has three. However, since the residual function successfully characterized these structures at lower S/N, we suspect that the high-s/n structures are extensions of this phenomenon, and that they have similar origins even though each pulse shows structures that occur on different timescales with differing amplitudes. We have unsuccessfully attempted to extend the residual model to these pulses beyond the Bessel function s third zeros. The additional wave features present at high S/N cannot be adequately fitted by simply extending the Bessel function limits and changing the t 0, a, Ω, and s values. The function incorrectly models these additional waves because the residual model was fine-tuned to work with moderate S/N GRB pulse structures; it cannot be extrapolated to the unexpected structures we encounter at high S/N. The wavelike structures at the beginning of each pulse are observed to have similar characteristics to those at the end of the pulse, except that they appear to be temporally inverted. The residual function described in Equation 7 partially explains this, as it was developed with an implicit recognition that GRB pulse residuals were best fit by temporally-inverted structures. The reason why this model worked was not understood. However, closer examination shows that the many maxima and minima of the pulse residual structures occur in a temporally-mirrored order, with near-reversible shapes. The structures at the beginning of each pulse are decent reversed facsimiles of the structures at the end of the pulse. Furthermore, the wavelike features at the beginning of pulses are compressed relative to those at the ends of pulses, apparently coordinated with the pulse s asymmetry. When an odd number of waves is found,

14 14 Hakkila et al. this indicates that the time of reflection occurs near a wave maximum. When an even number of waves is found, this indicates that the time of reflection lies at a wave minimum. Since the residual function is an incomplete representation of the residual structure, we choose to more accurately characterize the mirrored residual structure in our sample. We do this by folding the residuals over at a pulse s temporal symmetry point (the time of reflection, which we call t 0;mirror ), then stretching the folded reversed residuals by an amount s mirror until the folded wave prior to the reflection time (generally on the rising part of the pulse light curve) matches the wave following the reflection time (generally on the decaying part of the light curve). We find these values iteratively, by cycling t 0;mirror through a range of possible values, folding and reversing the residuals at each possible time of reflection, then cycling s mirror through a range of possible values so as to stretch the reversed-residuals to match the forward- and reversed stretched residual light curves. The CCF calculated at each step is used to compare the quality of the matches, with the optimal CCF value identifying the best characterization of t 0;mirror and s mirror. We note that t 0;mirror is similar to t 0 (defined in Equation 7) except that t 0;mirror is no longer necessarily measured from the time when the residual function is a maximum (t 0 is), and that s mirror is similar to s except that it is measured from t 0;mirror rather than from t 0. Our folding and stretching approach also uses all residuals rather than just those fitted by Equation 7. In other words, t 0;mirror and s mirror are less model-dependent definitions than t 0 and s, since the former require only a pulse fit to be made prior to their measurement whereas the latter require both a pulse fit and a residual fit. The results of mirroring and stretching the residuals of our six pulses are shown in Figures The times given in these figures are referenced to t 0;mirror (seventh column of Table 2) rather than to the trigger time (t = 0 seconds), and these times only describe the decay portion of the residuals (dotted line), because the rise times (solid line) have been time-inverted and stretched. Despite the bright, non-randomly distributed residual structures, the fits are good. The optimum values of the CCF exceed 0.6 for these residuals, indicating a strong correlation between the structures preceding the time of reflection and those following it. The CCFs of three of the pulses (143p2, 5614, and 7301p2) exceed 0.8. We note that the fit we chose for 143p1, with a CCF of 0.60, was slightly worse than the optimal fit of 0.64 we chose this fit because it pairs all but one of the residual peaks whereas the optimal fit misses three peaks in order to pair the two deepest troughs. The fitted characteristics of the mirrored residual parameters s mirror and t 0;mirror are compared to their counterparts s and t 0 in Table 2. The pulse asymmetry κ is included for comparison with the s value for each pulse, as these parameters have previously been shown to anti-correlate; this anti-correlation is only loosely observed from these six pulses. Also indicated are the best-fit CCF values used to obtain the time-reversible characteristics. The values are generally in good agreement, and differ primarily for the reasons described above. We estimate the significance of the mirroring by comparing the time-inverted, stretched residuals preceding the time of reflection to those following it. We use the p values obtained from a Spearman Rank-Order test to compare the two distributions, after first interpolating the time-reversed and stretched residuals so that they contain the same number of data points as the better-sampled post-t 0;mirror residuals. This approach gives us a likelihood that the residuals are related, without providing information about how closely the intensities match one another. The p values found for all pulses are significant (p < 10 8 ), indicating that residuals following t 0;mirror are time-reversed and stretched versions of those preceding it. To ensure that our fitting procedure has not inadvertently biased the results, we repeat the analysis after excluding all intensities less than 3σ away from the mean background rate (this data truncation eliminates any temporal correlations that might result from incompletely-fitted monotonic pulse components and/or backgrounds). We define p bright to be the Spearman Rank-Order correlation p value obtained from comparing these truncated data. The results of the Spearman Rank-Order test are shown in the final column of Table 2; all six pulses have p bright < The strongest matches are found for 7301p2, 249, and 143p2 (small p bright values indicate support for the stretched temporal inversion hypothesis). The least significant match is found for 5614 due to the small number of bright data points this pulse contains. We conclude that the time-reversed and stretched residuals are not an artifact of our data analysis procedure. Furthermore, we will show in Section 4 that the increased significance of p relative to p bright, when coupled with hardness evolution, indicates additional evidence for faint time-reversed and stretched structures extending farther from the time of reflection as opposed to inadequate monotonic pulse-fitting. We note that half the measured t 0;mirror values (143p1, 5614, and 7301p1) occur at times where the residuals are at a minimum rather than at a maximum. Since t 0 by definition measures reflection from a maximum value of the residuals, t 0;mirror and t 0 can differ significantly (see Table 2), as can s and s mirror. This result demonstrates that