-
Simple single-field inflation models with arbitrarily small tensor/scalar ratio
Authors:
Nina K. Stein,
William H. Kinney
Abstract:
We construct a family of simple single-field inflation models consistent with Planck / BICEP Keck bounds which have a parametrically small tensor amplitude and no running of the scalar spectral index. The construction consists of a constant-roll hilltop inflaton potential with the end of inflation left as a free parameter induced by higher-order operators which become dominant late in inflation. T…
▽ More
We construct a family of simple single-field inflation models consistent with Planck / BICEP Keck bounds which have a parametrically small tensor amplitude and no running of the scalar spectral index. The construction consists of a constant-roll hilltop inflaton potential with the end of inflation left as a free parameter induced by higher-order operators which become dominant late in inflation. This construction directly demonstrates that there is no lower bound on the tensor/scalar ratio for simple single-field inflation models.
△ Less
Submitted 27 December, 2022; v1 submitted 11 October, 2022;
originally announced October 2022.
-
Cyclic Cosmology and Geodesic Completeness
Authors:
William H. Kinney,
Nina K. Stein
Abstract:
We consider recently proposed bouncing cosmological models for which the Hubble parameter is periodic in time, but the scale factor grows from one cycle to the next as a mechanism for shedding entropy. Since the scale factor for a flat universe is equivalent to an overall conformal factor, it has been argued that this growth corresponds to a physically irrelevant rescaling, and such bouncing unive…
▽ More
We consider recently proposed bouncing cosmological models for which the Hubble parameter is periodic in time, but the scale factor grows from one cycle to the next as a mechanism for shedding entropy. Since the scale factor for a flat universe is equivalent to an overall conformal factor, it has been argued that this growth corresponds to a physically irrelevant rescaling, and such bouncing universes can be made perfectly cyclic, extending infinitely into the past and future. We show that any bouncing universe which uses growth of the scale factor to dissipate entropy must necessarily be geodesically past-incomplete, and therefore cannot be truly cyclic in time.
△ Less
Submitted 3 May, 2022; v1 submitted 28 October, 2021;
originally announced October 2021.
-
Natural Inflation After Planck 2018
Authors:
Nina K. Stein,
William H. Kinney
Abstract:
We calculate high-precision constraints on Natural Inflation relative to current observational constraints from Planck 2018 + BICEP/Keck(BK15) Polarization + BAO on $r$ and $n_S$, including post-inflationary history of the universe. We find that, for conventional post-inflationary dynamics, Natural Inflation with a cosine potential is disfavored at greater than 95\% confidence out by current data.…
▽ More
We calculate high-precision constraints on Natural Inflation relative to current observational constraints from Planck 2018 + BICEP/Keck(BK15) Polarization + BAO on $r$ and $n_S$, including post-inflationary history of the universe. We find that, for conventional post-inflationary dynamics, Natural Inflation with a cosine potential is disfavored at greater than 95\% confidence out by current data. If we assume protracted reheating characterized by $\overline{w}>1/3,$ Natural Inflation can be brought into agreement with current observational constraints. However, bringing unmodified Natural Inflation into the 68\% confidence region requires values of $T_{\mathrm{re}}$ below the scale of electroweak symmetry breaking. The addition of a SHOES prior on the Hubble Constant $H_0$ only worsens the fit.
△ Less
Submitted 14 June, 2021; v1 submitted 3 June, 2021;
originally announced June 2021.
-
A probabilistic approach to determination of Ceres' average surface composition from Dawn VIR and GRaND data
Authors:
H. Kurokawa,
B. L. Ehlmann,
M. C. De Sanctis,
M. G. A. Lapôtre,
T. Usui,
N. T. Stein,
T. H. Prettyman,
A. Raponi,
M. Ciarniello
Abstract:
The Visible-Infrared Mapping Spectrometer (VIR) on board the Dawn spacecraft revealed that aqueous secondary minerals -- Mg-phyllosilicates, NH4-bearing phases, and Mg/Ca carbonates -- are ubiquitous on Ceres. Ceres' low reflectance requires dark phases, which were assumed to be amorphous carbon and/or magnetite (~80 wt.%). In contrast, the Gamma Ray and Neutron Detector (GRaND) constrained the ab…
▽ More
The Visible-Infrared Mapping Spectrometer (VIR) on board the Dawn spacecraft revealed that aqueous secondary minerals -- Mg-phyllosilicates, NH4-bearing phases, and Mg/Ca carbonates -- are ubiquitous on Ceres. Ceres' low reflectance requires dark phases, which were assumed to be amorphous carbon and/or magnetite (~80 wt.%). In contrast, the Gamma Ray and Neutron Detector (GRaND) constrained the abundances of C (8-14 wt.%) and Fe (15-17 wt.%). Here, we reconcile the VIR-derived mineral composition with the GRaND-derived elemental composition. First, we model mineral abundances from VIR data, including either meteorite-derived insoluble organic matter, amorphous carbon, magnetite, or combination as the darkening agent and provide statistically rigorous error bars from a Bayesian algorithm combined with a radiative-transfer model. Elemental abundances of C and Fe are much higher than is suggested by the GRaND observations for all models satisfying VIR data. We then show that radiative transfer modeling predicts higher reflectance from a carbonaceous chondrite of known composition than its measured reflectance. Consequently, our second models use multiple carbonaceous chondrite endmembers, allowing for the possibility that their specific textures or minerals other than carbon or magnetite act as darkening agents, including sulfides and tochilinite. Unmixing models with carbonaceous chondrites eliminate the discrepancy in elemental abundances of C and Fe. Ceres' average reflectance spectrum and elemental abundances are best reproduced by carbonaceous-chondrite-like materials (40-70 wt.%), IOM or amorphous carbon (10 wt.%), magnetite (3-8 wt.%), serpentine (10-25 wt.%), carbonates (4-12 wt.%), and NH4-bearing phyllosilicates (1-11 wt.%).
△ Less
Submitted 30 October, 2020;
originally announced November 2020.
-
Incorporating Uncertainties in Atomic Data Into the Analysis of Solar and Stellar Observations: A Case Study in Fe XIII
Authors:
Xixi Yu,
Giulio Del Zanna,
David C. Stenning,
Jessi Cisewski-Kehe,
Vinay L. Kashyap,
Nathan Stein,
David A. van Dyk,
Harry P. Warren,
Mark A. Weber
Abstract:
Information about the physical properties of astrophysical objects cannot be measured directly but is inferred by interpreting spectroscopic observations in the context of atomic physics calculations. Ratios of emission lines, for example, can be used to infer the electron density of the emitting plasma. Similarly, the relative intensities of emission lines formed over a wide range of temperatures…
▽ More
Information about the physical properties of astrophysical objects cannot be measured directly but is inferred by interpreting spectroscopic observations in the context of atomic physics calculations. Ratios of emission lines, for example, can be used to infer the electron density of the emitting plasma. Similarly, the relative intensities of emission lines formed over a wide range of temperatures yield information on the temperature structure. A critical component of this analysis is understanding how uncertainties in the underlying atomic physics propagates to the uncertainties in the inferred plasma parameters. At present, however, atomic physics databases do not include uncertainties on the atomic parameters and there is no established methodology for using them even if they did. In this paper we develop simple models for the uncertainties in the collision strengths and decay rates for Fe XIII and apply them to the interpretation of density sensitive lines observed with the EUV Imagining spectrometer (EIS) on Hinode. We incorporate these uncertainties in a Bayesian framework. We consider both a pragmatic Bayesian method where the atomic physics information is unaffected by the observed data, and a fully Bayesian method where the data can be used to probe the physics. The former generally increases the uncertainty in the inferred density by about a factor of 5 compared with models that incorporate only statistical uncertainties. The latter reduces the uncertainties on the inferred densities, but identifies areas of possible systematic problems with either the atomic physics or the observed intensities.
△ Less
Submitted 17 September, 2018;
originally announced September 2018.
-
The ACS Survey of Galactic Globular Clusters XIV: Bayesian Single-Population Analysis of 69 Globular Clusters
Authors:
R. Wagner-Kaiser,
A. Sarajedini,
T. von Hippel,
D. C. Stenning,
D. A. van Dyk,
E. Jeffery,
E. Robinson,
N. Stein,
J. Anderson,
W. H. Jefferys
Abstract:
We use Hubble Space Telescope (HST) imaging from the ACS Treasury Survey to determine fits for single population isochrones of 69 Galactic globular clusters. Using robust Bayesian analysis techniques, we simultaneously determine ages, distances, absorptions, and helium values for each cluster under the scenario of a "single" stellar population on model grids with solar ratio heavy element abundanc…
▽ More
We use Hubble Space Telescope (HST) imaging from the ACS Treasury Survey to determine fits for single population isochrones of 69 Galactic globular clusters. Using robust Bayesian analysis techniques, we simultaneously determine ages, distances, absorptions, and helium values for each cluster under the scenario of a "single" stellar population on model grids with solar ratio heavy element abundances. The set of cluster parameters is determined in a consistent and reproducible manner for all clusters using the Bayesian analysis suite BASE-9. Our results are used to re-visit the age-metallicity relation. We find correlations with helium and several other parameters such as metallicity, binary fraction, and proxies for cluster mass. The helium abundances of the clusters are also considered in the context of CNO abundances and the multiple population scenario.
△ Less
Submitted 28 February, 2017;
originally announced February 2017.
-
A Bayesian Analysis of the Ages of Four Open Clusters
Authors:
Elizabeth J. Jeffery,
Ted von Hippel,
David A. van Dyk,
David C. Stenning,
Elliot Robinson,
Nathan Stein,
W. H. Jefferys
Abstract:
In this paper we apply a Bayesian technique to determine the best fit of stellar evolution models to find the main sequence turn off age and other cluster parameters of four intermediate-age open clusters: NGC 2360, NGC 2477, NGC 2660, and NGC 3960. Our algorithm utilizes a Markov chain Monte Carlo technique to fit these various parameters, objectively finding the best-fit isochrone for each clust…
▽ More
In this paper we apply a Bayesian technique to determine the best fit of stellar evolution models to find the main sequence turn off age and other cluster parameters of four intermediate-age open clusters: NGC 2360, NGC 2477, NGC 2660, and NGC 3960. Our algorithm utilizes a Markov chain Monte Carlo technique to fit these various parameters, objectively finding the best-fit isochrone for each cluster. The result is a high-precision isochrone fit. We compare these results with the those of traditional "by-eye" isochrone fitting methods. By applying this Bayesian technique to NGC 2360, NGC 2477, NGC 2660, and NGC 3960, we determine the ages of these clusters to be 1.35 +/- 0.05, 1.02 +/- 0.02, 1.64 +/- 0.04, and 0.860 +/- 0.04 Gyr, respectively. The results of this paper continue our effort to determine cluster ages to higher precision than that offered by these traditional methods of isochrone fitting.
△ Less
Submitted 2 November, 2016;
originally announced November 2016.
-
Planck Limits on Non-canonical Generalizations of Large-field Inflation Models
Authors:
Nina K. Stein,
William H. Kinney
Abstract:
In this paper, we consider two case examples of Dirac-Born-Infeld (DBI) generalizations of canonical large-field inflation models, characterized by a reduced sound speed, $c_{S} < 1$. The reduced speed of sound lowers the tensor-scalar ratio, improving the fit of the models to the data, but increases the equilateral-mode non-Gaussianity, $f^\mathrm{equil.}_\mathrm{NL}$, which the latest results fr…
▽ More
In this paper, we consider two case examples of Dirac-Born-Infeld (DBI) generalizations of canonical large-field inflation models, characterized by a reduced sound speed, $c_{S} < 1$. The reduced speed of sound lowers the tensor-scalar ratio, improving the fit of the models to the data, but increases the equilateral-mode non-Gaussianity, $f^\mathrm{equil.}_\mathrm{NL}$, which the latest results from the Planck satellite constrain by a new upper bound. We examine constraints on these models in light of the most recent Planck and BICEP/Keck results, and find that they have a greatly decreased window of viability. The upper bound on $f^\mathrm{equil.}_\mathrm{NL}$ corresponds to a lower bound on the sound speed and a corresponding lower bound on the tensor-scalar ratio of $r \sim 0.01$, so that near-future Cosmic Microwave Background observations may be capable of ruling out entire classes of DBI inflation models. The result is, however, not universal: infrared-type DBI inflation models, where the speed of sound increases with time, are not subject to the bound.
△ Less
Submitted 27 March, 2017; v1 submitted 28 September, 2016;
originally announced September 2016.
-
Detecting Relativistic X-ray Jets in High-Redshift Quasars
Authors:
Kathryn McKeough,
Aneta Siemiginowska,
C. C. Cheung,
Lukasz Stawarz,
Vinay L. Kashyap,
Nathan Stein,
Vasileios Stampoulis,
David A. van Dyk,
J. F. C. Wardle,
N. P. Lee,
D. E. Harris,
D. A. Schwartz,
Davide Donato,
Laura Maraschi,
Fabrizio Tavecchio
Abstract:
We analyze Chandra X-ray images of a sample of 11 quasars that are known to contain kiloparsec scale radio jets. The sample consists of five high-redshift (z >= 3.6) flat-spectrum radio quasars, and six intermediate redshift (2.1 < z < 2.9) quasars. The dataset includes four sources with integrated steep radio spectra and seven with flat radio spectra. A total of 25 radio jet features are present…
▽ More
We analyze Chandra X-ray images of a sample of 11 quasars that are known to contain kiloparsec scale radio jets. The sample consists of five high-redshift (z >= 3.6) flat-spectrum radio quasars, and six intermediate redshift (2.1 < z < 2.9) quasars. The dataset includes four sources with integrated steep radio spectra and seven with flat radio spectra. A total of 25 radio jet features are present in this sample. We apply a Bayesian multi-scale image reconstruction method to detect and measure the X-ray emission from the jets. We compute deviations from a baseline model that does not include the jet, and compare observed X-ray images with those computed with simulated images where no jet features exist. This allows us to compute p-value upper bounds on the significance that an X- ray jet is detected in a pre-determined region of interest. We detected 12 of the features unambiguously, and an additional 6 marginally. We also find residual emission in the cores of 3 quasars and in the background of 1 quasar that suggest the existence of unresolved X-ray jets. The dependence of the X-ray to radio luminosity ratio on redshift is a potential diagnostic of the emission mechanism, since the inverse Compton scattering of cosmic microwave background photons (IC/CMB) is thought to be redshift dependent, whereas in synchrotron models no clear redshift dependence is expected. We find that the high-redshift jets have X-ray to radio flux ratios that are marginally inconsistent with those from lower redshifts, suggesting that either the X-ray emissions is due to the IC/CMB rather than the synchrotron process, or that high redshift jets are qualitatively different.
△ Less
Submitted 12 September, 2016;
originally announced September 2016.
-
Bayesian Analysis of Two Stellar Populations in Galactic Globular Clusters III: Analysis of 30 Clusters
Authors:
R. Wagner-Kaiser,
D. C. Stenning,
A. Sarajedini,
T. von Hippel,
D. A. van Dyk,
E. Robinson,
N. Stein,
W. H. Jefferys
Abstract:
We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of 30 Galactic Globular Clusters to characterize two distinct stellar populations. A sophisticated Bayesian technique is employed to simultaneously sample the joint posterior distribution of age, distance, and extinction for each cluster, as well as unique helium values for two populations within e…
▽ More
We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of 30 Galactic Globular Clusters to characterize two distinct stellar populations. A sophisticated Bayesian technique is employed to simultaneously sample the joint posterior distribution of age, distance, and extinction for each cluster, as well as unique helium values for two populations within each cluster and the relative proportion of those populations. We find the helium differences among the two populations in the clusters fall in the range of ~0.04 to 0.11. Because adequate models varying in CNO are not presently available, we view these spreads as upper limits and present them with statistical rather than observational uncertainties. Evidence supports previous studies suggesting an increase in helium content concurrent with increasing mass of the cluster and also find that the proportion of the first population of stars increases with mass as well. Our results are examined in the context of proposed globular cluster formation scenarios. Additionally, we leverage our Bayesian technique to shed light on inconsistencies between the theoretical models and the observed data.
△ Less
Submitted 6 September, 2016;
originally announced September 2016.
-
The Power of Principled Bayesian Methods in the Study of Stellar Evolution
Authors:
Ted von Hippel,
David A. van Dyk,
David C. Stenning,
Elliot Robinson,
Elizabeth Jeffery,
Nathan Stein,
William H. Jefferys,
Erin O'Malley
Abstract:
It takes years of effort employing the best telescopes and instruments to obtain high-quality stellar photometry, astrometry, and spectroscopy. Stellar evolution models contain the experience of lifetimes of theoretical calculations and testing. Yet most astronomers fit these valuable models to these precious datasets by eye. We show that a principled Bayesian approach to fitting models to stellar…
▽ More
It takes years of effort employing the best telescopes and instruments to obtain high-quality stellar photometry, astrometry, and spectroscopy. Stellar evolution models contain the experience of lifetimes of theoretical calculations and testing. Yet most astronomers fit these valuable models to these precious datasets by eye. We show that a principled Bayesian approach to fitting models to stellar data yields substantially more information over a range of stellar astrophysics. We highlight advances in determining the ages of star clusters, mass ratios of binary stars, limitations in the accuracy of stellar models, post-main-sequence mass loss, and the ages of individual white dwarfs. We also outline a number of unsolved problems that would benefit from principled Bayesian analyses.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
Bayesian Analysis of Two Stellar Populations in Galactic Globular Clusters II: NGC 5024, NGC 5272, and NGC 6352
Authors:
R. Wagner-Kaiser,
D. C. Stenning,
E. Robinson,
T. von Hippel,
A. Sarajedini,
D. A. van Dyk,
N. Stein,
W. H. Jefferys
Abstract:
We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of Galactic Globular Clusters to find and characterize two stellar populations in NGC 5024 (M53), NGC 5272 (M3), and NGC 6352. For these three clusters, both single and double-population analyses are used to determine a best fit isochrone(s). We employ a sophisticated Bayesian analysis technique to…
▽ More
We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of Galactic Globular Clusters to find and characterize two stellar populations in NGC 5024 (M53), NGC 5272 (M3), and NGC 6352. For these three clusters, both single and double-population analyses are used to determine a best fit isochrone(s). We employ a sophisticated Bayesian analysis technique to simultaneously fit the cluster parameters (age, distance, absorption, and metallicity) that characterize each cluster. For the two-population analysis, unique population level helium values are also fit to each distinct population of the cluster and the relative proportions of the populations are determined. We find differences in helium ranging from $\sim$0.05 to 0.11 for these three clusters. Model grids with solar $α$-element abundances ([$α$/Fe] =0.0) and enhanced $α$-elements ([$α$/Fe]=0.4) are adopted.
△ Less
Submitted 20 April, 2016;
originally announced April 2016.
-
Bayesian Analysis of Two Stellar Populations in Galactic Globular Clusters I: Statistical and Computational Methods
Authors:
D. C. Stenning,
R. Wagner-Kaiser,
E. Robinson,
D. A. van Dyk,
T. von Hippel,
A. Sarajedini,
N. Stein
Abstract:
We develop a Bayesian model for globular clusters composed of multiple stellar populations, extending earlier statistical models for open clusters composed of simple (single) stellar populations (vanDyk et al. 2009, Stein et al. 2013). Specifically, we model globular clusters with two populations that differ in helium abundance. Our model assumes a hierarchical structuring of the parameters in whi…
▽ More
We develop a Bayesian model for globular clusters composed of multiple stellar populations, extending earlier statistical models for open clusters composed of simple (single) stellar populations (vanDyk et al. 2009, Stein et al. 2013). Specifically, we model globular clusters with two populations that differ in helium abundance. Our model assumes a hierarchical structuring of the parameters in which physical properties---age, metallicity, helium abundance, distance, absorption, and initial mass---are common to (i) the cluster as a whole or to (ii) individual populations within a cluster, or are unique to (iii) individual stars. An adaptive Markov chain Monte Carlo (MCMC) algorithm is devised for model fitting that greatly improves convergence relative to its precursor non-adaptive MCMC algorithm. Our model and computational tools are incorporated into an open-source software suite known as BASE-9. We use numerical studies to demonstrate that our method can recover parameters of two-population clusters, and also show model misspecification can potentially be identified. As a proof of concept, we analyze the two stellar populations of globular cluster NGC 5272 using our model and methods. (BASE-9 is available from GitHub: https://github.com/argiopetech/base/releases).
△ Less
Submitted 21 April, 2016; v1 submitted 20 April, 2016;
originally announced April 2016.
-
Preprocessing Solar Images while Preserving their Latent Structure
Authors:
Nathan M Stein,
David A van Dyk,
Vinay L Kashyap
Abstract:
Telescopes such as the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory, a NASA satellite, collect massive streams of high resolution images of the Sun through multiple wavelength filters. Reconstructing pixel-by-pixel thermal properties based on these images can be framed as an ill-posed inverse problem with Poisson noise, but this reconstruction is computationally expensive and…
▽ More
Telescopes such as the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory, a NASA satellite, collect massive streams of high resolution images of the Sun through multiple wavelength filters. Reconstructing pixel-by-pixel thermal properties based on these images can be framed as an ill-posed inverse problem with Poisson noise, but this reconstruction is computationally expensive and there is disagreement among researchers about what regularization or prior assumptions are most appropriate. This article presents an image segmentation framework for preprocessing such images in order to reduce the data volume while preserving as much thermal information as possible for later downstream analyses. The resulting segmented images reflect thermal properties but do not depend on solving the ill-posed inverse problem. This allows users to avoid the Poisson inverse problem altogether or to tackle it on each of $\sim$10 segments rather than on each of $\sim$10$^7$ pixels, reducing computing time by a factor of $\sim$10$^6$. We employ a parametric class of dissimilarities that can be expressed as cosine dissimilarity functions or Hellinger distances between nonlinearly transformed vectors of multi-passband observations in each pixel. We develop a decision theoretic framework for choosing the dissimilarity that minimizes the expected loss that arises when estimating identifiable thermal properties based on segmented images rather than on a pixel-by-pixel basis. We also examine the efficacy of different dissimilarities for recovering clusters in the underlying thermal properties. The expected losses are computed under scientifically motivated prior distributions. Two simulation studies guide our choices of dissimilarity function. We illustrate our method by segmenting images of a coronal hole observed on 26 February 2015.
△ Less
Submitted 14 December, 2015;
originally announced December 2015.
-
Detecting Unspecified Structure in Low-Count Images
Authors:
Nathan M. Stein,
David A. van Dyk,
Vinay L. Kashyap,
Aneta Siemiginowska
Abstract:
Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the ima…
▽ More
Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the image. To infer image structure, we conduct a Bayesian analysis of a full model that uses a multiscale component to allow flexible departures from the posited null model. As a test statistic, we use a tail probability of the posterior distribution under the full model. This choice of test statistic allows us to estimate a computationally efficient upper bound on a p-value that enables us to draw strong conclusions even when there are limited computational resources that can be devoted to simulations under the null model. We demonstrate the statistical performance of our method on simulated images. Applying our method to an X-ray image of the quasar 0730+257, we find significant evidence against the null model of a single point source and uniform background, lending support to the claim of an X-ray jet.
△ Less
Submitted 15 October, 2015;
originally announced October 2015.
-
Bayesian Analysis for Stellar Evolution with Nine Parameters (BASE-9): User's Manual
Authors:
Ted von Hippel,
Elliot Robinson,
Elizabeth Jeffery,
Rachel Wagner-Kaiser,
Steven DeGennaro,
Nathan Stein,
David Stenning,
William H. Jefferys,
David van Dyk
Abstract:
BASE-9 is a Bayesian software suite that recovers star cluster and stellar parameters from photometry. BASE-9 is useful for analyzing single-age, single-metallicity star clusters, binaries, or single stars, and for simulating such systems. BASE-9 uses Markov chain Monte Carlo and brute-force numerical integration techniques to estimate the posterior probability distributions for the age, metallici…
▽ More
BASE-9 is a Bayesian software suite that recovers star cluster and stellar parameters from photometry. BASE-9 is useful for analyzing single-age, single-metallicity star clusters, binaries, or single stars, and for simulating such systems. BASE-9 uses Markov chain Monte Carlo and brute-force numerical integration techniques to estimate the posterior probability distributions for the age, metallicity, helium abundance, distance modulus, and line-of-sight absorption for a cluster, and the mass, binary mass ratio, and cluster membership probability for every stellar object. BASE-9 is provided as open source code on a version-controlled web server. The executables are also available as Amazon Elastic Compute Cloud images. This manual provides potential users with an overview of BASE-9, including instructions for installation and use.
△ Less
Submitted 13 November, 2014;
originally announced November 2014.
-
The White Dwarf Age of NGC 2477
Authors:
Elizabeth J. Jeffery,
Ted von Hippel,
Steven DeGennaro,
David van Dyk,
Nathan Stein,
William H. Jefferys
Abstract:
We present deep photometric observations of the open cluster NGC 2477 using HST/WFPC2. By identifying seven cluster white dwarf candidates, we present an analysis of the white dwarf age of this cluster, using both the traditional method of fitting isochrones to the white dwarf cooling sequence, and by employing a new Bayesian statistical technique that has been developed by our group. This new met…
▽ More
We present deep photometric observations of the open cluster NGC 2477 using HST/WFPC2. By identifying seven cluster white dwarf candidates, we present an analysis of the white dwarf age of this cluster, using both the traditional method of fitting isochrones to the white dwarf cooling sequence, and by employing a new Bayesian statistical technique that has been developed by our group. This new method performs an objective, simultaneous model fit of the cluster and stellar parameters (namely age, metallicity, distance, reddening, as well as individual stellar masses, mass ratios, and cluster membership) to the photometry. Based on this analysis, we measure a white dwarf age of 1.035 +/- 0.054 +/- 0.087 Gyr (uncertainties represent the goodness of model fits and discrepancy among models, respectively), in good agreement with the cluster's main sequence turnoff age. This work is part of our ongoing work to calibrate main sequence turnoff and white dwarf ages using open clusters, and to improve the precision of cluster ages to the ~5% level.
△ Less
Submitted 16 February, 2011;
originally announced February 2011.
-
New Techniques to Determine Ages of Open Clusters Using White Dwarfs
Authors:
E. J. Jeffery,
T. von Hippel,
W. H. Jefferys,
D. E. Winget,
N. Stein,
S. DeGennaro
Abstract:
Currently there are two main techniques for independently determining the ages of stellar populations: main sequence evolution theory (via cluster isochrones) and white dwarf cooling theory. Open clusters provide the ideal environment for the calibration of these two clocks. Because current techniques to derive cluster ages from white dwarfs are observationally challenging, we discuss the feasib…
▽ More
Currently there are two main techniques for independently determining the ages of stellar populations: main sequence evolution theory (via cluster isochrones) and white dwarf cooling theory. Open clusters provide the ideal environment for the calibration of these two clocks. Because current techniques to derive cluster ages from white dwarfs are observationally challenging, we discuss the feasibility of determining white dwarf ages from the brighter white dwarfs alone. This would eliminate the requirement of observing the coolest (i.e., faintest) white dwarfs. We discuss our method for testing this new idea, as well as the required photometric precision and prior constraints on metallicity, distance, and reddening. We employ a new Bayesian statistical technique to obtain and interpret results.
△ Less
Submitted 27 November, 2006;
originally announced November 2006.
-
Inverting Color-Magnitude Diagrams to Access Precise Star Cluster Parameters: A Bayesian Approach
Authors:
Ted von Hippel,
William H. Jefferys,
James Scott,
Nathan Stein,
D. E. Winget,
Steven DeGennaro,
Albert Dam,
Elizabeth Jeffery
Abstract:
We demonstrate a new Bayesian technique to invert color-magnitude diagrams of main sequence and white dwarf stars to reveal the underlying cluster properties of age, distance, metallicity, and line-of-sight absorption, as well as individual stellar masses. The advantages our technique has over traditional analyses of color-magnitude diagrams are objectivity, precision, and explicit dependence on…
▽ More
We demonstrate a new Bayesian technique to invert color-magnitude diagrams of main sequence and white dwarf stars to reveal the underlying cluster properties of age, distance, metallicity, and line-of-sight absorption, as well as individual stellar masses. The advantages our technique has over traditional analyses of color-magnitude diagrams are objectivity, precision, and explicit dependence on prior knowledge of cluster parameters. Within the confines of a given set of often-used models of stellar evolution, the initial-final mass relation, and white dwarf cooling, and assuming photometric errors that one could reasonably achieve with the Hubble Space Telescope, our technique yields exceptional precision for even modest numbers of cluster stars. For clusters with 50 to 400 members and one to a few dozen white dwarfs, we find typical internal errors of sigma([Fe/H]) <= 0.03 dex, sigma(m-M) <= 0.02 mag, and sigma(Av) <= 0.01 mag. We derive cluster white dwarf ages with internal errors of typically only 10% for clusters with only three white dwarfs and almost always <= 5% with ten white dwarfs. These exceptional precisions will allow us to test white dwarf cooling models and standard stellar evolution models through observations of white dwarfs in open and globular clusters.
△ Less
Submitted 17 March, 2006;
originally announced March 2006.