We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

We present a classification of galaxies in the Pan-STARRS1 (PS1) 3$\pi$ survey based on their recent star formation history and morphology. Specifically, we train and test two Random Forest (RF) classifiers using photometric features (colors and moments) from the PS1 data release 2. Labels for the morphological classification are taken from Huertas-Company+2011, while labels for the star formation fraction (SFF) are from the Blanton+2005 catalog. We find that colors provide more predictive accuracy than photometric moments. We morphologically classify galaxies as either early- or late-type, and our RF model achieves a 78\% classification accuracy. Our second model classifies galaxies as having either a low-to-moderate or high SFF. This model achieves an 89\% classification accuracy. We apply both RF classifiers to the entire PS1 $3\pi$ dataset, allowing us to assign two scores to each PS1 source: $P_\mathrm{HSFF}$, which quantifies the probability of having a high SFF, and $P_\mathrm{spiral}$, which quantifies the probability of having a late-type morphology. Finally, as a proof of concept, we apply our classification framework to supernova (SN) host-galaxies from the Zwicky Transient Factory and the Lick Observatory Supernova Search samples. We show that by selecting on $P_\mathrm{HSFF}$ or $P_\mathrm{spiral}$ it is possible to significantly enhance or suppress the fraction of core-collapse SNe (or thermonuclear SNe) in the sample with respect to random guessing. This result demonstrates how contextual information can aid transient classifications at the time of first detection. In the current era of spectroscopically-starved time-domain astronomy, prompt automated classification is paramount.

Detections of SiO masers from the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey more tightly define the region where Oxygen-rich (O) Asymptotic Giant Branch (AGB) stars reside in multiple infrared (IR) color-color diagrams. Using MSX and 2MASS data along with radio spectra from the BAaDE survey we find that three main populations were observed in the BAaDE survey: O-rich AGB stars of which about $73\%$ host SiO masers, Carbon-rich (C) AGB stars which do not host these masers, and a small contaminating set of possible Young Stellar Objects (YSOs). The distinction between YSOs and AGB stars can be drawn using only MSX data, specifically the $[D]-[E]$ color, while the shorter wavelengths provided by 2MASS are necessary to divide potential C- and O-rich AGB stars. Divisions similar to these have been seen in multiple earlier IR-studies, but BAaDE currently provides a sample of $\sim$ 15000 sources which far exceeds previous studies in sample size, and therefore provides much more distinct divisions. With these IR distinctions in place, we discuss the sources which are exceptions in either their molecular detections or IR colors, as well as the distribution of the three populations in Galactic coordinates.

Approximately one-third of the gamma-ray sources in the third Fermi-LAT catalog are unidentified or unassociated with objects at other wavelengths. Observations with Swift-XRT have yielded possible counterparts in $\sim$30% of these source regions. The objective of this work is to identify the nature of these possible counterparts, utilizing their gamma ray properties coupled with the Swift derived X-ray properties. The majority of the known sources in the Fermi catalogs are blazars, which constitute the bulk of the extragalactic gamma-ray source population. The galactic population on the other hand is dominated by pulsars. Blazars and pulsars occupy different parameter space when X-ray fluxes are compared with various gamma-ray properties. In this work, we utilize the X-ray observations performed with the Swift-XRT for the unknown Fermi sources and compare their X-ray and gamma-ray properties to differentiate between the two source classes. We employ two machine learning algorithms, decision tree and random forest classifier, to our high signal-to-noise ratio sample of 217 sources, each of which correspond to Fermi unassociated regions. The accuracy score for both methods were found to be 97% and 99%, respectively. The random forest classifier, which is based on the application of a multitude of decision trees, associated a probability value (P$_{bzr}$) for each source to be a blazar. This yielded 173 blazar candidates with P$_{bzr}$ $\geq$ 90% for each of these sources, and 134 of these possible blazar source associations had P$_{bzr}$ $\geq$ 99%. The results yielded 13 sources with P$_{bzr}$ $\leq$ 10%, which we deemed as reasonable candidates for pulsars, 7 of which result with P$_{bzr}$ $\leq$ 1%. There were 31 sources that exhibited intermediate probabilities and were termed ambiguous due to their unclear characterization as a pulsar or a blazar.

