Engineering topics

Hogg, David W

Publications and source records attributed to Hogg, David W.

15 records

DOE OSTI2018

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey II: Report of a Community Workshop on the Scientific Synergies Between the SPHEREx Survey and Other Astronomy Observatories

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R > 35 and a 5-$\sigma$ sensitivity AB > 19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

STAR-GALAXY CLASSIFICATION IN MULTI-BAND OPTICAL IMAGING

Ground-based optical surveys such as PanSTARRS, DES, and LSST will produce large catalogs to limiting magnitudes of r {approx}> 24. Star-galaxy separation poses a major challenge to such surveys because galaxies-even very compact galaxies-outnumber halo stars at these depths. We investigate photometric classification techniques on stars and galaxies with intrinsic FWHM <0.2 arcsec. We consider unsupervised spectral energy distribution template fitting and supervised, data-driven support vector machines (SVMs). For template fitting, we use a maximum likelihood (ML) method and a new hierarchical Bayesian (HB) method, which learns the prior distribution of template probabilities from the data. SVM requires training data to classify unknown sources; ML and HB do not. We consider (1) a best-case scenario (SVM{sub best}) where the training data are (unrealistically) a random sampling of the data in both signal-to-noise and demographics and (2) a more realistic scenario where training is done on higher signal-to-noise data (SVM{sub real}) at brighter apparent magnitudes. Testing with COSMOS ugriz data, we find that HB outperforms ML, delivering {approx}80% completeness, with purity of {approx}60%-90% for both stars and galaxies. We find that no algorithm delivers perfect performance and that studies of metal-poor main-sequence turnoff stars may be challenged by poor star-galaxy separation. Using the Receiver Operating Characteristic curve, we find a best-to-worst ranking of SVM{sub best}, HB, ML, and SVM{sub real}. We conclude, therefore, that a well-trained SVM will outperform template-fitting methods. However, a normally trained SVM performs worse. Thus, HB template fitting may prove to be the optimal classification method in future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THE MILKY WAY'S CIRCULAR-VELOCITY CURVE BETWEEN 4 AND 14 kpc FROM APOGEE DATA

We measure the Milky Way's rotation curve over the Galactocentric range 4 kpc {approx}< R {approx}< 14 kpc from the first year of data from the Apache Point Observatory Galactic Evolution Experiment. We model the line-of-sight velocities of 3365 stars in 14 fields with b = 0 Degree-Sign between 30 Degree-Sign {<=} l {<=} 210 Degree-Sign out to distances of 10 kpc using an axisymmetric kinematical model that includes a correction for the asymmetric drift of the warm tracer population ({sigma} {sub R} Almost-Equal-To 35 km s{sup -1}). We determine the local value of the circular velocity to be V{sub c} (R {sub 0}) = 218 {+-} 6 km s{sup -1} and find that the rotation curve is approximately flat with a local derivative between -3.0 km s{sup -1} kpc{sup -1} and 0.4 km s{sup -1} kpc{sup -1}. We also measure the Sun's position and velocity in the Galactocentric rest frame, finding the distance to the Galactic center to be 8 kpc < R {sub 0} < 9 kpc, radial velocity V {sub R, Sun} = -10 {+-} 1 km s{sup -1}, and rotational velocity V {sub {phi}, Sun} = 242{sup +10} {sub -3} km s{sup -1}, in good agreement with local measurements of the Sun's radial velocity and with the observed proper motion of Sgr A*. We investigate various systematic uncertainties and find that these are limited to offsets at the percent level, {approx}2 km s{sup -1} in V{sub c} . Marginalizing over all the systematics that we consider, we find that V{sub c} (R {sub 0}) < 235 km s{sup -1} at >99 % confidence. We find an offset between the Sun's rotational velocity and the local circular velocity of 26 {+-} 3 km s{sup -1}, which is larger than the locally measured solar motion of 12 km s{sup -1}. This larger offset reconciles our value for V{sub c} with recent claims that V{sub c} {approx}> 240 km s{sup -1}. Combining our results with other data, we find that the Milky Way's dark-halo mass within the virial radius is {approx}8 Multiplication-Sign 10{sup 11} M {sub Sun }.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

