We report on the observation and measurement of astrometry, photometry, morphology, and activity of the interstellar object 3I/ATLAS, also designated C/2025 N1 (ATLAS) with the NSF-DOE Vera C. Rubin Observatory. Comet 3I/ATLAS, the third known interstellar object, was discovered on UT 2025 July 1. Rubin Observatory had coincidentally collected images of the object's region of the sky during routine commissioning. Facilitated by Rubin's high resolution and large aperture, we successfully recovered object detections from Rubin observations spanning UT 2025 June 21 (10 days before discovery, when 3I/ATLAS was 4.5 au from the Sun) through the date of discovery, and we acquired additional images through UT 2025 July 20 as part of commissioning. We measure on-sky locations of 3I/ATLAS in Rubin ugrizy bands, with a typical precision of similar to 70 mas, and briefly describe the reason this is coarser than our measured static source astrometric precision of similar to 3 mas in Rubin images. We measure grizy magnitudes of 3I/ATLAS photometry at similar to 0.01 mag precision, detecting no short-term photometric variability above 0.01 mag. We derive an estimated near-nucleus dust-to-nucleus scattering cross-sectional ratio of eta greater than or similar to 13 on UT 2025 July 2 based on Rubin photometry and an upper limit nucleus size computed from Hubble Space Telescope observations. We find Rubin colors of g - r = (0.657 +/- 0.013) mag, r - i = (0.235 +/- 0.018) mag, i - z = (0.147 +/- 0.042) mag, and z - y = (0.047 +/- 0.052) mag. These data represent the earliest observations of this object by a large (greater than or similar to 8 m class) telescope and illustrate the type of measurements (and discoveries) Rubin's Legacy Survey of Space and Time will provide after it begins in early 2026.

Gamma-ray bursts are the most luminous electromagnetic events in the Universe. Their prompt gamma-ray emission has typical durations between a fraction of a second and several minutes. A rare subset of these events have durations in excess of a thousand seconds, referred to as ultra-long gamma-ray bursts. Here, we report the discovery of the longest gamma-ray burst ever seen with a similar to 25 000 s gamma-ray duration, GRB 250702B, and characterize this event using data from four instruments in the InterPlanetary Network and the Monitor of All-sky X-ray Image. We find a hard spectrum, subsecond variability, and high total energy, which are only known to arise from ultrarelativistic jets powered by a rapidly spinning stellar-mass central engine. These properties and the extreme duration are together incompatible with all confirmed gamma-ray burst progenitors and nearly all models in the literature. This burst is naturally explained with the helium merger model, where a field binary ends when a black hole falls into a stripped star and proceeds to consume and explode it from within. Under this paradigm, GRB 250702B adds to the growing evidence that helium stars expand and that some ultra-long GRBs have similar evolutionary pathways as collapsars, stellar-mass gravitational wave sources, and potentially rare types of supernovae.

We present a complete analysis of Fermi Large Area Telescope (LAT) data of GRB 221009A, the brightest gamma-ray burst (GRB) ever detected. The burst emission above 30 MeV detected by the LAT preceded, by 1 s, the low-energy (<10 MeV) pulse that triggered the Fermi Gamma-Ray Burst Monitor (GBM), as has been observed in other GRBs. The prompt phase of GRB 221009A lasted a few hundred seconds. It was so bright that we identify a bad time interval of 64 s caused by the extremely high flux of hard X-rays and soft gamma rays, during which the event reconstruction efficiency was poor and the dead time fraction quite high. The late-time emission decayed as a power law, but the extrapolation of the late-time emission during the first 450 s suggests that the afterglow started during the prompt emission. We also found that high-energy events observed by the LAT are incompatible with synchrotron origin, and, during the prompt emission, are more likely related to an extra component identified as synchrotron self-Compton (SSC). A remarkable 400 GeV photon, detected by the LAT 33 ks after the GBM trigger and directionally consistent with the location of GRB 221009A, is hard to explain as a product of SSC or TeV electromagnetic cascades, and the process responsible for its origin is uncertain. Because of its proximity and energetic nature, GRB 221009A is an extremely rare event.