We report on the first 1,432 sources observed using the Atacama Large Millimeter/submillimeter Array (ALMA), from the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey, which aims to obtain tens of thousands of line-of-sight velocities from SiO masers in Asymptotic Giant Branch (AGB) stars in the Milky Way. A 71% detection rate of 86 GHz SiO masers is obtained from the infrared color-selected sample, and increases to 80% when considering the likely oxygen-rich stars using Midcourse Space Experiment (MSX) colors isolated in a region where [D]-[E] <= 1.38. Based on Galactic distributions, the presence of extended CS emission, and likely kinematic associations, the population of sources with [D]-[E] > 1.38 probably consists of young stellar objects, or alternatively, planetary nebulae. For the SiO detections we examined whether individual SiO transitions provide comparable stellar line-of-sight velocities, and found that any SiO transition is suitable for determining a stellar AGB line-of-sight velocity. Finally, we discuss the relative SiO detection rates and line strengths in the context of current pumping models.

We consider interesting results on the kinematics of the bulge, beginning with Mould's (1983) measurement of the velocity dispersion of the late M giants in the bulge, and carrying the discussion on to the present day, to discuss the largest survey ever undertaken of luminous SiO maser AGB stars. This SiO maser population of AGB stars is too luminous to have evolved from the oldest globular cluster age stellar population, yet has a velocity dispersion consistent with the oldest populations (the red clump stars) in the bulge. We identify two groups of maser stars by their kinematics; the subgroup with kinematics similar to the bulge has the highest bolometric luminosities with some giants reaching $M_{bol}<-6$. In a new result, we find that the most luminous maser stars with bulge kinematics are concentrated to the Galactic plane with the most luminous maers being found at $b<\mid 1.5^\circ \mid$; while the disk maser population shows no such difference. This mirrors a similar concentration to the Galactic Center found for long period Miras, noted in earlier studies. We conclude that the maser population represents another facet of the growing tension between evidence supporting a bulge population dominated by old stars, and a significant fraction of intermediate age stars.

We observed the 43 GHz v=1, 2, and 3 and 86 GHz v=1 SiO maser transitions quasi-simultaneously for a Mira-variable-dominated sample of over 80 sources from the Bulge Asymmetries and Dynamical Evolution (BAaDE) project, using ATCA, and statistically compared the relative line strengths. On average, the 43 GHz v=1 line is brighter than the 86 GHz v=1 line by a factor of 1.36+/-0.15. As a result, an 86 GHz v=1 observed sample can be observed to 85.9% of the distance of a 43 GHz v=1 observed sample using the same sensitivity. We discuss what impact this may have on the BAaDE Galactic plane survey using the VLA and ALMA. Despite fewer v=3 detections, specific trends are discerned or strengthened when the 43 GHz v=3 line is detected. In particular the 43 and 86 GHz v=1 lines are on average equal for sources with no detectable 43 GHz v=3 emission, but the 43 GHz v=1 line strength is on average about twice as bright as the 86 GHz v=1 line for sources with detectable 43 GHz v=3 emission. Some weak correlations are found between line strengths and Midcourse Space Experiment flux densities and colors, which are tightened when considering only sources with detectable 43 GHz v=3 emission. We discuss these trends in the context of a radiative pumping model to highlight how the 43 GHz v=3 line, when coupled with the v=1 and v=2 lines, can further our understanding of variable conditions like density in the circumstellar envelopes.

We present the first results from the BAaDE (Bulge Asymmetries and Dynamic Evolution) survey. Though only a subset of the complete survey ($\sim$2700 out of $\sim$20000 final sources), our data comprise the largest radio kinematic survey to date of stellar SiO masers observed toward the Galactic bulge and plane from $-15^\circ < l < +12^\circ$ and $-6^\circ < b < +6^\circ$. Our sources include a substantial number of line-of-sight (LoS) velocities in high extinction regions within $\pm 1^\circ$ of the Galactic plane. When matched with 2MASS photometry, our radio-detected sample lies significantly brighter than and red-ward of the first red giant branch tip, reaching extremes of $(J-Ks)_0 > 8$, colors consistent with Mira variables and mass losing AGB stars. We see a clean division into two kinematic populations: a kinematically cold ($\sigma \sim 50\rm\ km\ sec^{-1}$) population that we propose is in the foreground disk, consisting of giants with 2MASS $Ks<5.5$, and a kinematically hot ($\sigma \sim 100\rm\ km\ sec^{-1}$) candidate bulge/bar population for most giants with $Ks>5.5$. Only the kinematically hot giants with $Ks>5.5$ include the reddest stars. Adopting 8.3 kpc to the Galactic Center, and correcting for foreground extinction, we find that most of the sources have $M_{bol} \sim -5$, consistent with their being luminous, and possibly intermediate age, AGB stars. We note some tension between the possibly intermediate age of the kinematically hot population, and its high velocity dispersion compared to the disk.

