The Virgin Island Robotic Telescope is located at the Etelman Observatory, St Thomas, since 2002. We will present its evolution since that date with the changes we have performed in order to modify an automated instrument, needing human supervision, to a fully robotic observatory. The system is based on ROS (Robotic Observatory Software) developed for TAROT and now installed on various observatories across the world (Calern, La Silla, Zadko, Les Markes, Etelman Observatory).

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

Isolated Swift gamma-ray burst (GRB) pulses, like their higher-energy BATSE counterparts, emit the bulk of their pulsed emission as a hard-to-soft component that can be fitted by the Norris et al. (2005) empirical pulse model. This signal is overlaid by a fainter, three-peaked signal that can be modeled by an empirical wave-like function (Hakkila and Preece, 2014): the two fits combine to reproduce GRB pulses with distinctive three-peaked shapes. The precursor peak appears on or before the pulse rise and is often the hardest component, the central peak is the brightest, and the decay peak converts exponentially decaying emission into a long, soft, power-law tail. Accounting for systematic instrumental differences, the general characteristics of the fitted pulses are remarkably similar. Isolated GRB pulses are dominated by hard-to-soft evolution; this is more pronounced for asymmetric pulses than for symmetric ones. Isolated GRB pulses can also exhibit intensity tracking behaviors that, when observed, are tied to the timing of the three peaks: pulses with the largest maximum hardnesses are hardest during the precursor, those with smaller maximum hardnesses are hardest during the central peak, and all pulses can re-harden during the central peak and/or during the decay peak. Since these behaviors are essentially seen in all isolated pulses, the distinction between "hard-to-soft" and "intensity-tracking" pulses really no longer applies. Additionally, the triple-peaked nature of isolated GRB pulses seems to indicate that energy is injected on three separate occasions during the pulse duration: theoretical pulse models need to account for this.

Correlations among pulse properties in the prompt emission of long GRBs can potentially be used as cosmological distance indicators to estimate redshifts of GRBs to which these pulses belong. We demonstrate application of this technique to a sample of GRBs for which redshifts are not known. We also study the scatter of predicted redshifts of pulses found within individual bursts. We explore the characteristics of this scatter in hopes of identifying systematic corrections and/or pulse subsets that can be used to increase the technique's reliability.

We derive a new peak lag vs. peak luminosity relation in gamma-ray burst (GRB) pulses. We demonstrate conclusively that GRB spectral lags are pulse rather than burst properties and show how the lag vs. luminosity relation determined from CCF measurements of burst properties is essentially just a rough measure of this newly derived relation for individual pulses. We further show that most GRB pulses have correlated properties: short-lag pulses have shorter durations, are more luminous, and are harder within a burst than long-lag pulses. We also uncover a new pulse duration vs. pulse peak luminosity relation, and indicate that long-lag pulses often precede short-lag pulses. Although most pulse behaviors are supportive of internal shocks (including long-lag pulses), we identify some pulse shapes that could result from external shocks.

A simple scheme delineates Long GRBs with similar time histories using the Internal Luminosity Function power-law index and the spectral lag. Several generalizations are made about time history morphologies: 1) GRBs with long spectral lags contain fewer pulses that are broader than those found in bursts with short spectral lags, 2) short-lag bursts with small ILF power-law indices have many narrow pulses and are highly variable, while long-lag bursts with small ILF power-law indices are characterized by broad, smooth pulses and have low variability. GRB time history morphologies primarily identify intrinsic rather than extrinsic characteristics based on correlations with gamma-ray luminosity, afterglow luminosity, and numbers of pulses. These characteristics result because internal relativistic effects due to bulk Lorentz factor are larger than cosmological effects, and because the numbers and shapes of pulses indicate different efficiencies and forms of GRB energy release. Single-pulsed GRBs are characterized either by large ILF power indices (indicating a range of jet opening angles and Lorentz factors, with a FRED pulse shape), or they have long lags (large jet opening angles with low Lorentz factors, with either a FRED pulse shape or an unpeaked, smooth pulse shape). They also have lower-luminosity afterglows than multi-pulsed GRBs. GRBs with simple time histories are often associated with Type Ibc supernovae. This suggests that some single-pulsed GRBs contain single, beamed blast waves that are similar to and have characteristics that overlap with those of many supernovae. Such a connection may not exist between multi-pulsed GRBs and supernovae.

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.

We have identified two quiescent GRBs (bursts having two or more widely-separated distinct emission episodes) in which the post-quiescent emission exhibits distinctly different characteristics than the pre-quiescent emission. In these two cases (BATSE GRBs 960530 and 980125), the second emission episode has a longer lag, a smoother morphology, and softer spectral evolution than the first episode. Although the pre-quiescent emission satisfies the standard internal shock paradigm, we demonstrate that the post-quiescent emission is more consistent with external shocks. We infer that some observed soft, faint, long-lag GRBs are external shocks in which the internal shock signature is not observed. We further note that the peak luminosity ratio between quiescent episodes is not in agreement with the ratio predicated by the lag vs. peak luminosity relationship. We briefly discuss these observations in terms of current collapsar jet models.

