Search

The Legacy Survey of Space and Time, to be conducted with the Vera C. Rubin Observatory, is poised to revolutionize our understanding of the Solar System by providing an unprecedented wealth of data on various objects, including the elusive interstellar objects (ISOs). Detecting and classifying ISOs is crucial for studying the composition and diversity of materials from other planetary systems. However, the rarity and brief observation windows of ISOs, coupled with the vast quantities of data to be generated by LSST, create significant challenges for their identification and classification. This study aims to address these challenges by exploring the application of machine learning algorithms to the automated classification of ISO tracklets in simulated LSST data. We employed various machine learning algorithms, including random forests (RFs), stochastic gradient descent (SGD), gradient boosting machines (GBMs), and neural networks (NNs), to classify ISO tracklets in simulated LSST data. We demonstrate that GBM and RF algorithms outperform SGD and NN algorithms in accurately distinguishing ISOs from other Solar System objects. RF analysis shows that many derived Digest2 values are more important than direct observables in classifying ISOs from the LSST tracklets. The GBM model achieves the highest precision, recall, and F1 score, with values of 0.9987, 0.9986, and 0.9987, respectively. These findings lay the foundation for the development of an efficient and robust automated system for ISO discovery using LSST data, paving the way for a deeper understanding of the materials and processes that shape planetary systems beyond our own. The integration of our proposed machine learning approach into the LSST data processing pipeline will optimize the survey's potential for identifying these rare and valuable objects, enabling timely follow-up observations and further characterization., Comment: 11 pages, 4 figures, 6 tables

Nongravitational accelerations in the absence of observed activity have recently been identified on NEOs, opening the question of the prevalence of anisotropic mass-loss in the near-Earth environment. Motivated by the necessity of nongravitational accelerations to identify 2010 VL$_{65}$ and 2021 UA$_{12}$ as a single object, we investigate the problem of linking separate apparitions in the presence of nongravitational perturbations. We find that nongravitational accelerations on the order of $10^{-9}$ au/d$^2$ can lead to a change in plane-of-sky positions of $\sim10^3$ arcsec between apparitions. Moreover, we inject synthetic tracklets of hypothetical nongravitationally-accelerating NEOs into the Minor Planet Center orbit identification algorithms. We find that at large nongravitational accelerations ($|A_i|\geq10^{-8}$ au/d$^2$) these algorithms fail to link a significant fraction of these tracklets. We further show that if orbits can be determined for both apparitions, the tracklets will be linked regardless of nongravitational accelerations, although they may be linked to multiple objects. In order to aid in the identification and linkage of nongravitationally accelerating objects, we propose and test a new methodology to search for unlinked pairs. When applied to the current census of NEOs, we recover the previously identified case but identify no new linkages. We conclude that current linking algorithms are generally robust to nongravitational accelerations, but objects with large nongravitational accelerations may potentially be missed. While current algorithms are well-positioned for the anticipated increase in the census population from future survey missions, it may be possible to find objects with large nongravitational accelerations hidden in isolated tracklet pairs., Comment: 16 pages, 5 figures and 5 appendix figures. Accepted to ApJ

The prompt spectra of gamma-ray bursts are known to follow broadband continuum behavior over decades in energy. GRB 221009A, given the moniker the brightest of all time (BOAT), is the brightest gamma-ray burst identified in half a century of observations, and was first identified by the Fermi Gamma-ray Burst Monitor (GBM). On behalf of the Fermi-GBM Team, Lesage et al. (2023) described the initial GBM analysis. Ravasio et al. (2024) report the identification of a spectral line in part of the prompt emission of this burst, which they describe as evolving over 80 s from $\sim$12 MeV to 6 MeV. We report a GBM Team analysis on the Ravasio Line: 1) We cannot identify an instrumental effect that could have produced this signal, and 2) our method of calculating the statistical significance of the line shows it easily exceeds the 5$\sigma$ discovery threshold. We additionally comment on the claim of the line beginning at earlier time intervals, up to 37 MeV, as reported in Zhang et al. (2024). We find that it is reasonable to utilize these measurements for characterization of the line evolution, with caution. We encourage theoretical studies exploring this newly discovered gamma-ray burst spectral feature, unless any rigorous alternative explanation unrelated to the emission from GRB 221009A is identified.

We present the detection and analysis of GRB 231115A, a candidate extragalactic magnetar giant flare (MGF) observed by Fermi/GBM and localized by INTEGRAL to the starburst galaxy M82. This burst exhibits distinctive temporal and spectral characteristics that align with known MGFs, including a short duration and a high peak energy. Gamma-ray analyses reveal significant insights into this burst, supporting conclusions already established in the literature: our time-resolved spectral studies provide further evidence that GRB 231115A is indeed a MGF. Significance calculations also suggest a robust association with M82, further supported by a high Bayes factor that minimizes the probability of chance alignment with a neutron star merger. Despite extensive follow-up efforts, no contemporaneous gravitational wave or radio emissions were detected. The lack of radio emission sets stringent upper limits on possible radio luminosity. Constraints from our analysis show no fast radio bursts (FRBs) associated with two MGFs. X-ray observations conducted post-burst by Swift/XRT and XMM/Newton provided additional data, though no persistent counterparts were identified. Our study underscores the importance of coordinated multi-wavelength follow-up and highlights the potential of MGFs to enhance our understanding of short GRBs and magnetar activities in the cosmos. Current MGF identification and follow-up implementation are insufficient for detecting expected counterparts; however, improvements in these areas may allow for the recovery of follow-up signals with existing instruments. Future advancements in observational technologies and methodologies will be crucial in furthering these studies.

Photometric measurements allow the determination of an asteroid's absolute magnitude, which often represents the sole means to infer its size. Photometric observations can be obtained in a variety of filters that can be unique to a specific observatory. Those observations are then calibrated into specific bands with respect to reference star catalogs. In order to combine all the different measurements for evaluation, photometric observations need to be converted to a common band, typically V-band. Current band-correction schemes in use by IAU's Minor Planet Center, JPL's Center for Near Earth Object Studies and ESA's NEO Coordination Centre use average correction values for the apparent magnitude derived from photometry of asteroids as the corrections are dependent on the typically unknown spectrum of the object to be corrected. By statistically analyzing the photometric residuals of asteroids, we develop a new photometric correction scheme that does not only consider the band, but also accounts for reference catalog and observatory. We describe a new statistical photometry correction scheme for asteroid observations with debiased corrections. Testing this scheme on a reference group of asteroids, we see a 36% reduction in the photometric residuals. Moreover, the new scheme leads to a more accurate and debiased determination of the H-G magnitude system and, in turn, to more reliable inferred sizes. We discuss the significant shift in the corrections with this "DePhOCUS" debiasing system, its limitations, and the impact for photometric and physical properties of all asteroids, especially Near-Earth Objects.