Blazars are a class of active galactic nuclei (AGNs) known for their very rapid variabilty in the high energy regions of the electromagnetic spectrum. Despite this known fast variability, X-ray observations have generally not revealed variability in blazars with rate doubling or halving timescales less than approximately 15 min. Since its launch, the Swift X-ray Telescope has obtained 0.2-10 keV X-ray data on 143 AGNs, including blazars, through intense target of opportunity observations that can be analyzed in a multiwavelength context and used to model jet parameters, particularly during flare states. We have analyzed this broad Swift data set in a search for short timescale variability in blazars that could limit the size of the emission region in the blazar jet. While we do find several low-significance possible flares with potential indications of rapid variability, we find no strong evidence for rapid ($<$15 minutes) doubling or halving times in flares in the soft X-ray energy band for the AGNs analyzed.

HESS J0632+057 is the only gamma-ray binary known so far whose position in the sky allows observations with ground-based observatories both in the northern and southern hemispheres. Here we report on long-term observations of HESS J0632+057 conducted with the VERITAS and H.E.S.S. Cherenkov Telescopes and the X-ray Satellite Swift, spanning a time range from 2004 to 2012 and covering most of the system's orbit. The VHE emission is found to be variable, and is correlated with that at X-ray energies. An orbital period of $315 ^{+6}_{-4}$ days is derived from the X-ray data set, which is compatible with previous results, $P = (321 \pm 5$) days. The VHE light curve shows a distinct maximum at orbital phases close to 0.3, or about 100 days after periastron passage, which coincides with the periodic enhancement of the X-ray emission. Furthermore, the analysis of the TeV data shows for the first time a statistically significant ($> 6.5 \sigma$) detection at orbital phases 0.6--0.9. The obtained gamma-ray and X-ray light curves and the correlation of the source emission at these two energy bands are discussed in the context of the recent ephemeris obtained for the system. Our results are compared to those reported for other gamma-ray binaries.

We describe a long-term Swift monitoring program of Fermi gamma-ray sources, particularly the 23 gamma-ray "sources of interest." We present a systematic analysis of the Swift X-ray Telescope light curves and hardness ratios of these sources, and we calculate excess variability. We present data for the time interval of 2004 December 22 through 2012 August 31. We describe the analysis methods used to produce these data products, and we discuss the availability of these data in an online repository, which continues to grow from more data on these sources and from a growing list of additional sources. This database should be of use to the broad astronomical community for long term studies of the variability of these objects and for inclusion in multi-wavelength studies.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between January 2008 and February 2011, resulting in 56.2 hours of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as 5.2 +/- 1.1_stat +/- 2.6_sys * 10^-12 photons cm^-2 s^-1 above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an photon index of Gamma 3.4 +/- 0.5_stat +/- 0.3_sys and a flux normalization of 1.6 +/- 0.3_stat +/- 0.8_sys * 10^-11 photons cm^-2 s^-1 at 300 GeV. We also present multiwavelength results taken in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.

HESS J0632+057 is a variable, point-like source of Very High Energy ($>100$ GeV) gamma-rays located in the Galactic plane. It is positionally coincident with a Be star, it is a variable radio and X-ray source, has a hard X-ray spectrum, and has low radio flux. These properties suggest that the object may be a member of the rare class of TeV/X-ray binary systems. The definitive confirmation of this would be the detection of a periodic orbital modulation of the flux at any wavelength. We have obtained {\it Swift} X-ray telescope observations of the source from MJD 54857 to 55647 (Jan. 2009 - Mar. 2011) to test the hypothesis that HESS J0632+057 is an X-ray/TeV binary. We show that these data exhibit flux modulation with a period of $321 \pm 5$ days and we evaluate the significance of this period by calculating the null hypothesis probability, allowing for stochastic flaring. This periodicity establishes the binary nature of HESS J0632+057.

During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.