GALAXY GROWTH BY MERGING IN THE NEARBY UNIVERSE

We measure the mass growth rate by merging for a wide range of galaxy types. We present the small-scale (0.014 h {sup -1} {sub 70} Mpc < r < 11 h {sub 70} {sup -1} Mpc) projected cross-correlation functions w(r {sub p}) of galaxy subsamples from the spectroscopic sample of the NYU Value-Added Galaxy Catalog (5 Multiplication-Sign 10{sup 5} galaxies of redshifts 0.03 < z < 0.15) with galaxy subsamples from the Sloan Digital Sky Survey imaging (4 Multiplication-Sign 10{sup 7} galaxies). We use smooth fits to de-project the two-dimensional functions w(r {sub p}) to obtain smooth three-dimensional real-space cross-correlation functions {xi}(r) for each of several spectroscopic subsamples with each of several imaging subsamples. Because close pairs are expected to merge, the three-space functions and dynamical evolution time estimates provide galaxy accretion rates. We find that the accretion onto massive blue galaxies and onto red galaxies is dominated by red companions, and that onto small-mass blue galaxies, red and blue galaxies make comparable contributions. We integrate over all types of companions and find that at fixed stellar mass, the total fractional accretion rates onto red galaxies ({approx}3 h {sub 70} percent per Gyr) are greater than that onto blue galaxies ({approx}1 h {sub 70} percent per Gyr). These rates are almost certainly overestimates because we have assumed that all close pairs merge as quickly as the merger time that we used. One conclusion of this work is that if the total growth of red galaxies from z = 1 to z = 0 is mainly due to merging, the merger rates must have been higher in the past.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THE VERTICAL MOTIONS OF MONO-ABUNDANCE SUB-POPULATIONS IN THE MILKY WAY DISK

We present the vertical kinematics of stars in the Milky Way's stellar disk inferred from Sloan Digital Sky Survey/Sloan Extension for Galactic Understanding and Exploration (SDSS/SEGUE) G-dwarf data, deriving the vertical velocity dispersion, {sigma}{sub z}, as a function of vertical height |z| and Galactocentric radius R for a set of 'mono-abundance' sub-populations of stars with very similar elemental abundances [{alpha}/Fe] and [Fe/H]. We find that all mono-abundance components exhibit nearly isothermal kinematics in |z|, and a slow outward decrease of the vertical velocity dispersion: {sigma}{sub z}(z, R | [{alpha}/Fe], [Fe/H]) Almost-Equal-To {sigma}{sub z}([{alpha}/Fe], [Fe/H]) Multiplication-Sign exp (- (R - R{sub 0})/7 kpc). The characteristic velocity dispersions of these components vary from {approx}15 km s{sup -1} for chemically young, metal-rich stars with solar [{alpha}/Fe], to {approx}> 50 km s{sup -1} for metal-poor stars that are strongly [{alpha}/Fe]-enhanced, and hence presumably very old. The mean {sigma}{sub z} gradient (d{sigma}{sub z}/dz) away from the mid-plane is only 0.3 {+-} 0.2 km s{sup -1} kpc{sup -1}. This kinematic simplicity of the mono-abundance components mirrors their geometric simplicity; we have recently found their density distribution to be simple exponentials in both the z- and R-directions. We find a continuum of vertical kinetic temperatures ({proportional_to}{sigma}{sup 2}{sub z}) as a function of ([{alpha}/Fe], [Fe/H]), which contribute to the total stellar surface-mass density approximately as {Sigma}{sub R{sub 0}}({sigma}{sup 2}{sub z}){proportional_to} exp(-{sigma}{sup 2}{sub z}). This and the existence of isothermal mono-abundance populations with intermediate dispersions (30-40 km s{sup -1}) reject the notion of a thin-thick-disk dichotomy. This continuum of disk components, ranging from old, 'hot', and centrally concentrated ones to younger, cooler, and radially extended ones, argues against models where the thicker disk portions arise from massive satellite infall or heating; scenarios where either the oldest disk portion was born hot, or where internal evolution plays a major role, seem the most viable. In addition, the wide range of {sigma}{sub z}([{alpha}/Fe], [Fe/H]) combined with a constant {sigma}{sub z}(z) for each abundance bin provides an independent check on the precision of the SEGUE-derived abundances: {delta}{sub [{alpha}/Fe]} Almost-Equal-To 0.07 dex and {delta}{sub [Fe/H]} Almost-Equal-To 0.15 dex. The slow radial decline of the vertical dispersion presumably reflects the decrease in disk surface-mass density. This measurement constitutes a first step toward a purely dynamical estimate of the mass profile of the stellar and gaseous disk in our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THE SPATIAL STRUCTURE OF MONO-ABUNDANCE SUB-POPULATIONS OF THE MILKY WAY DISK

