Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gamma-ray bursts”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

AT2025ulz and S250818k: Leveraging DESI Spectroscopy in the Hunt for a Kilonova Associated with a Subsolar-mass Gravitational-wave Candidate

On 2025 August 18, the LIGO–Virgo–KAGRA collaboration reported a subthreshold gravitational-wave candidate detection consistent with a subsolar-mass neutron star merger, denoted S250818k. An optical transient, AT2025ulz, was discovered within the localization region. AT2025ulz initially appeared to meet the expected behavior of kilonova emission, the telltale signature of a binary neutron star merger. The transient subsequently rebrightened after ∼5 days and was classified as a Type IIb supernova. In this work, we analyze the observations of its host galaxy obtained by the Dark Energy Spectroscopic Instrument (DESI). From the DESI spectrum, we obtain a secure redshift of z = 0.084840 ± 0.000006. If S250818k has an astrophysical origin, this places the transient within 2σ of the gravitational-wave distance and results in an integral overlap between the gravitational-wave alert and the transient location of log 10 $\mathcal{I} ≈ 3.9 - 4.2$. Our analysis of the host galaxy’s spectral energy distribution reveals a star-forming, dusty galaxy with stellar mass ∼10 10 M ⊙ , broadly consistent with the population of both short gamma-ray bursts and core-collapse supernova host galaxies. We also present our follow-up of DESI-selected candidate host galaxies using the Fraunhofer Telescope at the Wendelstein Observatory, and show the promise of DESI for associating or rejecting candidate electromagnetic counterparts to gravitational-wave alerts. These results emphasize the value of DESI’s extensive spectroscopic dataset in rapidly characterizing host galaxies, enabling spectroscopic host subtraction, and guiding targeted follow-up.

