A broad and widely used class of stationary, linear, additive time series models can have statistical properties which many authors have asserted imply that the underlying process must be non-linear, non-stationary, multiplicative, or inconsistent with shot noise. This result is demonstrated with exact and numerical evaluation of the model flux distribution function and dependence of flux standard deviation on mean flux (here and in the literature called the \emph{rms-flux relation}). These models can: (1) exhibit normal, log-normal or other flux distributions; (2) show linear or slightly non-linear rms-mean flux dependencies; as well as (3) match arbitrary second order statistics of the time series data. Accordingly the above assertions cannot be made on the basis of statistical time series analysis alone. Also discussed are ambiguities in the meaning of terms relevant to this study -- \emph{linear}, \emph{stationary} and \emph{multiplicative} -- and functions that can transform observed fluxes to a normal distribution as well or better than the logarithm.

Almost 10 yr of $\gamma$-ray observations with the Fermi Large Area Telescope (LAT) have revealed extreme $\gamma$-ray outbursts from flat spectrum radio quasars (FSRQs), temporarily making these objects the brightest $\gamma$-ray emitters in the sky. Yet, the location and mechanisms of the $\gamma$-ray emission remain elusive. We characterize long-term $\gamma$-ray variability and the brightest $\gamma$-ray flares of six FSRQs. Consecutively zooming in on the brightest flares, which we identify in an objective way through Bayesian blocks and a hill-climbing algorithm, we find variability on subhour time scales and as short as minutes for two sources in our sample (3C279, CTA102) and weak evidence for variability at time scales less than the Fermi satellite's orbit of 95 minutes for PKS1510-089 and 3C454.3. This suggests extremely compact emission regions in the jet. We do not find any signs for $\gamma$-ray absorption in the broad-line region (BLR), which indicates that $\gamma$-rays are produced at distances greater than hundreds of gravitational radii from the central black hole. This is further supported by a cross-correlation analysis between $\gamma$-ray and radio/millimeter light curves, which is consistent with $\gamma$-ray production at the same location as the millimeter core for 3C273, CTA102, and 3C454.3. The inferred locations of the $\gamma$-ray production zones are still consistent with the observed decay times of the brightest flares if the decay is caused by external Compton scattering with BLR photons. However, the minute-scale variability is challenging to explain in such scenarios.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

We apply classical machine vision and machine deep learning methods to prototype signal classifiers for the search for extraterrestrial intelligence. Our novel approach uses two-dimensional spectrograms of measured and simulated radio signals bearing the imprint of a technological origin. The studies are performed using archived narrow-band signal data captured from real-time SETI observations with the Allen Telescope Array and a set of digitally simulated signals designed to mimic real observed signals. By treating the 2D spectrogram as an image, we show that high quality parametric and non-parametric classifiers based on automated visual analysis can achieve high levels of discrimination and accuracy, as well as low false-positive rates. The (real) archived data were subjected to numerous feature-extraction algorithms based on the vertical and horizontal image moments and Huff transforms to simulate feature rotation. The most successful algorithm used a two-step process where the image was first filtered with a rotation, scale and shift-invariant affine transform followed by a simple correlation with a previously defined set of labeled prototype examples. The real data often contained multiple signals and signal ghosts, so we performed our non-parametric evaluation using a simpler and more controlled dataset produced by simulation of complex-valued voltage data with properties similar to the observed prototypes. The most successful non-parametric classifier employed a wide residual (convolutional) neural network based on pre-existing classifiers in current use for object detection in ordinary photographs. These results are relevant to a wide variety of research domains that already employ spectrogram analysis from time-domain astronomy to observations of earthquakes to animal vocalization analysis.

Self-organization is a property of dissipative nonlinear processes that are governed by an internal driver and a positive feedback mechanism, which creates regular geometric and/or temporal patterns and decreases the entropy, in contrast to random processes. Here we investigate for the first time a comprehensive number of 16 self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous {\sl order out of chaos}, during the evolution from an initially disordered system to an ordered stationary system, via quasi-periodic limit-cycle dynamics, harmonic mechanical resonances, or gyromagnetic resonances. The internal driver can be gravity, rotation, thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational instability, the Rayleigh-B\'enard convection instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or loss-cone instability. Physical models of astrophysical self-organization processes involve hydrodynamic, MHD, and N-body formulations of Lotka-Volterra equation systems.

We demonstrate the effectiveness of a relatively straightforward analysis of the complex 3D Fourier transform of galaxy coordinates derived from redshift surveys. Numerical demonstrations of this approach are carried out on a volume-limited sample of the Sloan Digital Sky Survey redshift survey. The direct unbinned transform yields a complex 3D data cube quite similar to that from the Fast Fourier Transform (FFT) of finely binned galaxy positions. In both cases deconvolution of the sampling window function yields estimates of the true transform. Simple power spectrum estimates from these transforms are roughly consistent with those using more elaborate methods. However we concentrate on the less often studied Fourier phase spectrum, a simple and general framework for characterizing non-Gaussianity, more easily interpretable than the tangled, incomplete multi-point methods conventionally used. No significant signature of non-Gaussianity has been found in the relatively small data set analyzed, but we identify some threads of modern large scale inference methodology that will presumably yield detections in new wider and deeper surveys.

We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for char- acterizing the upper limits across a large area of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, dif- ferences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914.