The spatial, kinematic, and elemental-abundance structure of the Milky Way's stellar disk is complex, and has been difficult to dissect with local spectroscopic or global photometric data. Here, we develop and apply a rigorous density modeling approach for Galactic spectroscopic surveys that enables investigation of the global spatial structure of stellar sub-populations in narrow bins of [{alpha}/Fe] and [Fe/H], using 23,767 G-type dwarfs from SDSS/SEGUE, which effectively sample 5 kpc < R{sub GC} < 12 kpc and 0.3 kpc {approx}< |Z| {approx}< 3 kpc. We fit models for the number density of each such ([{alpha}/Fe] and [Fe/H]) mono-abundance component, properly accounting for the complex spectroscopic SEGUE sampling of the underlying stellar population, as well as for the metallicity and color distributions of the samples. We find that each mono-abundance sub-population has a simple spatial structure that can be described by a single exponential in both the vertical and radial directions, with continuously increasing scale heights ( Almost-Equal-To 200 pc to 1 kpc) and decreasing scale lengths (>4.5 kpc to 2 kpc) for increasingly older sub-populations, as indicated by their lower metallicities and [{alpha}/Fe] enhancements. That the abundance-selected sub-components with the largest scale heights have the shortest scale lengths is in sharp contrast with purely geometric 'thick-thin disk' decompositions. To the extent that [{alpha}/Fe] is an adequate proxy for age, our results directly show that older disk sub-populations are more centrally concentrated, which implies inside-out formation of galactic disks. The fact that the largest scale-height sub-components are most centrally concentrated in the Milky Way is an almost inevitable consequence of explaining the vertical structure of the disk through internal evolution. Whether the simple spatial structure of the mono-abundance sub-components and the striking correlations between age, scale length, and scale height can be plausibly explained by satellite accretion or other external heating remains to be seen.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

PHOTOMETRIC REDSHIFTS AND QUASAR PROBABILITIES FROM A SINGLE, DATA-DRIVEN GENERATIVE MODEL

We describe a technique for simultaneously classifying and estimating the redshift of quasars. It can separate quasars from stars in arbitrary redshift ranges, estimate full posterior distribution functions for the redshift, and naturally incorporate flux uncertainties, missing data, and multi-wavelength photometry. We build models of quasars in flux-redshift space by applying the extreme deconvolution technique to estimate the underlying density. By integrating this density over redshift, one can obtain quasar flux densities in different redshift ranges. This approach allows for efficient, consistent, and fast classification and photometric redshift estimation. This is achieved by combining the speed obtained by choosing simple analytical forms as the basis of our density model with the flexibility of non-parametric models through the use of many simple components with many parameters. We show that this technique is competitive with the best photometric quasar classification techniques-which are limited to fixed, broad redshift ranges and high signal-to-noise ratio data-and with the best photometric redshift techniques when applied to broadband optical data. We demonstrate that the inclusion of UV and NIR data significantly improves photometric quasar-star separation and essentially resolves all of the redshift degeneracies for quasars inherent to the ugriz filter system, even when included data have a low signal-to-noise ratio. For quasars spectroscopically confirmed by the SDSS 84% and 97% of the objects with Galaxy Evolution Explorer UV and UKIDSS NIR data have photometric redshifts within 0.1 and 0.3, respectively, of the spectroscopic redshift; this amounts to about a factor of three improvement over ugriz-only photometric redshifts. Our code to calculate quasar probabilities and redshift probability distributions is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THE EXTREME SMALL SCALES: DO SATELLITE GALAXIES TRACE DARK MATTER?

