Dominion Astrophysical Observatory

The tasks of the DAOPHOT program, developed to exploit the capability of photometrically linear image detectors to perform stellar photometry in crowded fields, are discussed. Raw CCD images are prepared prior to analysis, and following the obtaining of an initial star list with the FIND program, synthetic aperture photometry is performed on the detected objects with the PHOT routine. A local sky brightness and a magnitude are computed for each star in each of the specified stellar apertures, and for crowded fields, the empirical point-spread function must then be obtained for each data frame. The GROUP routine divides the star list for a given frame into optimum subgroups, and then the NSTAR routine is used to obtain photometry for all the stars in the frame by means of least- squares profile fits. The process is illustrated with images of stars in a crowded field, and shortcomings and possible improvements of the program are considered.

We present here the final results of the Hubble Space Telescope (HST) Key Project to measure the Hubble constant. We summarize our method, the results, and the uncertainties, tabulate our revised distances, and give the implications of these results for cosmology. Our results are based on a Cepheid calibration of several secondary distance methods applied over the range of about 60–400 Mpc. The analysis presented here benefits from a number of recent improvements and refinements, including (1) a larger LMC Cepheid sample to define the fiducial period-luminosity (PL) relations, (2) a more recent HST Wide Field and Planetary Camera 2 (WFPC2) photometric calibration, (3) a correction for Cepheid metallicity, and (4) a correction for incompleteness bias in the observed Cepheid PL samples. We adopt a distance modulus to the LMC (relative to which the more distant galaxies are measured) of μ0 = 18.50 ± 0.10 mag, or 50 kpc. New, revised distances are given for the 18 spiral galaxies for which Cepheids have been discovered as part of the Key Project, as well as for 13 additional galaxies with published Cepheid data. The new calibration results in a Cepheid distance to NGC 4258 in better agreement with the maser distance to this galaxy. Based on these revised Cepheid distances, we find values (in km s-1 Mpc-1) of H0 = 71 ± 2 ± 6 (systematic) (Type Ia supernovae), H0 = 71 ± 3 ± 7 (Tully-Fisher relation), H0 = 70 ± 5 ± 6 (surface brightness fluctuations), H0 = 72 ± 9 ± 7 (Type II supernovae), and H0 = 82 ± 6 ± 9 (fundamental plane). We combine these results for the different methods with three different weighting schemes, and find good agreement and consistency with H0 = 72 ± 8 km s-1 Mpc-1. Finally, we compare these results with other, global methods for measuring H0.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The comoving luminosity density of the Universe is estimated from the CFRS faint galaxy sample in three wavebands (2800A, 4400A and 1 micron) over the redshift range 0 < z < 1. In all three wavebands, the comoving luminosity density increases markedly with redshift. For a (q_0 = 0.5, Omega = 1.0) cosmological model, the comoving luminosity density increases as $(1+z)^{2.1 \pm 0.5}$ at 1 micron, as $(1+z)^{2.7 \pm 0.5}$ at 4400A and as $(1+z)^{3.9 \pm 0.75}$ at 2800A, these exponents being reduced by 0.43 and 1.12 for (0.05,0.1) and (-0.85,0.1) cosmological models respectively. The variation of the luminosity density with epoch can be reasonably well modelled by an actively evolving stellar population with a Salpeter initial mass function (IMF) extending to 125 M_sun, a star-formation rate declining with a power 2.5, and a turn-on of star-formation at early epochs. A Scalo (1986) IMF extending to the same mass limit produces too many long-lived low mass stars. This rapid evolution of the star-formation rate and comoving luminosity density of the Universe is in good agreement with the conclusions of Pei and Fall (1995) from their analysis of the evolving metallicity of the Universe. One consequence of this evolution is that the physical luminosity density at short wavelengths has probably declined by two orders of magnitude since z ~ 1.