New generation TeV gamma-ray telescopes have discovered many new sources, including several enigmatic unidentified TeV objects. HESS J0632+057 is a particularly interesting unidentified TeV source since: it is a point source, it has a possible hard-spectrum X-ray counterpart and a positionally consistent Be star, it has evidence of long-term VHE flux variability, and it is postulated to be a newly detected TeV/X-ray binary. We have obtained Swift X-ray telescope observations of this source from MJD 54857 to 54965, in an attempt to ascertain its nature and to investigate the hypothesis that it's a previously unknown X-ray/TeV binary. Variability and spectral properties similar to those of the other 3 known X-ray/TeV binaries have been observed, with measured flux increases by factors of approximately 3. X-ray variability is present on multiple timescales including days to months; however, no clear signature of periodicity is present on the timescales probed by these data. If binary modulation is present and dominating the measured variability, then the period of the orbit is likely to be more than 54 days (half of this campaign), or it has a shorter period with a variable degree of flux modulation on successive high states. If the two high states measured to date are due to binary modulation, then the favored period is approximately 35-40 days. More observations are required to determine if this object is truly a binary system and to determine the extent that the measured variability is due to inter-orbit flaring effects or periodic binary modulation.

We present the catalog of sources detected in the first 22 months of data from the hard X-ray survey (14--195 keV) conducted with the BAT coded mask imager on the \swift satellite. The catalog contains 461 sources detected above the 4.8 sigma level with BAT. High angular resolution X-ray data for every source from Swift XRT or archival data have allowed associations to be made with known counterparts in other wavelength bands for over 97% of the detections, including the discovery of ~30 galaxies previously unknown as AGN and several new Galactic sources. A total of 266 of the sources are associated with Seyfert galaxies (median redshift z ~ 0.03) or blazars, with the majority of the remaining sources associated with X-ray binaries in our Galaxy. This ongoing survey is the first uniform all sky hard X-ray survey since HEAO-1 in 1977. Since the publication of the 9-month BAT survey we have increased the number of energy channels from 4 to 8 and have substantially increased the number of sources with accurate average spectra. The BAT 22-month catalog is the product of the most sensitive all-sky survey in the hard X-ray band, with a detection sensitivity (4.8 sigma) of 2.2e-11 erg/cm2/s (1 mCrab) over most of the sky in the 14--195 keV band.

We present results from a long-term monitoring campaign on the TeV binary LSI +61 303 with VERITAS at energies above 500 GeV, and in the 2-10 keV hard X-ray bands with RXTE and Swift, sampling nine 26.5 day orbital cycles between September 2006 and February 2008. The binary was observed by VERITAS to be variable, with all integrated observations resulting in a detection at the 8.8 sigma (2006/2007) and 7.3 sigma (2007/2008) significance level for emission above 500 GeV. The source was detected during active periods with flux values ranging from 5 to 20% of the Crab Nebula, varying over the course of a single orbital cycle. Additionally, the observations conducted in the 2007-2008 observing season show marginal evidence (at the 3.6 sigma significance level) for TeV emission outside of the apastron passage of the compact object around the Be star. Contemporaneous hard X-ray observations with RXTE and Swift show large variability with flux values typically varying between 0.5 and 3.0*10^-11 ergs cm^-2 s^-1 over a single orbital cycle. The contemporaneous X-ray and TeV data are examined and it is shown that the TeV sampling is not dense enough to detect a correlation between the two bands.

We report on the results of a long term X-ray monitoring campaign of the galactic binary LS I +61 303 performed by the Rossi X-ray Timing Explorer. This dataset consists of 1 ks pointings taken every other day between 2007 August 28 until 2008 February 2. The observations covered six full cycles of the 26.496 day binary period and constitute the largest continuous X-ray monitoring dataset on LS I +61 303 to date with this sensitivity. There is no statistically strong detection of modulation of flux or photon index with orbital phase; however, we do find a strong correlation between flux and photon index, with the spectrum becoming harder at higher fluxes. The dataset contains three large flaring episodes, the largest of these reaching a flux level of 7.2 (+0.1,-0.2)*10^-11 erg cm^-2 s^-1 in the 3-10 keV band, which is a factor of three times larger than flux levels typically seen in the system. Analysis of these flares shows the X-ray emission from LS I +61 303 changing by up to a factor of six over timescales of several hundred seconds as well as doubling times as fast as 2 seconds. This is the fastest variability ever observed from LS I +61 303 at this wavelength and places constraints on the size of the X-ray emitting region.