GW230529 is the first compact binary coalescence detected by the LIGO-Virgo-KAGRA collaboration with at least one component mass confidently in the lower mass-gap, corresponding to the range 3-5$M_{\odot}$. If interpreted as a neutron star-black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM) emission. However, no EM counterpart has been reported. At the merger time $t_0$, Swift-BAT and Fermi-GBM together covered 100$\%$ of the sky. Performing a targeted search in a time window $[t_0-20 \text{s},t_0+20 \text{s}]$, we report no detection by the Swift-BAT and the Fermi-GBM instruments. Combining the position-dependent $\gamma-$ray flux upper limits and the gravitational-wave posterior distribution of luminosity distance, sky localization and inclination angle of the binary, we derive constraints on the characteristic luminosity and structure of the jet possibly launched during the merger. Assuming a top-hat jet structure, we exclude at 90$\%$ credibility the presence of a jet which has at the same time an on-axis isotropic luminosity $\gtrsim 10^{48}$ erg s$^{-1}$, in the bolometric band 1 keV-10 MeV, and a jet opening angle $\gtrsim 15$ deg. Similar constraints are derived testing other assumptions about the jet structure profile. Excluding GRB 170817A, the luminosity upper limits derived here are below the luminosity of any GRB observed so far., Comment: 18 pages, 1 table, 11 figures