We measure spectral indexes for 1823 galaxies in the Canadian Network for Observational Cosmology 1 (CNOC1) sample of 15 X-ray luminous clusters at 0.18 < z < 0.55 to investigate the mechanisms responsible for differential evolution between galaxy cluster and field environments. The radial trends of D4000, W0(Hδ), and W0(O ) are all consistent with an age sequence, in the sense that the last episode of star formation occurred more recently in galaxies farthest from the cluster center. Throughout the cluster environment, galaxies show evidence for older stellar populations than field galaxies; they have weaker W0(O ) and W0(Hδ) lines and stronger D4000 indexes. From our primary sample of 1413 galaxies, statistically corrected for incompleteness and selection effects, we identify a sample of K+A galaxies, which have strong Hδ absorption lines [W0(Hδ) > 5 Å] but no [O II] emission [W0(O ) < 5 Å], perhaps indicative of recently terminated star formation. The observed fraction of 4.4% ± 0.7% in the cluster sample is an overestimate due to a systematic effect that results from the large uncertainties on individual spectral index measurements. Corrected for this bias, we estimate that K+A galaxies make up only 2.1% ± 0.7% of the cluster sample and 0.1% ± 0.7% of the field. From the subsample of galaxies more luminous than Mr = -18.8 + 5 log h, which is statistically representative of a complete sample to this limit, the corrected fraction of K+A galaxies is 1.5% ± 0.8% in the cluster and 1.2% ± 0.8% in the field. Compared with the z ≈ 0.1 fraction of 0.30%, the fraction of K+A galaxies in the CNOC1 field sample is greater by perhaps a factor of 4, but with only 1 σ significance; no further evolution of this fraction is detectable over our redshift range. We compare our data with the results of PEGASE and GISSEL96 spectrophotometric models and conclude, from the relative fractions of red and blue galaxies with no [O II] λ3727 emission and strong Hδ absorption, that up to 1.9% ± 0.8% of the cluster population may have had its star formation recently truncated without a starburst. However, this is still not significantly greater than the fraction of such galaxies in the field, 3.1% ± 1.0%. Furthermore, we do not detect an excess of cluster galaxies that have unambiguously undergone starbursts within the last 1 Gyr. In fact, at 6.3% ± 2.1%, the A+em galaxies that Poggianti et al. have recently suggested are dusty starbursts are twice as common in the field as in the cluster environment. Our results imply that these cluster environments are not responsible for inducing starbursts; thus, the increase in cluster blue galaxy fraction with redshift may not be a strictly cluster-specific phenomenon. We suggest that the truncation of star formation in clusters may largely be a gradual process, perhaps due to the exhaustion of gas in the galactic disks over fairly long timescales; in this case differential evolution may result because field galaxies can refuel their disks with gas from extended halos, thus regenerating star formation, while cluster galaxies may not have such halos and so continue to evolve passively.

We examine the origin of clustercentric gradients in the star formation rates and colors of rich cluster galaxies within the context of a simple model where clusters are built through the ongoing accretion of field galaxies. The model assumes that after galaxies enter the cluster their star formation rates decline on a timescale of a few Gyrs, the typical gas consumption timescale of disk galaxies in the field. Such behaviour might be expected if tides and ram pressure strip off the gaseous envelopes that normally fuel star formation in spirals over a Hubble time. Combining these timescales with mass accretion histories derived from N-body simulations of cluster formation in a Lambda-CDM universe, we reproduce the systematic differences observed in the color distribution of cluster and field galaxies, as well as the strong suppression of star formation in cluster galaxies and its dependence on clustercentric radius. The simulations also indicate that a significant fraction of galaxies beyond the virial radius of the cluster may have been within the main body of the cluster in the past, a result that explains naturally why star formation in the outskirts of clusters (and as far out as two virial radii) is systematically suppressed relative to the field. The agreement with the data beyond the cluster virial radius is also improved if we assume that stripping happens within lower mass systems, before the galaxy is accreted into the main body of the cluster. We conclude that the star formation rates of cluster galaxies depend primarily on the time elapsed since their accretion onto massive virialized systems, and that the cessation of star formation may have taken place gradually over a few Gyrs.

