Search

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 (BTI) of 64 seconds 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 seconds 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., Comment: 60 pages, 38 figures, 9 tables

After the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Here, we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. To this end, we apply a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment with a unique star formation history and hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolution CMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favourable to the creation of black holes, low-mass galaxies contain few massive clusters and therefore have a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations., Comment: 14 pages, 9 figures, 3 tables

We present an incremental version (4FGL-DR4, for Data Release 4) of the fourth Fermi-LAT catalog containing 7194 gamma-ray sources. Based on the first 14 years of science data in the energy range from 50 MeV to 1 TeV, it uses the same analysis methods as the 4FGL-DR3 catalog did for 12 years of data, with only a few improvements. The spectral parameters, spectral energy distributions, light curves and associations are updated for all sources. We add four new extended sources and modify two existing ones. Among the 6658 4FGL-DR3 sources, we delete 14 and change the localization of 10, while 32 are newly associated, eleven associations are changed and three associations are discarded. We add 546 point sources, among which 8 are considered identified and 229 have a plausible counterpart at other wavelengths. Most are just above the detection threshold, and 14 are transient sources below the detection threshold that can affect the light curves of nearby sources., Comment: Data files at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/14yr_catalog/. Refereed paper is DOI 10.3847/1538-4365/ac6751

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results., Comment: 142 pages. Accepted by the Astrophysical Journal Supplement

Structured, or tabular, data is the most common format in data science. While deep learning models have proven formidable in learning from unstructured data such as images or speech, they are less accurate than simpler approaches when learning from tabular data. In contrast, modern tree-based Machine Learning (ML) models shine in extracting relevant information from structured data. An essential requirement in data science is to reduce model inference latency in cases where, for example, models are used in a closed loop with simulation to accelerate scientific discovery. However, the hardware acceleration community has mostly focused on deep neural networks and largely ignored other forms of machine learning. Previous work has described the use of an analog content addressable memory (CAM) component for efficiently mapping random forests. In this work, we focus on an overall analog-digital architecture implementing a novel increased precision analog CAM and a programmable network on chip allowing the inference of state-of-the-art tree-based ML models, such as XGBoost and CatBoost. Results evaluated in a single chip at 16nm technology show 119x lower latency at 9740x higher throughput compared with a state-of-the-art GPU, with a 19W peak power consumption.

We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing analysis techniques we probe the spectral and temporal evolution of GRB 221009A. We find no emission prior to the GBM trigger time (t0; 2022 October 9 at 13:16:59.99 UTC), indicating that this is the time of prompt emission onset. The triggering pulse exhibits distinct spectral and temporal properties suggestive of the thermal, photospheric emission of shock-breakout, with significant emission up to $\sim$15 MeV. We characterize the onset of external shock at t0+600 s and find evidence of a plateau region in the early-afterglow phase which transitions to a slope consistent with Swift-XRT afterglow measurements. We place the total energetics of GRB 221009A in context with the rest of the GBM sample and find that this GRB has the highest total isotropic-equivalent energy ($\textrm{E}_{\gamma,\textrm{iso}}=1.0\times10^{55}$ erg) and second highest isotropic-equivalent luminosity ($\textrm{L}_{\gamma,\textrm{iso}}=9.9\times10^{53}$ erg/s) based on redshift of z = 0.151. These extreme energetics are what allowed us to observe the continuously emitting central engine of GBM from the beginning of the prompt emission phase through the onset of early afterglow., Comment: 26 pages 7 figures - accepted for publication in ApJL

With nearly a hundred gravitational wave detections, the origin of black hole mergers has become a key question. Here, we focus on understanding the typical galactic environment in which binary black hole mergers arise. To this end, we synthesize progenitors of binary black hole mergers as a function of the redshift of progenitor formation, present-day formation galaxy mass, and progenitor stellar metallicity for $240$ star formation and binary evolution models. We provide guidelines to infer the formation galaxy properties and time of formation, highlighting the interplay between the star formation rate and the efficiency of forming merging binary black holes from binary stars, both of which strongly depend on metallicity. We find that across models, over 50% of BBH mergers have a progenitor metallicity of a few tenths of Solar metallicity, however, inferring formation galaxy properties strongly depends on both the binary evolution model and global metallicity evolution. The numerous, low-mass black holes ($\mathrm{\lesssim 15\,M_{\odot}}$) trace the bulk of the star formation in galaxies heavier than the Milky Way ($M_\mathrm{Gal}$ $\mathrm{\gtrsim 10^{10.5}\,M_{\odot}}$). In contrast, heavier BBH mergers typically stem from larger black holes forming in lower metallicity dwarf galaxies ($M_\mathrm{Gal}$ $\mathrm{\lesssim 10^{9}\,M_{\odot}}$). We find that the progenitors of detectable binary black holes tend to arise from dwarf galaxies at a lower formation redshift ($\lesssim \, 1$). We also produce a posterior probability of the progenitor environment for any detected gravitational wave signal. For the massive GW150914 merger, we show that it likely came from a very low metallicity ($Z$ $\mathrm{\lesssim}\,0.025\,\mathrm{Z_{\odot}}$) environment., Comment: 17 pages, 15 figures, 1 table, Accpeted for publication by MNRAS

