Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5\sigma$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3\sigma$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science.

We introduce a new capability of the Neil Gehrels Swift Observatory, to provide event data from the Burst Alert Telescope (BAT) on demand in response to transients detected by other instruments. These event data are not continuously available due to the large telemetry load, but are critical to recovering weak or sub-threshold GRBs that are not triggered onboard, such as the likely counterparts to GW-detected off-axis binary neutron star mergers. We show that the availability of event data can effectively increase the rate of detections, and arcminute localizations, of GRB 170817-like bursts by >400%. We describe a spacecraft commanding pipeline purpose built to enable this science; to our knowledge the first fully autonomous extremely-low-latency commanding of a space telescope for scientific purposes. This pipeline has been successfully run in its complete form since early 2020, and has resulted in the recovery of BAT event data for >700 externally triggered events to date (GWs, neutrinos, GRBs triggered by other facilities, FRBs), now running with a success rate of ~90%. We exemplify the utility of this new capability by using the resultant data to (1) set the most sensitive (8 sigma) upper limits of $8.1\times10^{-8}$ erg cm$^{-2}$ s$^{-1}$ (14-195 keV) on prompt 1s duration short GRB-like emission within $\pm$ 15s of the unmodelled GW burst candidate S200114f, and (2) provide arcminute localizations for short GRB 200325A and other bursts. We show that using data from GUANO to localize GRBs discovered elsewhere, we can increase the net rate of arcminute localized GRBs by 10-20% per year. Along with the scientific yield of more sensitive searches for sub-threshold GRBs, the new capabilities designed for this project will enable further rapid response Target of Opportunity capabilities for Swift, and have implications for the design of future rapid-response space telescopes.

We present the Gravitational Wave Treasure Map, a tool to coordinate, visualize, and assess the electromagnetic follow-up of gravitational wave (GW) events. With typical GW localization regions of hundreds to thousands of square degrees and dozens of active follow-up groups, the pursuit of electromagnetic (EM) counterparts is a challenging endeavor, but the scientific payoff for early discovery of any counterpart is clear. With this tool, we provide a website and API interface that allows users to easily see where other groups have searched and better inform their own follow-up search efforts. A strong community of Treasure Map users will increase the overall efficiency of EM counterpart searches and will play a fundamental role in the future of multi-messenger astronomy.

The IceCube report of a $\sim 3.5\sigma$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 05056+056 in 2014-2015 and the 2017 detection of a high-energy neutrino, IceCube-170922A, during a gamma-ray flare from the same blazar, have revived the interest in scenarios for neutrino production in blazars. We perform comprehensive analyses on the long-term electromagnetic emission of TXS 05056+056 using optical, X-ray, and gamma-ray data from the All-Sky Automated Survey for Supernovae (ASAS-SN), the Neil Gehrels Swift Observatory (Swift), the Monitor of All-sky X-ray Image (MAXI), and the Fermi Large Area Telescope (Fermi-LAT). We also perform numerical modeling of the spectral energy distributions (SEDs) in four epochs prior to 2017 with contemporaneous gamma-ray and lower energy (optical and/or X-ray) data. We find that the multi-epoch SEDs are consistent with a hybrid leptonic scenario, where the gamma-rays are produced in the blazar zone via external inverse Compton scattering of accelerated electrons, and high-energy neutrinos are produced via the photomeson production process of co-accelerated protons. The multi-epoch SEDs can be satisfactorily explained with the same jet parameters and variable external photon density and electron luminosity. Using the maximal neutrino flux derived for each epoch, we put an upper limit of $\sim0.4-2$ on the muon neutrino number in ten years of IceCube observations. Our results are consistent with the IceCube-170922A detection, which can be explained as an upper fluctuation from the average neutrino rate expected from the source, but in strong tension with the 2014-2015 neutrino flare.

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.

The Neil Gehrels Swift Observatory carried out prompt searches for gravitational wave (GW) events detected by the LIGO/Virgo Collaboration (LVC) during the second observing run ("O2"). Swift performed extensive tiling of eight LVC triggers, two of which had very low false-alarm rates (GW 170814 and the epochal GW 170817), indicating a high confidence of being astrophysical in origin; the latter was the first GW event to have an electromagnetic counterpart detected. In this paper we describe the follow-up performed during O2 and the results of our searches. No GW electromagnetic counterparts were detected; this result is expected, as GW 170817 remained the only astrophysical event containing at least one neutron star after LVC's later retraction of some events. A number of X-ray sources were detected, with the majority of identified sources being active galactic nuclei. We discuss the detection rate of transient X-ray sources and their implications in the O2 tiling searches. Finally, we describe the lessons learned during O2, and how these are being used to improve the \swift\ follow-up of GW events. In particular, we simulate a population of GRB afterglows to evaluate our source ranking system's ability to differentiate them from unrelated and uncatalogued X-ray sources. We find that $\approx$60-70% of afterglows whose jets are oriented towards Earth will be given high rank (i.e., "interesting" designation) by the completion of our second follow-up phase (assuming their location in the sky was observed), but that this fraction can be increased to nearly 100% by performing a third follow-up observation of sources exhibiting fading behavior.

