skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on January 1, 2026

Title: ScatteringOptics.jl: An Interstellar Scattering Framework in the Julia Programming Language
ScatteringOptics.jl is an astronomy software package developed in the Julia programming language (Bezanson et al., 2017). It implements physical models for the anisotropic scattering of radio waves, which arise from turbulence in the ionized interstellar medium. This toolkit excels in simulating and modeling the temporal, spatial, and spectral effects of interstellar scintillation in the strong scattering regime, taking advantage of Julia’s speed and composability. The package provides essential functionalities for modeling, analyzing, and interpreting the images of the Galactic Center’s supermassive black hole, Sagittarius A*, especially with the Event Horizon Telescope (EHT), as well as the images of extremely high brightness temperature emissions in active galactic nuclei using space very long baseline interferometry.  more » « less
Award ID(s):
2107681
PAR ID:
10618363
Author(s) / Creator(s):
; ;
Editor(s):
Foreman-Mackey, Dan
Publisher / Repository:
Journal of Open Source Software
Date Published:
Journal Name:
Journal of Open Source Software
Volume:
10
Issue:
105
ISSN:
2475-9066
Page Range / eLocation ID:
6354
Subject(s) / Keyword(s):
astronomy VLBI black holes radio astronomy Julia
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present the first event-horizon-scale images and spatiotemporal analysis of Sgr A* taken with the Event Horizon Telescope in 2017 April at a wavelength of 1.3 mm. Imaging of Sgr A* has been conducted through surveys over a wide range of imaging assumptions using the classical CLEAN algorithm, regularized maximum likelihood methods, and a Bayesian posterior sampling method. Different prescriptions have been used to account for scattering effects by the interstellar medium toward the Galactic center. Mitigation of the rapid intraday variability that characterizes Sgr A* has been carried out through the addition of a “variability noise budget” in the observed visibilities, facilitating the reconstruction of static full-track images. Our static reconstructions of Sgr A* can be clustered into four representative morphologies that correspond to ring images with three different azimuthal brightness distributions and a small cluster that contains diverse nonring morphologies. Based on our extensive analysis of the effects of sparse ( u , v )-coverage, source variability, and interstellar scattering, as well as studies of simulated visibility data, we conclude that the Event Horizon Telescope Sgr A* data show compelling evidence for an image that is dominated by a bright ring of emission with a ring diameter of ∼50 μ as, consistent with the expected “shadow” of a 4 × 10 6 M ⊙ black hole in the Galactic center located at a distance of 8 kpc. 
    more » « less
  2. Abstract The Galactic electron density model NE2001 describes the multicomponent ionized structure of the Milky Way interstellar medium. NE2001 forward models the dispersion and scattering of compact radio sources, including pulsars, fast radio bursts, active galactic nuclei, and masers, and the model is routinely used to predict the distances of radio sources lacking independent distance measures. Here we present the open-source package NE2001p, a fully Python implementation of NE2001. The model parameters are identical to NE2001 but the computational architecture is optimized for Python, yielding small (<1%) numerical differences between NE2001p and the Fortran code. NE2001p can be used on the command-line and through Python scripts available on PyPI. Future package releases will include modular extensions aimed at providing short-term improvements to model accuracy, including a modified thick disk scale height and additional clumps and voids. This implementation of NE2001 is a springboard to a next-generation Galactic electron density model now in development. 
    more » « less
  3. Abstract Cosmic-ray transport on galactic scales depends on the detailed properties of the magnetized, multiphase interstellar medium (ISM). In this work, we postprocess a high-resolution TIGRESS magnetohydrodynamic simulation modeling a local galactic disk patch with a two-moment fluid algorithm for cosmic-ray transport. We consider a variety of prescriptions for the cosmic rays, from a simple, purely diffusive formalism with constant scattering coefficient, to a physically motivated model in which the scattering coefficient is set by the critical balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We separately focus on cosmic rays with kinetic energies of ∼1 GeV (high-energy) and ∼30 MeV (low energy), respectively important for ISM dynamics and chemistry. We find that simultaneously accounting for advection, streaming, and diffusion of cosmic rays is crucial for properly modeling their transport. Advection dominates in the high-velocity, low-density hot phase, while diffusion and streaming are more important in higher-density, cooler phases. Our physically motivated model shows that there is no single diffusivity for cosmic-ray transport: the scattering coefficient varies by four or more orders of magnitude, maximal at density n H ∼ 0.01 cm −3 . The ion-neutral damping of Alfvén waves results in strong diffusion and nearly uniform cosmic-ray pressure within most of the mass of the ISM. However, cosmic rays are trapped near the disk midplane by the higher scattering rate in the surrounding lower-density, higher-ionization gas. The transport of high-energy cosmic rays differs from that of low-energy cosmic rays, with less effective diffusion and greater energy losses for the latter. 
    more » « less
  4. Abstract Radio wave scattering can cause severe reductions in detection sensitivity for surveys of Galactic and extragalactic fast (∼ms duration) transients. While Galactic sources like pulsars undergo scattering in the Milky Way interstellar medium (ISM), extragalactic fast radio bursts (FRBs) can also experience scattering in their host galaxies and other galaxies intervening in their lines of sight. We assess Galactic and extragalactic scattering horizons for fast radio transients using a combination of NE2001 to model the dispersion measure and scattering time (τ) contributed by the Galactic disk, and independently constructed electron density models for the Galactic halo and other galaxies’ ISMs and halos that account for different galaxy morphologies, masses, densities, and strengths of turbulence. For source redshifts 0.5 ≤zs≤ 1, an all-sky, isotropic FRB population has simulated values ofτ(1 GHz) ranging from ∼1μs to ∼2 ms (90% confidence, observer frame) that are dominated by host galaxies, althoughτcan be ≫2 ms at low Galactic latitudes. A population atzs= 5 has 0.01 ≲τ≲ 300 ms at 1 GHz (90% confidence), dominated by intervening galaxies. About 20% of these high-redshift FRBs are predicted to haveτ> 5 ms at 1 GHz (observer frame), and ≳40% of FRBs betweenzs∼ 0.5–5 haveτ≳ 1 ms forν≤ 800 MHz. Our scattering predictions may be conservative if scattering from circumsource environments is significant, which is possible under specific conditions. The percentage of FRBs selected against from scattering could also be substantially larger than we predict if circumgalactic turbulence causes more small-scale (≪1 au) density fluctuations than observed from nearby halos. 
    more » « less
  5. Abstract To date, the search for radio technosignatures has focused on sky location as a primary discriminant between technosignature candidates and anthropogenic radio frequency interference (RFI). In this work, we investigate the possibility of searching for technosignatures by identifying the presence and nature of intensity scintillations arising from the turbulent, ionized plasma of the interstellar medium. Past works have detailed how interstellar scattering can both enhance and diminish the detectability of narrowband radio signals. We use the NE2001 Galactic free electron density model to estimate scintillation timescales to which narrowband signal searches would be sensitive, and discuss ways in which we might practically detect strong intensity scintillations in detected signals. We further analyze the RFI environment of the Robert C. Byrd Green Bank Telescope with the proposed methodology and comment on the feasibility of using scintillation as a filter for technosignature candidates. 
    more » « less