Reliable neutron star mass measurements are key to determining the equation-of-state of cold nuclear matter, but these are rare. "Black Widows" and "Redbacks" are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. While inclinations can be inferred from subtle features in optical light curves, such estimates may be systematically biased due to incomplete heating models and poorly-understood variability. Using data from the Fermi Large Area Telescope, we have searched for gamma-ray eclipses from 49 spider systems, discovering significant eclipses in 7 systems, including the prototypical black widow PSR B1957$+$20. Gamma-ray eclipses require direct occultation of the pulsar by the companion, and so the detection, or significant exclusion, of a gamma-ray eclipse strictly limits the binary inclination angle, providing new robust, model-independent pulsar mass constraints. For PSR B1957$+$20, the eclipse implies a much lighter pulsar ($M_{\rm psr} = 1.81 \pm 0.07\,M_{\odot}$) than inferred from optical light curve modelling., Comment: 31 pages, 4 figures, includes supplementary tables; published in Nature Astronomy

The next Galactic core-collapse supernova (CCSN) will be a unique opportunity to study within a fully multi-messenger approach the explosion mechanism responsible for the formation of neutron stars and stellar-mass black holes. State-of-the-art numerical simulations of those events reveal the complexity of the gravitational-wave emission which is highly stochastic. This challenges the possibility to infer the properties of the compact remnant and of its progenitor using the information encoded in the waveforms. In this paper we take further steps in a program we recently initiated to overcome those difficulties. In particular we show how oscillation modes of the proto-neutron star, highly visible in the gravitational-wave signal, can be used to reconstruct the time evolution of their physical properties. Extending our previous work where only the information from a single detector was used we here describe a new data-analysis pipeline that coherently combines gravitational-wave detectors' data and infers the time evolution of a combination of the mass and radius of the compact remnant. The performance of the method is estimated employing waveforms from 2D and 3D CCSN simulations covering a progenitor mass range between 11$\mathrm{M_{\odot}}$\, and 40$\mathrm{M_{\odot}}$\, and different equations of state for both a network of up to five second-generation detectors and the proposed third-generation detectors Einstein Telescope and Cosmic Explorer. Our study shows that it will be possible to infer PNS properties for CCSN events occurring in the vicinity of the Milky Way, up to the Large Magellanic Cloud, with the current generation of gravitational-wave detectors.

The Fermi Large Area Telescope (LAT) light curve repository (LCR) is a publicly available, continually updated library of gamma-ray light curves of variable Fermi-LAT sources generated over multiple timescales. The Fermi-LAT LCR aims to provide publication-quality light curves binned on timescales of 3 days, 7 days, and 30 days for 1525 sources deemed variable in the source catalog of the first 10 years of Fermi-LAT observations. The repository consists of light curves generated through full likelihood analyses that model the sources and the surrounding region, providing fluxes and photon indices for each time bin. The LCR is intended as a resource for the time-domain and multi-messenger communities by allowing users to quickly search LAT data to identify correlated variability and flaring emission episodes from gamma-ray sources. We describe the sample selection and analysis employed by the LCR and provide an overview of the associated data access portal., Comment: Accepted for publication in ApJ Supplement Series

An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all sources. The new reported AGNs include 587 blazar candidates and four radio galaxies. We describe the properties of the new sample and outline changes affecting the previously published one. We also introduce two new parameters in this release, namely the peak energy of the SED high-energy component and the corresponding flux. These parameters allow an assessment of the Compton dominance, the ratio of the Inverse-Compton to the synchrotron peak luminosities, without relying on X-ray data., Comment: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal Supplement Series. Fits files are available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/4LACDR3/ or alternatively at https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-h.fits and https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-l.fits

A paradox of requirements specifications as dominantly practiced in the industry is that they often claim to be object-oriented (OO) but largely rely on procedural (non-OO) techniques. Use cases and user stories describe functional flows, not object types. To gain the benefits provided by object technology (such as extendibility, reusability, reliability), requirements should instead take advantage of the same data abstraction concepts - classes, inheritance, information hiding - as OO design and OO programs. Many people find use cases and user stories appealing because of the simplicity and practicality of the concepts. Can we reconcile requirements with object-oriented principles and get the best of both worlds? This article proposes a unified framework. It shows that the concept of class is general enough to describe not only "objects" in a narrow sense but also scenarios such as use cases and user stories and other important artifacts such as test cases and oracles. Having a single framework opens the way to requirements that enjoy the benefits of both approaches: like use cases and user stories, they reflect the practical views of stakeholders; like object-oriented requirements, they lend themselves to evolution and reuse.

