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.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: NE2001p: A Native Python Implementation of the NE2001 Galactic Electron Density Model
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
Award ID(s):
2020265
PAR ID:
10501152
Author(s) / Creator(s):
;
Publisher / Repository:
Research Notes of the AAS
Date Published:
Journal Name:
Research Notes of the AAS
Volume:
8
Issue:
1
ISSN:
2515-5172
Page Range / eLocation ID:
17
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; (, The Astrophysical Journal)
    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
  2. ; ; ; (, The Astrophysical Journal)
    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
  3. ; ; ; ; ; ; ; ; ; ; et al (, The Astrophysical Journal)
    Abstract A pulsar’s scintillation bandwidth is inversely proportional to the scattering delay, making accurate measurements of scintillation bandwidth critical to characterize unmitigated delays in efforts to measure low-frequency gravitational waves with pulsar timing arrays. In this pilot work, we searched for a subset of known pulsars within ∼97% of the data taken with the Puerto Rico Ultimate Pulsar Processing Instrument for the AO327 survey with the Arecibo telescope, attempting to measure the scintillation bandwidths in the data set by fitting to the 2D autocorrelation function of their dynamic spectra. We successfully measured 38 bandwidths from 23 pulsars (six without prior literature values), finding that: almost all of the measurements are larger than the predictions from NE2001 and YMW16 (two popular galactic models); NE2001 is more consistent with our measurements than YMW16; Gaussian fits to the bandwidth are more consistent with both electron density models than Lorentzian ones; and for the 17 pulsars with prior literature values, the measurements between various sources often vary by a few factors. The success of Gaussian fits may be due to the use of Gaussian fits to train models in previous work. The variance of literature values over time could relate to the scaling factor used to compare measurements, but also seems consistent with time-varying interstellar medium parameters. This work can be extended to the rest of AO327 to further investigate these trends, highlighting the continuing importance of large archival data sets for projects beyond their initial conception. 
    more » « less
  4. ; ; (, International Journal of Quantum Chemistry)
    Abstract We report the implementation of a hierarchical equations of motion (HEOM) module within the open‐source Libra software. It includes the standard and scaled HEOM algorithms for computing the dynamics of open quantum systems interacting with a harmonic bath. The module allows the computing of the evolution of the reduced density matrix, as well as spectral lineshapes. The truncation, filtering, and “update list” schemes, as well as OpenMP parallelization, allow for further computational saving. The package is written in a mix of C++ and Python languages, delivering the best compromise between user friendliness and efficiency. The Python layer of the package takes advantage of standard Python libraries, such as h5py, which allows efficient storage and retrieval of the generated results. The package can be seamlessly used within Jupyter notebooks; its careful design shall provide the maximal convenience and intuitiveness to its users. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al (, Bioinformatics)
    Marschall, Tobias (Ed.)
    Abstract MotivationJBrowse Jupyter is a package that aims to close the gap between Python programming and genomic visualization. Web-based genome browsers are routinely used for publishing and inspecting genome annotations. Historically they have been deployed at the end of bioinformatics pipelines, typically decoupled from the analysis itself. However, emerging technologies such as Jupyter notebooks enable a more rapid iterative cycle of development, analysis and visualization. ResultsWe have developed a package that provides a Python interface to JBrowse 2’s suite of embeddable components, including the primary Linear Genome View. The package enables users to quickly set up, launch and customize JBrowse views from Jupyter notebooks. In addition, users can share their data via Google’s Colab notebooks, providing reproducible interactive views. Availability and implementationJBrowse Jupyter is released under the Apache License and is available for download on PyPI. Source code and demos are available on GitHub at https://github.com/GMOD/jbrowse-jupyter. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email