We investigate the radial distribution of galaxies within their host dark matter halos as measured in the Sloan Digital Sky Survey by modeling their small-scale clustering. Specifically, we model the Jiang et al. measurements of the galaxy two-point correlation function down to very small projected separations (10 h{sup -1} kpc {<=} r {<=} 400 h{sup -1} kpc), in a wide range of luminosity threshold samples (absolute r-band magnitudes of -18 up to -23). We use a halo occupation distribution framework with free parameters that specify both the number and spatial distribution of galaxies within their host dark matter halos. We assume one galaxy resides in the halo center and additional galaxies are considered satellites that follow a radial density profile similar to the dark matter Navarro-Frenk-White (NFW) profile, except that the concentration and inner slope are allowed to vary. We find that in low luminosity samples (M{sub r} < -19.5 and lower), satellite galaxies have radial profiles that are consistent with NFW. M{sub r} < -20 and brighter satellite galaxies have radial profiles with significantly steeper inner slopes than NFW (we find inner logarithmic slopes ranging from -1.6 to -2.1, as opposed to -1 for NFW). We define a useful metric of concentration, M{sub 1/10}, which is the fraction of satellite galaxies (or mass) that are enclosed within one-tenth of the virial radius of a halo. We find that M{sub 1/10} for low-luminosity satellite galaxies agrees with NFW, whereas for luminous galaxies it is 2.5-4 times higher, demonstrating that these galaxies are substantially more centrally concentrated within their dark matter halos than the dark matter itself. Our results therefore suggest that the processes that govern the spatial distribution of galaxies, once they have merged into larger halos, must be luminosity dependent, such that luminous galaxies become poor tracers of the underlying dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THE SDSS-III BARYON OSCILLATION SPECTROSCOPIC SURVEY: QUASAR TARGET SELECTION FOR DATA RELEASE NINE

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg{sup 2}, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations (BAOs) in the distribution of Ly{alpha} absorption from the spectra of a sample of {approx}150,000 z > 2.2 quasars. Along with measuring the angular diameter distance at z Almost-Equal-To 2.5, BOSS will provide the first direct measurement of the expansion rate of the universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars in the redshift range 2.2 < z < 3.5, where their colors tend to overlap those of the far more numerous stars. During the first year of the BOSS survey, quasar target selection (QTS) methods were developed and tested to meet the requirement of delivering at least 15 quasars deg{sup -2} in this redshift range, with a goal of 20 out of 40 targets deg{sup -2} allocated to the quasar survey. To achieve these surface densities, the magnitude limit of the quasar targets was set at g {<=} 22.0 or r {<=} 21.85. While detection of the BAO signature in the distribution of Ly{alpha} absorption in quasar spectra does not require a uniform target selection algorithm, many other astrophysical studies do. We have therefore defined a uniformly selected subsample of 20 targets deg{sup -2}, for which the selection efficiency is just over 50% ({approx}10 z > 2.20 quasars deg{sup -2}). This 'CORE' subsample will be fixed for Years Two through Five of the survey. For the remaining 20 targets deg{sup -2}, we will continue to develop improved selection techniques, including the use of additional data sets beyond the Sloan Digital Sky Survey (SDSS) imaging data. In this paper, we describe the evolution and implementation of the BOSS QTS algorithms during the first two years of BOSS operations (through 2011 July), in support of the science investigations based on these data, and we analyze the spectra obtained during the first year. During this year, 11,263 new z > 2.20 quasars were spectroscopically confirmed by BOSS, roughly double the number of previously known quasars with z > 2.20. Our current algorithms select an average of 15 z > 2.20 quasars deg{sup -2} from 40 targets deg{sup -2} using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS QTS.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