The 2020s will be the most data-rich decade of astronomy in history. As the scale and complexity of our surveys increase, the problem of scheduling becomes more critical. We must develop high-quality scheduling approaches, implement them as open-source software, and begin linking the typically separate stages of observation and data analysis.

Electromagnetic observations of the sky have been the basis for our study of the Universe for millennia, cosmic ray studies are now entering their second century, the first neutrinos from an astrophysical source were identified three decades ago, and gravitational waves were directly detected only four years ago. Detections of these messengers are now common. Astrophysics will undergo a revolution in the 2020s as multimessenger detections become routine. The 8th Astro2020 Thematic Area is Multimessenger Astronomy and Astrophysics, which includes the identification of the sources of gravitational waves, astrophysical and cosmogenic neutrinos, cosmic rays, and gamma-rays, and the coordinated multimessenger and multiwavelength follow-ups. Identifying and characterizing multimessenger sources enables science throughout and beyond astrophysics. Success in the multimessenger era requires: (i) sensitive coverage of the non-electromagnetic messengers, (ii) full coverage of the electromagnetic spectrum, with either fast-response observations or broad and deep high-cadence surveys, and (iii) improved collaboration, communication, and notification platforms.

Neutron star mergers, referring to both binary neutron star and neutron star black hole mergers, are the canonical multimessenger events. They have been detected across the electromagnetic spectrum, have recently been detected in gravitational waves, and are likely to produce neutrinos over several decades in energy. The non-thermal prompt and afterglow emission of short gamma-ray bursts and the quasi-thermal emission from the radioactively powered kilonovae provide distinct insights into the physics of neutron star mergers. When combined with direct information on coalescence from gravitational waves and neutrinos these sources may become the best understood astrophysical transients. Multimessenger observations of these cataclysmic events will determine sources of gravitational waves and astrophysical neutrinos, enable precision cosmology, and unique tests of fundamental physics, the origin of heavy elements, the behavior of relativistic jets, and the equation of state of supranuclear matter. In this white paper we present a summary of the science discoveries possible with multimessenger observations of neutron star mergers and provide recommendations to enable them in the new era of time-domain, multimessenger astronomy.

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.

The bright transient AT2018cow has been unlike any other known type of transient. Its high brightness, rapid rise and decay and initially nearly featureless spectrum are unprecedented and difficult to explain using models for similar burst sources. We present evidence for faint gamma-ray emission continuing for at least 8 days, and featureless spectra in the ultraviolet bands -- both unusual for eruptive sources. The X-ray variability of the source has a burst-like character. The UV-optical spectrum does not show any CNO line but is well described by a blackbody. We demonstrate that a model invoking the tidal disruption of a 0.1 - 0.4 Msun Helium White Dwarf (WD) by a 100,000 to one million solar mass Black Hole (BH) located in the outskirts of galaxy Z~137-068 could provide an explanation for most of the characteristics shown in the multi-wavelength observations. A blackbody-like emission is emitted from an opaque photosphere, formed by the debris of the WD disruption. Broad features showing up in the optical/infrared spectra in the early stage are probably velocity broadened lines produced in a transient high-velocity outward moving cocoon. The asymmetric optical/infrared lines that appeared at a later stage are emission from an atmospheric layer when it detached from thermal equilibrium with the photosphere, which undergoes more rapid cooling. The photosphere shrinks when its temperature drops, and the subsequent infall of the atmosphere produced asymmetric line profiles. Additionally, a non-thermal jet might be present, emitting X-rays in the 10-150 keV band.

Gamma-ray bursts (GRBs) occurring in the local Universe constitute an interesting sub-class of the GRB family, since their luminosity is on average lower than that of their cosmological analogs. We aim to contribute to the study of local bursts by reporting the case of GRB 171205A. This source was discovered by Swift Burst Alert Telescope (BAT) on 2017, December 5 and soon associated with a low redshift host galaxy (z=0.037), and an emerging SN (SN 2017iuk). We analyzed the full Swift, dataset, comprising the UV-Optical Telescope (UVOT), X-ray Telescope (XRT) and BAT data. In addition, we employed the Konus-Wind high energy data as a valuable extension at gamma-ray energies. The photometric SN signature is clearly visible in the UVOT u, b and v filters. The maximum emission is reached at ~ 13 (rest frame) days, and the whole bump resembles that of SN 2006aj, but lower in magnitude and with a shift in time of +2 d. A prebump in the v-band is also clearly visible, and this is the first time that such a feature is not observed achromatically in GRB-SNe. Its physical origin cannot be easily explained. The X-ray spectrum shows an intrinsic Hydrogen column density N_H,int = 7.4(+4.1 -3.6) X 10^20 / cm^2$, which is at the low end of the N_H,int, even considering just low redshift GRBs. The spectrum also features a thermal component, which is quite common in GRBs associated with SNe, but whose origin is still a matter of debate. Finally, the isotropic energy in the gamma-ray band, E_iso = 2.18(+0.63 -0.50) X 10^49 erg, is lower than those of cosmological GRBs. Combining this value with the peak energy in the same band, E_p=125(+141 -37) keV, implies that GRB 171205A is an outlier of the Amati relation, as are some other low redshift GRBs, and its emission mechanism should be different from that of canonical, farther away GRBs.