Hall, Xander J. [Carnegie Mellon University, Pitts↗
DOE OSTI2026

Exploring the connection between compact object mergers and fast X-ray transients

Context. The connection between compact object mergers and some extragalactic fast X-ray transients (FXRTs) has long been hypothesised but never ultimately established. Aims. In this work, we investigate two FXRTs, the LEIA X-ray Transient LXT 240402A and the Einstein Probe EP 250207b, whose precise positions lie close to nearby (z ≲ 0.1) quiescent galaxies with a negligible probability of chance coincidence, identifying them as particularly promising cases of merger-driven explosions in the local Universe. Methods. We used Chandra to derive accurate localisations for both events and secure otherwise ambiguous associations with their optical counterparts. Deep optical and near-infrared observations with VLT, GTC, and LBT were performed to characterise the surrounding environment and search for kilonova emission, the hallmark of neutron star mergers. Complementary early-time X-ray monitoring with Swift and Einstein Probe was used to constrain the non-thermal afterglow. Results. We find that both FXRTs remain compatible with a compact binary merger progenitor, which produced low-mass ejecta and kilonova emission subdominant to the afterglow. However, alternative explanations such as a distant (z ≳ 1) core-collapse supernova cannot be conclusively ruled out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hunting for gamma-ray emission from fast radio bursts

Context. Fast radio bursts (FRBs) are a recently discovered class of GHz-band, ms-duration, Jansky-level-flux astrophysical transients. Although hundreds of models have been proposed so far for FRB progenitors (the most popular ones involve magnetars), their physical origin and emission mechanism are still a mystery, making them one of the most compelling problems in astrophysics. Aims. FRBs are caused by astrophysical processes that are not yet understood. Exploring their high-energy counterpart is crucial for constraining their origin and emission mechanism. Methods. Thanks to more than 13 years of gamma-ray data collected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, and to more than 1000 FRB events (from 561 non-repeating and 22 repeating sources), one of the largest samples created thus far, we performed the largest and deepest search for high-energy emission from FRB sources to date (between 100 MeV and 1 TeV). In addition to the analysis involving individual FRB events on different timescales (from a few seconds up to several years), we performed, for the first time, a stacking analysis on the full sample of FRB events as well as a search for triplet photons in coincidence with the radio event. Results. We do not detect significant emission, reporting the most stringent constraints, on short timescales, for the FRB-like emission from SGR 1935+2154 with E γ <10 41 erg, corresponding to a factor η < 10 7 with respect to the emitted radio energy. Similarly, for the stacked signal of steady emission from all repeaters, the obtained upper limit (UL) on the FRBs luminosity (L γ < 1.6 × 10 43 erg s -1 is more than two orders of magnitude lower than those derived from the individual sources. Finally, no individual or triplet photons have been significantly associated with FRB events. We derived the LAT ms-sensitivity to be ~0.3 ph cm -2 s -1 and constrained the gamma-ray energy E γ,δT = 1 ms ≲ 10 47 (D L /150 Mpc) 2 erg, ruling out a gamma-ray-to-radio energy ratio greater than 10 9 on ms timescales. Conclusions. The results reported here represent the most stringent UL reported so far on the high-energy emission from FRBs on short and long time scales, as well as on cumulative emission and individual photon searches. While the origin of FRBs is still unclear, our work provides important constraints for FRB modelling, which might shed light on their emission mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Jetted and Turbulent Stellar Deaths: New LVK-detectable Gravitational-wave Sources

Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ~10–100 Hz, over a characteristic jet activity timescale ~10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ~10 events yr –1 . The detection rate drops to ~1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Stellar population and metal production in AGN discs

ABSTRACT As gravitational wave detections increase the number of observed compact binaries (consisting of neutron stars or blacks), we begin to probe the different conditions producing these binaries. Most studies of compact remnant formation focus either on stellar collapse from the evolution of field binary stars in gas-free environments or on the formation of stars in clusters where dynamical interactions capture the compact objects, forming binaries. But a third scenario exists. In this paper, we study the fate of massive stars formed, accrete gas, and evolve in the dense discs surrounding supermassive black holes. We calculate the explosions produced and compact objects formed by the collapse of these massive stars. Nucleosynthetic yields may provide an ideal, directly observable, diagnostic of the formation and fate of these stars in active galactic nuclei. We present a first study of the explosive yields from these stars, comparing these yields with the observed nucleosynthetic signatures in the discs around supermassive stars with quasars. We show that, even though these stars tend to form black holes, their rapid rotation leads to discs that can eject a considerable amount of iron during the collapse of the star. The nucleosynthetic yields from these stars can produce constraints on the number of systems formed in this manner, but further work is needed to exploit variations from the initial models presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

4th TDAMM Workshop White Paper

Time-Domain and Multi-Messenger Astrophysics (TDAMM) is entering a fundamentally new phase characterized by an unprecedented increase in the rate and diversity of astrophysical transient detections. The community is transitioning from a discovery-limited to a follow-up-limited era, driven by major investments across electromagnetic, gravitational-wave, and neutrino observatories. Upcoming facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and wide-field survey instruments will produce a deluge of time-domain alerts, reaching millions of events per night. Simultaneously, upgrades to the gravitational-wave network (LVK O5 and beyond) and neutrino observatories (IceCube Gen2) will significantly increase the detection rates of non-electromagnetic messengers. New high-energy missions and expansions of the InterPlanetary Network (IPN) will further enhance discovery capabilities across the gamma-ray and X-ray regimes. This convergence of capabilities represents a transformative opportunity: for the first time, the community will routinely detect rare and high-impact events across multiple messengers. However, the scientific return from these discoveries will depend critically on the ability to rapidly identify, prioritize, and coordinate follow-up observations across a heterogeneous and globally distributed set of facilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Collapsar Black Holes Are Likely Born Slowly Spinning

Collapsing stars constitute the main black hole (BH) formation channel, and are occasionally associated with the launch of relativistic jets that power γ-ray bursts (GRBs). Thus, collapsars offer an opportunity to infer the natal (before spin-up/down by accretion) BH spin directly from observations. We show that once the BH saturates with a large-scale magnetic flux, the jet power is dictated by the BH spin and mass accretion rate. Core-collapse simulations by Halevi et al. and GRB observations favor stellar density profiles that yield an accretion rate of $\dot m\approx 10^{-2} M_{\odot}s^{-1}$, weakly dependent on time. This leaves the spin as the main factor that governs the jet power. By comparing the jet power to characteristic GRB luminosities, we find that the majority of BHs associated with jets are likely born slowly spinning with a dimensionless spin of a ≃ 0.2, or a ≃ 0.5 for wobbling jets, with the main uncertainty originating in the unknown γ-ray radiative efficiency. This result could be applied to the entire core-collapse BH population, unless an anticorrelation between the stellar magnetic field and angular momentum is present. In a companion paper, Jacquemin-Ide et al., we show that regardless of the natal spin, the extraction of BH rotational energy leads to spin-down to a ≲ 0.2, consistent with gravitational-wave observations. We verify our results by performing the first 3D general-relativistic magnetohydrodynamic simulations of collapsar jets with characteristic GRB energies, powered by slowly spinning BHs. We find that jets of typical GRB power struggle to escape from the star, providing the first numerical indication that many jets fail to generate a GRB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole–Neutron Star Mergers

The ongoing LIGO–Virgo–KAGRA observing run O4 provides an opportunity to discover new multimessenger events, including binary neutron star (BNS) mergers such as GW170817 and the highly anticipated first detection of a multimessenger black hole–neutron star (BH–NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH–NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH–NS mergers. Adopting optimal binary properties, a mass ratio of q = 2, and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from premerger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of ~–15 a few hours after the merger, longer than previous estimates, which did not consider the first principles–based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of ≳1 Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within ~200 Mpc. The BH–NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multimessengers from the Radioactive Decay of r -process Nuclei

The radioactive β-decay of nuclei synthesized in the rapid neutron capture process (r-process) releases a variety of particles, including electrons, γ-rays, neutrinos, and neutrons. These particles provide a rich set of multimessenger signals that carry information about the astrophysical environments where neutron-rich nucleosynthesis occurs. In this work, we calculate from first principles the emission spectra resulting from the β-decay of r-process nuclei. Our approach incorporates detailed nuclear structure and decay data to model the energy distributions of each particle species. We couple the spectra with a nuclear reaction network simulation to obtain the temporal evolution of these distributions. We find that the emission distributions vary significantly in time and are nonthermal, with substantial average energies. We investigate these nuclear signals as a direct probe of heavy element formation and show that they are complementary observables to kilonova.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Ion Acceleration in Fermi-LAT Behind-the-Limb Solar Flares: The Role of Coronal Shock Waves

Here, we investigate the relationship between the gamma-ray emission measured with Fermi-LAT and radio signatures of coronal shock waves in four behind-the-limb (BTL) solar flares. All events were associated with metric type II radio burst. Both start and end times of the radio bursts were synchronized with the gamma-ray emission. The type II bursts associated with the BTL gamma-ray flares had higher speeds and lower formation heights than those of an average sample. These findings support the notion that the highly relativistic ions that produce the gamma-rays in BTL flares are accelerated at CME-driven propagating coronal shock waves rather than in large-scale coronal loops.

Corona↗
DOE OSTI2024

Supernova Axions Convert to Gamma Rays in Magnetic Fields of Progenitor Stars

It has long been established that axions could have been produced within the nascent proto-neutron star formed following the type II supernova SN1987A, escaped the star due to their weak interactions, and then converted to gamma rays in the Galactic magnetic fields; the nonobservation of a gamma-ray flash coincident with the neutrino burst leads to strong constraints on the axion-photon coupling for axion masses m a ≲10 –10 eV. In this Letter, we use SN1987A to constrain higher mass axions, all the way to m a ~10 –3 eV, by accounting for axion production from the Primakoff process, nucleon bremsstrahlung, and pion conversion along with axion-photon conversion on the still-intact magnetic fields of the progenitor star. Moreover, we show that gamma-ray observations of the next Galactic supernova, leveraging the magnetic fields of the progenitor star, could detect quantum chromodynamics axions for masses above roughly 50 μ⁢eV, depending on the supernova. We propose a new full-sky gamma-ray satellite constellation that we call the GALactic AXion Instrument for Supernova (GALAXIS) to search for such future signals along with related signals from extragalactic neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Radiation characterization summary for the WSMR fast burst reactor environment at the 6-inch location

The characterization of the neutron, prompt gamma-ray, and delayed gamma-ray radiation fields for the White Sands Missile Range (WSMR) Fast Burst Reactor, also known as molybdenum-alloy Godiva (Molly-G) has been assessed at the 6-inch irradiation location. The neutron energy spectra, uncertainties, and common radiation metrics are presented. Code-dependent recommended constants are given to facilitate the conversion of various dosimetry readings into radiation metrics desired by experimenters. The Molly-G core was designed and configured similarly to Godiva II, as an unreflected, unmoderated, cylindrical annulus of uranium-molybdenum-alloy fuel with molybdenum loading of 10%. At the 6-inch position, the axial fluence maximum is about 2.4×10 13 n/cm 2 per MJ of reactor energy; about 0.1% of the neutron fluence is below 1 keV and 96% is above 100 keV. The 1-MeV Damage-Equivalent Silicon (DES) fluence is estimated at 2.2×10 13 n/cm 2 per MJ of reactor energy. The prompt gamma-ray dose is roughly 2.5E+03 rad(Si) per MJ and the delayed gamma-ray dose is about 1.3E+03 rad(Si) per MJ.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Possible Counterpart Signal of the Fermi Bubbles at the Cosmic-Ray Positrons

The inner Galaxy has hosted cosmic-ray burst events, including those responsible for the gamma-ray Fermi bubbles and the eROSITA bubbles in X-rays. In this work, we study the Alpha Magnetic Spectrometer positron fraction and find three features around 12, 21, and 48 GeV, of which the lowest energy has a 1.4–4.9σ significance, depending on astrophysical background assumptions. Using background simulations that explain the cosmic-ray positron fraction, positron flux, and electron plus positron flux by primary and secondary cosmic rays and cosmic rays from local pulsars, we test these spectral features as originating from electron/positron burst events from the inner Galaxy. We find the 12 GeV feature to be explained by an event of age τ ≃ 3–10 Myr, in agreement with the proposed age of the Fermi bubbles. Furthermore, the energy in cosmic-ray electrons and positrons propagating along the Galactic disk and not within the Fermi bubbles volume is estimated to be 10 51.5 –10 57.5 erg, or O(10 –4 ) – O(1) the cosmic-ray energy causing the Fermi bubbles. We advocate that these positron fraction features are the counterpart signals of the Fermi bubbles, or of substructures within them, or of the eROSITA bubbles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

BatAnalysis - A Comprehensive Python Pipeline for Swift BAT Survey Analysis

The Swift Burst Alert Telescope (BAT) is a coded-aperture gamma-ray instrument with a large field of view that primarily operates in survey mode when it is not triggering on transient events. The survey data consist of 80- channel detector plane histograms that accumulate photon counts over periods of at least 5 minutes. These histograms are processed on the ground and are used to produce the survey data set between 14 and 195 keV. Survey data comprise >90% of all BAT data by volume and allow for the tracking of long-term light curves and spectral properties of cataloged and uncataloged hard X-ray sources. Until now, the survey data set has not been used to its full potential due to the complexity associated with its analysis and the lack of easily usable pipelines. Here, we introduce the BatAnalysis Python package, a wrapper for HEASoftpy, which provides a modern, opensource pipeline to process and analyze BAT survey data. BatAnalysis allows members of the community to use BAT survey data in more advanced analyses of astrophysical sources, including pulsars, pulsar wind nebula, active galactic nuclei, and other known/unknown transient events that may be detected in the hard X-ray band. We outline the steps taken by the Python code and exemplify its usefulness and accuracy by analyzing survey data of the Crab Nebula, NGC 2992, and a previously uncataloged MAXI transient. The BatAnalysis package allows for ~18 yr of BAT survey data to be used in a systematic way to study a large variety of astrophysical sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

BatAnalysis - A Comprehensive Python Pipeline for Swift BAT Survey Analysis

The Swift Burst Alert Telescope (BAT) is a coded-aperture gamma-ray instrument with a large field of view that primarily operates in survey mode when it is not triggering on transient events. The survey data consist of 80-channel detector plane histograms that accumulate photon counts over periods of at least 5 minutes. These histograms are processed on the ground and are used to produce the survey data set between 14 and 195 keV. Survey data comprise >90% of all BAT data by volume and allow for the tracking of long-term light curves and spectral properties of cataloged and uncataloged hard X-ray sources. Until now, the survey data set has not been used to its full potential due to the complexity associated with its analysis and the lack of easily usable pipelines. Here, we introduce the BatAnalysis Python package, a wrapper for HEASoftpy, which provides a modern, open-source pipeline to process and analyze BAT survey data. BatAnalysis allows members of the community to use BAT survey data in more advanced analyses of astrophysical sources, including pulsars, pulsar wind nebula, active galactic nuclei, and other known/unknown transient events that may be detected in the hard X-ray band. We outline the steps taken by the Python code and exemplify its usefulness and accuracy by analyzing survey data of the Crab Nebula, NGC 2992, and a previously uncataloged MAXI transient. The BatAnalysis package allows for ~18 yr of BAT survey data to be used in a systematic way to study a large variety of astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

High Elevation Radiation Array (HERA) detectors for airborne thunderstorm investigations

A high-energy atmospheric physics phenomenon, referred to as a terrestrial gamma ray flash (TGF), is associated with lightning and produces large bursts of energetic photon radiation. TGFs will be investigated using a suite of gamma-ray instruments designed and constructed to fly on ten United States Air Force (USAF) WC-130J Hurricane Hunter aircraft as part of an aircrew ionization study led by the Air Force Institute of Technology (AFIT) and the United States Air Force School of Aerospace Medicine (USAFSAM), in cooperation with the 53rd Weather Reconnaissance Squadron (WRS). Each instrument consists of one NaI and one plastic detector, a GPS timing device, and an instrument computer that performs data acquisition. High Elevation Radiation Array (HERA) detectors will be employed to maximize the chances of observing TGFs near their source and to gain a better understanding of their origin, mechanism, ubiquity, and to assess potential hazards posed to military and commercial aircrew and passengers. The HERA program, deployed on 10 separate Air Force aircraft over a multi-year campaign, will result in thousands of observational flight hours and be the largest concerted effort to date to observe TGFs in situ through aircraft observations. In this paper, we give an overview of the scientific goals of this campaign and how the HERA instruments have been designed to meet those goals. Here, we include a detailed description of the HERA instrument, along with mass model and signal processing simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Quick Look Report of Godiva U233 shot March 2022 Data Collection Activities

This report discussing the March 2022 irradiation of U233 on the Godiva critical assembly using burst mode and subsequent analysis of the fission products starting at 65 minutes post irradiation, and ending 12 days later. The Short-Lived Fission Product Yield (SLFPY) project is one of many efforts under the Nuclear Physics Multi-Lab and related Venture Project. The objective of the SLFPY project is to provide improved measurements of fission product yields for select nuclides. The SLFPY project consists of measurement and analysis campaigns. To produce fission data, samples of actinide material are irradiated in a high-flux neutron environment. These samples are quickly placed in front of a detection system where both gamma singles and coincidence spectra are collected as a function of time. Gamma spectroscopy is performed on the irradiated material on site near the irradiation facility, starting as soon as feasible after the irradiation and continuing for 7 days or longer. Following the data collection, time-dependent fission yield analyses are performed using gamma-ray spectroscopy. This document summarizes the measurement campaign performed March 21 - April 4, 2022 where sub-gram quantities of 233 U nitrate (UO 2 (NO 3 ) 2 ) were irradiated in the Godiva reactor. The Godiva burst #2076 conducted under experiment IER 504 had a temperature change of +156 degrees Celsius, lasted 13.32 microseconds, and had a reactivity of $\$$1.089. The targets consisted of uranium oxide solution deposited in a quartz tube and packed with quartz wool. The flame-sealed quartz ampoules were used to contain the uranium samples for the irradiation and count activities. On March 22, 2022, at 10:27:00 a.m. the three uranium ampoules, along with a witness foil pack, were irradiated in the Godiva critical assembly using burst mode operations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Periodic Gamma-Ray Modulation of the Blazar PG 1553+113 Confirmed by Fermi-LAT and Multiwavelength Observations

A 2.1 yr periodic oscillation of the gamma-ray flux from the blazar PG 1553+113 has previously been tentatively identified in ∼7 yr of data from the Fermi Large Area Telescope. After 15 yr of Fermi sky-survey observations, doubling the total time range, we report >7 cycle gamma-ray modulation with an estimated significance of 4σ against stochastic red noise. Independent determinations of oscillation period and phase in the earlier and the new data are in close agreement (chance probability <0.01). Pulse timing over the full light curve is also consistent with a coherent periodicity. Multiwavelength new data from Swift X-Ray Telescope, Burst Alert Telescope, and UVOT, and from KAIT, Catalina Sky Survey, All-Sky Automated Survey for Supernovae, and Owens Valley Radio Observatory ground-based observatories as well as archival Rossi X-Ray Timing Explorer satellite-All Sky Monitor data, published optical data of Tuorla, and optical historical Harvard plates data are included in our work. Optical and radio light curves show clear correlations with the gamma-ray modulation, possibly with a nonconstant time lag for the radio flux. We interpret the gamma-ray periodicity as possibly arising from a pulsational accretion flow in a sub-parsec binary supermassive black hole system of elevated mass ratio, with orbital modulation of the supplied material and energy in the jet. Other astrophysical scenarios introduced include instabilities, disk and jet precession, rotation or nutation, and perturbations by massive stars or intermediate-mass black holes in polar orbit.

79 ASTRONOMY AND ASTROPHYSICS↗