AN AFFINE-INVARIANT SAMPLER FOR EXOPLANET FITTING AND DISCOVERY IN RADIAL VELOCITY DATA

Markov chain Monte Carlo (MCMC) proves to be powerful for Bayesian inference and in particular for exoplanet radial velocity fitting because MCMC provides more statistical information and makes better use of data than common approaches like chi-square fitting. However, the nonlinear density functions encountered in these problems can make MCMC time-consuming. In this paper, we apply an ensemble sampler respecting affine invariance to orbital parameter extraction from radial velocity data. This new sampler has only one free parameter, and does not require much tuning for good performance, which is important for automatization. The autocorrelation time of this sampler is approximately the same for all parameters and far smaller than Metropolis-Hastings, which means it requires many fewer function calls to produce the same number of independent samples. The affine-invariant sampler speeds up MCMC by hundreds of times compared with Metropolis-Hastings in the same computing situation. This novel sampler would be ideal for projects involving large data sets such as statistical investigations of planet distribution. The biggest obstacle to ensemble samplers is the existence of multiple local optima; we present a clustering technique to deal with local optima by clustering based on the likelihood of the walkers in the ensemble. We demonstrate the effectiveness of the sampler on real radial velocity data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

THE PRISM MULTI-OBJECT SURVEY (PRIMUS). I. SURVEY OVERVIEW AND CHARACTERISTICS

We present the PRIsm MUlti-object Survey (PRIMUS), a spectroscopic faint galaxy redshift survey to z {approx} 1. PRIMUS uses a low-dispersion prism and slitmasks to observe {approx}2500 objects at once in a 0.18 deg{sup 2} field of view, using the Inamori Magellan Areal Camera and Spectrograph camera on the Magellan I Baade 6.5 m telescope at Las Campanas Observatory. PRIMUS covers a total of 9.1 deg{sup 2} of sky to a depth of i{sub AB} {approx} 23.5 in seven different deep, multi-wavelength fields that have coverage from the Galaxy Evolution Explorer, Spitzer, and either XMM or Chandra, as well as multiple-band optical and near-IR coverage. PRIMUS includes {approx}130,000 robust redshifts of unique objects with a redshift precision of {sigma}{sub z}/(1 + z) {approx} 0.005. The redshift distribution peaks at z {approx} 0.6 and extends to z = 1.2 for galaxies and z = 5 for broad-line active galactic nuclei. The motivation, observational techniques, fields, target selection, slitmask design, and observations are presented here, with a brief summary of the redshift precision; a forthcoming paper presents the data reduction, redshift fitting, redshift confidence, and survey completeness. PRIMUS is the largest faint galaxy survey undertaken to date. The high targeting fraction ({approx}80%) and large survey size will allow for precise measures of galaxy properties and large-scale structure to z {approx} 1.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

A SYSTEMATIC SEARCH FOR MASSIVE BLACK HOLE BINARIES IN THE SLOAN DIGITAL SKY SURVEY SPECTROSCOPIC SAMPLE

We present the results of a systematic search for massive black hole binaries in the Sloan Digital Sky Survey (SDSS) spectroscopic database. We focus on bound binaries, under the assumption that one of the black holes is active. In this framework, the broad lines associated with the accreting black hole are expected to show systematic velocity shifts with respect to the narrow lines, which trace the rest frame of the galaxy. For a sample of 54,586 quasars and 3929 galaxies at redshifts 0.1 < z < 1.5, we brute-force model each spectrum as a mixture of two quasars at two different redshifts. The spectral model is a data-driven dimensionality reduction of the SDSS quasar spectra based on a matrix factorization. We identified 32 objects with peculiar spectra. Nine of them can be interpreted as black hole binaries. This doubles the number of known black hole binary candidates. We also report on the discovery of a new class of extreme double-peaked emitters with exceptionally broad and faint Balmer lines. For all the interesting sources, we present detailed analysis of the spectra and discuss possible interpretations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