We present a comprehensive multiwavelength analysis of GRB 240825A, a bright gamma-ray burst (GRB) detected by Fermi and Swift, with a prompt duration (T_ 90 4 sec in 50-300 keV) near the boundary separating short and long GRBs, prompting a detailed investigation into its classification and progenitor. Using classical prompt metrics (duration, minimum variability timescale (MVT), lag, and spectral hardness) and modern classification techniques (machine-learning (ML) based t-SNE, support vector machine, energy-hardness-duration, and ε≡ E_γ,iso,52 / E_p,z,2^5/3), we find GRB 240825A exhibits hybrid characteristics. The short MVT (13.830 ± 1.574 ms), rest-frame duration, and ML-based classification indicate a merger-like or ambiguous nature, while its energetics and position on the Amati relation favor a collapsar origin. We conducted deep optical and NIR photometric and spectroscopic late-time search for an associated supernova (SN)/kilonova (KN) and the host galaxy using 10.4 m GTC and 8.4 m binocular LBT telescopes. No bright SN (like SN 1998bw) is detected down to stringent limits (e.g., m_r > 26.1 mag at 17.59 days), despite a redshift of z = 0.659 measured from GTC spectroscopy. Host galaxy SED modeling with Prospector indicates a massive, dusty, and star-forming galaxy-typical of collapsar GRB hosts, though with low sSFR and large offset. We compare these findings with hybrid events like GRB 211211A, GRB 230307A, GRB 200826A, including SNe-GRBs, and conclude that GRB 240825A likely originated from a massive star collapse, with the associated supernova obscured by a dusty host environment or low luminosity SN with absolute magnitude M_V fainter than -18.0. This study emphasizes the need for multiwavelength follow-up and a multi-layered classification to determine GRB progenitors.

Gamma-ray bursts (GRBs) are cosmic events occurring at large distances beyond our galaxy. They provide a unique opportunity to study electromagnetic patterns not seen elsewhere. When the collimated GRB outflow interacts with the outer layers of a star or the wind generated by a binary neutron star merger, it releases energy, forming a quasi-spherical outflow around it. This broad outflow begins to radiate once it has transferred enough energy to the surrounding medium. We have developed a new analytical model that describes the synchrotron afterglow scenario of the quasi-spherical outflow, including factors such as stratified density, self-absorption regime, and the fraction of electrons accelerated by the shock front. We also successfully describe the multiwavelength observations of a sample of llGRB afterglows (GRB 980425, 031203, 060218, 100316D, 130603B, 150101B and 171205A) that exhibited a late component, analysed in both stellar wind and constant-density environments. Our analysis shows that a constant-density environment is favoured. Additionally, we consider the multiwavelength upper limits of the short bursts reported in the Swift-BAT data base.

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 (MPC), JPL’s Center for Near Earth Object Studies (CNEOS) and ESA’s NEO Coordination Centre (NEOCC) 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 analyzed nearly 500000 observations submitted to the MPC from 468 asteroids with published and independently determined high confidence H and G values. 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.

Imminent impactors may be detected only a few hours before their impact with Earth, providing a brief opportunity to characterize them before impact. We describe the characterization of imminent impactor 2024 RW _1 , which was discovered by the Catalina Sky Survey on 2024 September 4 at 05:43 UTC, before it entered the atmosphere near the northern Philippines at 16:39 UTC. We observed 2024 RW _1 with the Astrophysical Research Consortium Telescope Imaging Camera on the Apache Point Astrophysical Research Consortium’s 3.5 m telescope on 2024 September 4 at 10:16 UTC. We obtained g , r , i , and z photometry of 2024 RW _1 , yielding color indices of g – r = 0.47 ± 0.04, r – i = 0.13 ± 0.04, i – z = –0.11 ± 0.07, and g – i = 0.60 ± 0.04, corresponding to a spectral slope of 0.67 ± 0.40 %/100 nm. The closest match to an asteroid spectral type is with B-type asteroids from the C complex. We detect variations in the time-series photometry of the asteroid with an amplitude of ∼0.75, and a double-peaked rotation period of ∼1900 s. Assuming a visible albedo of 0.07 ± 0.03, a density of ∼1500 kg m ^−3 , and a calculated absolute magnitude of 30.92 ± 0.05, we estimate that the asteroid has a diameter of 3.3 ± 0.7 m and a total mass of ∼28,000 kg. Comparing our astrometric orbital solutions to NEOMOD3, the most likely source of 2024 RW _1 is the 3:1 main-belt mean-motion resonance (77% probability) followed by the ν _6 resonance (13% probability), consistent with its organic B-type nature.