With its last observing run, the LIGO, Virgo, and KAGRA collaboration has detected almost one hundred gravitational waves from compact binary coalescences. A common approach to studying the population properties of the observed binaries is to use phenomenological models to describe the spin, mass, and redshift distributions. More recently, with the aim of providing a clearer link to astrophysical processes forming the observed compact binaries coalescences, several authors have proposed to employ synthetic catalogs for population studies. In this paper, we review how to employ and interpret synthetic binary catalogs for gravitational-wave progenitors studies. We describe how to build multi-channel merger rates and describe their associated probabilities focusing on stellar progenitor properties. We introduce a method to quantify the match between the phenomenological reconstruction of merger rates with synthetic catalogs. We detail the implementation of synthetic catalogs for multi-channel hierarchical Bayesian inference, highlighting computational aspects and issues related to hyper-prior choice. We find that when inferring stellar progenitors' properties from gravitational-wave observations, the relative efficiency in compact objects production should be taken into account. Finally, by simulating binary black hole detections with LIGO and Virgo sensitivity expected for the O4 observing run, we present two case studies related to the inference of the common envelope efficiency and progenitor metallicity of the binary black holes. We finally discuss how progenitors' properties can be linked to binary black hole properties., Comment: 13 pages, 12 figures, submitted to MNRAS

Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Supernova Remnants (SNRs), IC 443, W44, W49B and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 years of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS~I +61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites., Comment: Accepted for publication in The Astrophysical Journal

After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models., Comment: 3 figures in the main text. 3 figures and 8 tables are in the supplementary material

The 559 Hz black-widow pulsar PSR J1555-2908, originally discovered in radio, is also a bright gamma-ray pulsar. Timing its pulsations using 12 yr of Fermi-LAT gamma-ray data reveals long-term variations in its spin frequency that are much larger than is observed from other millisecond pulsars. While this variability in the pulsar rotation rate could be intrinsic "timing noise", here we consider an alternative explanation: the variations arise from the presence of a very-low-mass third object in a wide multi-year orbit around the neutron star and its low-mass companion. With current data, this hierarchical-triple-system model describes the pulsar's rotation slightly more accurately than the best-fitting timing-noise model. Future observations will show if this alternative explanation is correct., Comment: 6 pages, 2 figures, submitted to ApJL

We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the radio and gamma-ray pulsation discovery and timing, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ~10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system., Comment: 15 pages, 6 figures, accepted by ApJ

Purpose: To develop and evaluate a free-breathing respiratory motion compensated 4D (3D+respiration) $T_2$-weighted turbo spin echo sequence with application to radiology and MR-guided radiotherapy. Methods: k-space data are continuously acquired using a rewound Cartesian acquisition with spiral profile ordering (rCASPR) to provide matching contrast to the conventional linear phase encode ordering and to sort data into multiple respiratory phases. Low-resolution respiratory-correlated 4D images were reconstructed with compressed sensing and used to estimate non-rigid deformation vector fields, which were subsequently used for a motion compensated image reconstruction. rCASPR sampling was compared to linear and CASPR sampling in terms of point-spread-function (PSF) and image contrast with in silico, phantom and in vivo experiments. Reconstruction parameters for low-resolution 4D-MRI (spatial resolution and temporal regularization) were determined using a grid search. The proposed motion compensated rCASPR was evaluated in eight healthy volunteers and compared to free-breathing scans with linear sampling. Image quality was compared based on visual inspection and quantitatively by means of the gradient entropy. Results: rCASPR provided a superior PSF (similar in ky and narrower in kz) and showed no considerable differences in images contrast compared to linear sampling. The optimal 4D-MRI reconstruction parameters were spatial resolution=$4.5 mm^3$ and $\lambda_t=10^{-4}$. The groupwise average gradient entropy was 22.31 for linear, 22.20 for rCASPR, 22.14 for soft-gated rCASPR and 22.02 for motion compensated rCASPR. Conclusion: The proposed motion compensated rCASPR enables high quality free-breathing T2-TSE with minimal changes in image contrast and scan time. The proposed method therefore enables direct transfer of clinically used 3D TSE sequences to free-breathing., Comment: 19 pages, 11 figures

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources., Comment: accepted in ApJS; follow-up paper to 1902.10045

Serological tests are important for understanding the physiopathology and following the evolution of the Covid-19 pandemic. Assays based on flow cytometry (FACS) of tissue culture cells expressing the spike (S) protein of SARS-CoV-2 have repeatedly proven to perform slightly better than the plate-based assays ELISA and CLIA (chemiluminescent immuno-assay), and markedly better than lateral flow immuno-assays (LFIA). Here, we describe an optimized and very simple FACS assay based on staining a mix of two Jurkat cell lines, expressing either high levels of the S protein (Jurkat-S) or a fluorescent protein (Jurkat-R expressing m-Cherry, or Jurkat-G, expressing GFP, which serve as an internal negative control). We show that the Jurkat-S\&R-flow test has a much broader dynamic range than a commercial ELISA test and performs at least as well in terms of sensitivity and specificity. Also, it is more sensitive and quantitative than the hemagglutination-based test HAT, which we described recently. The Jurkat-flow test requires only a few microliters of blood; thus, it can be used to quantify various Ig isotypes in capillary blood collected from a finger prick. It can be used also to evaluate serological responses in mice, hamsters, cats and dogs. FACS tests offer a very attractive solution for laboratories with access to tissue culture and flow cytometry who want to monitor serological responses in humans or in animals, and how these relate to susceptibility to infection, or re-infection, by the virus, and to protection against Covid-19.