The Ninth Catalogue of Spectroscopic Binary Orbits (http://sb9.astro.ulb.ac.be) continues the series of compilations of spectroscopic orbits carried out over the past 35 years by Batten and collaborators. As of 2004 May 1st, the new Catalogue holds orbits for 2386 systems. Some essential differences between this catalogue and its predecessors are outlined and three straightforward applications are presented: (1) completeness assessment: period distribution of SB1s and SB2s; (2) shortest periods across the H-R diagram; (3) period-eccentricity relation.

The 10.7 cm solar radio flux, or F 10.7 is, along with sunspot number, one of the most widely used indices of solar activity. This paper describes the equipment and procedures used to make the measurements and to calibrate them, and discusses some of the “most‐asked” questions about the data.

Stars appearing in CCD images obtained over 224 nights during the course of 69 observing runs have been calibrated to the Johnson/Kron-Cousins BVRI photometric system defined by the equatorial standards of Landolt (1992, AJ, 104, 340). More than 15,000 stars suitable for use as photometric standards have been identified, where "suitable" means that the star has been observed five or more times during photometric conditions and has a standard error of the mean magnitude less than 0.02 mag in at least two of the four bandpasses, and shows no significant evidence of intrinsic variability. Many of these stars are in the same fields as Landolt's equatorial standards or Graham's (1982, PASP, 94, 244) southern E-region standards, but are considerably fainter. This enhances the value of those fields for the calibration of photometry obtained with large telescopes. Other standards have been defined in fields containing popular objects of astrophysical interest, such as star clusters and famous galaxies, extending Landolt-system calibrators to declinations far from the equator and to stars of sub-Solar chemical abundances. I intend to continue to improve and enlarge this set of photometric standard stars as more observing runs are reduced. The full current database of photometric indices is being made freely available via a site on the World-Wide Web, or by direct request to the author. Although the contents of the database will evolve in detail, at any given time it should represent the largest sample of precise BVRI broad-band photometric standards available anywhere.

Abstract The derivation of a closed expression is presented to calculate the vertical component of the gravitational attraction of a right rectangular prism, with sides parallel to the coordinate axis. As any configuration can be expressed as the sum of prisms of various sizes and densities, the computation of the total gravitational effect of bodies of arbitrary shapes at any point outside of or on the boundary of the bodies is straightforward. To calculate the gravitational effect of the 'unit' building element a subroutine called Prism has been developed, tested, and incorporated, in one program to calculate terrain corrections, and in another program for three-dimensional analysis of a gravity field.

We are obtaining spectra, spectral types and basic physical parameters for the nearly 3600 dwarf and giant stars earlier than M0 in the Hipparcos catalog within 40pc of the Sun. Here we report on results for 1676 stars in the southern hemisphere observed at Cerro Tololo Interamerican Observatory and Steward Observatory. These results include new, precise, homogeneous spectral types, basic physical parameters (including the effective temperature, surface gravity, and metallicity, [M/H]) and measures of the chromospheric activity of our program stars. We include notes on astrophysically interesting stars in this sample, the metallicity distribution of the solar neighborhood and a table of solar analogues. We also demonstrate that the bimodal nature of the distribution of the chromospheric activity parameter log R'HK depends strongly on the metallicity, and we explore the nature of the ``low-metallicity'' chromospherically active K-type dwarfs.

(abridged) We measure spectral indices for 1823 galaxies in the CNOC1 sample of fifteen X-ray luminous clusters at 0.18

Since the launch of Hubble Space Telescope (HST) 9 yr ago, Cepheid distances to 25 galaxies have been determined for the purpose of calibrating secondary distance indicators. Eighteen of these have been measured by the HST Key Project team, six by the Supernova Calibration Project, and one independently by Tanvir. Collectively, this work sets out an array of survey markers over the region within 25 Mpc of the Milky Way. A variety of secondary distance indicators can now be calibrated, and the accompanying four papers employ the full set of 25 galaxies to consider the Tully-Fisher relation, the fundamental plane of elliptical galaxies, Type Ia supernovae, and surface brightness fluctuations. When calibrated with Cepheid distances, each of these methods yields a measurement of the Hubble constant and a corresponding measurement uncertainty. We combine these measurements in this paper, together with a model of the velocity field, to yield the best available estimate of the value of H0 within the range of these secondary distance indicators and its uncertainty. The uncertainty in the result is modeled in an extensive simulation we call the 'virtual Key Project.' The velocity-field model includes the influence of the Virgo cluster, the Great Attractor, and the Shapley supercluster, but does not play a significant part in determining the result. The result is H0 = 71 ± 6 km s-1 Mpc-1. The largest contributor to the uncertainty of this 67% confidence level result is the distance of the Large Magellanic Cloud, which has been assumed to be 50 ± 3 kpc. This takes up the first 6.5% of our 9% error budget. Other contributors are the photometric calibration of the WFPC2 instrument, which takes up 4.5%, deviations from uniform Hubble flow in the volume sampled (lesssim2%), the composition sensitivity of the Cepheid period-luminosity relation (4%), and departures from a universal reddening law (~1%). These are the major components that , when combined in quadrature, make up the 9% total uncertainty. If the LMC distance modulus were systematically smaller by 1 σ than that adopted here, the derived value of the Hubble constant would increase by 4 km s-1 Mpc-1. Most of the significant systematic errors are capable of amelioration in future work. These include the uncertainty in the photometric calibration of WFPC2, the LMC distance, and the reddening correction. A NICMOS study is in its preliminary reduction phase, addressing the last of these concerns. Various empirical analyses have suggested that Cepheid distance moduli are affected by metallicity differences. If we adopted the composition sensitivity obtained in the Key Project's study of M101, and employed the oxygen abundances measured spectroscopically in each of the Cepheid fields we have studied, the value of the Hubble constant would be reduced by 4% ± 2% to 68 ± 6 km s-1 Mpc-1.

The morphological properties of galaxies in the range 21 < I < 25 mag in the Hubble Deep Field are investigated using a quantitative classification system based on measurements of the central concentration and asymmetry of galaxian light. The class distribution of objects in the Hubble Deep Field is strongly skewed towards highly asymmetric objects, relative to distributions from both the HST Medium Deep Survey at I<22 mag and an artificially redshifted sample of local galaxies. The steeply rising number count-magnitude relation for irregular/peculiar/merging systems at I < 22 mag reported by Glazebrook et al. continues to at least I = 25 mag. Although these peculiar systems are predominantly blue at optical wavelengths, a significant fraction also exhibit red U–B colours, which may indicate that they are at high redshift. Beyond Glazebrook et al.'s magnitude limit, the spiral counts appear to rise more steeply than high-normalization no-evolution predictions, whereas those of elliptical/S0 galaxies only slightly exceed such predictions and may turn over beyond I ∼ 24 mag. These results are compared with those from previous investigations of faint galaxy morphology with HST, and the possible implications are briefly discussed. The large fraction of peculiar/irregular/merging systems in the Hubble Deep Field suggests that by I ∼ 25 mag the conventional Hubble system no longer provides an adequate description of the morphological characteristics of a high fraction of field galaxies.

We have embarked on a project, under the aegis of the Nearby Stars (NStars)/ Space Interferometry Mission Preparatory Science Program to obtain spectra, spectral types, and, where feasible, basic physical parameters for the 3600 dwarf and giant stars earlier than M0 within 40 parsecs of the sun. In this paper we report on the results of this project for the first 664 stars in the northern hemisphere. These results include precise, homogeneous spectral types, basic physical parameters (including the effective temperature, surface gravity and the overall metallicity, [M/H]) and measures of the chromospheric activity of our program stars. Observed and derived data presented in this paper are also available on the project's website at http://stellar.phys.appstate.edu/ .

Understanding diffuse Galactic radio emission is interesting both in its own right and for minimizing foreground contamination of cosmological measurements. cosmic microwave background experiments have focused on frequencies 10 GHz, whereas 21-cm tomography of the high-redshift universe will mainly focus on 0.2 GHz, for which less is currently known about Galactic emission. Motivated by this, we present a global sky model derived from all publicly available total power large-area radio surveys, digitized with optical character recognition when necessary and compiled into a uniform format, as well as the new Villa Elisa data extending the 1.42-GHz map to the entire sky. We quantify statistical and systematic uncertainties in these surveys by comparing them with various global multifrequency model fits. We find that a principal component based model with only three components can fit the 11 most accurate data sets (at 10, 22, 45 and 408 MHz and 1.42, 2.326, 23, 33, 41, 61, 94 GHz) to an accuracy around 1-10 per cent depending on frequency and sky region. Both our data compilation and our software returning a predicted all-sky map at any frequency from 10 MHz to 100 GHz are publicly available at http://space.mit.edu/home/angelica/gsm.

We describe observations and abundance analysis of a high-resolution, high-S/N survey of 168 stars, most of which are metal-poor dwarfs. We follow a self-consistent LTE analysis technique to determine the stellar parameters and abundances, and estimate the effects of random and systematic uncertainties on the resulting abundances. Element-to-iron ratios are derived for key alpha, odd, Fe-peak, r- and s-process elements. Effects of Non-LTE on the analysis of Fe I lines are shown to be very small on the average. Spectroscopically determined surface gravities are derived that are quite close to those obtained from Hipparcos parallaxes. 1.

The quantitative morphological classification of distant galaxies is essential to the understanding of the evolution of galaxies over the history of the Universe. This paper presents Hubble Space Telescope WFPC2 F606W and F814W photometric structural parameters for 7450 galaxies in the ``Groth Strip.'' These parameters are based on a two-dimensional bulge+disk surface brightness model and were obtained using an automated reduction and analysis pipeline described in detail here. A first set of fits was performed separately in each bandpass, and a second set of fits was performed simultaneously on both bandpasses. The information produced by these two types of fits can be used to explore different science goals. Systematic and random fitting errors in all structural parameters as well as bulge and disk colors are carefully characterized through extensive sets of simulations. The results of these simulations are given in catalogs similar to the real science catalogs so that both real and simulated measurements can be sampled according to the same selection criteria to show biases and errors in the science data subset of interest. The effects of asymmetric structures on the recovered bulge+disk fitting parameters are also explored through simulations. The full multidimensional photometric survey selection function of the Groth Strip is also computed. This selection function, coupled to bias maps from simulations, provides a complete and objective reproduction of the observational limits, and these limits can be applied to theoretical predictions from galaxy evolution models for direct comparisons with the data.

In this paper we present a new statistic for quantifying galaxy morphology based on measurements of the Gini coefficient of galaxy light distributions. This statistic is easy to measure and is commonly used in econometrics to measure how wealth is distributed in human populations. When applied to galaxy images, the Gini coefficient provides a quantitative measure of the inequality with which a galaxy's light is distributed amongst its constituent pixels. We measure the Gini coefficient of local galaxies in the Early Data Release of the Sloan Digital Sky Survey and demonstrate that this quantity is closely correlated with measurements of central concentration, but with significant scatter. This scatter is almost entirely due to variations in the mean surface brightness of galaxies. By exploring the distribution of galaxies in the three-dimensional parameter space defined by the Gini coefficient, central concentration, and mean surface brightness, we show that all nearby galaxies lie on a well-defined two-dimensional surface (a slightly warped plane) embedded within a three-dimensional parameter space. By associating each galaxy sample with the equation of this plane, we can encode the morphological composition of the entire SDSS g-band sample using the following three numbers: 22.451, 5.366, 7.010. The i-band sample is encoded as: 22.149, 5.373, and 7.627.

The bimodal globular cluster (GC) metallicity distributions of many giant elliptical galaxies are often cited as evidence for the formation of such galaxies through mergers involving gas-rich spirals. In such models, the metal- rich GCs are assumed to have formed during the merger process. We explore an alternative possibility: that these metal-rich clusters represent the galaxy's intrinsic GC population and that the metal-poor component of the observed GC metallicity distribution arises from the capture of GCs from other galaxies, either through mergers or through tidal stripping. Starting with plausible assumptions for the initial galaxy luminosity function and for the dependence of GC metallicity on parent galaxy luminosity, we show that the growth of a pre-existing seed galaxy through mergers and tidal stripping is accompanied by the capture of metal-poor GCs whose properties are similar to those which are observed to surround giant ellipticals. We describe a method of using the observed number of metal-poor and metal-rich GCs to infer the merger histories of individual elliptical galaxies, and use this technique to derive limits on the number of galaxies and total luminosity accreted to date by M49. We argue that although GC specific frequency is conserved in galaxy mergers, the same may not be true of tidal stripping by the mean field of the host galaxy cluster. Comparisons of model GC metallicity distributions and specific frequencies to those observed for the well-studied galaxies M49 and M87 show that it is possible to explain their bimodal GC metallicity distributions and discordant specific frequencies without resorting to the formation of new GCs in mergers or by invoking multiple bursts of GC formation.