An official website of the United States government
Abstract We introduce a new, open-source, Python-based package,extrabol, for inferring the bolometric light curve evolution of extragalactic thermal transients.extraboluses non-parametric Gaussian Process regression for light curve estimation that requires minimal user interaction.extrabolis available via GitHub.
more »
« less
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
- PAR ID:
- 10367726
- 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
-
-
Abstract It has become increasingly useful to answer questions in gravitational-wave astronomy usingtransdimensionalmodels, where the number of free parameters can be varied depending on the complexity required to fit the data. Given the growing interest in transdimensional inference, we introduce a new package for the Bayesian inference Library (Bilby), calledtBilby. ThetBilbypackage allows users to set up transdimensional inference calculations using the existingBilbyarchitecture with off-the-shelf nested samplers and/or Markov Chain Monte Carlo algorithms. Transdimensional models are particularly helpful when seeking to test theoretically uncertain predictions described by phenomenological models. For example, bursts of gravitational waves can be modeled using a superposition ofNwavelets, whereNis itself a free parameter. Short pulses are modeled with small values ofN, whereas longer, more complicated signals are represented with a large number of wavelets stitched together. Other transdimensional models have been used to describe instrumental noise and the population properties of gravitational-wave sources. We provide a few demonstrations oftBilby, including fitting the gravitational-wave signal GW150914 with a superposition ofNsine-Gaussian wavelets. We outline our plans to further develop thetBilbycode suite for a broader range of transdimensional problems.more » « less
-
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 103M⊙to 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
-
Obtaining accurate and deterministic predictions of the risks associated with the presence of contaminants in aquifers is an illusive goal given the presence of heterogeneity in hydrological properties and limited site characterization data. For such reasons, a probabilistic framework is needed to quantify the risks in groundwater systems. In this work, we present a computational toolboxVisU-HydRAthat aims to statistically characterize and visualize metrics that are relevant in risk analysis with the ultimate goal of supporting decision making. TheVisU-HydRAcomputational toolbox is an open-source Python package that can be linked to a series of existing codes such as MODFLOW and PAR2, a GPU-accelerated transport simulator. To illustrate the capabilities of the computational toolbox, we simulate flow and transport in a heterogeneous aquifer within a Monte Carlo framework. The computational toolbox allows to compute the probability of a contaminant’s concentration exceeding a safe threshold value as well as the uncertainty associated with the loss of resilience of the aquifer. To ensure consistency and a reproducible workflow, a step-by-step tutorial is provided and available on a GitHub repository.more » « less
-
In this article, we introduce the packagebinsreg, which implements the binscatter methods developed by Cattaneo et al. (2024a, arXiv:2407.15276 [stat.EM]; 2024b,American Economic Review114: 1488–1514). The package comprises seven commands:binsreg, binslogit, binsprobit, binsqreg, binstest binspwc, andbinsregselect. The first four commands implement binscatter plotting, point estimation, and uncertainty quantification (confidence intervals and confidence bands) for least-squares linear binscatter regression (binsreg) and for nonlinear binscatter regression (binslogitfor logit regression,binsprobitfor. probit regression, andbinsqregfor quantile regression). The next two commands focus on pointwise and uniform inference:binstestimplements hypothesis testing procedures for parametric specifications and for nonparametric shape restrictions of the unknown regression function, whilebinspwcimplements multigroup pairwise statistical comparisons. The last command,binsregselect, implements. data-driven number-of-bins selectors. The commands offer binned scatterplots and allow for covariate adjustment, weighting, clustering, and multisample analysis, which is useful when studying treatment-effect heterogeneity in randomizec and observational studies, among many other features.more » « less
