We present a study of interstellar comet 2I/2019 Q4 (Borisov) using both preperihelion and postperihelion observations spanning late September 2019 through late January 2020. The intrinsic brightness of the comet was observed to continuously decline throughout the timespan, likely due to the decreasing effective scattering cross-section as a result of volatile sublimation with a slope of $-0.43 \pm 0.02$ km$^{2}$ d$^{-1}$. We witnessed no significant change in the slightly reddish colour of the comet, with mean values of $\left \langle g - r \right \rangle = 0.68 \pm 0.04$, $\left \langle r - i \right \rangle = 0.23 \pm 0.03$, and the normalised reflectivity gradient across the $g$ and $i$ bands $\overline{S'} \left(g,i\right) = \left(10.6 \pm 1.4\right)$ % per $10^3$ \AA, all unremarkable in the context of solar system comets. Using the available astrometric observations, we confidently detect the nongravitational acceleration of the comet following a shallow heliocentric distance dependency of $r_{\rm H}^{-1 \pm 1}$. Accordingly, we estimate that the nucleus is most likely $\lesssim$0.4 km in radius, and that a fraction of $\gtrsim$0.2% of the total nucleus in mass has been eroded due to the sublimation activity since the earliest observation of the comet in 2018 December by the time of perihelion. Our morphology simulation suggests that the dust ejection speed increased from $\sim$4 m s$^{-1}$ in September 2019 to $\sim$7 m s$^{-1}$ around perihelion for the optically dominant dust grains of $\beta \sim 0.01$, and that the observable dust grains are no smaller than micron size.

We present a study of C/2019 Y4 (ATLAS) using Sloan $gri$ observations from mid-January to early April 2020. During this timespan, the comet brightened with a growth in the effective cross-section of $\left(2.0 \pm 0.1 \right) \times 10^{2}$ m$^{2}$ s$^{-1}$ from the beginning to $\sim$70 d preperihelion in late March 2020, followed by a brightness fade and the comet gradually losing the central condensation. Meanwhile, the comet became progressively bluer, and was even bluer than the Sun ($g - r \approx 0.2$) when the brightness peaked, likely due to activation of subterranean fresh volatiles exposed to sunlight. With the tailward-bias corrected astrometry we found an enormous radial nongravitational parameter, $A_{1} = \left(+2.25 \pm 0.13\right) \times 10^{-7}$ au d$^{-2}$ in the heliocentric motion of the comet. Taking all of these finds into consideration, we conclude that the comet has disintegrated since mid-March 2020. By no means was the split new to the comet, as we quantified that the comet had undergone another split event around last perihelion $\sim$5 kyr ago, during which its sibling C/1844 Y1 (Great Comet) was produced, with the in-plane component of the separation velocity $\gtrsim$1 m s$^{-1}$. We constrained that the nucleus of C/2019 Y4 before disintegration was $\gtrsim$60 m in radius, and has been protractedly ejecting dust grains of $\sim$10-40 $\mu$m (assuming dust bulk density 0.5 g cm$^{-3}$) with ejection speed $\sim$30 m s$^{-1}$ in early March 2020 and increased to $\sim$80 m s$^{-1}$ towards the end of the month for grains of $\sim$10 $\mu$m.

We present Hubble Space Telescope observations of a photometric outburst and splitting event in interstellar comet 2I/Borisov. The outburst, first reported with the comet outbound at 2.8 AU (Drahus et al.~2020), was caused by the expulsion of solid particles having a combined cross-section about 100 sq. km and a mass in 0.1 mm sized particles about 2e7 kg. The latter corresponds to 1e-4 of the mass of the nucleus, taken as a sphere of radius 500 m. A transient ``double nucleus'' was observed on UT 2020 March 30 (about three weeks after the outburst), having a cross-section about 0.6 sq. km and corresponding dust mass 1e5 kg. The secondary was absent in images taken on and before March 28, and in images taken on and after April 03. The unexpectedly delayed appearance and rapid disappearance of the secondary are consistent with an origin through rotational bursting of one or more large (meter-sized) boulders under the action of outgassing torques, following their ejection from the main nucleus. Overall, our observations reveal that the outburst and splitting of the nucleus are minor events involving a negligible fraction of the total mass: 2I/Borisov will survive its passage through the planetary region largely unscathed.

Hubble Space Telescope observations of interstellar comet 2I/Borisov near perihelion show the ejection of large (>~100 um) particles at <~9 m/s speeds, with estimated mass-loss rates of ~35 kg/s. The total mass loss from comet Borisov corresponds to loss of a surface shell on the nucleus only ~0.4 m thick. This shell is thin enough to be susceptible to past chemical processing in the interstellar medium by cosmic rays, meaning that the ejected materials cannot necessarily be considered as pristine. Our high-resolution images reveal persistent asymmetry in the dust coma, best explained by a thermal lag on the rotating nucleus causing peak mass loss to occur in the comet nucleus afternoon. In this interpretation, the nucleus rotates with an obliquity of 30 deg (pole direction RA = 205 deg and Dec. = 52 deg). The subsolar latitude varied from -35 deg (southern solstice) at the time of discovery to 0 deg (equinox) in 2020 January, suggesting the importance of seasonal effects. Subsequent activity likely results from regions freshly activated as the northern hemisphere is illuminated for the first time.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

We present high resolution imaging observations of interstellar comet 2I/Borisov (formerly C/2019 Q4) obtained using the Hubble Space Telescope. Scattering from the comet is dominated by a coma of large particles (characteristic size 0.1 mm) ejected anisotropically. Convolution modeling of the coma surface brightness profile sets a robust limit to the spherical-equivalent nucleus radius r_n < 0.5 km (geometric albedo 0.04 assumed). We obtain an independent constraint based on the non-gravitational acceleration of the nucleus, finding r_n > 0.2 km (nucleus density 500 kg/m3 assumed). The profile and the non-gravitational constraints cannot be simultaneously satisfied if density < 25 kg/m3; the nucleus of comet Borisov cannot be a low density fractal assemblage of the type proposed elsewhere for the nucleus of 1I/'Oumuamua. We show that the spin-up timescale to outgassing torques, even at the measured low production rates, is comparable to or shorter than the residence time in the Sun's water sublimation zone. The spin angular momentum of the nucleus should be changed significantly during the current solar fly-by. Lastly, we find that the differential interstellar size distribution in the 0.5 mm to 100 m size range can be represented by power laws with indices < 4 and that interstellar bodies of 100 m size scale strike Earth every one to two hundred million years.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Object 1I/2017 U1 (`Oumuamua) is the first interstellar small body ever discovered in the solar system. By the time of discovery, it had already passed perihelion. To investigate the behavior of `Oumuamua around perihelion, we searched for it in Solar and Heliospheric Observatory (SOHO) and Solar TErrestrial RElations Observatory (STEREO) images from early 2017 September (preperihelion), but did not detect it. The nondetection of `Oumuamua by STEREO renders more stringent constraints on its physical properties thanks to the extreme forward-scattering observing geometry. Assuming geometric albedo $p_V = 0.1$, the effective scattering cross-section of any dust coma was $\lesssim \left(2.1 \pm 0.2 \right) \times 10^{4}$ m$^{2}$. Assuming it behaved like a typical solar-system comet this would correspond to a total mass of $\lesssim 20 \pm 2$ kg, and a water production rate of $\lesssim \left(6.1 \pm 0.5 \right) \times 10^{25}$ s$^{-1}$ at heliocentric distance $r_{\rm H} = 0.375$ au. If scaled to post-discovery $r_{\rm H}$, the water production rate would be smaller than any of the previously reported upper limits by at least an order of magnitude. To exhibit the reported nongravitational motion with our default assumptions requires a nucleus bulk density $\lesssim$40 kg m$^{-3}$; higher bulk densities are possible for other assumptions. Alternatively, we show that thermal fracturing could have plausibly removed an inert surface layer between these observations and discovery, thus initiating activity after `Oumuamua left the field of view of STEREO.

The first simultaneous detection of a short gamma-ray burst (SGRB) with a gravitational-wave (GW) signal ushered in a new era of multi-messenger astronomy. In order to increase the number of SGRB-GW simultaneous detections, we need full sky coverage in the gamma-ray regime. BurstCube, a CubeSat for Gravitational Wave Counterparts, aims to expand sky coverage in order to detect and localize gamma-ray bursts (GRBs). BurstCube will be comprised of 4 Cesium Iodide scintillators coupled to arrays of Silicon photo-multipliers on a 6U CubeSat bus (a single U corresponds to cubic unit $\sim$10 cm $\times$ 10 cm $\times$ 10 cm) and will be sensitive to gamma-rays between 50 keV and 1 MeV, the ideal energy range for GRB prompt emission. BurstCube will assist current observatories, such as $Swift$ and $Fermi$, in the detection of GRBs as well as provide astronomical context to gravitational wave events detected by Advanced LIGO, Advanced Virgo, and KAGRA. BurstCube is currently in its development and testing phase to prepare for launch readiness in the fall of 2021. We present the mission concept, preliminary performance, and status.

We report on a multi-observatory campaign to examine asteroid 3200 Phaethon during its December 2017 close approach to Earth, in order to improve our measurements of its fundamental parameters, and to search for surface variations, cometary activity and fragmentation. The mean colors of Phaethon are B-V = 0.702 +/- 0.004, V-R = 0.309 +/- 0.003, R-I = 0.266 +/- 0.004, neutral to slightly blue, consistent with previous classifications of Phaethon as a F-type or B-type asteroid. Variations in Phaethon's B-V colors (but not V-R or R-I) with observer sub-latitude are seen and may be associated with craters observed by the Arecibo radar. High cadence photometry over phases from 20 to 100 degrees allows a fit to the values of the HG photometric parameters; H = 14.57 +/- 0.02, 13.63 +/- 0.02, 13.28 +/- 0.02, 13.07 +/- 0.02; G = 0.00 +/- 0.01, -0.09 +/- 0.01, -0.10 +/- 0.01, -0.08 +/- 0.01 in the BVRI filters respectively; the negative G values are consistent with other observations of F type asteroids. Light curve variations were seen that are also consistent with concavities reported by Arecibo, indicative of large craters on Phaethon's surface whose ejecta may be the source of the Geminid meteoroid stream. A search for gas/dust production set an upper limit of 0.06 +/- 0.02 kg/s when Phaethon was 1.449 AU from the Sun, and 0.2 +/- 0.1 kg/s at 1.067 AU. A search for meter-class fragments accompanying Phaethon did not find any whose on-sky motion was not also consistent with background main belt asteroids.

Main-belt asteroid (6478) Gault was observed to show cometary features in early 2019. To investigate the cause, we conducted {\it BVR} observations at Xingming Observatory, China, from 2019 January to April. The two tails were formed around 2018 October 26--November 08, and 2018 December 29--2019 January 08, respectively, and consisted of dust grains of $\gtrsim$20 $\mu$m to 3 mm in radius ejected at a speed of $0.15 \pm 0.05$ m s$^{-1}$ and following a broken power-law size distribution bending at grain radius $\sim$70 $\mu$m (bulk density 1 g cm$^{-3}$ assumed). The total mass of dust within a $10^4$ km-radius aperture around Gault declined from $\sim$$9 \times 10^6$ kg since 2019 January at a rate of $2.28 \pm 0.07$ kg s$^{-1}$, but temporarily surged around 2019 March 25, because Earth thence crossed the orbital plane of Gault, within which the ejected dust was mainly distributed. No statistically significant colour or short-term lightcurve variation was seen. Nonetheless we argue that Gault is currently subjected to rotational instability. Using the available astrometry, we did not detect any nongravitational acceleration in the orbital motion of Gault.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

We present a photometric and dynamical study of comet C/2010 U3 (Boattini), which was seen active in prediscovery data as early as 2005 November at a new inbound record heliocentric distance $r_{\rm H} = 25.8$ au. Two outburst events around 2009 and 2017 were observed. The coma and tail of the comet consist of dust grains of $\sim$10 $\mu$m in radius, ejected protractedly at speeds $\lesssim$50 m s$^{-1}$ near the subsolar point, and are subjected to the Lorentz force, solar gravitation and radiation pressure force altogether. The prolonged activity indicates that sublimation of supervolatiles (e.g., CO, CO$_2$) is at play, causing a net mass-loss rate $\gtrsim$1 kg s$^{-1}$. To sustain the mass loss, the nucleus radius has to be $\gtrsim$0.1 km. The color of the cometary dust, similar to other long-period comets, is redder than the solar colors, but we also observed potential color variations when the comet was at $10 < r_{\rm H} < 15$ au, concurrent with the onset of crystallisation of amorphous water ice, if at all. Using publicly available and our refined astrometric measurements, we estimated the precise trajectory of the comet and propagated it backward to its previous perihelion. We found that the comet has visited the planetary region $1.96 \pm 0.04$ Myr ago, with barycentric perihelion distance $q = 8.364 \pm 0.004$ au. Thus, C/2010 U3 (Boattini) is almost certainly a dynamically old comet from the Oort cloud, and the observed activity cannot be caused by retained heat from the previous apparition.

Distant long-period comet C/2017 K2 has been outside the planetary region of the solar system for 3 Myr, negating the possibility that heat retained from the previous perihelion could be responsible for its activity. This inbound comet is also too cold for water ice to sublimate and too cold for amorphous water ice, if present, to crystallize. C/2017 K2 thus presents an ideal target in which to investigate the mechanisms responsible for activity in distant comets. We have used Hubble Space Telescope to study the comet in the pre-perihelion distance range 13.8 to 15.9 AU. The coma maintains a logarithmic surface brightness gradient $m = -1.010\pm$0.004, consistent with steady-state mass loss. The absence of a radiation pressure swept tail indicates that the effective particle size is large (0.1 mm) and the mass loss rate is $\sim$200 kg s$^{-1}$, remarkable for a comet still beyond the orbit of Saturn. Extrapolation of the photometry indicates that activity began in 2012.1, at 25.9$\pm$0.9 AU, where the blackbody temperature is only 55 K. This large distance and low temperature suggest that cometary activity is driven by the sublimation of a super-volatile ice (e.g.~CO), presumably preserved by K2's long-term residence in the Oort cloud. The mass loss rate can be sustained by CO sublimation from an area $\lesssim 2$ km$^2$, if located near the hot sub-solar point on the nucleus. However, while the drag force from sublimated CO is sufficient to lift millimeter sized particles against the gravity of the cometary nucleus, it is 10$^2$ to 10$^3$ times too small to eject these particles against inter-particle cohesion. Our observations thus require either a new understanding of the physics of inter-particle cohesion or the introduction of another mechanism to drive distant cometary mass loss. We suggest thermal fracture and electrostatic supercharging in this context.

Near-Earth asteroid (3200) Phaethon is notable for its association to a strong annual meteor shower, the Geminids, indicative of one or more episodes of mass ejection in the past. The mechanism of Phaethon's past activity is not yet understood. Here we present a Hubble Space Telescope (HST) search of meter-sized fragments in the vicinity of Phaethon, carried out during Phaethon's historic approach to the Earth in mid-December of 2017. Numerical simulation conducted to guide HST's pointing also show that the dynamical evolution of Phaethon-originated particles is quick, as ejected materials take no longer than $\sim250$ yr to spread to the entire orbit of Phaethon. Our search was completed down to 4-meter-class limit (assuming Phaethon-like albedo) and was expected to detect 0.035% particles ejected by Phaethon in the last several decades. The negative result of our search capped the total mass loss of Phaethon over the past few dozen orbits to be $10^{12}$ kg at $3\sigma$ level, taking the best estimates of size power-law from meteor observations and spacecraft data. Our result also implies a millimeter-sized dust flux of $<10^{-12} \mathrm{m^{-2} s^{-1}}$ within 0.1 au of Phaethon, suggesting that any Phaethon-bound mission is unlikely to encounter dense dust clouds.

It depends. Our experiment reveals that, given an optically thin coma, generally, the smaller the signal ratio of nucleus to coma, the less reliable is the cometary nucleus-extraction technique. We strongly suggest the technique only be applied to cases where the nucleus signal occupies $\gtrsim$10% of the total signal wherein the bias is no more than a few percent. Otherwise there is probably no way to debias results from this technique in reality, since its reliability is highly sensitive to entangling complications, including the coma profile, and the point-spread function (PSF).

We present a study of two newly discovered heliocentric hyperbolic objects -- 2017 U7 and 2018 C2. Both are possibly thermally evolved comets. Observations of the latter in 2018 March from Xingming Observatory revealed that it has a color similar to those of the long-period comets, Trojans, and D-type asteroids: $m_B - m_V = 0.75 \pm 0.03$, $m_V - m_R = 0.41 \pm 0.02$, and $m_R - m_I = 0.37 \pm 0.03$. A possible extremely faint coma of $\sim$9$"$ across was observed. The radial profile of 2018 C2 in comparison with those of the field stars helps confirm its cometary nature. Based on our convolution model, its mass-loss rate is estimated to be $0.7 \pm 0.2$ kg s$^{-1}$. Assuming geometric albedo $p_R = 0.04$, its effective radius is $4.4 \pm 0.5$ km, which means that the fraction of active area is merely $\sim$10$^{-5}$-10$^{-4}$. Our N-body dynamical simulations show that both objects are most likely dynamically old members from the Oort cloud. 2017 U7 has a $\sim$60% chance to escape the solar system in the 1.5 Myr following its current perihelion passage due to a moderate close encounter with Jupiter in 2020 May. On the condition that no disintegration occurs, 2018 C2 will revisit the inner solar system 0.13-0.14 Myr later, with perihelion distance $\sim$2 AU.

Jupiter-family Comet D/1770 L1 (Lexell) was the first discovered Near-Earth Object (NEO), and passed the Earth on 1770 Jul 1 at a recorded distance of 0.015 au. The comet was subsequently lost due to unfavorable observing circumstances during its next apparition followed by a close encounter with Jupiter in 1779. Since then, the fate of D/Lexell has attracted interest from the scientific community, and now we revisit this long-standing question. We investigate the dynamical evolution of D/Lexell based on a set of orbits recalculated using the observations made by Charles Messier, the comet's discoverer, and find that there is a $98\%$ chance that D/Lexell remains in the Solar System by the year of 2000. This finding remains valid even if a moderate non-gravitational effect is imposed. Messier's observations also suggest that the comet is one of the largest known near-Earth comets, with a nucleus of $\gtrsim 10$ km in diameter. This implies that the comet should have been detected by contemporary NEO surveys regardless of its activity level if it has remained in the inner Solar System. We identify asteroid 2010 JL$_{33}$ as a possible descendant of D/Lexell, with a $0.8\%$ probability of chance alignment, but a direct orbital linkage of the two bodies has not been successfully accomplished. We also use the recalculated orbit to investigate the meteors potentially originating from D/Lexell. While no associated meteors have been unambiguously detected, we show that meteor observations can be used to better constrain the orbit of D/Lexell despite the comet being long lost.

We present a study of comet C/2017 K2 (PANSTARRS) using prediscovery archival data taken from 2013 to 2017. Our measurements show that the comet has been marginally increasing in activity since at least 2013 May (heliocentric distance of $r_{\mathrm{H}} = 23.7$ AU pre-perihelion). We estimate the mass-loss rate during the period 2013--2017 as $\overline{\dot{M}} \approx \left(2.4 \pm 1.1 \right) \times 10^{2}$ kg s$^{-1}$, which requires a minimum active surface area of $\sim$10--10$^2$ km$^{2}$ for sublimation of supervolatiles such as CO and CO$_2$, by assuming a nominal cometary albedo $p_V = 0.04 \pm 0.02$. The corresponding lower limit to the nucleus radius is a few kilometers. Our Monte Carlo dust simulations show that dust grains in the coma are $\gtrsim0.5$ mm in radius, with ejection speeds from $\sim$1--3 m s$^{-1}$, and have been emitted in a protracted manner since 2013, confirming estimates by Jewitt et al. (2017). The current heliocentric orbit is hyperbolic. Our N-body backward dynamical integration of the orbit suggests that the comet is most likely (with a probability of $\sim$98\%) from the Oort spike. The calculated median reciprocal of the semimajor axis 1 Myr ago was $a_{\mathrm{med}}^{-1} = \left( 3.61 \pm 1.71 \right) \times 10^{-5}$ AU$^{-1}$ (in a reference system of the solar-system barycentre).

We present observations showing in-bound long-period comet C/2017 K2 (PANSTARRS) to be active at record heliocentric distance. Nucleus temperatures are too low (60 K to 70 K) either for water ice to sublimate or for amorphous ice to crystallize, requiring another source for the observed activity. Using the Hubble Space Telescope we find a sharply-bounded, circularly symmetric dust coma 10$^5$ km in radius, with a total scattering cross section of $\sim$10$^5$ km$^2$. The coma has a logarithmic surface brightness gradient -1 over much of its surface, indicating sustained, steady-state dust production. A lack of clear evidence for the action of solar radiation pressure suggests that the dust particles are large, with a mean size $\gtrsim$ 0.1 mm. Using a coma convolution model, we find a limit to the apparent magnitude of the nucleus $V >$ 25.2 (absolute magnitude $H >$ 12.9). With assumed geometric albedo $p_V$ = 0.04, the limit to the nucleus circular equivalent radius is $<$ 9 km. Pre-discovery observations from 2013 show that the comet was also active at 23.7 AU heliocentric distance. While neither water ice sublimation nor exothermic crystallization can account for the observed distant activity, the measured properties are consistent with activity driven by sublimating supervolatile ices such as CO$_2$, CO, O$_2$ and N$_2$. Survival of supervolatiles at the nucleus surface is likely a result of the comet's recent arrival from the frigid Oort cloud.

We review the evidence for buried ice in the asteroid belt; specifically the questions around the so-called Main Belt Comets (MBCs). We summarise the evidence for water throughout the Solar System, and describe the various methods for detecting it, including remote sensing from ultraviolet to radio wavelengths. We review progress in the first decade of study of MBCs, including observations, modelling of ice survival, and discussion on their origins. We then look at which methods will likely be most effective for further progress, including the key challenge of direct detection of (escaping) water in these bodies.

BurstCube will detect long GRBs, attributed to the collapse of massive stars, short GRBs (sGRBs), resulting from binary neutron star mergers, as well as other gamma-ray transients in the energy range 10-1000 keV. sGRBs are of particular interest because they are predicted to be the counterparts of gravitational wave (GW) sources soon to be detectable by LIGO/Virgo. BurstCube contains 4 CsI scintillators coupled with arrays of compact low-power Silicon photomultipliers (SiPMs) on a 6U Dellingr bus, a flagship modular platform that is easily modifiable for a variety of 6U CubeSat architectures. BurstCube will complement existing facilities such as Swift and Fermi in the short term, and provide a means for GRB detection, localization, and characterization in the interim time before the next generation future gamma-ray mission flies, as well as space-qualify SiPMs and test technologies for future use on larger gamma-ray missions. The ultimate configuration of BurstCube is to have a set of $\sim10$ BurstCubes to provide all-sky coverage to GRBs for substantially lower cost than a full-scale mission.

We present a photometric and astrometric study of the split active asteroid P/2016 J1 (PANSTARRS). The two components (hereafter J1-A and J1-B) separated either $\sim$1500 days (2012 May to June) or 2300 days (2010 April) prior to the current epoch, with a separation speed $V_{\mathrm{sep}} = 0.70 \pm 0.02$ m s$^{-1}$ for the former scenario, or $0.83 \pm 0.06$ m s$^{-1}$ for the latter. Keck photometry reveals that the two fragments have similar, Sun-like colors which are comparable to the colors of primitive C- and G-type asteroids. With a nominal comet-like albedo, $p_{R} = 0.04$, the effective, dust-contaminated cross sections are estimated to be 2.4 km$^{2}$ for J1-A, and 0.5 km$^{2}$ for J1-B. We estimate that the nucleus radii lie in the range $140 \lesssim R_{\mathrm{N}} \lesssim 900$ m for J1-A and $40 \lesssim R_{\mathrm{N}} \lesssim 400$ m, for J1-B. A syndyne-synchrone simulation shows that both components have been active for 3 to 6 months, by ejecting dust grains at speeds $\sim$0.5 m s$^{-1}$ with rates $\sim$1 kg s$^{-1}$ for J1-A and 0.1 kg s$^{-1}$ for J1-B. In its present orbit, the rotational spin-up and devolatilization times of 2016 J1 are very small compared to the age of the solar system, raising the question of why this object still exists. We suggest that ice that was formerly buried within this asteroid became exposed at the surface, perhaps via a small impact, and that sublimation torques then rapidly drove it to break-up. Further disintegration events are anticipated due to the rotational instability.

Comets can exhibit non-gravitational accelerations caused by recoil forces due to anisotropic mass loss. So might active asteroids. We present an astrometric investigation of 18 active asteroids in search of non-gravitational acceleration. Statistically significant (signal-to-noise ratio (SNR) $> 3$) detections are obtained in three objects: 313P/Gibbs, 324P/La Sagra and (3200) Phaethon. The strongest and most convincing detection ($>$7$\sigma$ in each of three orthogonal components of the acceleration), is for the $\sim$1 km diameter nucleus of 324P/La Sagra. A 4.5$\sigma$ detection of the transverse component of the acceleration of 313P/Gibbs (also $\sim$1 km in diameter) is likely genuine too, as evidenced by the stability of the solution to the rejection or inclusion of specific astrometric datasets. We also find a 3.4$\sigma$ radial-component detection for $\sim$5 km diameter (3200) Phaethon, but this detection is more sensitive to the inclusion of specific datasets, suggesting that it is likely spurious in origin. The other 15 active asteroids in our sample all show non-gravitational accelerations consistent with zero. We explore different physical mechanisms which may give rise to the observed non-gravitational effects, and estimate mass-loss rates from the non-gravitational accelerations. We present a revised momentum-transfer law based on a physically realistic sublimation model for future work on non-gravitational forces, but note that it has little effect on the derived orbital elements.

We present a study of the active asteroid (3200) Phaethon in the 2016 apparition using the STEREO spacecraft and compare the results with data from the previous two perihelia in 2009 and 2012. Once again, Phaethon brightened by $\sim$2 mag soon after its perihelion passage, contradicting expectations from the phase function of a macroscopic monolithic body. Subsequently, a short antisolar tail of $\sim$0{\deg}.1 in length was formed within $\sim$1 day and quickly disappeared. No trail was seen. Our syndyne-synchrone analysis indicates that the tail was comprised of submicron to micron particles and can be approximated by a synchrone coinciding with the outburst. We estimate that the outburst has released a mass of $\sim$10$^{4}$--10$^{5}$ kg, comparable to the two mass ejections in 2009 and 2012, and that the average mass-loss rate is $\sim$0.1--1 kg s$^{-1}$. The forward-scattering effect hinted at low level activity of Phaethon prior to the outburst, which increased the effective cross section by merely $\lesssim$1 km$^2$. Without the forward-scattering enhancement, detecting such activity at side-scattering phase angles is very difficult. The forward-scattering effect also reinforces that the ejected dust grains rather than gas emissions were responsible for the activity of Phaethon. Despite Phaethon's reactivation, it is highly unlikely that the Geminid meteoroid stream can be sustained by similar perihelion mass-loss events.

Encke-type comet 332P/Ikeya-Murakami is experiencing cascading fragmentation events during its 2016 apparition. It is likely the first splitting Encke-type comet ever being observed. A nongravitational solution to the astrometry reveals a statistical detection of the radial and transverse nongravitational parameters, $A_{1} = \left(1.54 \pm 0.39\right) \times 10^{-8}$ AU day$^{-2}$, and $A_{2} = \left(7.19 \pm 1.92\right) \times 10^{-9}$ AU day$^{-2}$, respectively, which implies a nucleus erosion rate of $\left(0.91 \pm 0.17\right)$% per orbital revolution. The mass-loss rate likely has to be supported by a much larger fraction of an active surface area than known cases of short-period comets; it may be relevant to the ongoing fragmentation. We failed to detect any serendipitous pre-discovery observations of the comet in archival data from major sky surveys, whereby we infer that 332P used to be largely inactive, and is perhaps among the few short-period comets which have been reactivated from weakly active or dormant states. We therefore constrain an upper limit to the nucleus size as $2.0 \pm 0.2$ km in radius. A search for small bodies in similar orbits to that of 332P reveals comet P/2010 B2 (WISE) as the best candidate. From an empirical generalised Jupiter-family (Encke-type included) comet population model, we estimate the likelihood of chance alignment of the 332P--P/2010 B2 pair to be 1 in 33, a small number indicative of a genetic linkage between the two comets on a statistical basis. The pair possibly originated from a common progenitor which underwent a disintegration event well before the twentieth century.

We present initial time-resolved observations of the split comet 332P/Ikeya-Murakami taken using the Hubble Space Telescope. Our images reveal a dust-bathed cluster of fragments receding from their parent nucleus at projected speeds in the range 0.06 to 3.5 m s$^{-1}$ from which we estimate ejection times from October to December 2015. The number of fragments with effective radii $\gtrsim$20 m follows a differential power law with index $\gamma$ = -3.6$\pm$0.6, while smaller fragments are less abundant than expected from an extrapolation of this power-law. We argue that, in addition to losses due to observational selection, torques from anisotropic outgassing are capable of destroying the small fragments by driving them quickly to rotational instability. Specifically, the spin-up times of fragments $\lesssim$20 m in radius are shorter than the time elapsed since ejection from the parent nucleus. The effective radius of the parent nucleus is $r_e \le$ 275 m (geometric albedo 0.04 assumed). This is about seven times smaller than previous estimates and results in a nucleus mass at least 300 times smaller than previously thought. The mass in solid pieces, $2\times10^9$ kg, is about 4% of the mass of the parent nucleus. As a result of its small size, the parent nucleus also has a short spin-up time. Brightness variations in time-resolved nucleus photometry are consistent with rotational instability playing a role in the release of fragments.

We describe 2016 January to April observations of the fragments of 332P/Ikeya-Murakami, a comet earlier observed in a 2010 October outburst (Ishiguro et al 2014). We present photometry of the fragments, and perform simulations to infer the time of breakup. We argue that the eastern-most rapidly brightening fragment ($F4$) best corresponds to the original nucleus, rather than the initial bright fragment $F1$. We compute radial and tangential non-gravitational parameters, $A_1 = (1.5 \pm 0.4) \times 10^{-8}$ AU day$^{-2}$ and $(7.2 \pm 1.9) \times 10^{-9}$ AU day$^{-2}$; both are consistent with zero at the $4\sigma$ level. Monte Carlo simulations indicate that the fragments were emitted on the outbound journey well after the 2010 outburst, with bright fragment $F1$ splitting in mid--2013 and the fainter fragments within months of the 2016 January recovery. Western fragment $F7$ is the oldest, dating from 2011. We suggest that the delayed onset of the splitting is consistent with a self-propagating crystallization of water ice.

LOFT-P is a concept for a NASA Astrophysics Probe-Class (<$1B) X-ray timing mission, based on the LOFT concept originally proposed to ESAs M3 and M4 calls. LOFT-P requires very large collecting area (>6 m^2, >10x RXTE), high time resolution, good spectral resolution, broad-band spectral coverage (2-30 keV), highly flexible scheduling, and an ability to detect and respond promptly to time-critical targets of opportunity. It addresses science questions such as: What is the equation of state of ultra dense matter? What are the effects of strong gravity on matter spiraling into black holes? It would be optimized for sub-millisecond timing to study phenomena at the natural timescales of neutron star surfaces and black hole event horizons and to measure mass and spin of black holes. These measurements are synergistic to imaging and high-resolution spectroscopy instruments, addressing much smaller distance scales than are possible without very long baseline X-ray interferometry, and using complementary techniques to address the geometry and dynamics of emission regions. A sky monitor (2-50 keV) acts as a trigger for pointed observations, providing high duty cycle, high time resolution monitoring of the X-ray sky with ~20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger studies. A probe-class mission concept would employ lightweight collimator technology and large-area solid-state detectors, technologies which have been recently greatly advanced during the ESA M3 study. Given the large community interested in LOFT (>800 supporters, the scientific productivity of this mission is expected to be very high, similar to or greater than RXTE (~2000 refereed publications). We describe the results of a study, recently completed by the MSFC Advanced Concepts Office, that demonstrates that LOFT-P is feasible within a NASA probe-class mission budget.

V404 Cygni was discovered in 1989 by the $Ginga$ X-ray satellite during its only previously observed X-ray outburst and soon after confirmed as a black hole binary. On June 15, 2015, the Gamma Ray Burst Monitor (GBM) triggered on a new outburst of V404 Cygni. We present 13 days of GBM observations of this outburst including Earth occultation flux measurements, spectral and temporal analysis. The Earth occultation fluxes reached 30 Crab with detected emission to 100 keV and determined, via hardness ratios, that the source was in a hard state. At high luminosity, spectral analysis between 8 and 300 keV showed that the electron temperature decreased with increasing luminosity. This is expected if the protons and electrons are in thermal equilibrium during an outburst with the electrons cooled by the Compton scattering of softer seed photons from the disk. However, the implied seed photon temperatures are unusually high, suggesting a contribution from another source, such as the jet. No evidence of state transitions is seen during this time period. The temporal analysis reveals power spectra that can be modeled with two or three strong, broad Lorentzians, similar to the power spectra of black hole binaries in their hard state.

Jupiter-family comet 15P/Finlay has been reportedly quiet in activity for over a century but has harbored two outbursts during its 2014/2015 perihelion passage. Here we present an analysis of these two outbursts using a set of cometary observations. The outbursts took place between 2014 Dec. 15.4--16.0 UT and 2015 Jan. 15.5--16.0 UT as constrained by ground-based and spacecraft observations. We find a characteristic ejection speed of $V_0=300$ to $650 \mathrm{m \cdot s^{-1}}$ for the ejecta of the first outburst and $V_0=550$ to $750 \mathrm{m \cdot s^{-1}}$ for that of the second outburst using a Monte Carlo dust model. The mass of the ejecta is calculated to be $M_\mathrm{d}=2$ to $3\times10^5 \mathrm{kg}$ for the first outburst and $M_\mathrm{d}=4$ to $5\times10^5 \mathrm{kg}$ for the second outburst, corresponds to less than $10^{-7}$ of the nucleus mass. The specific energy of the two outbursts is found to be $0.3$ to $2\times10^5 \mathrm{J \cdot kg^{-1}}$. We also revisit the long-standing puzzle of the non-detection of the hypothetical Finlayid meteor shower by performing a cued search using the 13-year data from the Canadian Meteor Orbit Radar, which does not reveal any positives. The Earth will pass the 2014/2015 outburst ejecta around 2021 Oct. 6 at 22 h UT to Oct. 7 at 1 h UT, with a chance for some significant meteor activity in the radio range, which may provide further clues to the Finlayid puzzle. A southerly radiant in the constellation of Ara will favor the observers in the southern tip of Africa.

We present studies of C/2015 D1 (SOHO), the first sunskirting comet ever seen from ground stations over the past half century. The Solar and Heliospheric Observatory (SOHO) witnessed its peculiar light curve with a huge dip followed by a flareup around perihelion: the dip was likely caused by sublimation of olivines, directly evidenced by a coincident temporary disappearance of the tail. The flareup likely reflects a disintegration event, which we suggest was triggered by intense thermal stress established within the nucleus interior. Photometric data reveal an increasingly dusty coma, indicative of volatile depletion. A catastrophic mass loss rate of $\sim$10$^{5}$ kg s$^{-1}$ around perihelion was seen. Ground-based Xingming Observatory spotted the post-perihelion debris cloud. Our morphological simulations of post-perihelion images find newly released dust grains of size $a \gtrsim 10$ $\mu$m in radius, however, a temporal increase in $a_{\min}$ was also witnessed, possibly due to swift dispersions of smaller grains swept away by radiation forces without replenishment. Together with the fading profile of the light curve, a power law dust size distribution with index $\gamma = 3.2 \pm 0.1$ is derived. We detected no active remaining cometary nuclei over $\sim$0.1 km in radius in post-perihelion images acquired at Lowell Observatory. Applying radial non-gravitational parameter, $\mathcal{A}_{1} = \left(1.209 \pm 0.118 \right) \times 10^{-6}$ AU day$^{-2}$, from an isothermal water-ice sublimation model to the SOHO astrometry significantly reduces residuals and sinusoidal trends in the orbit determination. The nucleus mass $\sim$10$^{8}$--10$^{9}$ kg, and the radius $\sim$50--150 m (bulk density $\rho_{\mathrm{d}} = 0.4$ g cm$^{-3}$ assumed) before the disintegration are deduced from the photometric data; consistent results were determined from the non-gravitational effects.

It is speculated that some weakly active comets may be transitional objects between active and dormant comets. These objects are at a unique stage of the evolution of cometary nuclei, as they are still identifiable as active comets, in contrast to inactive comets that are observationally indistinguishable from low albedo asteroids. In this paper, we present a synthesis of comet and meteor observations of Jupiter-family comet 209P/LINEAR, one of the most weakly active comets recorded to-date. Images taken by the Xingming 0.35-m telescope and the Gemini Flamingo-2 camera are modeled by a Monte Carlo dust model, which yields a low dust ejection speed ($1/10$ of that of moderately active comets), dominance of large dust grains, and a low dust production of $0.4~\mathrm{kg \cdot s^{-1}}$ at 19~d after the 2014 perihelion passage. We also find a reddish nucleus of 209P/LINEAR that is similar to D-type asteroids and most Trojan asteroids. Meteor observations with the Canadian Meteor Orbit Radar (CMOR), coupled with meteoroid stream modeling, suggest a low dust production of the parent over the past few hundred orbits, although there are hints of a some temporary increase in activity in the 18th century. Dynamical simulations indicate 209P/LINEAR may have resided in a stable near-Earth orbit for $\sim 10^4$~yr, which is significantly longer than typical JFCs. All these lines of evidence imply that 209P/LINEAR as an aging comet quietly exhausting its remaining near surface volatiles. We also compare 209P/LINEAR to other low activity comets, where evidence for a diversity of the origin of low activity is seen.

The High Altitude Water Cherenkov (HAWC) Observatory continuously observes gamma-rays between 100 GeV to 100 TeV in an instantaneous field of view of about 2 steradians above the array. The large amount of raw data, the importance of small number statistics, the large dynamic range of gamma-ray signals in time (1 - $10^8$ sec) and angular extent (0.1 - 100 degrees), and the growing need to directly compare results from different observatories pose some special challenges for the analysis of HAWC data. To address these needs, we have designed and implemented a modular analysis framework based on the method of maximum likelihood. The framework facilitates the calculation of a binned Poisson Log-likelihood value for a given physics model (i.e., source model), data set, and detector response. The parameters of the physics model (sky position, spectrum, angular extent, etc.) can be optimized through a likelihood maximization routine to obtain a best match to the data. In a similar way, the parameters of the detector response (absolute pointing, angular resolution, etc.) can be optimized using a well-known source such as the Crab Nebula. The framework was designed concurrently with the Multi-Mission Maximum Likelihood (3ML) architecture, and allows for the definition of a general collection of sources with individually varying spectral and spatial morphologies. Compatibility with the 3ML architecture allows to easily perform powerful joint fits with other observatories. In this contribution, we overview the design and capabilities of the HAWC analysis framework, stressing the overarching design points that have applicability to other astronomical and cosmic-ray observatories.

Astrophysical sources are now observed by many different instruments at different wavelengths, from radio to high-energy gamma-rays, with an unprecedented quality. Putting all these data together to form a coherent view, however, is a very difficult task. Each instrument has its own data format, software and analysis procedure, which are difficult to combine. It is for example very challenging to perform a broadband fit of the energy spectrum of the source. The Multi-Mission Maximum Likelihood framework (3ML) aims to solve this issue, providing a common framework which allows for a coherent modeling of sources using all the available data, independent of their origin. At the same time, thanks to its architecture based on plug-ins, 3ML uses the existing official software of each instrument for the corresponding data in a way which is transparent to the user. 3ML is based on the likelihood formalism, in which a model summarizing our knowledge about a particular region of the sky is convolved with the instrument response and compared to the corresponding data. The user can choose between a frequentist analysis, and a Bayesian analysis. In the former, parameters of the model are optimized in order to obtain the best match to the data (i.e., the maximum of the likelihood). In the latter, the priors specified by the user are used to build the posterior distribution, which is then sampled with Markov Chain Monte Carlo or Multinest. Our implementation of this idea is very flexible, allowing the study of point sources as well as extended sources with arbitrary spectra. We will review the problem we aim to solve, the 3ML concepts and its innovative potential.

We present initial observations of the newly-discovered active asteroid 313P/Gibbs (formerly P/2014 S4), taken to characterize its nucleus and comet-like activity. The central object has a radius $\sim$0.5 km (geometric albedo 0.05 assumed). We find no evidence for secondary nuclei and set (with qualifications) an upper limit to the radii of such objects near 25 m, assuming the same albedo. Both aperture photometry and a morphological analysis of the ejected dust show that mass-loss is continuous at rates $\sim$0.2 to 0.4 kg s$^{-1}$, inconsistent with an impact origin. Large dust particles, with radii $\sim$50 to 100 $\mu$m, dominate the optical appearance. At 2.4 AU from the Sun, the surface equilibrium temperatures are too low for thermal or desiccation stresses to be responsible for the ejection of dust. No gas is spectroscopically detected (limiting the gas mass loss rate to $<$1.8 kg s$^{-1}$). However, the protracted emission of dust seen in our data and the detection of another episode of dust release near perihelion, in archival observations from 2003, are highly suggestive of an origin by the sublimation of ice. Coincidentally, the orbit of 313P/Gibbs is similar to those of several active asteroids independently suspected to be ice sublimators, including P/2012 T1, 238P/Read and 133P/Elst-Pizarro, suggesting that ice is abundant in the outer asteroid belt.

The Kreutz family of sungrazing comets contains over 2~000 known members, many of which are believed to be under $\sim 100$~m sizes (\textit{mini~comets}) and have only been studied at small heliocentric distances ($r_\mathrm{H}$) with space-based SOHO/STEREO spacecraft. To understand the brightening process of mini Kreutz comets, we conducted a survey using CFHT/MegaCam at moderate $r_\mathrm{H}$ guided by SOHO/STEREO observations. We identify two comets that should be in our search area but are not detected, indicating that the comets have either followed a steeper brightening rate within the previously-reported rapid brighten stage (the \textit{brightening~burst}), or the brightening burst starts earlier than expected. We present a composite analysis of pre-perihelion light-curves of five Kreutz comets that covered to $\sim 1$~AU. We observe a significant diversity in the light-curves that can be used to grossly classify them into two types: C/Ikeya-Seki and C/SWAN follow the canonical $r_\mathrm{H}^{-4}$ while the others follow $r_\mathrm{H}^{-7}$. In particular, C/SWAN seems to have undergone an outburst ($\Delta m>5$~mag) or a rapid brightening ($n\gtrsim11$) between $r_\mathrm{H}=1.06$~AU to $0.52$~AU, and shows hints of structural/compositional differences compared to other bright Kreutz comets. We also find evidence that the Kreutz comets as a population lose their mass less efficiently than the dynamically new comet, C/ISON, and are relatively devoid of species that drive C/ISON's activity at large $r_\mathrm{H}$. Concurrent observations of C/STEREO in different wavelengths also suggest that a blue-ward species such as CN may be the main driver for brightening burst instead of previously-thought sodium.

The dynamically new comet, C/2013 A1 (Siding Spring), is to make a close approach to Mars on 2014 October 19 at 18:30 UT at a distance of 40+/-1 Martian radius. Such extremely rare event offers a precious opportunity for the spacecrafts on Mars to closely study a dynamically new comet itself as well as the planet-comet interaction. Meanwhile, the high speed meteoroids released from C/Siding Spring also pose a threat to physically damage the spacecrafts. Here we present our observations and modeling results of C/Siding Spring to characterize the comet and assess the risk posed to the spacecrafts on Mars. We find that the optical tail of C/Siding Spring is dominated by larger particles at the time of the observation. Synchrone simulation suggests that the comet was already active in late 2012 when it was more than 7 AU from the Sun. By parameterizing the dust activity with a semi-analytic model, we find that the ejection speed of C/Siding Spring is comparable to comets such as the target of the Rosetta mission, 67P/Churyumov-Gerasimenko. Under nominal situation, the simulated dust cone will miss the planet by about 20 Martian radius. At the extreme ends of uncertainties, the simulated dust cone will engulf Mars, but the meteoric influx at Mars is still comparable to the nominal sporadic influx, seemly indicating that intense and enduring meteoroid bombardment due to C/Siding Spring is unlikely. Further simulation also suggests that gravitational disruption of the dust tail may be significant enough to be observable at Earth.

The large majority of extragalactic very high energy (VHE; E>100 GeV) sources belongs to the class of active galactic nuclei (AGN), in particular the BL Lac sub-class. AGNs are characterized by an extremely bright and compact emission region, powered by a super-massive black hole (SMBH) and an accretion disk, and relativistic outflows (jets) detected all across the electro-magnetic spectrum. In BL Lac sources the jet axis is oriented close to the line of sight, giving rise to a relativistic boosting of the emission. In radio galaxies, on the other hand, the jet makes a larger angle to the line of sight allowing to resolve the central core and the jet in great details. The giant radio galaxy M 87 with its proximity (1 6Mpc) and its very massive black hole ((3-6) x 10^9 M_solar) provides a unique laboratory to investigate VHE emission in such objects and thereby probe particle acceleration to relativistic energies near SMBH and in jets. M 87 has been established as a VHE emitter since 2005. The VHE emission displays strong variability on time-scales as short as a day. It has been subject of a large joint VHE and multi-wavelength (MWL) monitoring campaign in 2008, where a rise in the 43 GHz VLBA radio emission of the innermost region (core) was found to coincide with a flaring activity at VHE. This had been interpreted as a strong indication that the VHE emission is produced in the direct vicinity of the SMBH black hole. In 2010 again a flare at VHE was detected triggering further MWL observations with the VLBA, Chandra, and other instruments. At the same time M 87 was also observed with the Fermi-LAT telescope at GeV energies and the European VLBI Network (EVN). In this contribution preliminary results from the campaign will be presented.

We describe a Chandra X-ray ToO project designed to isolate the site of TeV flaring in the radio galaxy M87. To date, we have triggered the Chandra observations only once (2010 April) and by the time of the first of our 9 observations, the TeV flare had ended. However, we found that the X-ray intensity of the unresolved nucleus was at an elevated level for our first observation. Of the more than 60 Chandra observations we have made of the M87 jet covering 9 years, the nucleus was measured at a comparably high level only 3 times. Two of these occasions can be associated with TeV flaring, and at the time of the third event, there were no TeV monitoring activities. From the rapidity of the intensity drop of the nucleus, we infer that the size of the emitting region is of order a few light days times the unknown beaming factor; comparable to the same sort of estimate for the TeV emitting region. We also find evidence of spectral evolution in the X-ray band which seems consistent with radiative losses affecting the non-thermal population of the emitting electrons within the unresolved nucleus.

The multiwavelength observation of the nearby radio galaxy M87 provides a unique opportunity to study in detail processes occurring in Active Galactic Nuclei from radio waves to TeV gamma-rays. Here we report the detection of gamma-ray emission above 250 GeV from M87 in spring 2007 with the VERITAS atmospheric Cherenkov telescope array and discuss its correlation with the X-ray emission. The gamma-ray emission is measured to be point-like with an intrinsic source radius less than 4.5 arcmin. The differential energy spectrum is fitted well by a power-law function: dPhi/dE=(7.4+-1.3_{stat}+-1.5_{sys})(E/TeV)^{-2.31+-0.17_{stat}+-0.2_{sys}} 10^{-9}m^{-2}s^{-1}TeV^{-1}. We show strong evidence for a year-scale correlation between the gamma-ray flux reported by TeV experiments and the X-ray emission measured by the ASM/RXTE observatory, and discuss the possible short-time-scale variability. These results imply that the gamma-ray emission from M87 is more likely associated with the core of the galaxy than with other bright X-ray features in the jet.