Swift-XRT observations of the X-ray emission from gamma ray bursts (GRBs) and during the GRB afterglow have led to many new results during the past two years. One of these exciting results is that approximately 1/3-1/2 of GRBs contain detectable X-ray flares. The mean fluence of the X-ray flares is ~10 times less than that of the initial prompt emission, but in some cases the flare is as energetic as the prompt emission itself. The flares display fast rises and decays, and they sometimes occur at very late times relative to the prompt emission (sometimes as late as 10^5 s after T_0) with very high peak fluxes relative to the underlying afterglow decay that has clearly begun prior to some flares. The temporal and spectral properties of the flares are found to favor models in which flares arise due to the same GRB internal engine processes that spawned the prompt GRB emission. Therefore, both long and short GRB internal engine models must be capable of producing high fluences in the X-ray band at very late times.

Observations of gamma ray bursts (GRBs) with Swift produced the initially surprising result that many bursts have large X-ray flares superimposed on the underlying afterglow. The flares were sometimes intense, had rapid rise and decay phases, and occurred late relative to the ``prompt'' phase. Some remarkable flares are observed with fluence comparable to the prompt GRB fluence. Many GRBs have several flares, which are sometimes overlapping. Short, intense, repetitive, and late flaring can be most easily understood within the context of the standard fireball model with the internal engine that powers the prompt GRB emission in an active state at late times. However, other models for flares have been proposed. Flare origin can be investigated by comparing the flare spectra to that of the afterglow and the initial prompt emission. In this work, we have analyzed all significant X-ray flares from the first 110 GRBs observed by Swift. From this sample 33 GRBs were found to have significant X-ray flares, with 77 flares that were detected above the 3$\sigma$ level. In addition to temporal analysis presented in a companion paper, a variety of spectral models have been fit to each flare. In some cases, we find that the spectral fits favor a Band function model, which is more akin to the prompt emission than to that of an afterglow. We find that the average fluence of the flares is 2.4e-7 erg/cm^2/s in the 0.2-10 keV energy band, which is approximately a factor of ten below the average prompt GRB fluence. These results, when combined with those presented in the companion paper on temporal properties of flares, supports the hypothesis that most X-ray flares are late-time activity of the internal engine that spawned the initial GRB; not an afterglow related effect.

We present the first systematic investigation of the morphological and timing properties of flares in GRBs observed by Swift/XRT. We consider a large sample drawn from all GRBs detected by Swift, INTEGRAL and HETE-2 prior to 2006 Jan 31, which had an XRT follow-up and which showed significant flaring. Our sample of 33 GRBs includes long and short, at low and high redshift, and a total of 69 flares. The strongest flares occur in the early phases, with a clear anti-correlation between the flare peak intensity and the flare time of occurrence. Fitting each X-ray flare with a Gaussian model, we find that the mean ratio of the width and peak time is <Delta t / t > = 0.13+/-0.10, albeit with a large scatter. Late flares at times > 2000 seconds have long durations, Delta t>300 s, and can be very energetic compared to the underlying continuum. We further investigated if there is a clear link between the number of pulses detected in the prompt phase by BAT and the number of X-ray flares detected by XRT, finding no correlation. However, we find that the distribution of intensity ratios between successive BAT prompt pulses and that between successive XRT flares is the same, an indication of a common origin for gamma-ray pulses and X-ray flares. All evidence indicates that flares are indeed related to the workings of the central engine and, within the standard fireball scenario, originate from internal shocks rather than external shocks. While all flares can be explained by long-lasting engine activity, 29/69 flares may also be explained by refreshed shocks. However, 10 can only be explained by prolonged activity of the central engine.

The Swift X-ray Telescope (XRT) has discovered that flares are quite common in early X-ray afterglows of Gamma-Ray Bursts (GRBs), being observed in roughly 50% of afterglows with prompt followup observations. The flares range in fluence from a few percent to ~ 100% of the fluence of the prompt emission (the GRB). Repetitive flares are seen, with more than 4 successive flares detected by the XRT in some afterglows. The rise and fall times of the flares are typically considerably smaller than the time since the burst. These characteristics suggest that the flares are related to the prompt emission mechanism, but at lower photon energies. We conclude that the most likely cause of these flares is late-time activity of the GRB central engine.