Since their first discovery in the late 1960s, Gamma-ray bursts have attracted an exponentially growing interest from the international community due to their central role in the most highly debated open questions of the modern research of astronomy, astrophysics, cosmology, and fundamental physics. These range from the intimate nuclear composition of high density material within the core of ultra-dense neuron stars, to stellar evolution via the collapse of massive stars, the production and propagation of gravitational waves, as well as the exploration of the early Universe by unveiling first stars and galaxies (assessing also their evolution and cosmic re-ionization). GRBs have stimulated in the past $\sim$50 years the development of cutting-edge technological instruments for observations of high energy celestial sources from space, leading to the launch and successful operations of many different scientific missions (several of them still in data taking mode nowadays). In this review, we provide a brief description of the GRB-dedicated missions from space being designed and developed for the future. The list of these projects, not meant to be exhaustive, shall serve as a reference to interested readers to understand what is likely to come next to lead the further development of GRB research and associated phenomenology., Comment: Accepted for publication on Universe. Invited review, contribution to the Universe Special Issue "Recent Advances in Gamma Ray Astrophysics and Future Perspectives", P. Romano eds. (https://www.mdpi.com/journal/universe/special_issues/7299902Z97)

We explore large-scale cosmic structure using the spatial distribution of 542 gamma-ray bursts (GRBs) having accurately measured positions and spectroscopic redshifts. Prominent cosmological clusters are identified in both the northern and southern galactic hemispheres (avoiding extinction effects in the plane of the Milky Way) using the Bootstrap Point-Radius method. The Northern Galactic hemisphere contains a significant group of four GRBs in the redshift range 0.59 < z < 0.62 (with a Bootstrap probability of p = 0.012) along with the previously-identified Hercules-Corona Borealis GreatWall (in the revised redshift range 0.9 < z < 2.1, p = 0.017). The Southern Galactic hemisphere contains the previously-identified Giant GRB Ring (p = 0.022) along with another possible cluster of 7 - 9 GRBs at 1.179 < z < 1.444 (p = 0.031). Additionally, both the Hercules-Corona Borealis Great Wall and the Giant GRB Ring have become more prominent as the GRB sample size has grown. The approach used here underscores the potential value of GRB clustering as a probe of large-scale cosmic structure, complementary to galaxy and quasar clustering. Because of the vast scale on which GRB clustering provides valuable insights, it is important that optical GRB monitoring continue so that additional spectroscopic redshift measurements should be obtained., Comment: 11 pages, 15 figures, submitted to MNRAS

We explore fitting gamma-ray burst spectra with three physically-motivated models, and thus revisit the viability of synchrotron radiation as the primary source of GRB prompt emission. We pick a sample of 100 bright GRBs observed by the Fermi Gamma-ray Burst Monitor (GBM), based on their energy flux values. In addition to the standard empirical spectral models used in previous GBM spectroscopy catalogs, we also consider three physically-motivated models; (a) a Thermal Synchrotron model, (b) a Band model with a High-energy Cutoff, and (c) a Smoothly Broken Power Law (SBPL) model with a Multiplicative Broken Power Law (MBPL). We then adopt the Bayesian information criterion (BIC) to compare the fits obtained and choose the best model. We find that 42% of the GRBs from the fluence spectra and 23% of GRBs from the peak-flux spectra have one of the three physically-motivated models as their preferred one. From the peak-flux spectral fits, we find that the low-energy index distributions from the empirical model fits for long GRBs peak around the synchrotron value of -2/3, while the two low-energy indices from the SBPL+MBPL fits of long GRBs peak close to the -2/3 and -3/2 values expected for a synchrotron spectrum below and above the cooling frequency., Comment: arXiv admin note: text overlap with arXiv:2103.13528

Magnetars are slowly rotating neutron stars that possess the strongest magnetic fields ($10^{14}-10^{15} \mathrm{G}$) known in the cosmos. They display a range of transient high-energy electromagnetic activity. The brightest and most energetic of these events are the gamma-ray bursts (GRBs) known as magnetar giant flares (MGFs), with isotropic energy $E\approx10^{44}-10^{46} \mathrm{erg}$. There are only seven detections identified as MGFs to date: three unambiguous events occurred in our Galaxy and the Magellanic Clouds, and the other four MGF candidates are associated with nearby star-forming galaxies. As all seven identified MGFs are bright at Earth, additional weaker events remain unidentified in archival data. We conducted a search of the Fermi Gamma-ray Burst Monitor (GBM) database for candidate extragalactic MGFs and, when possible, collected localization data from the Interplanetary Network (IPN) satellites. Our search yielded one convincing event, GRB 180128A. IPN localizes this burst with NGC 253, commonly known as the Sculptor Galaxy. This event is the second MGF in modern astronomy to be associated with this galaxy and the first time two bursts are associated with a single galaxy outside our own. Here, we detail the archival search criteria that uncovered this event and its spectral and temporal properties, which are consistent with expectations for a MGF. We also discuss the theoretical implications and finer burst structures resolved from various binning methods. Our analysis provides observational evidence for an eighth identified MGF., Comment: 15 pages, 8 figures

The second {\itshape Fermi}/LAT gamma-ray burst (GRB) catalog (2FLGC) spanning the first decade of operations by the LAT collaboration was recently released. The closure relations of the synchrotron forward shock (FS) model are not able to reproduce a sizeable portion of the afterglow-phase light curves in this collection, indicating that there may be a large contribution from some other mechanism. Recently, synchrotron self-Compton (SSC) light curves from the reverse shock (RS) regions were derived in the thick- and thin-shell regime for a constant-density medium, and it was demonstrated that analytical light curves could explain the~GeV flare observed in several bursts from 2FLGC, including GRB 160509A. Here, we generalise the SSC RS scenario from the constant density to a stratified medium, and show that this contribution helps to describe the early light curves exhibited in some {\itshape Fermi}/LAT-detected bursts. As a particular case, we model a sample of eight bursts that exhibited a short-lasting emission with the synchrotron and SSC model from FS and RS regions, evolving in a stellar-wind environment, constraining the microphysical parameters, the circumburst density, the bulk Lorentz factor, and the fraction of shock-accelerated electrons. We demonstrate that the highest-energy photons can only be described by the SSC from the forward-shock region., Comment: 35 pages, 12 figures, 4 tables, accepted in mnras

We describe enhancements to the digest2 software, a short-arc orbit classifier for heliocentric orbits. Digest2 is primarily used by the Near-Earth Object (NEO) community to flag newly discovered objects for a immediate follow-up and has been a part of NEO discovery process for more than 15 years. We have updated the solar system population model used to weight the digest2 score according to the 2023 catalog of known solar system orbits and extended the list of mean uncertainties for 140 observatory codes. Moreover, we have added Astrometry Data Exchange Standard (ADES) input format support to digest2, which provides additional information for the astrometry, such as positional uncertainties for each detection. The digest2 code was also extended to read the roving observer astrometric format as well as the ability to compute a new parameter from the provided astrometric uncertainties ($RMS'$) that can serve as an indicator of in-tracklet curvature when compared with tracklet's great-circle fit RMS. Comparison with the previous version of digest2 confirmed the improvement in accuracy of NEO identification and found that using ADES XML input significantly reduces the computation time of the digest2., Comment: 14 pages, 8 tables, 6 figures, appendix

We present the discovery of the radio afterglow of the short $\gamma$-ray burst (GRB) 210726A, localized to a galaxy at a photometric redshift of $z\sim 2.4$. While radio observations commenced $\lesssim 1~$day after the burst, no radio emission was detected until $\sim11~$days. The radio afterglow subsequently brightened by a factor of $\sim 3$ in the span of a week, followed by a rapid decay (a "radio flare"). We find that a forward shock afterglow model cannot self-consistently describe the multi-wavelength X-ray and radio data, and underpredicts the flux of the radio flare by a factor of $\approx 5$. We find that the addition of substantial energy injection, which increases the isotropic kinetic energy of the burst by a factor of $\approx 4$, or a reverse shock from a shell collision are viable solutions to match the broad-band behavior. At $z\sim 2.4$, GRB 210726A is among the highest redshift short GRBs discovered to date as well as the most luminous in radio and X-rays. Combining and comparing all previous radio afterglow observations of short GRBs, we find that the majority of published radio searches conclude by $\lesssim 10~$days after the burst, potentially missing these late rising, luminous radio afterglows., Comment: 29 pages, 10 figures, accepted to ApJ

The Interplanetary Network (IPN) is a detection, localization and alert system that utilizes the arrival time of transient signals in gamma-ray detectors on spacecraft separated by planetary baselines to geometrically locate the origin of these transients. Due to the changing astrophysical landscape and the new emphasis on time domain and multi-messenger astrophysics (TDAMM) from the Pathways to Discovery in Astronomy and Astrophysics for the 2020s, this Gamma-ray Transient Network Science Analysis Group was tasked to understand the role of the IPN and high-energy monitors in this new era. The charge includes describing the science made possible with these facilities, tracing the corresponding requirements and capabilities, and highlighting where improved operations of existing instruments and the IPN would enhance TDAMM science. While this study considers the full multiwavelength and multimessenger context, the findings are specific to space-based high-energy monitors. These facilities are important both for full characterization of these transients as well as facilitating follow-up observations through discovery and localization. The full document reports a brief history of this field, followed by our detailed analyses and findings in some 68 pages, providing a holistic overview of the role of the IPN and high-energy monitors in the coming decades., Comment: Terms of Reference and additional information on the Science Analysis Group are available at https://pcos.gsfc.nasa.gov/sags/gtn-sag.php

The presence of delayed GeV emission after a strong transient, such as a GRB (Gamma-Ray Burst), in the VHE (Very-High Energy, $E>100$ GeV) band can be the signature of a non-zero magnetic field in the intergalactic medium. We used a synchrotron self-Compton multiwavelength model to infer an analytical description of the intrinsic VHE spectrum (corrected for absorption by the Extragalactic Background Light, EBL) of GRB$\,$190114C to predict the lightcurves and SEDs of the delayed emission with Monte Carlo simulations for different IGMF (Intergalactic Magnetic Field) configurations (strengths $B=8\times10^{-21}$ G, $10^{-20}$ G, $3\times 10^{-20}$G and correlation length $\lambda>1$ Mpc), and compared them with the Fermi-LAT (Fermi Large Area Telescope) limits computed for several exposure times. We found that Fermi LAT is not sensitive enough to constrain any IGMF strengths using GRB$\,$190114C., Comment: 8 pages, 3 figures, accepted for publication to Physical Review D

GRB 221009A has been referred to as the Brightest Of All Time (the BOAT). We investigate the veracity of this statement by comparing it with a half century of prompt gamma-ray burst observations. This burst is the brightest ever detected by the measures of peak flux and fluence. Unexpectedly, GRB 221009A has the highest isotropic-equivalent total energy ever identified, while the peak luminosity is at the $\sim99$th percentile of the known distribution. We explore how such a burst can be powered and discuss potential implications for ultra-long and high-redshift gamma-ray bursts. By geometric extrapolation of the total fluence and peak flux distributions GRB 221009A appears to be a once in 10,000 year event. Thus, while it almost certainly not the BOAT over all of cosmic history, it may be the brightest gamma-ray burst since human civilization began., Comment: Version accepted to ApJL. Also adds proper acknowledgements

We present the IXPE observation of GRB 221009A which includes upper limits on the linear polarization degree of both prompt and afterglow emission in the soft X-ray energy band. GRB 221009A is an exceptionally bright gamma-ray burst (GRB) that reached Earth on 2022 October 9 after travelling through the dust of the Milky Way. The Imaging X-ray Polarimetry Explorer (IXPE) pointed at GRB 221009A on October 11 to observe, for the first time, the 2-8 keV X-ray polarization of a GRB afterglow. We set an upper limit to the polarization degree of the afterglow emission of 13.8% at a 99% confidence level. This result provides constraints on the jet opening angle and the viewing angle of the GRB, or alternatively, other properties of the emission region. Additionally, IXPE captured halo-rings of dust-scattered photons which are echoes of the GRB prompt emission. The 99% confidence level upper limit to the prompt polarization degree depends on the background model assumption and it ranges between ~55% to ~82%. This single IXPE pointing provides both the first assessment of X-ray polarization of a GRB afterglow and the first GRB study with polarization observations of both the prompt and afterglow phases., Comment: Accepted by Astrophysical Journal Letters

We report statistically significant detections of non-radial nongravitational accelerations based on astrometric data in the photometrically inactive objects 1998 KY$_{26}$, 2005 VL$_1$, 2016 NJ$_{33}$, 2010 VL$_{65}$, 2016 RH$_{120}$, and 2010 RF$_{12}$. The magnitudes of the nongravitational accelerations are greater than those typically induced by the Yarkovsky effect and there is no radiation-based, non-radial effect that can be so large. Therefore, we hypothesize that the accelerations are driven by outgassing, and calculate implied H$_2$O production rates for each object. We attempt to reconcile outgassing induced acceleration with the lack of visible comae or photometric activity via the absence of surface dust and low levels of gas production. Although these objects are small and some are rapidly rotating, surface cohesive forces are stronger than the rotational forces and rapid rotation alone cannot explain the lack of surface debris. It is possible that surface dust was removed previously, perhaps via outgassing activity that increased the rotation rates to their present day value. We calculate dust production rates of order $\sim10^{-4}$ g s$^{-1}$ in each object assuming that the nuclei are bare, within the upper limits of dust production from a sample stacked image of 1998 KY$_{26}$ of $\dot{M}_{\rm Dust}<0.2$ g s$^{-1}$. This production corresponds to brightness variations of order $\sim0.0025\%$, which are undetectable in extant photometric data. We assess the future observability of each of these targets, and find that the orbit of 1998 KY$_{26}$ -- which is also the target for the extended Hayabusa2 mission -- exhibits favorable viewing geometry before 2025., Comment: 21 pages, 9 figures, 3 tables, accepted for publication in the Planetary Science Journal

Gamma-ray bursts (GRBs) are the most powerful explosions in the universe. How efficiently the jet converts its energy to radiation is a long-standing problem and it is poorly constrained. The standard model invokes a relativistic fireball with a bright photosphere emission component. A definitive diagnosis of GRB radiation components and measurement of GRB radiative efficiency require prompt emission and afterglow data with high-resolution and wide-band coverage in time and energy. Here we report a comprehensive temporal and spectral analysis of the TeV-emitting bright GRB 190114C. Its fluence is one of the highest of all GRBs detected so far, which allows us to perform a high-resolution study of the prompt emission spectral properties and their temporal evolution down to a timescale of about 0.1 s. We observe that each of the initial pulses has a thermal component contributing $\sim20\%$ of the total energy, the corresponding temperature and the inferred Lorentz factor of the photosphere evolve following broken power-law shapes. From the observation of the non-thermal spectra and the light curve, the onset of afterglow corresponding to the deceleration of the fireball is considered at $\sim 6$~s. By incorporating the thermal and the non-thermal observations, as well as the photosphere and the synchrotron radiative mechanisms, we can directly derive the fireball energy budget with little dependence on hypothetical parameters and to measure a $\sim 16\%$ radiative efficiency for this GRB. With the fireball energy budget derived, the afterglow microphysics parameters can also be constrained directly from the data., Comment: 27 pages, 8 figures (including 14 panels), 6 tables, accepted for publication in The Astrophysical Journal Letters

2009 MS9 is a trans-Neptunian object (TNO) whose perihelion brings it close to the distance where some long period comets are seen to become active. Knowing this, and the fact that this object appears to brighten in excess of it's predicted nucleus brightness suggests that 2009 MS9 has a delayed onset of activity brought on by the sublimation of a species more volatile than water. In this paper we characterize 2009 MS9's physical properties and investigate potential outgassing through composite images, sublimation models, and measurements of spectral reflectivity. We find that deep composite images of the object at various epochs along its orbit show no evidence of dust yet place sensitive limits to the dust production. We estimate the nucleus radius to be 11.5 km $\pm 3.5$ km using thermal IR modeling from NEOWISE data and use this and data pre-perihelion to estimate a geometric albedo of 0.25. We compare a CO-sublimation activity model to its post perihelion heliocentric light curve and find this data supports an active fractional area of $5 \times 10^{-6}$ assuming 2 $\mu$m sized grains and other typical comet parameters. The spectral reflectivity of the surface materials obtained with the Gemini Observatory and CFHT at different epochs shows a reddening spectral slope. We compare the physical properties of 2009 MS9 to both TNO and comet populations, and speculate that 2009 MS9's reddening may be due to the buildup of a dust mantle on the surface and could be an explanation of why TNOs exhibit a color bimodality., Comment: 25 pages, 8 figures, Accepted to PSJ

We present the discovery of the first millimeter afterglow of a short-duration $\gamma$-ray burst (SGRB) and the first confirmed afterglow of an SGRB localized by the GUANO system on Swift. Our Atacama Large Millimeter/Sub-millimeter Array (ALMA) detection of SGRB 211106A establishes an origin in a faint host galaxy detected in Hubble Space Telescope (HST) imaging at $0.7\lesssim z\lesssim1.4$. From the lack of a detectable optical afterglow, coupled with the bright millimeter counterpart, we infer a high extinction, $A_{\rm V}\gtrsim2.6$ mag along the line of sight, making this the one of the most highly dust-extincted SGRBs known to date. The millimeter-band light curve captures the passage of the synchrotron peak from the afterglow forward shock and reveals a jet break at $t_{\rm jet}=29.2^{+4.5}_{-4.0}$~days. For a presumed redshift of $z=1$, we infer an opening angle, $\theta_{\rm jet}=(15.5\pm1.4)$~degrees, and beaming-corrected kinetic energy of $\log(E_{\rm K}/{\rm erg})=51.8\pm0.3$, making this one of the widest and most energetic SGRB jets known to date. Combining all published millimeter-band upper limits in conjunction with the energetics for a large sample of SGRBs, we find that energetic outflows in high density environments are more likely to have detectable millimeter counterparts. Concerted afterglow searches with ALMA should yield detection fractions of 24-40% on timescales of $\gtrsim2$~days at rates $\approx0.8$-1.6 per year, outpacing the historical discovery rate of SGRB centimeter-band afterglows., Comment: 20 pages, 9 figures, 4 tables. Version accepted for publication in ApJ Letters

Gamma-rays, the most energetic photons, carry information from the far reaches of extragalactic space with minimal interaction or loss of information. They bring messages about particle acceleration in environments so extreme they cannot be reproduced on earth for a closer look. Gamma-ray astrophysics is so complementary with collider work that particle physicists and astroparticle physicists are often one in the same. Gamma-ray instruments, especially the Fermi Gamma-ray Space Telescope, have been pivotal in major multi-messenger discoveries over the past decade. There is presently a great deal of interest and scientific expertise available to push forward new technologies, to plan and build space- and ground-based gamma-ray facilities, and to build multi-messenger networks with gamma rays at their core. It is therefore concerning that before the community comes together for planning exercises again, much of that infrastructure could be lost to a lack of long-term planning for support of gamma-ray astrophysics. Gamma-rays with energies from the MeV to the EeV band are therefore central to multiwavelength and multi-messenger studies to everything from astroparticle physics with compact objects, to dark matter studies with diffuse large scale structure. These goals and new discoveries have generated a wave of new gamma-ray facility proposals and programs. This paper highlights new and proposed gamma-ray technologies and facilities that have each been designed to address specific needs in the measurement of extreme astrophysical sources that probe some of the most pressing questions in fundamental physics for the next decade. The proposed instrumentation would also address the priorities laid out in the recent Astro2020 Decadal Survey, a complementary study by the astrophysics community that provides opportunities also relevant to Snowmass., Comment: Contribution to Snowmass 2021

The TESS exoplanet-hunting mission detected the rising and decaying optical afterglow of GRB 191016A, a long Gamma-Ray Burst (GRB) detected by Swift-BAT but without prompt XRT or UVOT follow-up due to proximity to the moon. The afterglow has a late peak at least 1000 seconds after the BAT trigger, with a brightest-detected TESS datapoint at 2589.7 s post-trigger. The burst was not detected by Fermi-LAT, but was detected by Fermi-GBM without triggering, possibly due to the gradual nature of rising light curve. Using ground-based photometry, we estimate a photometric redshift of $z_\mathrm{phot} = 3.29\pm{0.40}$. Combined with the high-energy emission and optical peak time derived from TESS, estimates of the bulk Lorentz factor $\Gamma_\mathrm{BL}$ range from $90-133$. The burst is relatively bright, with a peak optical magnitude in ground-based follow-up of $R=15.1$ mag. Using published distributions of GRB afterglows and considering the TESS sensitivity and sampling, we estimate that TESS is likely to detect $\sim1$ GRB afterglow per year above its magnitude limit., Comment: 10 pages, 7 figures. Accepted for publication in the Astrophysical Journal

We present X-ray and multi-band optical observations of the afterglow and host galaxy of GRB 180418A, discovered by ${\it Swift}$/BAT and ${\it Fermi}$/GBM. We present a reanalysis of the GBM and BAT data deriving durations of the prompt emission of $T_{90}\approx$2.56s and $\approx$1.90s, respectively. Modeling the ${\it Fermi}$/GBM catalog of 1405 bursts (2008-2014) in the Hardness-$T_{90}$ plane, we obtain a probability of $\approx$60% that GRB 180418A is a short-hard burst. From a combination of ${\it Swift}$/XRT and ${\it Chandra}$ observations, the X-ray afterglow is detected to $\approx$38.5 days after the burst, and exhibits a single power-law decline with $F_{\rm X} \propto t^{-0.98}$. Late-time Gemini observations reveal a faint r$\approx$25.69 mag host galaxy at an angular offset of $\approx$0.16''. At the likely redshift range of z$\approx$1-2.25, we find that the X-ray afterglow luminosity of GRB 180418A is intermediate between short and long GRBs at all epochs during which there is contemporaneous data, and that GRB 180418A lies closer to the $E_{\gamma,{\rm peak}}-E_{\gamma,{\rm iso}}$ correlation for short GRBs. Modeling the multi-wavelength afterglow with the standard synchrotron model, we derive the burst explosion properties and find a jet opening angle of $\theta_{\rm j} \gtrsim 9-14^{\circ}$. If GRB 180418A is a short GRB that originated from a neutron star merger, it has one of the brightest and longest-lived afterglows along with an extremely faint host galaxy. If instead the event is a long GRB that originated from a massive star collapse, it has among the lowest luminosity afterglows, and lies in a peculiar space in terms of the Hardness-$T_{90}$ and $E_{\gamma,{\rm peak}}-E_{\gamma,{\rm iso}}$ planes., Comment: 26 pages, 9 figures, 3 tables (Accepted, ApJ)

The detection of TeV gamma-ray bursts (GRBs) brought new opportunities for studying the physics of particle acceleration at relativistic shocks. The \hess telescopes recently observed very-high-energy (VHE) emission from a nearby low-luminosity GRB, GRB 190829A. Follow-up observations with, e.g., Swift-XRT, revealed unusual flare activities at $\sim 10^3\rm~s$, which can be caused by a long-lasting central engine. We show that the VHE emission during the H.E.S.S. observation time is naturally produced in the external inverse-Compton (EIC) scenario, where seed photons supplied by the flares or other late-time dissipation are upscattered to VHE energies by the non-thermal electrons accelerated at the external forward shock. Our calculations show that the EIC flare nearly coincides with the late-prompt flare, but extends $\sim3-4$ times longer than the duration of the late-prompt flare. The preferred kinetic energy and initial Lorentz factor used in our model are $\sim 10^{52}\rm~erg$ and $\sim 20$, respectively. Understanding the mechanisms of the VHE emission from low-luminosity GRBs will help us constrain the properties of the outflow and the central engine activities, as well as the particle acceleration mechanism., Comment: 12 pages, 6 figures, accepted for publication in ApJ. HESS data are added, figures are updated, appendix is extended, main conclusions unchanged

The capability of the Fermi Gamma-ray Burst Monitor (GBM) to localize gamma-ray bursts (GRBs) is evaluated for two different automated algorithms: the GBM Team's RoboBA algorithm and the independently developed BALROG algorithm. Through a systematic study utilizing over 500 GRBs with known locations from instruments like Swift and the Fermi LAT, we directly compare the effectiveness of, and accurately estimate the systematic uncertainty for, both algorithms. We show simple adjustments to the GBM Team's RoboBA, in operation since early 2016, yields significant improvement in the systematic uncertainty, removing the long tail identified in the systematic, and improves the overall accuracy. The systematic uncertainty for the updated RoboBA localizations is $1.8^\circ$ for 52% of GRBs and $4.1^\circ$ for the remaining 48%. Both from public reporting by BALROG and our systematic study, we find the systematic uncertainty of $1-2^\circ$ quoted in GCN circulars for bright GRBs localized by BALROG is an underestimate of the true magnitude of the systematic, which we find to be $2.7^\circ$ for 74% of GRBs and $33^\circ$ for the remaining 26%. We show that, once the systematic uncertainty is considered, the RoboBA 90% localization confidence regions can be more than an order of magnitude smaller in area than those produced by BALROG., Comment: Accepted for publication in ApJ. Added Appendix showing localization comparisons for each GRB in the near-realtime public reporting sample

We describe the digest2 software package, a fast, short-arc orbit classifier for small Solar System bodies. The digest2 algorithm has been serving the community for more than 13 years. The code provides a score, D2, which represents a pseudo-probability that a tracklet belongs to a given Solar System orbit type. digest2 is primarily used as a classifier for Near-Earth Object (NEO) candidates, to identify those to be prioritized for follow-up observation. We describe the historical development of digest2 and demonstrate its use on real and synthetic data. We find that digest2 can accurately and precisely distinguish NEOs from non- NEOs. At the time of detection, 14% of NEO tracklets and 98.5% of non-NEOs tracklets have D2 below the critical value of D2 = 65. 94% of our simulated NEOs achieved the maximum D2 = 100 and 99.6% of NEOs achieved $D2 \ge 65$ at least once during the simulated 10-year timeframe. We demonstrate that D2 varies as a function of time, rate of motion, magnitude and sky-plane location, and show that NEOs tend to have lower D2 at low Solar elongations close to the ecliptic. We use our findings to recommend future development directions for the digest2 code., Comment: 19 figures, 7 tables

Recently, Rakshit et al. (2018) reported the discovery of SDSS J103024.95$+$551622.7, a radio-loud narrow-line Seyfert 1 galaxy having a $\sim 100$ kpc scale double-lobed radio structure. Here we analyse archival radio interferometric imaging data taken with the Very Large Array (VLA) at 5 GHz, and with the Very Long Baseline Array (VLBA) at 4.3 and 7.6 GHz. Two hotspots and a compact core are detected with the VLA at arcsec scale, while a single milliarcsec-scale compact radio core is seen with the highest resolution VLBA observations. The Fermi Large Area Telescope did not detect $\gamma$-ray emission at the position of this source. In the mid-infrared, the Wide-field Infrared Survey Explorer satellite light curve, covering more than 7 years and including the most recent data points, hints on flux density variability at 3.4 $\mu$m. Our findings support the notion that this source is a young version of Fanaroff-Riley type II radio galaxies., Comment: 9 pages, 2 figures; accepted for publication in Astrophysics and Space Science

We present observations of the possible short GRB 180418A in $\gamma$-rays, X-rays, and in the optical. Early optical photometry with the TAROT and RATIR instruments show a bright peak ($\approx$ 14.2 AB mag) between $T+28$ and $T+90$ seconds that we interpret as the signature of a reversal shock. Later observations can be modeled by a standard forward shock model and show no evidence of jet break, allowing us to constrain the jet collimation to $\theta_j> 7^\circ$. Using deep late-time optical observations we place an upper limit of $r>24$ AB mag on any underlying host galaxy. The detection of the afterglow in the \textit{Swift} UV filters constrains the GRB redshift to $z<1.3$ and places an upper bound on the $\gamma$-ray isotropic equivalent energy $E_{\rm{\gamma,iso}} < 3 \times 10^{51}$ erg. The properties of this GRB (e.g. duration, hardness ratio, energetic, and environment) lie at the intersection between short and long bursts, and we can not conclusively identify its type. We estimate that the probability that it is drawn from the population of short GRBs is 10\%-30\%., Comment: Accepted por publication in ApJ

In this document, we detail the improvements made to the Fermi GBM targeted sub-threshold search for counterparts to LIGO/Virgo gravitational-wave triggers. We describe the implemented changes and compare the sensitivity of the O3 search to that of the version of the search that operated during O2. Overall, we have improved both the sensitivity and speed of the targeted search. Further improvements to the search have been made for the O3b observing run, including automated upperlimits estimation and incorporating the updated localization systematic with the new version of the search., Comment: Updates and improvements for operations during LIGO/Virgo O3b. New sections on upper limits calculation, localization systematic uncertainty in the search, and a case study on the sub-threshold candidate during O3a: Fermi GBM-190816

The conventional view of stellar-mass binary black hole (sBBH) mergers is that there should not be enough matter present to produce a detectable electromagnetic transient. However, there ARE a number of mechanisms for producing such a counterpart through accretion of matter from various reserves, charged black holes, or interactions with magnetic or exotic fields. After reviewing these mechanisms, we describe what we can learn from multi-messenger observations of sBBH mergers in the areas of stellar evolution and compact binary formation, cosmological measurements, searches for charged black holes, and tests of general relativity and fundamental physics. We end with some remarks about the need to support both observing capabilities and modeling in the next decade and beyond., Comment: A science white paper submitted to Astro2020; this version with additional endorsers and a few small additions to the text

We studied the Near-Earth Asteroid (NEA) candidates posted on the Minor Planet Center's Near-Earth Object Confirmation Page (NEOCP) between years 2013 and 2016. Out of more than 17,000 NEA candidates, while the majority became either new discoveries or were associated with previously known objects, about 11% were unable to be followed-up or confirmed. We further demonstrate that of the unconfirmed candidates, 926+/-50 are likely to be NEAs, representing 18% of discovered NEAs in that period. Only 11% (~93) of the unconfirmed NEA candidates were large (having absolute magnitude H<22). To identify the reasons why these NEAs were not recovered, we analyzed those from the most prolific asteroid surveys: Pan-STARRS, the Catalina Sky Survey, the Dark Energy Survey, and the Space Surveillance Telescope. We examined the influence of plane-of-sky positions and rates of motion, brightnesses, submission delays, and computed absolute magnitudes, as well as correlations with the phase of the moon and seasonal effects. We find that delayed submission of newly discovered NEA candidate to the NEOCP drove a large fraction of the unconfirmed NEA candidates. A high rate of motion was another significant contributing factor. We suggest that prompt submission of suspected NEA discoveries and rapid response to fast moving targets and targets with fast growing ephemeris uncertainty would allow better coordination among dedicated follow-up observers, decrease the number of unconfirmed NEA candidates, and increase the discovery rate of NEAs., Comment: 35 pages, 19 figures, 8 tables

We report the discovery of a $H_r = 3.4\pm0.1$ dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO$_{179}$ is red with $(g-r)=0.88 \pm 0.21$, roughly round, and slowly rotating, with a period of $30.6$ hr. Estimates of its albedo imply a diameter of 600--900~km. Observations sampling the span between 2005--2016 provide an exceptionally well-determined orbit for 2010 JO$_{179}$, with a semi-major axis of $78.307\pm0.009$ au, distant orbits known to this precision are rare. We find that 2010 JO$_{179}$ librates securely within the 21:5 mean-motion resonance with Neptune on hundred-megayear time scales, joining the small but growing set of known distant dwarf planets on metastable resonant orbits. These imply a substantial trans-Neptunian population that shifts between stability in high-order resonances, the detached population, and the eroding population of the scattering disk.

We perform high fidelity simulations of a wide-field telescopic survey searching for Near-Earth Objects (NEO) larger than 140\,m, focusing on the observational and detection model, detection efficiency and accuracy. As a test survey we select the Large Synoptic Survey Telescope. We use its proposed pointings for a 10-year mission and model the detection of near-Earth objects in the fields. We discuss individual model parameters for magnitude losses, vignetting, fading, asteroid rotation and colors, fill factor, limiting magnitude, rate of motion, field shape and rotation and survey patterns and we assess results in terms of the cumulative completeness of the detected population as a function of size and time. Additionally, we examine the sources of modeling uncertainty and derive the overall NEO population completeness for the baseline LSST survey to be $55\pm5$% for NEOs with absolute magnitude brighter than 22. Including already discovered objects and ongoing surveys, the population completeness at the end of the LSST baseline survey should reach $\sim 77\%$., Comment: 44 pages, 20 figures, 11 tables. arXiv admin note: substantial text overlap with arXiv:1705.06209

We have conducted a detailed simulation of LSST's ability to link near-Earth and main belt asteroid detections into orbits. The key elements of the study were a high-fidelity detection model and the presence of false detections in the form of both statistical noise and difference image artifacts. We employed the Moving Object Processing System (MOPS) to generate tracklets, tracks and orbits with a realistic detection density for one month of the LSST survey. The main goals of the study were to understand whether a) the linking of Near-Earth Objects (NEOs) into orbits can succeed in a realistic survey, b) the number of false tracks and orbits will be manageable, and c) the accuracy of linked orbits would be sufficient for automated processing of discoveries and attributions. We found that the overall density of asteroids was more than 5000 per LSST field near opposition on the ecliptic, plus up to 3000 false detections per field in good seeing. We achieved 93.6% NEO linking efficiency for H<22 on tracks composed of tracklets from at least three distinct nights within a 12-day interval. The derived NEO catalog was comprised of 96% correct linkages. Less than 0.1% of orbits included false detections, and the remainder of false linkages stemmed from main belt confusion, which was an artifact of the short time span of the simulation. The MOPS linking efficiency can be improved by refined attribution of detections to known objects and by improved tuning of the internal kd-tree linking algorithms., Comment: 41 pages, 23 figures, 6 tables. arXiv admin note: substantial text overlap with arXiv:1705.06209

We present the \emph{Fermi} Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger (BBH) event GW170104. No candidate electromagnetic counterparts was detected by either GBM or LAT. A detailed analysis of the GBM and LAT data over timescales from seconds to days covering the LIGO localization region is presented. The resulting flux upper bound from the GBM is (5.2--9.4)$\times$10$^{-7}$ erg cm$^{-2}$ s$^{-1}$ in the 10-1000 keV range and from the LAT is (0.2--13)$\times$10$^{-9}$ erg cm$^{-2}$ s$^{-1}$ in the 0.1--1 GeV range. We also describe the improvements to our automated pipelines and analysis techniques for searching for and characterizing the potential electromagnetic counterparts for future gravitational wave events from Advanced LIGO/VIRGO.

This report describes the methodology and results of an assessment study of the performance of the Large Synoptic Survey Telescope (LSST) in its planned efforts to detect and catalog near-Earth objects (NEOs). The baseline LSST survey approach is designed to make only two visits to a given field in a given night, leading to two possible NEO detections per night. These nightly pairs must be linked across nights to derive orbits of moving objects. However, the presence of false detections in the data stream leads to the possibility of high rates of false tracklets, and the ensuing risk that the resulting orbit catalog may be contaminated by false orbits. To understand the challenges stemming from a two-visit-per-night cadence, we conducted high-fidelity linkage tests on a full-density simulated LSST detection stream. We also sought to quantify the overall performance of LSST as an NEO discovery system under the hypothesis that detections from the baseline cadence can be successfully linked. Our simulations revealed that in 10 years LSST would catalog ~60% of NEOs with absolute magnitude H<22, which is a proxy for 140m and larger objects. This results neglects linking losses and the contribution of any other NEO surveys. Including our worst-case linking efficiency we reach a overall performance assessment of 55% completeness of NEOs with H<22. We estimate that survey mis-modeling could account for systematic errors of up to 5%. We find that restricting the evaluation metric to so-called Potentially Hazardous Asteroids (PHAs) increases the completeness by 3-4%, and that including the benefits of past and expected future NEO survey activity increases completeness at the end of the baseline LSST survey by 15-20%. Assembling these results leads to a projection that by the end of the baseline LSST survey the NEO catalog will be 80+/-5% complete for PHAs with H < 22.

We performed a statistical analysis of the astrometric errors for the major asteroid surveys. We analyzed the astrometric residuals as a function of observation epoch, observed brightness and rate of motion, finding that astrometric errors are larger for faint observations and some stations improved their astrometric quality over time. Based on this statistical analysis we develop a new weighting scheme to be used when performing asteroid orbit determination. The proposed weights result in ephemeris predictions that can be conservative by a factor as large as 1.5. However, the new scheme is more robust with respect to outliers and better handles faint detections., Comment: 25 pages, 11 figures, 5 tables

We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely $\sim$2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust tail models show that those two components became active simultaneously $\sim$250 days before the current perihelion, with comparable maximum loss rates of $\sim$0.7 kg s$^{-1}$ and $\sim$0.5 kg s$^{-1}$, and total ejected masses of 8$\times$10$^{6}$ kg and 6$\times$10$^{6}$ kg for fragments J1-A and J1-B, respectively. In consequence, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6 to 9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process., Comment: Accepted by ApJ Letters, Feb. 17, 2017

We derived 90% confidence limits (CL) on the interstellar number density ($\rho_{IS}^{CL}$) of interstellar objects (ISO; comets and asteroids) as a function of the slope of their size-frequency distribution and limiting absolute magnitude. To account for gravitational focusing, we first generated a quasi-realistic ISO population to ~750 au from the Sun and propagated it forward in time to generate a steady state population of ISOs with heliocentric distance <50 au. We then simulated the detection of the synthetic ISOs using pointing data for each image and average detection efficiencies for each of three contemporary solar system surveys --- PS1, the Mt. Lemmon Survey, and the Catalina Sky Survey. These simulations allowed us to determine the surveys' combined ISO detection efficiency under several different but realistic modes of identifying ISOs in the survey data. Some of the synthetic detected ISOs had eccentricities as small as 1.01 --- in the range of the largest eccentricities of several known comets. Our best CL of $\rho_{IS}^{CL} = 1.4 \times 10^{-4}$ au$^{-3}$ implies that the expectation that extra-solar systems form like our solar system, eject planetesimals in the same way, and then distribute them throughout the galaxy, is too simplistic, or that the SFD or behavior of ISOs as they pass through our solar system is far from expectations., Comment: accepted for publication in AJ

The merger of a binary neutron star system is suggested as the central engine of short gamma-ray bursts (sGRBs). For the merger process, simulations predict that magnetic field is amplified beyond magnetar field strength by Kelvin-Helmholtz instabilities. With the Large Area Telescope (LAT), bursts have been detected that show a temporally extended component in coincidence with a short-lasting peak at the end of the prompt phase. We show that the presence of these LAT components in a sGRB could provide evidence of magnetic field amplification in the neutron star merger., Comment: 6 pages and 2 figures, revised version. To appear in the proceedings of the conference "High-Energy Phenomena in Relativistic Outflows (HEPRO V)", held in La Plata, Argentina (5-8 October 2015)

We determine the absolute magnitude (H) distribution (or size-frequency distribution, SFD; $N(H) \propto 10^{\alpha H}$ where $\alpha$ is the slope of the distribution) for near-Earth objects (NEO) with $1322$. There is also another change in slope from steep to shallow around H=27. The three ARM target candidates detected by Pan-STARRS1 in one year of surveying have a corrected SFD with slope $\alpha = 0.40^{+0.33}_{-0.45}$. We also show that the window for follow up observations of small (H$\gtrsim$22) NEOs with the NASA IRTF telescope and Arecibo and Goldstone radars are extremely short - on order of days, and procedures for fast response must be implemented in order to measure physical characteristics of small Earth-approaching objects. CFHT's MegaCam and Pan-STARRS1 have longer observing windows and are capable of following-up more NEOs due to their deeper limiting magnitudes and wider fields of view., Comment: 48 pages, 11 figures, 2 tables

The bright long gamma-ray burst GRB 141207A was observed by the {\it Fermi Gamma-ray Space Telescope} and detected by both instruments onboard. The observations show that the spectrum in the prompt phase is not well described by the canonical empirical Band function alone, and that an additional power-law component is needed. In the early phase of the prompt emission, a modified blackbody with a hard low-energy photon index ($\alpha$ = +0.2 -- +0.4) is detected, which suggests a photospheric origin. In a finely time-resolved analysis, the spectra are also well fitted by the modified blackbody combined with a power-law function. We discuss the physical parameters of the photosphere such as the bulk Lorentz factor of the relativistic flow and the radius. We also discuss the physical origin of the extra power-law component observed during the prompt phase in the context of different models such as leptonic and hadronic scenarios in the internal shock regime and synchrotron emission in the external forward shock. In the afterglow phase, the temporal and spectral behaviors of the temporally extended high-energy emission and the fading X-ray emission detected by XRT on-board {\it Swift} are consistent with synchrotron emission in a radiative external forward shock., Comment: 21 pages, 10 figures, 3 tables, accepted for the publication in ApJ