skip to main content

Title: Pyspeckit: A Spectroscopic Analysis and Plotting Package
Abstract

pyspeckitis a toolkit and library for spectroscopic analysis in Python. We describe thepyspeckitpackage and highlight some of its capabilities, such as interactively fitting a model to data, akin to the historically widely-usedsplotfunction inIRAF.pyspeckitemploys the Levenberg–Marquardt optimization method via thempfitandlmfitimplementations, and important assumptions regarding error estimation are described here. Wrappers to usepymcandemceeas optimizers are provided. A parallelized wrapper to fit lines in spectral cubes is included. As part of theastropyaffiliated package ecosystem,pyspeckitis open source and open development, and welcomes input and collaboration from the community.

Authors:
; ; ; ; ; ;
Award ID(s):
2008101
Publication Date:
NSF-PAR ID:
10367726
Journal Name:
The Astronomical Journal
Volume:
163
Issue:
6
Page Range or eLocation-ID:
Article No. 291
ISSN:
0004-6256
Publisher:
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The Transiting Exoplanet Survey Satellite (TESS) has an exceptionally large plate scale of 21″ px−1, causing most TESS light curves to record the blended light of multiple stars. This creates a danger of misattributing variability observed by TESS to the wrong source, which would invalidate any analysis. We developed a method that can localize the origin of variability on the sky to better than one fifth of a pixel. Given measured frequencies of variability (e.g., from periodogram analysis), we show that the best-fit sinusoid amplitudes to raw light curves extracted from each pixel are distributed in the same way as light from the variable source. The primary assumption of this method is that other nearby stars are not variable at the same frequencies. Essentially, we are using the high frequency resolution of TESS to overcome limitations from its low spatial resolution. We have implemented our method in an open-source Python package,TESS_localize(github.com/Higgins00/TESS-Localize), that determines the location of a variable source on the sky and the most likely Gaia source given TESS pixel data and a set of observed frequencies of variability. Our method utilizes models of the TESS pixel response function, and we characterize systematics in themore »residuals of fitting these models to data. We find that even stars more than three pixels outside a photometric aperture can produce significant contaminant signals in the extracted light curves. Given the ubiquity of source blending in TESS light curves, verifying the source of observed variability should be a standard step in TESS analyses.

    « less
  2. Abstract

    To better understand the orbital dynamics of exoplanets around close binary stars, i.e., circumbinary planets (CBPs), we applied techniques from dynamical systems theory to a physically motivated set of solutions in the Circular Restricted Three-Body Problem (CR3BP). We applied Floquet theory to characterize the linear dynamical behavior—static, oscillatory, or exponential—surrounding planar circumbinary periodic trajectories (limit cycles). We computed prograde and retrograde limit cycles and analyzed their geometries, stability bifurcations, and dynamical structures. Orbit and stability calculations are exact computations in the CR3BP and reproducible through the open-source Python packagepyraa. The periodic trajectories (doi.org/10.5281/zenodo.7532982) produce a set of noncrossing, dynamically cool circumbinary orbits conducive to planetesimal growth. For mass ratiosμ∈ [0.01, 0.50], we found recurring features in the prograde families. These features include (1) an innermost near-circular trajectory, inside which solutions have resonant geometries, (2) an innermost stable trajectory (ac≈ 1.61 − 1.85abin) characterized by a tangent bifurcating limit cycle, and (3) a region of dynamical instability (a≈ 2.1abin; Δa≈ 0.1abin), the exclusion zone, bounded by a pair of critically stable trajectories bifurcating limit cycles. The exterior boundary of the exclusion zone is consistent with prior determinations ofacaround a circular binary. We validate our analytic results withN-body simulations and applymore »them to the Pluto–Charon system. The absence of detected CBPs in the inner stable region, between the prograde exclusion zone andac, suggests that the exclusion zone may inhibit the inward migration of CBPs.

    « less
  3. Abstract

    Given its large plate scale of 21″ pixel−1, analyses of data from the Transiting Exoplanet Survey Satellite (TESS) space telescope must be wary of source confusion from blended light curves, which creates the potential to attribute observed photometric variability to the wrong astrophysical source. We explore the impact of light curve contamination on the detection of fast yellow pulsating supergiant (FYPS) stars as a case study to demonstrate the importance of confirming the source of detected signals in the TESS pixel data. While some of the FYPS signals have already been attributed to contamination from nearby eclipsing binaries, others are suggested to be intrinsic to the supergiant stars. In this work, we carry out a detailed analysis of the TESS pixel data to fit the source locations of the dominant signals reported for 17 FYPS stars with the Python packageTESS_localize. We are able to reproduce the detections of these signals for 14 of these sources, obtaining consistent source locations for four. Three of these originate from contaminants, while the signal reported for BZ Tuc is likely a spurious frequency introduced to the light curve of this 127 day Cepheid by the data processing pipeline. Other signals are not significantmore »enough to be localized with our methods, or have long periods that are difficult to analyze given other TESS systematics. Since no localizable signals hold up as intrinsic pulsation frequencies of the supergiant targets, we argue that unambiguous detection of pulsational variability should be obtained before FYPS are considered a new class of pulsator.

    « less
  4. Abstract

    A large fraction of Type Ia supernova (SN Ia) observations over the next decade will be in the near-infrared (NIR), at wavelengths beyond the reach of the current standard light-curve model for SN Ia cosmology, SALT3 (∼2800–8700 Å central filter wavelength). To harness this new SN Ia sample and reduce future light-curve standardization systematic uncertainties, we train SALT3 at NIR wavelengths (SALT3-NIR) up to 2μm with the open-source model-training softwareSALTshaker, which can easily accommodate future observations. Using simulated data, we show that the training process constrains the NIR model to ∼2%–3% across the phase range (−20 to 50 days). We find that Hubble residual (HR) scatter is smaller using the NIR alone or optical+NIR compared to optical alone, by up to ∼30% depending on filter choice (95% confidence). There is significant correlation between NIR light-curve stretch measurements and luminosity, with stretch and color corrections often improving HR scatter by up to ∼20%. For SN Ia observations expected from the Roman Space Telescope, SALT3-NIR increases the amount of usable data in the SALT framework by ∼20% at redshiftz≲ 0.4 and by ∼50% atz≲ 0.15. The SALT3-NIR model is part of the open-sourceSNCosmoandSNANASN Ia cosmology packages.

  5. Abstract

    Two magnetic braking models are implemented inMESAfor use in theMISTstellar model grids. Stars less than about 1.3 solar masses are observed to spin down over time through interaction with their magnetized stellar winds (i.e., magnetic braking). This is the basis for gyrochronology and is fundamental to the evolution of lower-mass stars. The detailed physics behind magnetic braking are uncertain, as are 1D stellar evolution models. Thus, we calibrate our models and compare to data from open clusters. Each braking model tested here is capable of reproducing aspects of the data, with important distinctions; neither fully accounts for the observations. The Matt et al. prescription matches the slowly rotating stars observed in open clusters but tends to overestimate the presence of rapidly rotating stars. The Garraffo et al. prescription often produces too much angular momentum loss to accurately match the observed slow sequence for lower-mass stars but reproduces the bimodal nature of slowly and rapidly rotating stars observed in open clusters fairly well. Our models additionally do not reproduce the observed solar lithium depletion, corroborating previous findings that effects other than rotation may be important. We find additional evidence that some level of mass dependency may be missing in thesemore »braking models to match the rotation periods observed in clusters older than 1 Gyr better.

    « less