Context. We present a comprehensive analysis of the digest2 parameters for candidates of the Near-Earth Object Confirmation Page (NEOCP) that were reported between 2019 and 2024. Our study proposes methods for significantly reducing the inclusion of non-NEO objects on the NEOCP. Despite the substantial increase in near-Earth object (NEO) discoveries in recent years, only about half of the NEOCP candidates are ultimately confirmed as NEOs. Therefore, much observing time is spent following up on non-NEOs. Furthermore, approximately 11% of the candidates remain unconfirmed because the follow-up observations are insufficient. These are nearly 600 cases per year. Aims. To reduce false positives and minimize wasted resources on non-NEOs, we refine the posting criteria for NEOCP based on a detailed analysis of all digest2 scores. Methods. We investigated 30 distinct digest2 parameter categories for candidates that were confirmed as NEOs and non-NEOs. From this analysis, we derived a filtering mechanism based on selected digest2 parameters that were able to exclude 20% of the non-NEOs from the NEOCP while maintaining a minimal loss of true NEOs. We also investigated the application of four machine-learning (ML) techniques, that is, the gradient-boosting machine (GBM), the random forest (RF) classifier, the stochastic gradient descent (SGD) classifier, and neural networks (NN) to classify NEOCP candidates as NEOs or non-NEOs. Based on digest2 parameters as input, our ML models achieved a precision of approximately 95% in distinguishing between NEOs and non-NEOs. Results. Combining the digest2 parameter filter with an ML-based classification model, we demonstrate a significant reduction in non-NEOs on the NEOCP that exceeds 80%, while limiting the loss of NEO discovery tracklets to 5.5%. Importantly, we show that most follow-up tracklets of initially misclassified NEOs are later correctly identified as NEOs. This effectively reduces the net loss of true NEOs to approximately 1%. Conclusions. A greater purity of NEO candidates on the NEOCP would allow follow-up observers to allocate more resources to confirming high-priority objects. This would improve the overall observational efficiency and the confirmation rate of NEO discoveries. We suggest that our methods are used as part of the NEOCP pipeline.

In the past few decades, large universal structures have been found that challenge the homogeneity and isotropy expected in standard cosmological models. The largest of these, identified as the Hercules–Corona Borealis Great Wall, was found in 2014 in the northern galactic hemisphere in the redshift range of 1.6≤z≤2.1. Subsequent studies used an increasing gamma-ray burst database to show that the cluster was unlikely to have been caused by statistical sampling uncertainties. This study re-examines burst clustering in the northern galactic hemisphere using a recently developed methodology. Evidence is provided that the Hercules–Corona Borealis Great Wall cluster is larger than previously thought, with members potentially spanning the redshift range of 0.33≤z≤2.43. The extension of this cluster’s size does not appear to have been due to statistical variations or sampling biases.

Magnetar giant flares (MGFs) are the extremely short, energetic transients originating from highly magnetized neutron stars. When observed in nearby galaxies, these rare events are nearly indistinguishable from cosmological short gamma-ray bursts. We present the analysis of GRB 231115A, a candidate extragalactic MGF observed by Fermi/GBM and localized by INTEGRAL to the starburst galaxy M82. This burst exhibits distinctive temporal and spectral characteristics, including a short duration and a high peak energy, consistent with known MGFs. Time-resolved analysis reveals rapid spectral evolution and a clear correlation between luminosity and spectral hardness, providing robust evidence of relativistic outflows. Archival Chandra data identified point sources within the GRB 231115A localization consistent with the theoretical maximum persistent emission luminosity, though no definitive counterpart was found. Simulations indicate that any transient emission associated with GRB 231115A would require energies exceeding those of typical magnetar bursts to be detectable by current instruments. While the tail of a MGF originating from outside of the Milky Way and its satellite galaxies has never been detected, analysis suggests that such emission could be observable at M82's distance with instruments like Swift/XRT or NICER, though no tail was identified for this event. These findings underscore the need for improved follow-up strategies and technological advancements to enhance MGF detection and characterization.

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.