THE AROMATIC FEATURES IN VERY FAINT DWARF GALAXIES

We present optical and mid-infrared photometry of a statistically complete sample of 29 dwarf galaxies (M{sub r} > - 15 mag) selected from the Sloan Digital Sky Survey (SDSS) spectroscopic sample and observed in the mid-infrared with Spitzer IRAC. This sample contains nearby (redshift {approx}<0.005) galaxies 3 mag fainter than previously studied samples. We compare our sample with other star-forming galaxies that have been observed with both IRAC and SDSS. We examine the relationship of the infrared color, [3.6]-[7.8], sensitive to polycyclic aromatic hydrocarbon (PAH) abundance and also hot dust and stellar continuum, with star formation rates (SFRs), oxygen abundances, and radiation hardness, all estimated by optical emission lines. Consistent with studies of more luminous dwarfs, we find that these dwarf galaxies show much redder [3.6]-[7.8] color than luminous galaxies with similar specific SFRs. Unlike luminous galaxies, we find that these dwarf galaxies show no significant dependence at all of the [3.6]-[7.8] color on SFR, oxygen abundance, or radiation hardness, despite the fact that the sample spans a significant range in all of these quantities. When the dwarfs in our sample are compared with more luminous dwarfs, we find that the [3.6]-[7.8] color, potentially tracing the PAH emission, depends on oxygen abundance and radiation hardness. However, these two parameters are correlated with one another as well; we break this degeneracy by looking at the PAH-oxygen abundance relation at a fixed radiation hardness and the PAH-hardness relation at a fixed oxygen abundance. This test shows that the [3.6]-[7.8] color in dwarf galaxies appears to depend more directly on oxygen abundance based on the data currently available.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

ARE THE ULTRA-FAINT DWARF GALAXIES JUST CUSPS?

We develop a technique to investigate the possibility that some of the recently discovered ultra-faint dwarf satellites of the Milky Way might be cusp caustics rather than gravitationally self-bound systems. Such cusps can form when a stream of stars folds, creating a region where the projected two-dimensional surface density is enhanced. In this work, we construct a Poisson maximum likelihood test to compare the cusp and exponential models of any substructure on an equal footing. We apply the test to the Hercules dwarf (d {approx} 113 kpc, M{sub V} {approx} -6.2, e {approx} 0.67). The flattened exponential model is strongly favored over the cusp model in the case of Hercules, ruling out at high confidence that Hercules is a cusp catastrophe. This test can be applied to any of the Milky Way dwarfs, and more generally to the entire stellar halo population, to search for the cusp catastrophes that might be expected in an accreted stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2010

INFERRING THE ECCENTRICITY DISTRIBUTION

Standard maximum-likelihood estimators for binary-star and exoplanet eccentricities are biased high, in the sense that the estimated eccentricity tends to be larger than the true eccentricity. As with most non-trivial observables, a simple histogram of estimated eccentricities is not a good estimate of the true eccentricity distribution. Here, we develop and test a hierarchical probabilistic method for performing the relevant meta-analysis, that is, inferring the true eccentricity distribution, taking as input the likelihood functions for the individual star eccentricities, or samplings of the posterior probability distributions for the eccentricities (under a given, uninformative prior). The method is a simple implementation of a hierarchical Bayesian model; it can also be seen as a kind of heteroscedastic deconvolution. It can be applied to any quantity measured with finite precision-other orbital parameters, or indeed any astronomical measurements of any kind, including magnitudes, distances, or photometric redshifts-so long as the measurements have been communicated as a likelihood function or a posterior sampling.

79 ASTRONOMY AND ASTROPHYSICS↗