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.

 
more » « less
Award ID(s):
2008101
NSF-PAR ID:
10367726
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astronomical Journal
Volume:
163
Issue:
6
ISSN:
0004-6256
Format(s):
Medium: X Size: Article No. 291
Size(s):
["Article No. 291"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    SkyPortalis an open-source software package designed to discover interesting transients efficiently, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, crossmatching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, aSkyPortal-based platform has been operating at scale for >2 yr for the Zwicky Transient Facility Phase II community, with hundreds of users, containing tens of millions of time-domain sources, interacting with dozens of telescopes, and enabling community reporting. WhileSkyPortalemphasizes rich user experiences across common front-end workflows, recognizing that scientific inquiry is increasingly performed programmatically,SkyPortalalso surfaces an extensive and well-documented application programming interface system. From back-end and front-end software to data science analysis tools and visualization frameworks, theSkyPortaldesign emphasizes the reuse and leveraging of best-in-class approaches, with a strong extensibility ethos. For instance,SkyPortalnow leverages ChatGPT large language models to generate and surface source-level human-readable summaries automatically. With the imminent restart of the next generation of gravitational-wave detectors,SkyPortalnow also includes dedicated multimessenger features addressing the requirements of rapid multimessenger follow-up: multitelescope management, team/group organizing interfaces, and crossmatching of multimessenger data streams with time-domain optical surveys, with interfaces sufficiently intuitive for newcomers to the field. This paper focuses on the detailed implementations, capabilities, and early science results that establishSkyPortalas a community software package ready to take on the data science challenges and opportunities presented by this next chapter in the multimessenger era.

     
    more » « less
  2. Abstract

    Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics codeFLASHto theN-body codeph4and the stellar evolution codeSeBausing the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisionalN-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds, and supernovae inFLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ∼1000 stars formed from clouds with a mass range of 103Mto 105M. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H iiregion shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.

     
    more » « less
  3. 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 the 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.

     
    more » « less
  4. Abstract

    We present the first integrated-light, TESS-based light curves for star clusters in the Milky Way, Small Magellanic Cloud, and Large Magellanic Cloud. We explore the information encoded in these light curves, with particular emphasis on variability. We describe our publicly available packageelk, which is designed to extract the light curves by applying principal component analysis to perform background light correction and incorporating corrections for TESS systematics, allowing us to detect variability on timescales shorter than ∼10 days. We perform a series of checks to ensure the quality of our light curves, removing observations where systematics are identified as dominant features, and deliver light curves for 348 previously cataloged open and globular clusters. Where TESS has observed a cluster in more than one observing sector, we provide separate light curves for each sector (for a total of 2204 light curves). We explore in detail the light curves of star clusters known to contain high-amplitude Cepheid and RR Lyrae variable stars, and we confirm that the variability of these known variables is still detectable when summed together with the light from thousands of other stars. We also demonstrate that even some low-amplitude stellar variability is preserved when integrating over a stellar population.

     
    more » « less
  5. 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 apply 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.

     
    more » « less