Detection of the IceCube-170922A neutrino coincident with the flaring blazar TXS 0506+056, the first and only 3-sigma high-energy neutrino source association to date, offers a potential breakthrough in our understanding of high-energy cosmic particles and blazar physics. We present a comprehensive analysis of TXS 0506+056 during its flaring state, using newly collected Swift, NuSTAR, and X-shooter data with Fermi observations and numerical models to constrain the blazar's particle acceleration processes and multimessenger (electromagnetic and high-energy neutrino) emissions. Accounting properly for electromagnetic cascades in the emission region, we find a physically-consistent picture only within a hybrid leptonic scenario, with gamma-rays produced by external inverse-Compton processes and high-energy neutrinos via a radiatively-subdominant hadronic component. We derive robust constraints on the blazar's neutrino and cosmic-ray emissions and demonstrate that, because of cascade effects, the 0.1-100keV emissions of TXS 0506+056 serve as a better probe of its hadronic acceleration and high-energy neutrino production processes than its GeV-TeV emissions. If the IceCube neutrino association holds, physical conditions in the TXS 0506+056 jet must be close to optimal for high-energy neutrino production, and are not favorable for ultra-high-energy cosmic-ray acceleration. Alternatively, the challenges we identify in generating a significant rate of IceCube neutrino detections from TXS 0506+056 may disfavor single-zone models. In concert with continued operations of the high-energy neutrino observatories, we advocate regular X-ray monitoring of TXS 0506+056 and other blazars in order to test single-zone blazar emission models, clarify the nature and extent of their hadronic acceleration processes, and carry out the most sensitive possible search for additional multimessenger sources.

The Swift Gamma-Ray Burst Explorer is a uniquely capable mission, with three on-board instruments and rapid slewing capabilities. It serves as a fast-response satellite observatory for everything from gravitational-wave counterpart searches to cometary science. Swift averages 125 different observations per day, and is consistently over-subscribed, responding to about one-hundred Target of Oportunity (ToO) requests per month from the general astrophysics community, as well as co-pointing and follow-up agreements with many other observatories. Since launch in 2004, the demands put on the spacecraft have grown consistently in terms of number and type of targets as well as schedule complexity. To facilitate this growth, various scheduling tools and helper technologies have been built by the Swift team to continue improving the scientific yield of the Swift mission. However, these tools have been used only to assist humans in exploring the local pareto surface and for fixing constraint violations. Because of the computational complexity of the scheduling task, no automation tool has been able to produce a plan of equal or higher quality than that produced by a well-trained human, given the necessary time constraints. In this proceeding we formalize the Swift Scheduling Problem as a dynamic fuzzy Constraint Satisfaction Problem (DF-CSP) and explore the global solution space. We detail here several approaches towards achieving the goal of surpassing human quality schedules using classical optimization and algorithmic techniques, as well as machine learning and recurrent neural network (RNN) methods. We then briefly discuss the increased scientific yield and benefit to the wider astrophysics community that would result from the further development and adoption of these technologies.

$\textit{Swift}$'s rapid slewing, flexible planning, and multi-wavelength instruments make it the most capable space-based follow-up engine for finding poorly localized sources. During O1 and O2 $\textit{Swift}$ successfully tiled hundreds of square-degrees of sky in the LVC localization regions, searching for, and identifying, possible X-ray and UV/O transients in the field. $\textit{Swift}$ made important contributions to the discovery and characterization of the kilonova AT 2017gfo, discovering the UV emission and providing the deepest X-ray upper limits in the first 24 hours after the trigger, strongly constraining the dynamics and geometry of the counterpart. $\textit{Swift}$ tiled 92% of the galaxy convolved error region down to average X-ray flux sensitivities of $10^{-12}$ erg $\mbox{cm}^{-2}$ $\mbox{s}^{-1}$, significantly increasing our confidence that AT 2017gfo is indeed the counterpart to GW 170817 and sGRB 170817. However, there remains significant room for improvement of $\textit{Swift}$'s follow-up in preparation for O3. This will take the form of both revised observation strategy based on detailed analysis of the results from O2, and significant changes to $\textit{Swift}$'s operational capabilities. These improvements are necessary both for maximizing the likelihood that $\textit{Swift}$ finds a counterpart, and minimizing the impact that follow-up activities have on other $\textit{Swift}$ science priorities. We outline areas of improvement to the observing strategy itself for optimal tiling of the LVC localization regions. We also discuss ongoing work on operational upgrades that will decrease latency in our response time, and minimize impact on pre-planned observations, while maintaining spacecraft health and safety.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).