Unsupervised pattern recognition algorithms support the existence of three gamma-ray burst classes; Class I (long, large fluence bursts of intermediate spectral hardness), Class II (short, small fluence, hard bursts), and Class III (soft bursts of intermediate durations and fluences). The algorithms surprisingly assign larger membership to Class III than to either of the other two classes. A known systematic bias has been previously used to explain the existence of Class III in terms of Class I; this bias allows the fluences and durations of some bursts to be underestimated (Hakkila et al., ApJ 538, 165, 2000). We show that this bias primarily affects only the longest bursts and cannot explain the bulk of the Class III properties. We resolve the question of Class III existence by demonstrating how samples obtained using standard trigger mechanisms fail to preserve the duration characteristics of small peak flux bursts. Sample incompleteness is thus primarily responsible for the existence of Class III. In order to avoid this incompleteness, we show how a new dual timescale peak flux can be defined in terms of peak flux and fluence. The dual timescale peak flux preserves the duration distribution of faint bursts and correlates better with spectral hardness (and presumably redshift) than either peak flux or fluence. The techniques presented here are generic and have applicability to the studies of other transient events. The results also indicate that pattern recognition algorithms are sensitive to sample completeness; this can influence the study of large astronomical databases such as those found in a Virtual Observatory.

The connection between Gamma-Ray Bursts (GRBs) and their afterglows is currently not well understood. Afterglow models of synchrotron emission generated by external shocks in the GRB fireball model predict emission detectable in the gamma-ray regime ($\gax 25$ keV). In this paper, we present a temporal and spectral analysis of a subset of BATSE GRBs with smooth extended emission tails to search for signatures of the ``early high-energy afterglow'', i.e., afterglow emission that initially begins in the gamma-ray phase and subsequently evolves into X-Ray, uv, optical, and radio emission as the blast wave is decelerated by the ambient medium. From a sample of 40 GRBs we find that the temporal decays are best described with a power-law $\sim t^{\beta}$, rather than an exponential, with a mean index $<\beta > \approx -2$. Spectral analysis shows that $\sim 20%$ of these events are consistent with fast-cooling synchrotron emission for an adiabatic blast wave; three of which are consistent with the blast wave evolution of a jet, with $F_{\nu} \sim t^{-p}$. This behavior suggests that, in some cases, the emission may originate from a narrow jet, possibly consisting of ``nuggets'' whose angular size are less than $1 / \Gamma$, where $\Gamma$ is the bulk Lorentz factor.

A possible relationship between the presence of a radio afterglow and gamma-ray burst spectral hardness is discussed. The correlation is marginally significant; the spectral hardness of the bursts with radio afterglows apparently results from a combination of the break energy Ebreak and the high-energy spectral index beta. If valid, this relationship would indicate that the afterglow does carry information pertaining to the GRB central engine.

We present evidence for burst emission from SGR 1900+14 with a power-law high energy spectrum extending beyond 500 keV. Unlike previous detections of high energy photons during bursts from SGRs, these emissions are not associated with high-luminosity burst intervals. Not only is the emission hard, but the spectra are better fit by Band's GRB function rather than by the traditional optically-thin thermal bremsstrahlung model. We find that the spectral evolution within these hard events obeys a hardness/intensity anti-correlation. Temporally, these events are distinct from typical SGR burst emissions in that they are longer (~ 1 s) and have relatively smooth profiles. Despite a difference in peak luminosity of > 1E+11 between these bursts from SGR 1900+14 and cosmological GRBs, there are striking temporal and spectral similarities between the two kinds of bursts, aside from spectral evolution. We outline an interpretation of these events in the context of the magnetar model.

In this Letter, we present the first evidence in the BATSE data for a prompt high-energy (25-300 keV) afterglow component from a gamma-ray burst (GRB), GRB980923. The event consists of rapid variabilty lasting ~40 s followed by a smooth power law emission tail lasting ~400 s. An abrupt change in spectral shape is found when the tail becomes noticeable. Our analysis reveals that the spectral evolution in the tail of the burst mimics that of a cooling synchrotron spectrum, similar to the spectral evolution of the low-energy afterglows for GRBs. This evidence for a separate emission component is consistent with the internal-external shock scenario in the relativistic fireball picture. In particular, it illustrates that the external shocks can be generated during the gamma-ray emission phase, as in the case of GRB990123.

Following a BeppoSAX alert (Piro 1999a) and the discovery of the OT at SAAO (Vreeswijk et al. 1999a), we observed GRB 990510 with the FORS instrument on ESO's VLT Unit 1 (`Antu'). The burst is unremarkable in gamma rays, but in optical is the first one to show good evidence for jet-like outflow (Stanek et al. 1999, Harrison et al. 1999). We report the detection of significant linear polarization in the afterglow: it is (1.6 +/- 0.2)% 0.86 days after trigger, and after 1.81 days is consistent with that same value, but much more uncertain. The polarization angle is constant on a time scale of hours, and may be constant over one day. We conclude that the polarization is intrinsic to the source and due to the synchrotron nature of the emission, and discuss the random and ordered field geometries that may be responsible for it.