There have recently been several reports of apparently periodic variations in the light curves of quasars, e.g. PG 1302-102 by Graham et al. (2015a). Any quasar showing periodic oscillations in brightness would be a strong candidate to be a close binary supermassive black hole and, in turn, a candidate for gravitational wave studies. However, normal quasars -- powered by accretion onto a single, supermassive black hole -- usually show stochastic variability over a wide range of timescales. It is therefore important to carefully assess the methods for identifying periodic candidates from among a population dominated by stochastic variability. Using a Bayesian analysis of the light curve of PG 1302-102, we find that a simple stochastic process is preferred over a sinusoidal variations. We then discuss some of the problems one encounters when searching for rare, strictly periodic signals among a large number of irregularly sampled, stochastic time series, and use simulations of quasar light curves to illustrate these points. From a few thousand simulations of steep spectrum (`red noise') stochastic processes, we find many simulations that display few-cycle periodicity like that seen in PG 1302-102. We emphasise the importance of calibrating the false positive rate when the number of targets in a search is very large.

We present results on multifrequency Very Long Baseline Array (VLBA) monitoring observations of the double-image gravitationally lensed blazar JVAS B0218+357. Multi-epoch observations started less than one month after the gamma-ray flare detected in 2012 by the Large Area Telescope on board Fermi, and spanned a 2-month interval. The radio light curves did not reveal any significant flux density variability, suggesting that no clear correlation between the high energy and low-energy emission is present. This behaviour was confirmed also by the long-term Owens Valley Radio Observatory monitoring data at 15 GHz. The milliarcsecond-scale resolution provided by the VLBA observations allowed us to resolve the two images of the lensed blazar, which have a core-jet structure. No significant morphological variation is found by the analysis of the multi-epoch data, suggesting that the region responsible for the gamma-ray variability is located in the core of the AGN, which is opaque up to the highest observing frequency of 22 GHz.

PSR B1259-63/LS 2883 is a gamma-ray binary system containing a radio pulsar in a highly elliptical ~3.4-year orbit around a Be star. In its 2010 periastron passage, multiwavelength emission from radio to TeV was observed, as well as an unexpected GeV flare measured by the Fermi Large Area Telescope (LAT). Here, we report the results of LAT monitoring of PSR B1259-63 during its most recent 2014 periastron passage. We compare the gamma-ray behavior in this periastron with the former in 2010 and find that PSR B1259-63 shows a recurrent GeV flare. The similarities and differences in the phenomenology of both periastron passages are discussed.

Analysis of solar radius measurements acquired by the Michelson Doppler Imager on the SOHO spacecraft supports previously reported evidence of solar internal r-mode oscillations in Mt Wilson radius data and in nuclear-decay data acquired at the Lomonosov Moscow State University. The frequencies of these oscillations are compatible with oscillations in a putative inner tachocline that separates a slowly rotating core from the radiative envelope.

The Fermi-LAT has observed new gamma-ray flares from the blazar B0218+357 during July 2014. While no significant change in the gamma-ray spectrum has been previously observed through the flaring phase in late-2012, during this recent high activity the source displayed an exceptionally hard spectrum. The latter led to the detection of very high energy (VHE, E > 100 GeV) gamma-rays from B0218+357 by the MAGIC telescopes, establishing this source as the most distant TeV emitter known to date. In addition to the detection of VHE emission, this blazar is of particular interest since it is known to be a double-image gravitationally lensed system with a lens delay of 11.46 +/- 0.16 days measured in gamma-rays. We present the Fermi-LAT study of the July 2014 flares and discuss them in the context of previous measurements.

We report on the interannual variability of the atmospheric ice/dust cycle in the Martian polar regions for Mars Years 28-30. We used CRISM emission phase function measurements to derive atmospheric dust optical depths and data from the MARCI instrument to derive atmospheric water ice optical depths. We have used autocorrelation and cross correlation functions in order to quantify the degree to which dust and ice are correlated throughout both polar regions during Mars Years 28-29. We find that in the south polar region, dust has the tendency to "self clear", demonstrated by negative autocorrelation around the central peak. This does not occur in the north polar region. In the south polar region, dust and ice are temporally and spatially anti correlated. In the north polar region, this relationship is reversed, however temporal correlation of northern dust and ice clouds is weak - 6 times weaker than the anticorrelation in the south polar region. Our latitudinal autocorrelation functions allow us to put average spatial sizes of event cores and halos. Dust events in the south are largest, affecting almost the entire pole, whereas dust storms are smaller in the north. Ice clouds in north are similar in latitudinal extent to those in the south (both have halos < 10{\deg}). Using cross-correlation functions of water ice and dust, we find that dust events temporally lag ice events by 35-80 degrees of solar longitude in the north and south poles, which is likely due to seasonality of dust and ice events.

The major uncertainties in studies of the multi-scale structure of the Universe arise not from observational errors but from the variety of legitimate definitions and detection methods for individual structures. To facilitate the study of these methodological dependencies we have carried out 12 different analyses defining structures in various ways. This has been done in a purely geometrical way by utilizing the HOP algorithm as a unique parameter-free method of assigning groups of galaxies to local density maxima or minima. From three density estimation techniques (smoothing kernels, Bayesian Blocks and self organizing maps) applied to three data sets (the Sloan Digital Sky Survey Data Release 7, the Millennium Simulation and randomly distributed points) we tabulate information that can be used to construct catalogs of structures connected to local density maxima and minima. The resulting sizes follow continuous multi-scale distributions with no indication of the presence of a discrete hierarchy. We also introduce a novel void finder that utilizes a method to assemble Delaunay tetrahedra into connected structures and characterizes regions very nearly empty of galaxies in the source catalog.

Using data from the Fermi Large Area Telescope (LAT), we report the first clear gamma-ray measurement of a delay between flares from the gravitationally lensed images of a blazar. The delay was detected in B0218+357, a known double-image lensed system, during a period of enhanced gamma-ray activity with peak fluxes consistently observed to reach >20-50 times its previous average flux. An auto-correlation function analysis identified a delay in the gamma-ray data of 11.46 +/- 0.16 days (1 sigma) that is ~1 day greater than previous radio measurements. Considering that it is beyond the capabilities of the LAT to spatially resolve the two images, we nevertheless decomposed individual sequences of superposing gamma-ray flares/delayed emissions. In three such ~8-10 day-long sequences within a ~4-month span, considering confusion due to overlapping flaring emission and flux measurement uncertainties, we found flux ratios consistent with ~1, thus systematically smaller than those from radio observations. During the first, best-defined flare, the delayed emission was detailed with a Fermi pointing, and we observed flux doubling timescales of ~3-6 hrs implying as well extremely compact gamma-ray emitting regions.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

This presentation describes the Bayesian Block algorithm in the context of its application to analysis of time series data from the Fermi Gamma Ray Space Telescope. More generally this algorithm performs optimal segmentation analysis on sequential data in any mode, with arbitrary sampling and in the presence of gaps and exposure variations. A new procedure for correcting for backgrounds is also described.

Claims by the DAMA (DArk MAtter) collaboration to have detected an annually varying signal consistent with models of dark matter appear to be at variance with results from other dark-matter searches. To further understand the DAMA results, we have carried out an independent analysis of DAMA data reconstructed from published figures. In addition to reexamining the Lomb-Scargle and chi-square analyses previously carried out by the DAMA collaboration, we carry out two new likelihood analyses and a new chi-square analysis, focusing attention on the treatment of experimental errors and binning. We confirm the existence of an annual oscillation, with a maximum in early June, but at a lower significance level than previously reported.

In this document, we describe the scientific potential of blazar observations with a X-ray polarimetry mission like GEMS (Gravity and Extreme Magnetism SMEX). We describe five blazar science investigations that such a mission would enable: (i) the structure and the role of magnetic fields in AGN jets, (ii) analysis of the polarization of the synchrotron X-ray emission from AGN jets, (iii) discrimination between synchrotron self-Compton and external Compton models for blazars with inverse Compton emission in the X-ray band, (iv) a precision study of the polarization properties of the X-ray emission from Cen-A, (v) tests of Lorentz Invariance based on X-ray polarimetric observations of blazars. We conclude with a discussion of a straw man observation program and recommended accompanying multiwavelength observations.

Analysis of over 36 years of time series data from the NSO/AFRL/Sac Peak K-line monitoring program elucidates five components of the variation of the seven measured chromospheric parameters: (a) the solar cycle (period ~ 11 years), (b) quasi-periodic variations (periods ~100 days), (c) a broad band stochastic process (wide range of periods), (d) rotational modulation, and (e) random observational errors, independent of (a)-(d). Correlation and power spectrum analyses elucidate periodic and aperiodic variation of these parameters. Time-frequency analysis illuminates periodic and quasi periodic signals, details of frequency modulation due to differential rotation, and in particular elucidates the rather complex harmonic structure (a) and (b) at time scales in the range ~0.1 - 10 years. These results using only full-disk data suggest that similar analyses will be useful at detecting and characterizing differential rotation in stars from stellar light-curves such as thosebeing produced by NASA's Kepler observatory. Component (c) consists of variations over a range of timescales, in the manner of a 1/f random process with a power-law slope index that varies in a systematic way. A time-dependent Wilson-Bappu effect appears to be present in the solar cycle variations (a), but not in the more rapid variations of the stochastic process (c). Component (d) characterizes differential rotation of the active regions. Component (e) is of course not characteristic of solar variability, but the fact that the observational errors are quite small greatly facilitates the analysis of the other components. The data analyzed in this paper can be found at the National Solar Observatory web site http://nsosp.nso.edu/cak_mon/, or by file transfer protocol at ftp://ftp.nso.edu/idl/cak.parameters

We present the first Kepler monitoring of a strongly variable BL Lac, W2R1926+42. The light curve covers 181 days with ~0.2% errors, 30 minute sampling and >90% duty cycle, showing numerous delta I/I > 25% flares over timescales as short as a day. The flux distribution is highly skewed and non-Gaussian. The variability shows a strong rms-flux correlation with the clearest evidence to date for non-linearity in this relation. We introduce a method to measure periodograms from the discrete autocorrelation function, an approach that may be well-suited to a wide range of Kepler data. The periodogram is not consistent with a simple power-law, but shows a flattening at frequencies below 7x10-5 Hz. Simple models of the power spectrum, such as a broken power law, do not produce acceptable fits, indicating that the Kepler blazar light curve requires more sophisticated mathematical and physical descriptions than currently in use.

We examine the expected X-ray polarization properties of neutron-star X-ray sources of various types, e.g., accretion and rotation powered pulsars, magnetars, and low-mass X-ray binaries. We summarize the model calculations leading to these expected properties. We describe how a comparison of these with their observed properties, as inferred from GEMS data, will probe the essential dynamical, electromagnetic, plasma, and emission processes in neutron-star binaries, discriminate between models of these processes, and constrain model parameters. An exciting goal is the first observational demonstration in this context of the existence of vacuum resonance, a fundamental quantum electrodynamical phenomenon first described in the 1930s.

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

We present here a summary of the scientific goals behind the Gravity and Extreme Magnetism SMEX (GEMS) X-ray polarimetry mission's black hole (BH) observing program. The primary targets can be divided into two classes: stellar-mass galactic BHs in accreting binaries, and super-massive BHs in the centers of active galactic nuclei (AGN). The stellar-mass BHs can in turn be divided into various X-ray spectral states: thermal-dominant (disk), hard (radio jet), and steep power-law (hot corona). These different spectral states are thought to be generated by different accretion geometries and emission mechanisms. X-ray polarization is an ideal tool for probing the geometry around these BHs and revealing the specific properties of the accreting gas.

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it - an improved and generalized version of Bayesian Blocks (Scargle 1998) - that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piecewise linear and piecewise exponential representations, multi-variate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by (Arias-Castro, Donoho and Huo 2003). In the spirit of Reproducible Research (Donoho et al. 2008) all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical wavelengths has a size of ~11 ly across, doubled its gamma-ray flux within eight hours. The peak photon flux was (186 +- 6) 10-7 cm-2 s-1 above 100 MeV, which corresponds to a 30-fold increase compared to the average value. During the flare, a new component emerged in the spectral energy distribution, which peaked at an energy of (375 +- 26) MeV at flare maximum. The observations imply that the emission region was likely relativistically beamed toward us and that variations in its motion are responsible for the observed spectral variability.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

We present our temporal and spectral analyses of 29 bursts from SGR J0501+4516, detected with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 13 days of the source activation in 2008 (August 22 to September 3). We find that the T90 durations of the bursts can be fit with a log-normal distribution with a mean value of ~ 123 ms. We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in photon space (i.e., using their deconvolved spectra) and find that these are very similar to the T90s estimated in count space (following a log-normal distribution with a mean value of ~ 124 ms). We fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest bursts with several models. We find that a single power law with an exponential cutoff model fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We expand on the physical interpretation of these two models and we compare their parameters and discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that Epeak decreases with energy flux (and fluence) to a minimum of ~30 keV at F=8.7e-6 erg/cm2/s, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs J1550-5418 and 1806-20. The isotropic luminosity corresponding to these flux values is roughly similar for all sources (0.4-1.5 e40 erg/s).

We analyze the Swift/BAT sample of short gamma-ray bursts, using an objective Bayesian Block procedure to extract temporal descriptors of the bursts' initial pulse complexes (IPCs). The sample comprises 12 and 41 bursts with and without extended emission (EE) components, respectively. IPCs of non-EE bursts are dominated by single pulse structures, while EE bursts tend to have two or more pulse structures. The medians of characteristic timescales - durations, pulse structure widths, and peak intervals - for EE bursts are factors of ~ 2-3 longer than for non-EE bursts. A trend previously reported by Hakkila and colleagues unifying long and short bursts - the anti-correlation of pulse intensity and width - continues in the two short burst groups, with non-EE bursts extending to more intense, narrower pulses. In addition we find that preceding and succeeding pulse intensities are anti-correlated with pulse interval. We also examine the short burst X-ray afterglows as observed by the Swift/XRT. The median flux of the initial XRT detections for EE bursts (~ 6 x 10^-10 erg cm^-2 s^-1) is ~> 20 x brighter than for non-EE bursts, and the median X-ray afterglow duration for EE bursts (~ 60,000 s) is ~ 30 x longer than for non-EE bursts. The tendency for EE bursts toward longer prompt-emission timescales and higher initial X-ray afterglow fluxes implies larger energy injections powering the afterglows. The longer-lasting X-ray afterglows of EE bursts may suggest that a significant fraction explode into more dense environments than non-EE bursts, or that the sometimes-dominant EE component efficiently powers the afterglow. Combined, these results favor different progenitors for EE and non-EE short bursts.

Three methods for detecting and characterizing structure in point data, such as that generated by redshift surveys, are described: classification using self-organizing maps, segmentation using Bayesian blocks, and density estimation using adaptive kernels. The first two methods are new, and allow detection and characterization of structures of arbitrary shape and at a wide range of spatial scales. These methods should elucidate not only clusters, but also the more distributed, wide-ranging filaments and sheets, and further allow the possibility of detecting and characterizing an even broader class of shapes. The methods are demonstrated and compared in application to three data sets: a carefully selected volume-limited sample from the Sloan Digital Sky Survey redshift data, a similarly selected sample from the Millennium Simulation, and a set of points independently drawn from a uniform probability distribution -- a so-called Poisson distribution. We demonstrate a few of the many ways in which these methods elucidate large scale structure in the distribution of galaxies in the nearby Universe.

A key goal of radio and $\gamma-$ray observations of active galactic nuclei is to characterize their time variability in order to elucidate physical processes responsible for the radiation. I describe algorithms for relevant time series analysis tools -- correlation functions, Fourier and wavelet amplitude and phase spectra, structure functions, and time-frequency distributions, all for arbitrary data modes and sampling schemes. For example radio measurements can be cross-analyzed with data streams consisting of time-tagged gamma-ray photons. Underlying these methods is the Bayesian block scheme, useful in its own right to characterize local structure in the light curves, and also prepare raw data for input to the other analysis algorithms. One goal of this presentation is to stimulate discussion of these methods during the workshop.

The term "false-alarm probability" denotes the probability that at least one out of M independent power values in a prescribed search band of a power spectrum computed from a white-noise time series is expected to be as large as or larger than a given value. The usual formula is based on the assumption that powers are distributed exponentially, as one expects for power measurements of normally distributed random noise. However, in practice one typically examines peaks in an over-sampled power spectrum. It is therefore more appropriate to compare the strength of a particular peak with the distribution of peaks in over-sampled power spectra derived from normally distributed random noise. We show that this leads to a formula for the false-alarm probability that is more conservative than the familiar formula. We also show how to combine these results with a Bayesian method for estimating the probability of the null hypothesis (that there is no oscillation in the time series), and we discuss as an example the application of these procedures to Super-Kamiokande solar neutrino data.

The initial pulse complex (IPC) in short gamma-ray bursts is sometimes accompanied by a softer, low-intensity extended emission (EE) component. In cases where such a component is not observed, it is not clear if it is present but below the detection threshold. Using Bayesian Block (BB) methods, we measure the EE component and show that it is present in one quarter of a Swift/BAT sample of 51 short bursts, as was found for the Compton/BATSE sample. We simulate bursts with EE to calibrate the BAT threshold for EE detection and show that this component would have been detected in nearly half of BAT short bursts if it were present, to intensities ~ 10^-2 counts cm^-2 s^-1, a factor of five lower than actually observed in short bursts. In the BAT sample the ratio of average EE intensity to IPC peak intensity, Rint, ranges over a factor of 25, Rint ~ 3 x 10^-3 to 8 x 10^-2. In comparison, for the average of the 39 bursts without an EE component, the 2-sigma upper limit is Rint < 8 x 10^-4. These results suggest that a physical threshold effect operates near Rint ~ few x 10^-3, below which the EE component is not manifest.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

A very small dispersion in the speed of light may be observable in Fermi time- and energy-tagged data on variable sources, such as gamma-ray bursts (GRB) and active galactic nuclei (AGN). We describe a method to compute the size of this effect by applying the Feynman sum-over-histories formalism for relativistic quantum electrodynamics to a discrete model of space-time called \emph{causal set theory}.

The well-studied VHE (E >100 GeV) blazar Mrk 501 was observed between March and May 2008 as part of an extensive multiwavelength observation campaign including radio, optical, X-ray and VHE gamma-ray instruments. Mrk 501 was in a low state of activity during the campaign, with a low VHE flux of about 20% the Crab Nebula flux. Nevertheless, significant flux variations could be observed in X-rays as well as $\gamma$-rays. Overall Mrk 501 showed increased variability when going from radio to gamma-ray energies. The broadband spectral energy distribution during the two different emission states of the campaign was well described by a homogeneous one-zone synchrotron self-Compton model. The high emission state was satisfactorily modeled by increasing the amount of high energy electrons with respect to the low emission state. This parameterization is consistent with the energy-dependent variability trend observed during the campaign.

The LAT instrument on the Fermi mission will reveal the rich spectral and temporal gamma-ray burst phenomena in the > 100 MeV band. The synergy with Fermi's GBM detectors will link these observations to those in the well explored 10-1000 keV range; the addition of the > 100 MeV band observations will resolve theoretical uncertainties about burst emission in both the prompt and afterglow phases. Trigger algorithms will be applied to the LAT data both onboard the spacecraft and on the ground. The sensitivity of these triggers will differ because of the available computing resources onboard and on the ground. Here we present the LAT's burst detection methodologies and the instrument's GRB capabilities.

The usual procedure for estimating the significance of a peak in a power spectrum is to calculate the probability of obtaining that value or a larger value by chance (known as the "p-value"), on the assumption that the time series contains only noise - typically that the measurements are derived from random samplings of a Gaussian distribution. We really need to know the probability that the time series is - or is not - compatible with the null hypothesis that the measurements are derived from noise. This probability can be calculated by Bayesian analysis, but this requires one to specify and evaluate a second hypothesis, that the time series does contain a contribution other than noise. We approach the problem of identifying this function in two ways. We first propose three simple conditions that it seems reasonable to impose on this function, and show that these conditions may be satisfied by a simple function with one free parameter. We then define two different ways of combining information derived from two independent power estimates. We find that this consistency condition may be satisfied, to good approximation, by a special case of the previously proposed likelihood function. We find that the resulting significance estimates are considerably more conservative than those usually associated with the p-values. As two examples, we apply the new procedure to two recent analyses of solar neutrino data: (a) power spectrum analysis of Super-Kamiokande data, and (b) the combined analysis of radiochemical neutrino data and irradiance data.

DisCan is a new algorithm implementing photon dispersion cancellation in order to measure energy-dependent delays in variable sources. This method finds the amount of reversed dispersion that optimally cancels any actual dispersion present. It applies to any time- and energy-tagged photon data, and can avoid binning in both time and energy. The primary motivation here is the search for quantum gravity based dispersion in future gamma ray burst data from the Gamma Ray Large Area Space Telescope (GLAST). Extrapolation of what is know about bursts at lower energies yields a reasonable prospect that photon dispersion effects consistent with some quantum gravity formalisms may be detected in sufficiently bright bursts. Short bursts have no or very small inherent lags, and are therefore better prospects than long ones, but even they suffer systematic error due to pulse asymmetry that may yield an irreducible uncertainty. We note that data at energies higher than about 0.1 TeV may not be useful for detecting dispersion in GRBs. Of several variants of the proposed algorithm, one based on Shannon information is consistently somewhat superior to all of the others we investigated.

The Compton Gamma Ray Observatory (CGRO) observed many types of data and one of them is the time-tagged photon events (TTE data). We use the Bayesian block analysis, using Bayesian statistics, analyses the TTE data. Our results; calculations of duration (T100), count rates (burst photon numbers in different channels) and count peaks (in 64, 16 and 4 ms). We present the duration, the peak duration and the distance between peaks distributions. Principal Component Analysis (PCA) has been also applied. The PCA shows interesting results, such as channel 4 (highest energy channel) probably is very important.

Currently, the best available probe of the early phase of gamma-ray burst (GRB) jet attributes is the prompt gamma-ray emission, in which several intrinsic and extrinsic variables determine GRB pulse evolution. Bright, usually complex bursts have many narrow pulses that are difficult to model due to overlap. However, the relatively simple, long spectral lag, wide-pulse bursts may have simpler physics and are easier to model. In this work we analyze the temporal and spectral behavior of wide pulses in 24 long-lag bursts, using a pulse model with two shape parameters -- width and asymmetry -- and the Band spectral model with three shape parameters. We find that pulses in long-lag bursts are distinguished both temporally and spectrally from those in bright bursts: the pulses in long spectral lag bursts are few in number, and ~ 100 times wider (10s of seconds), have systematically lower peaks in nu*F(nu), harder low-energy spectra and softer high-energy spectra. We find that these five pulse descriptors are essentially uncorrelated for our long-lag sample, suggesting that at least five parameters are needed to model burst temporal and spectral behavior. However, pulse width is strongly correlated with spectral lag; hence these two parameters may be viewed as mutual surrogates. We infer that accurate formulations for estimating GRB luminosity and total energy will depend on several gamma-ray attributes, at least for long-lag bursts. The prevalence of long-lag bursts near the BATSE trigger threshold, their predominantly low nF(n) spectral peaks, and relatively steep upper power-law spectral indices indicate that Swift will detect many such bursts.

In solar physics, especially in exploratory stages of research, it is often necessary to compare the power spectra of two or more time series. One may, for instance, wish to estimate what the power spectrum of the combined data sets might have been, or one may wish to estimate the significance of a particular peak that shows up in two or more power spectra. One may also on occasion need to search for a complex of peaks in a single power spectrum, such as a fundamental and one or more harmonics, or a fundamental plus sidebands, etc. Visual inspection can be revealing, but it can also be misleading. This leads one to look for one or more ways of forming statistics, which readily lend themselves to significance estimation, from two or more power spectra. We derive formulas for statistics formed from the sum, the minimum, and the product of two or more power spectra. A distinguishing feature of our formulae is that, if each power spectrum has an exponential distribution, each statistic also has an exponential distribution.

Many signal processing problems can be solved by maximizing the fitness of a segmented model over all possible partitions of the data interval. This letter describes a simple but powerful algorithm that searches the exponentially large space of partitions of $N$ data points in time $O(N^2)$. The algorithm is guaranteed to find the exact global optimum, automatically determines the model order (the number of segments), has a convenient real-time mode, can be extended to higher dimensional data spaces, and solves a surprising variety of problems in signal detection and characterization, density estimation, cluster analysis and classification.

Narayan and Heyl (2002) have developed a theoretical framework to convert suitable upper limits on type I X-ray bursts from accreting black hole candidates (BHCs) into evidence for an event horizon. However, no appropriate observational limit exists in the literature. In this paper we survey 2101.2 ks of data from the Unconventional Stellar Aspect (USA) X-ray timing experiment and 5142 ks of data from the Rossi X-ray Timing Explorer (RXTE) experiment to obtain a formal constraint of this type. 1122 ks of neutron star data yield a population averaged mean burst rate of 1.69 x 10-5 bursts s-1 while 6081 ks of BHC data yield a 95% confidence level upper limit of 4.9 x 10-7 bursts s-1. This is the first published limit of this type for Black Hole Candidates. Applying the theoretical framework of Narayan and Heyl (2002) we calculate regions of unstable luminosity where the neutron stars are expected to burst and the BHCs would be expected to burst if they had a surface. In this unstable luminosity region 464 ks of neutron star data yield an averaged mean burst rate of 4.1 x 10-5 bursts s-1 and 1512 ks of BHC data yield a 95% confidence level upper limit of 2.0 x 10-6 bursts s-1, and a limit of > 10\sigma that BHCs do not burst with a rate similar to the rate of neutron stars in these unstable regions. This gives further evidence that BHCs do not have surfaces unless there is some new physics occurring on their surface.

We present an analysis of long term X-ray monitoring observations of Circinus X-1 (Cir X-1) made with four different instruments: Vela 5B, Ariel V ASM, Ginga ASM, and RXTE ASM, over the course of more than 30 years. We use Lomb-Scargle periodograms to search for the ~16.5 day orbital period of Cir X-1 in each of these data sets and from this derive a new orbital ephemeris based solely on X-ray measurements, which we compare to the previous ephemerides obtained from radio observations. We also use the Phase Dispersion Minimization (PDM) technique, as well as FFT analysis, to verify the periods obtained from periodograms. we obtain dynamic periodograms (both Lomb-Scargle and PDM) of Cir X-1 during the RXTE era, showing the period evolution of Cir X-1, and also displaying some unexplained discrete jumps in the location of the peak power.

We present preliminary results on Gamma Ray Burst (GRB) triggers with the Gamma-ray Large Area Space Telescope (GLAST). After a brief summary of the detector layout, GLAST expected performances on GRB detection are recalled. Status report on the simulation software and preliminary triggers studies are then reported, already showing significant improvement on EGRET results.

This paper reports X-ray spectral observations of a relatively nearby z = 0.048 BL Lacertae (BL Lac) object 1ES1959+650, which is a potential TeV emitter. The observations include 31 short pointings made by the Unconventional Stellar Aspect (USA) Experiment on board the Advanced Research and Global Observation Satellite (ARGOS), and 17 pointings by the PCA on board the Rossi X-ray Timing Explorer (RXTE). Most of these observations were spaced by less than 1 day. \es was detected by the ARGOS USA detector in the range 1-16 keV, and by the PCA in the 2-16 keV range but at different times. During the closely spaced RXTE observations beginning on 2000 July 28, an ending of one flare and a start of another are visible, associated with spectral changes, where the photon index Gamma ranges between ~ 1.4 and 1.7, and the spectrum is harder when the source is brighter. This implies that 1ES1959 is an XBL-type blazar, with the X-ray emission likely to originate via the synchrotron process. The USA observations reveal another flare that peaked on 2000 November 14 and doubled the flux within a few days, again associated with spectral changes of the same form. The spectral variability correlated with the flux and timing characteristics of this object that are similar to those of other nearby BL Lacs, and suggest relativistic beaming with a Doppler factor delta > 1.6 and magnetic fields of the order of a few mG. We also suggest that the steady component of the X-ray emission -- present in this object as well as in other XBLs -- may be due to the large-scale relativistic jet (such as measured by Chandra in many radio-loud AGN), but pointing very closely to our line of sight.

We report on timing and spectral observations of the 2000 outburst of XTE J1550-564 made by the Unconventional Stellar Aspect (USA) Experiment on board the Advanced Research and Global Observation Satellite (ARGOS). We observe a low-frequency quasi-periodic oscillation (LFQPO) with a centroid frequency that tends to increase with increasing flux and a fractional rms amplitude which is correlated with the hardness ratio. The evolution of the hardness ratio (4--16 keV/1--4 keV) with time and source flux is examined. The hardness-intensity diagram (HID) shows a cyclical movement in the clockwise direction and possibly indicates the presence of two independent accretion flows. We observe a relationship between the USA 4--16 keV count rate and radio observations and discuss this in the context of previously observed correlations between X-ray, radio, optical and IR data. We examine our results in the context of models invoking two accretion flows: a thin disk and a hot sub-Keplerian flow.

We analyze BATSE time-tagged event (TTE) data for short gamma-ray bursts (T90 duration < 2.6 s), studying spectral lag vs. peak flux and duration, as well as the number of distinct pulse structures per burst. Performing the cross-correlation between two energy bands, we measure an average lag ~ 20-40 x shorter than for long bursts, and a lag distribution close to symmetric about zero - unlike long bursts. Using a "Bayesian Block" method to identify significantly distinct pulse peaks, we find an order of magnitude fewer pulses than found in studies of long bursts. The disparity in lag magnitude is discontinuous across the ~ 2-s valley between long and short bursts. Thus, short bursts do not appear to be representable as a continuation of long bursts' temporal characteristics.

We analyze BATSE 64-ms data for long gamma-ray bursts (T90 > 2.6 s), exploring the relationships between spectral lag, duration, and the number of distinct pulses per burst, Npulses. We measure durations using a brightness-independent technique. Within a similarly brightness-independent framework, we use a "Bayesian Block" method to find significant valleys and peaks, and thereby identify distinct pulses in bursts. Our results show that, across large dynamic ranges in peak flux and the Npulses measure, bursts have short lags and narrow pulses, while bursts with long lags tend to have just a few significant, wide pulses. There is a tendency for harder bursts to have few pulses, and these hardest bursts have short lags. Spectral lag and duration appear to be nearly independent - even for those bursts with relatively long lags: wider pulses tend to make up in duration for fewer pulses. Our brightness-independent analysis adds to the nascent picture of an intimate connection between pulse width, spectral lag, and peak luminosity for bursts with known redshifts: We infer that lower-luminosity bursts should have fewer episodes of organized emission - wider pulses with longer spectral lags.

We have studied the time evolution of the power spectrum of XTE J1550-564, using X-ray luminosity time series data obtained by the Rossi X-Ray Timing Explorer satellite. A number of important practical fundamental issues arise in the analysis of these data, including dealing with time-tagged event data, removal of noise from a highly non-stationary signal, and comparison of different time-frequency distributions. We present two new methods to understand the time frequency variations, and compare them to the dynamic power spectrum of Homan et al. All of the approaches provide evidence that the QPO frequency varies in a systematic way during the time evolution of the signal.

If the solar neutrino flux were constant, as is widely assumed, the histogram of flux measurements would be unimodal. On the other hand, sinusoidal or square-wave modulation (either periodic or stochastic) may lead to a bimodal histogram. We here present evidence that the neutrino flux histogram is in fact bimodal. We analyze all available data from gallium experiments, coordinating results from the GALLEX and GNO experiments into one data set, and adopting results from the SAGE experiment as another data set. The two histograms, from the two data sets, are consistent in showing peaks in the range 45-75 SNU and 90-120 SNU, with a valley in between. By combining the data into one data set, we may form more detailed histograms; these strengthen the case that the flux is bimodal. A preliminary statistical analysis indicates that the bimodal character of the solar neutrino flux is highly significant. Since the upper peak is close to the expected flux (120-140 SNU), we may infer that the neutrino deficit is due to time-varying attenuation of the flux produced in the core. We estimate the time scale of this variation to be in the range 10-60 days. Attenuation that varies on such a time scale is suggestive of the influence of solar rotation, and points towards a process involving the solar magnetic field in conjunction with a nonzero neutrino magnetic moment.

The USA experiment on ARGOS and RXTE have exensively observed the X-ray transient XTE J1118+480 during its recent outburst in 2000 April--June. We present detailed monitoring of the evolution of a low frequency QPO which drifts from 0.07 Hz to 0.15 Hz during the outburst. We examine possible correlations of the QPO frequency with the flux and spectral characteristics of the source, and compare this QPO to low frequency QPOs observed in other black hole candidates.

The fast temporal structures and cosmological distances of gamma-ray bursts (GRBs) afford a natural laboratory for testing theories of frequency-dependent propagation of high-energy photons, as predicted for quantum gravity (QG). We calibrate the sensitivity of the proposed Gamma-ray Large Area Space Telescope (GLAST) by performing simulations which include: the response of GLAST to a GRB fluence distribution; a distribution of spectral power-law indices similar to the EGRET sample; and consideration of gamma-gamma attenuation, significant above ~ 10 GeV for redshifts z > 3 - 5. We find that GLAST should detect > 200 GRBs per year, with sensitivity to a few tens of GeV for a few bursts. GLAST could detect the energy- and distance-dependent dispersion (10 ms / GeV / Gpc) predicted by QG with 1 - 2 years of observations. Attribution to QG would require correlation of GRB redshifts with the temporal and energetic signatures.

The USA Experiment is a new X-ray timing experiment with large collecting area and microsecond time resolution capable of conducting a broad program of studies of galactic X-ray binaries. USA is one of nine experiments aboard the Advanced Research and Global Observation Satellite which was launched February 23, 1999. USA is a collimated proportional counter X-ray telescope with about 1000 cm$^2$ of effective area per detector with two detectors sensitive to photons in the energy range 1-15 keV. A unique feature of USA is that photon events are time tagged by reference to an onboard GPS receiver allowing precise absolute time and location determination. We will present an overview of the USA instrument, capabilities, and scientific observing plan as well as the current status of the instrument.

Recent analysis of the Homestake data indicates that the solar neutrino flux contains a periodic variation that may be attributed to rotational modulation occurring deep in the solar interior, either in the tachocline or in the radiative zone. This paper presents an analysis of GALLEX data that yields supporting evidence of this rotational modulation at the 0.1% significance level. The depth of modulation inferred from the rotational signature is large enough to explain the neutrino deficit. The Rieger 157-day periodicity, first discovered in solar gamma-ray flares, is present also in Homestake data. A related oscillation with period 52 days is found in the GALLEX data. The relationship of these periods to the rotational period inferred from neutrino data suggests that they are due to r-mode oscillations.

We have conducted a search for aperiodic millisecond variability in the integrated 1 to 25 keV X-ray region of Cyg X-1. We have examined HEAO A-1 archival data and Rossi X-ray Timing Explorer (RXTE) guest observer data for evidence of excess power above the Poisson noise floor using the relative integral power analysis and the Fourier transform method. Our results are in disagreement with the results of Meekins et al. (1984). We attribute the discrepancy to an instrumental effect for which Meekins et al.(1984) did not apply a correction. With the correction we see no evidence for excess power above 25 Hz in the HEAO A-1 data. Our analysis of RXTE data is in agreement with previously published results of different data sets and shows no sign of excess power above 30 Hz.

I describe a new time-domain algorithm for detecting localized structures (bursts), revealing pulse shapes, and generally characterizing intensity variations. The input is raw counting data, in any of three forms: time-tagged photon events (TTE), binned counts, or time-to-spill (TTS) data. The output is the most likely segmentation of the observation into time intervals during which the photon arrival rate is perceptibly constant -- i.e. has a fixed intensity without statistically significant variations. Since the analysis is based on Bayesian statistics, I call the resulting structures Bayesian Blocks. Unlike most, this method does not stipulate time bins -- instead the data themselves determine a piecewise constant representation. Therefore the analysis procedure itself does not impose a lower limit to the time scale on which variability can be detected. Locations, amplitudes, and rise and decay times of pulses within a time series can be estimated, independent of any pulse-shape model -- but only if they do not overlap too much, as deconvolution is not incorporated. The Bayesian Blocks method is demonstrated by analyzing pulse structure in BATSE $\gamma$-ray data. The MatLab scripts and sample data can be found on the WWW at: http://george.arc.nasa.gov/~scargle/papers.html

Bayesian Blocks is a new time series algorithm for detecting localized structures (spikes or shots), revealing pulse shapes, and generally characterizing intensity variations. It maps raw counting data into a maximum likelihood piecewise constant representation of the underlying signal. This bin-free method imposes no lower limit on measurable time scales. Applied to BATSE TTE data, it reveals the shortest know burst structure -- a spike superimposed on the main burst in GRB 920229 = Trigger 1453, with rise and decay timescales ~ few 100 microseconds.

We investigate whether a simple cosmology can fit GRB results in both time dilation and Log N - Log P simultaneously. Simplifying assumptions include: all GRBs are spectrally identical to BATSE trigger 143, $\Omega=1$ universe, and no luminosity and number density evolution. Observational data used includes: the BATSE 3B peak brightness distribution (64-ms time scale), the Pioneer Venus Orbiter (PVO) brightness distribution, and the Norris et al. time dilation results for peak aligned profiles presented at this meeting. We find acceptable cosmological fits to the brightness distributions when placing BATSE trigger 143 at a redshift of 0.15 $\pm$ 0.10. This translates into a $(1 + z_{dim}) / (1 + z_{bright})$ factor of about 1.50 $\pm$ 0.50 between selected brightness extremes of the Norris et al. sample. Norris et al. estimate, however, that $(1 + z_{dim}) / (1 + z_{bright})$ $\approx$ 2.0 $\pm$ 0.5 when considering duration tests. The difference is marginal and could be accounted for by evolution. We therefore find that evolution of GRBs is preferred but not demanded.

Two tests for cosmological time dilation in $\gamma$-ray bursts -- the peak alignment and auto-correlation statistics -- involve averaging information near the times of peak intensity. Both tests require width corrections, assuming cosmological origin for bursts, since narrower temporal structure from higher energy would be redshifted into the band of observation, and since intervals between pulse structures are included in the averaging procedures. We analyze long ($>$ 2 s) BATSE bursts and estimate total width corrections for trial time-dilation factors (TDF = [1+$z_{\rm dim}$]/[1+$z_{\rm brt}$]) by time-dilating and redshifting bright bursts. Both tests reveal significant trends of increasing TDF with decreasing peak flux, but neither provides sufficient discriminatory power to distinguish between actual TDFs in the range 2--3.

If $\gamma$-ray bursts are at cosmological distances, then not only their constituent pulses but also the intervals between pulses should be time-dilated. Unlike time-dilation measures of pulse emission, intervals would appear to require negligible correction for redshift of narrower temporal structure from higher energy into the band of observation. However, stretching of pulse intervals is inherently difficult to measure without incurring a timescale-dependent bias since, as time profiles are stretched, more structure can appear near the limit of resolution. This problem is compounded in dimmer bursts because identification of significant structures becomes more problematic. We attempt to minimize brightness bias by equalizing signal-to-noise (s/n) level of all bursts. We analyze wavelet-denoised burst profiles binned to several resolutions, identifying significant fluctuations between pulse structures and interjacent valleys. When bursts are ranked by peak flux, an interval time-dilation signature is evident, but its magnitude and significance are dependent upon temporal resolution and s/n level.

We have measured the T90 and T50 durations of bright and dim gamma-ray bursts (GRBs) detected by the Compton Gamma Ray Observatory's Burst and Transient Source Experiment (BATSE). The T90 (T50) duration is defined as the interval over which 5\% (25\%) to 95\% (75\%) of the burst counts accumulate. Out of 775 bursts observed by BATSE, 159 bursts were analyzed; bursts with durations shorter than 1.5 s were excluded. A Kolmogorov-Smirnov test yields a probability of 6 x $10^{-5}$ that the T50 durations of the dim and bright samples are drawn from the same parent population. We find that the centroid and extent of the duration distribution for the dim sample are scaled by approximately a factor of two relative to those of the bright sample. The measured time dilation factor is not sensitive to choice of energy band. These results are quantitatively consistent with previous tests for time dilation in a smaller sample of BATSE bursts. The dimmer bursts, if cosmological, would lie at redshifts of order two.

If gamma-ray bursters are at cosmological distances - as suggested by their isotropic distribution on the sky and by their number-intensity relation - then the burst profiles will be stretched in time, by an amount proportional to the redshift, 1 + $z$. We have tested data from the {\it Compton} Gamma Ray Observatory's Burst and Transient Source Experiment (BATSE) for such time dilation. Our measures of time scale are constructed to avoid selection effects arising from intensity differences by rescale all bursts to fiducial levels of peak intensity and noise bias. The three tests involved total count rate above background, wavelet decomposition, and alignment of the highest peaks. In all three tests, the dim bursts are stretched by a factor of about two relative to the bright ones, over seven octaves of time scale. We calibrated the measurements by dilating synthetic bursts that approximate the temporal characteristics of bright BATSE bursts. Results are consistent with bursts of BATSE's peak-flux completeness limit being at cosmological distances corresponding to $z \sim 1$, and thus with independent cosmological interpretations of the BATSE number-intensity relation.