The study of gamma-ray burst (GRB) jets has focused predominantly on the gamma-ray portion of the spectral energy distribution (SED) to understand jet properties and their physics. Recent theoretical development has turned to the lower-energy side of the SED to test competing jet models. We considered the application of wide-field X-ray detectors to extend the observation of the SED and for better distinguishing spectral models, aimed at resolving theoretical features existing at or below the sensitivity of missions such as Fermi and Swift. A proposed SmallSat reference mission is introduced, and analysis is conducted on simulated GRBs to determine its improvement in understanding the SED compared with the Fermi-gamma-ray burst monitor (GBM). Detection rates of the reference mission are simulated using a GRB population model and convolved with the energy flux needed to resolve models to find estimated rates of GRBs that the reference mission can resolve better than Fermi-GBM. We discuss the methods and results along with the scientific context for this type of mission.

The second Fermi/LAT gamma-ray burst (GRB) Catalogue (2FLGC) spanning the first decade of operations by the Large Area Telescope (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 regimes 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 generalize 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 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 FS region.

The LargE Area burst Polarimeter (LEAP) is a proposed Compton scattering polarimeter that will, for the first time, measure the level of polarization for a significant number of GRBs with sufficient sensitivity to determine the magnetic field structure, composition, energy dissipation mechanism of GRB jets, and determine the prompt emission mechanism of GRBs. Once approved, LEAP will be deployed as an external payload on the International Space Station (ISS) where it will measure GRB polarization over the energy range from 50–1000 keV, perform GRB spectroscopy from 20 keV to 6 MeV, and self-sufficiently determine the source direction. LEAP is uniquely suited to fill a critical gap in our knowledge regarding GRBs, by exposing the underlying physics that governs astrophysical jets and the extreme environment surrounding newborn compact objects.

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.

Context. The Legacy Survey of Space and Time (LSST), 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. Aims. 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. Methods. 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. Results. Our results 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 (right ascension, declination, and magnitude) 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. Conclusions. 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.

Nongravitational accelerations in the absence of observed activity have recently been identified on near-Earth objects (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 VL65 and 2021 UA12 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 1 x 10-9 au day-2 can lead to a change in plane-of-sky positions of similar to 1 x 103 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 (divided by A i divided by >= 1 x 10-8 au day-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.

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 ∼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σ 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.

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 the first stars and galaxies (assessing also their evolution and cosmic re-ionization). GRBs in the past ∼50 years have stimulated 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 currently). 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 the associated phenomenology.

While the dominant radiation mechanism gamma-ray bursts (GRBs) remains a question of debate, synchrotron emission is one of the foremost candidates to describe the multi-wavelength afterglow observations. As such, it is expected that GRBs should present some degree of polarization across their evolution - presenting a feasible means of probing these bursts' energetic and angular properties. Although obtaining polarization data is difficult due to the inherent complexities regarding GRB observations, advances are being made, and theoretical modeling of synchrotron polarization is now more relevant than ever. In this manuscript, we present the polarization for a fiduciary model where the synchrotron forward-shock emission evolving in the radiative-adiabatic regime is described by a radially stratified off-axis outflow. This is parameterized with a power-law velocity distribution and decelerated in a constant-density and wind-like external environment. We apply this theoretical polarization model for selected bursts presenting evidence of off-axis afterglow emission, including the nearest orphan GRB candidates observed by the Neil Gehrels Swift Observatory and a few Gravitational Wave (GWs) events that could generate electromagnetic emission. In the case of GRB 170817A, we require the available polarimetric upper limits in radio wavelengths to constrain its magnetic field geometry.