Search

Over the past decades and even centuries, the astronomical community has accumulated a signif-icant heritage of recorded observations of a great many astronomical objects. Those records con-tain irreplaceable information about long-term evolutionary and non-evolutionary changes in our Universe, and their preservation and digitization is vital. Unfortunately, most of those data risk becoming degraded and thence totally lost. We hereby call upon the astronomical community and US funding agencies to recognize the gravity of the situation, and to commit to an interna-tional preservation and digitization efforts through comprehensive long-term planning supported by adequate resources, prioritizing where the expected scientific gains, vulnerability of the origi-nals and availability of relevant infrastructure so dictates. The importance and urgency of this issue has been recognized recently by General Assembly XXX of the International Astronomical Union (IAU) in its Resolution B3: "on preservation, digitization and scientific exploration of his-torical astronomical data". We outline the rationale of this promotion, provide examples of new science through successful recovery efforts, and review the potential losses to science if nothing it done., Comment: 8 pages, White Paper submitted to ASTRO2020 Decadal Survey

A complete census of Blue Horizontal Branch (BHB) and Blue Straggler Star (BSS) population within the 10$'$ radius from the center of the Globular Cluster, NGC 288 is presented, based on the images from the Ultraviolet Imaging Telescope (UVIT). The UV and UV$-$optical Colour-Magnitude Diagrams (CMDs) are constructed by combining the UVIT, HST-ACS and ground data and compared with the BaSTI isochrones generated for UVIT filters. We used stellar proper motions data from GAIA DR2 to select the cluster members. Our estimations of the temperature distribution of 110 BHB stars reveal two peaks with the main peak at $T_{eff}\sim$ 10,300 K with the distribution extending up to $T_{eff}\sim$ 18,000 K. We identify the well known photometric gaps including the G-jump in the BHB distribution which are located between the peaks. We detect a plateau in the FUV magnitude for stars hotter than $T_{eff}\sim$ 11,500 K (G-jump), which could be due to the effect of atomic diffusion. We detect 2 Extreme HB (EHB) candidates with temperatures ranging from 29,000 to 32,000 K. The radial distribution of 68 BSSs suggests that the bright BSSs are more centrally concentrated than the faint BSS and the BHB distribution. We find that the BSSs have a mass range of 0.86 - 1.25 M$_{\odot}$ and an age range of 2 - 10 Gyr with a peak at 1 M$_{\odot}$ and 4 Gyr respectively. This study showcases the importance of combining UVIT with HST, ground, and GAIA data in deriving HB and BSS properties., Comment: Accepted in MNRAS; 17 pages, 17 figures and 5 Tables. Full catalogue (Table 2) will be available in electronic format

The analysis of the early star formation history (SFH) of nearby galaxies, obtained from their resolved stellar populations is relevant as a test for cosmological models. However, the early time resolution of observationally derived SFHs is limited by several factors. Thus, direct comparison of observationally derived SFHs with those derived from theoretical models of galaxy formation is potentially biased. Here we investigate and quantify this effect. For this purpose, we analyze the duration of the early star formation activity in a sample of four Local Group dwarf galaxies and test whether they are consistent with being true fossils of the pre-reionization era; i.e., if the quenching of their star formation occurred before cosmic reionization by UV photons was completed. Two classical dSph (Cetus and Tucana) and two dTrans (LGS-3 and Phoenix) isolated galaxies with total stellar masses between $1.3\times 10^6$ to $7.2\times 10^6$ M$_\odot$ have been studied. Accounting for time resolution effects, the SFHs peak as much as 1.25 Gyr earlier than the optimal solutions. Thus, this effect is important for a proper comparison of model and observed SFHs. It is also shown that none of the analyzed galaxies can be considered a true-fossil of the pre-reionization era, although it is possible that the {\it outer regions} of Cetus and Tucana are consistent with quenching by reionization., Comment: To be published by the ApJ

It is widely accepted that stars do not form in isolation but result from the fragmentation of molecular clouds, which in turn leads to star cluster formation. Over time, clusters dissolve or are destroyed by interactions with molecular clouds or tidal stripping, and their members become part of the general field population. Star clusters are thus among the basic building blocks of galaxies. In turn, star cluster populations, from young associations and open clusters to old globulars, are powerful tracers of the formation, assembly, and evolutionary history of their parent galaxies. Although their importance had been recognised for decades, major progress in this area has only become possible in recent years, both for Galactic and extragalactic cluster populations. Star clusters are the observational foundation for stellar astrophysics and evolution, provide essential tracers of galactic structure, and are unique stellar dynamical environments. Star formation, stellar structure, stellar evolution, and stellar nucleosynthesis continue to benefit and improve tremendously from the study of these systems. Additionally, fundamental quantities such as the initial mass function can be successfully derived from modelling either the H-R diagrams or the integrated velocity structures of, respectively, resolved and unresolved clusters and cluster populations. Star cluster studies thus span the fields of Galactic and extragalactic astrophysics, while heavily affecting our detailed understanding of the process of star formation in dense environments.This report highlights science results of the last decade in the major fields covered by IAU Commission 37: Star clusters and associations., Comment: 17 pages, to be published in the Transactions of the IAU, by the International Astronomical Union. Abstract reduced to fit in

We have used a combination of broad-band near-infrared and optical Johnson-Cousins photometry to study the dust properties in the line of sight to the Galactic globular cluster M4. We have investigated the reddening effects in terms of absolute strength and variation across the cluster field, as well as the form of the reddening law defined by the type of dust. Here, we determine the ratio of absolute to selective extinction (R_V) in the line of sight towards M4, which is known to be a useful indicator for the type of dust and therefore characterizes the applicable reddening law. Our method is independent of age assumptions and appears to be significantly more precise and accurate than previous approaches. We obtain A_V/E(B-V)=3.76\pm0.07 (random error) for the dust in the line of sight to M4 for our set of filters. With this value, the distance to M4 is found to be 1.80\pm0.05 kpc (random error). A reddening map for M4 has been created, which reveals a spatial differential reddening of delta E(B-V)>0.2 mag across a field within 10' around the cluster centre and a total mean reddening of E(B-V)=0.37\pm0.01. In order to provide accurate zero points for the extinction coefficients of our photometric filters, we investigated the impact of stellar parameters such as temperature, surface gravity and metallicity on the extinction properties in different bandpasses. Using both synthetic ATLAS9 spectra and observed spectral energy distributions, we found similar sized effects for the range of temperature, surface gravity, and metallicity typical of globular cluster stars: each causes a change of about 3% in the necessary correction factor for each filter combination. From our calculations, we provide extinction zero points for Johnson-Cousins and 2MASS filters, spanning a wide range of stellar parameters and dust types, suited for accurate, object-specific extinction corrections., 24 pages, 19 figures, 13 tables

We present an analysis of the star formation history (SFH) of the transition-type (dIrr/dSph) Local Group galaxy \objectname[]{LGS-3} (Pisces) based on deep photometry obtained with the {\it Advanced Camera for Surveys} onboard the {\it Hubble Space Telescope}. Our analysis shows that the SFH of \objectname[]{LGS-3} is dominated by a main episode $\sim 11.7$ Gyr ago with a duration of $\sim$ 1.4 Gyr which formed $\sim 90%$ of the stars. Subsequently, \objectname[]{LGS-3} continued forming stars until the present, although at a much lower rate. The lack of early chemical enrichment is in contrast to that observed in the isolated dSph galaxies of comparable luminosity, implying that the dSphs were more massive and subjected to more tidal stripping. We compare the SFH of \objectname[]{LGS-3} with expectations from cosmological models. Most or all the star formation was produced in \objectname[]{LGS-3} after the reionization epoch, assumed to be completed at $z\sim6$ or $\sim 12.7$ Gyr ago. The total mass of the galaxy is estimated to be between 2 and $4\times 10^8$ M$_\odot$, corresponding to circular velocities between 28 km\ s$^{-1}$ to 36 km\ s$^{-1}$. These values are close to but somewhat above the limit of 30 km\ s$^{-1}$ below which the UV background is expected to prevent any star formation after reionization. Feedback from SNe associated with the initial episode of star formation (mechanical luminosity from SNe $L_w=5.3\times 10^{38}$ erg s$^{-1}$) is probably inadequate to completely blow away the gas. However, the combined effects of SN feedback and UV background heating might be expected to completely halt star formation at the reionization epoch for the low mass of \objectname[]{LGS-3}; this suggests that self-shielding is important to the early evolution of galaxies in this mass range., Comment: accepted to be published at ApJ

NGC 2808 is a unique globular cluster with not only a bimodal-horizontal branch (HB) but also with gaps on the blue horizontal branch. Adequate interpretation of the nature of the detected peculiarities in "bimodal" and "gap" clusters is of paramount importance for understanding the nature of the second parameter phenomenon and scenarios for the formation of the Galaxy. Although RR Lyrae variables are HB stars that can provide powerful constraints to models on the origin of bimodal HBs, unfortunately, until recently, only one RR Lyrae variable was known in this cluster. Here we present the first calibrated time-series CCD photometry for newly discovered fundamental mode RR Lyrae variables in NGC 2808, with observations over a range of twenty years. Investigations of RR Lyrae variable stars in this peculiar, bimodal-horizontal branch globular cluster are carried out to account for its formation, and the effects of helium enrichment and differential reddening. The Oosterhoff classification of NGC$\,$2808, which has recently been associated with a previously unknown dwarf galaxy in Canis Major, is also discussed., Comment: To appear in the Carnegie Observatories Astrophysics Series, Volume 5

(abridged) We present the first study of the variable star populations in the isolated dwarf spheroidal galaxies (dSph) Cetus and Tucana. Based on Hubble Space Telescope images obtained with the Advanced Camera for Surveys in the F475W and F814W bands, we identified 180 and 371 variables in Cetus and Tucana, respectively. The vast majority are RR Lyrae stars. In Cetus we also found three anomalous Cepheids, four candidate binaries and one candidate long-period variable (LPV), while six anomalous Cepheids and seven LPV candidates were found in Tucana. Of the RR Lyrae stars, 147 were identified as fundamental mode (RRab) and only eight as first-overtone mode (RRc) in Cetus, with mean periods of 0.614 and 0.363 day, respectively. In Tucana we found 216 RRab and 82 RRc giving mean periods of 0.604 and 0.353 day. These values place both galaxies in the so-called Oosterhoff Gap, as is generally the case for dSph. We calculated the distance modulus to both galaxies using different approaches based on the properties of RRab and RRc, namely the luminosity-metallicity and period-luminosity-metallicity relations, and found values in excellent agreement with previous estimates using independent methods: (m-M)_{0,Cet}=24.46+-0.12 and (m-M)_{0,Tuc}=24.74+-0.12, corresponding to 780+-40 kpc and 890+-50 kpc. We also found numerous RR Lyrae variables pulsating in both modes simultaneously (RRd): 17 in Cetus and 60 in Tucana. Tucana is, after Fornax, the second dSph in which such a large fraction of RRd (~17%) has been observed. We provide the photometry and pulsation parameters for all the variables, and compare the latter with values from the literature for well-studied dSph of the Local Group and Galactic globular clusters., 26 pages, 24 figures, in emulateapj format. To be published in ApJ. Some figures heavily degraded; See http://www.iac.es/project/LCID/?p=publications for a version with full resolution figures

We describe here the reduction methods that we developed to study the faintest red dwarfs and white dwarfs in an outer field of NGC6397, which was observed by \hst for 126 orbits in 2005. The particular challenge of this data set is that the faintest stars are not readily visible in individual exposures, so special care must be taken to combine the information in all the exposures in order to identify and measure them. Unfortunately, it is hard to find the faintest stars without also finding a large number of faint galaxies, so we developed specialized tools to distinguish between the point-like stars and the barely resolved galaxies. We found that artificial-star tests, while obviously necessary for completeness determination, can also play an important role in helping us optimize our finding and measuring algorithms. Although this paper focuses on this data set specifically, many of the techniques are new and might find application in other work, particularly when a large number of images is available for a single field., Comment: 35 pages, preprint style

We investigate in detail the white dwarf cooling sequence of the globular cluster Messier 4. In particular we study the influence of various systematic uncertainties, both observational and theoretical, on the determination of the cluster age from the white dwarf cooling sequence. These include uncertainties in the distance to the cluster and the extinction along the line of sight, as well as the white dwarf mass, envelope and core compositions and the white dwarf --main sequence mass relation. We find that fitting to the full two-dimensional colour-magnitude diagram offers a more robust method for age determination than the traditional method of fitting the one-dimensional white dwarf luminosity function. After taking into account the various uncertainties, we find a best fit age of 12.1 Gyr, with a 95% lower limit of 10.3 Gyr. We also perform fits using two other sets of cooling models from the literature. The models of Chabrier et al (2000) yield an encouragingly similar result, although the models of Salaris et al (2000) do not provide as good a fit. Our results support our previous determination of a delay between the formation of the Galactic halo and the onset of star formation in the Galactic disk., Comment: many pages, 43 postscript figures, submitted to ApJ

[Abridged] We present the results of a Hubble Space Telescope snapshot program to survey the inner region of the globular cluster NGC 6441 for its variable stars. A total of 57 variable stars was found including 38 RR Lyrae stars, 6 Population II Cepheids, and 12 long period variables. Of the RR Lyrae stars observed in this survey, 26 are pulsating in the fundamental mode with a mean period of 0.753d and 12 are first-overtone mode pulsators with a mean period of 0.365d. These values match up very well with those found in ground-based surveys. Combining all the available data for NGC 6441, we find mean periods of 0.759d and 0.375d for the RRab and RRc stars, respectively. We also find that the RR Lyrae in this survey are located in the same regions of a period-amplitude diagram as those found in ground-based surveys. Although NGC 6441 is a metal-rich globular cluster, its RR Lyrae more closely resemble those in Oosterhoff type II globular clusters. However, even compared to typical Oosterhoff type II systems, the mean period of its RRab stars is unusually long. We also derived I-band period-luminosity relations for the RR Lyrae stars. Of the six Population II Cepheids, five are of W Virginis type and one is a BL Herculis variable stars. This makes NGC 6441, along with NGC 6388, the most metal-rich globular cluster known to contain these types of variable stars. Another variable, V118, may also be a Population II Cepheid given its long period and its separation in magnitude from the RR Lyrae stars. We argue that there does not appear to be a change in the period-luminosity relation slope between the BL Herculis and W Virginis stars, but that a change of slope does occur when the RV Tauri stars are added to the period-luminosity relation., Comment: 28 pages, including 9 figures and 8 tables, emulateapj5/apjfonts style. Accepted by the Astronomical Journal. Approximate publication date September 2003. We recommend the interested reader to download the preprint with full-resolution figures, which can be found at http://www.astro.puc.cl/~mcatelan/Pritzl.zip

Using published photometry and spectroscopy, we construct the fundamental plane and D_n-Sigma relations in Leo I, Virgo and Fornax. The published Cepheid P-L relations to spirals in these clusters fixes the relation between angular size and metric distance for both the fundamental plane and D_n-Sigma relations. Using the locally calibrated fundamental plane, we infer distances to a sample of clusters with a mean redshift of cz \approx 6000 \kms, and derive a value of H_0=78+- 5+- 9 km/s/Mpc (random, systematic) for the local expansion rate. This value includes a correction for depth effects in the Cepheid distances to the nearby clusters, which decreased the deduced value of the expansion rate by 5% +- 5%. If one further adopts the metallicity correction to the Cepheid PL relation, as derived by the Key Project, the value of the Hubble constant would decrease by a further 6%+- 4%. These two sources of systematic error, when combined with a +- 6% error due to the uncertainty in the distance to the Large Magellanic Cloud, a +- 4% error due to uncertainties in the WFPC2 calibration, and several small sources of uncertainty in the fundamental plane analysis, combine to yield a total systematic uncertainty of +- 11%. We find that the values obtained using either the CMB, or a flow-field model, for the reference frame of the distant clusters, agree to within 1%. The Dn-Sigma relation also produces similar results, as expected from the correlated nature of the two scaling relations. A complete discussion of the sources of random and systematic error in this determination of the Hubble constant is also given, in order to facilitate comparison with the other secondary indicators being used by the Key Project., Comment: 21 pages, 3 figures, Accepted for publication in ApJ

We report on the detection of Cepheids and the first distance measurement to the spiral galaxy NGC 1326-A, a member of the Fornax cluster of galaxies. We have employed data obtained with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. Over a 49 day interval, a total of twelve V-band (F555W) and eight I-band (F814W) epochs of observation were obtained. Two photometric reduction packages, ALLFRAME and DoPHOT, have been employed to obtain photometry measures from the three Wide Field CCDs. Variability analysis yields a total of 17 Cepheids in common with both photometry datasets, with periods ranging between 10 and 50 days. Of these 14 Cepheids with high-quality lightcurves are used to fit the V and I period-luminosity relations and derive apparent distance moduli, assuming a Large Magellanic Cloud distance modulus (m-M) (LMC) = 18.50 +- 0.10 mag and color excess E(B-V) = 0.10 mag. Assuming A(V)/E(V-I) = 2.45, the DoPHOT data yield a true distance modulus to NGC 1326-A of (m-M)_0 = 31.36 +- 0.17 (random) +- 0.13 (systematic) mag, corresponding to a distance of 18.7 \pm 1.5 (random) \pm 1.2 (systematic) Mpc. The derived distance to NGC 1326-A is in good agreement with the distance derived previously to NGC 1365, another spiral galaxy member of the Fornax cluster. However the distances to both galaxies are significantly lower than to NGC 1425, a third Cepheid calibrator in the outer parts of the cluster., Comment: 33 pages A gzipped tar file containing 12 figures can be obtained from http://www.ipac.caltech.edu/H0kp/n1326a/n1326a.html

We report on the chemical abundances of stars in the Carina dwarf spheroidal galaxy (dSph) derived from low-resolution spectra. We have determined values of [Fe/H] for 52 stars from the reduced equivalent width of the Ca II infrared triplet lines. The Carina dSph has a mean metallicity of [Fe/H] = -1.99 +/- 0.08 and an intrinsic metallicity dispersion 0.25 dex (1 sigma). By directly determining the chemical abundances of Carina stars through spectroscopy, we can overcome the age-metallicity degeneracy inherent in color-magnitude diagrams (CMDs) and determine its star-formation history with unprecedented accuracy., Comment: 4 pages, 2 figures, to appear in in "Spectro-Photometric Dating of Stars and Galaxies", eds. I. Hubeny, S. Heap & R. Cornett, (ASP: San Francisco), ASP Conf. Series, in press

The final efforts of the HST Key Project on the Extragalactic Distance Scale are presented. Four distance indicators, the Surface Brightness Fluctuation method, the Fundamental Plane for early-type galaxies, the Tully-Fisher relation and the Type Ia Supernovae, are calibrated using Cepheid distances to galaxies within 25 Mpc. The calibration is then applied to distant samples reaching cz~10000 km/s and (in the case of SNIa) beyond. By combining the constraints imposed on the Hubble constant by the four distance indicators, we obtain H0 = 71+/-6 km/s/Mpc., Comment: To appear in the proceedings of the Cosmic Flows Workshop, Victoria, Canada, July 1999, eds. S. Courteau, M. Strauss & J. Willick, ASP series. 10 pages, 2 figures

Since the launch of the Hubble Space Telescope nine years ago, Cepheid distances to 25 galaxies have been determined for the purpose of calibrating secondary distance indicators. A variety of these can now be calibrated, and the accompanying papers by Sakai, Kelson, Ferrarese, and Gibson 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 H_0 within the range of these secondary distance indicators and its uncertainty. The result is H_0 = 71 +/- 6 km/sec/Mpc. 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 paper presents a detailed description of the acquisition and processing of a large body of imaging data for three fields in the globular cluster M4 taken with the Wide Field and Planetary Camera 2 aboard the Hubble Space Telescope. Analysis with the ALLFRAME package yielded the deepest photometry yet obtained for this cluster. The resulting data-set for 4708 stars (positions and calibrated photometry in V, I, and, in two fields, U) spanning approximately six cluster core radii is available on the AAS CD-ROM (or email a request to RAI). The scientific analysis is deferred to three companion papers, which investigate the significant white dwarf population discovered and the main sequence population., Comment: 14 pages, 7 figures, accepted for publication in ApJS

We report the detection of Cepheids and a new distance to the spiral galaxy NGC 2541, based on data obtained with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope (HST). A total of 25 exposures (divided into 13 epochs) are obtained using the F555W filter (transformed to Johnson V), and nine exposures (divided into five epochs) using the F814W filter (transformed to Cousins I). Photometric reduction of the data is performed using two independent packages, DoPHOT and DAOPHOT II/ALLFRAME, which give very good agreement in the measured magnitudes. A total of 34 bona fide Cepheids, with periods ranging from 12 to over 60 days, are identified based on both sets of photometry. By fitting V and I period-luminosity relations, apparent distance moduli are derived assuming a Large Magellanic Cloud distance modulus and mean color excess of 18.50 +/- 0.10 mag and E(B-V) = 0.10 mag respectively. Adopting A(V)/E(V-I)=2.45, we obtain a true distance modulus to NGC 2541 of 30.47 +/- 0.11 (random) +/- 0.12 (systematic) mag (D = 12.4 +/- 0.6 (random) +/- 0.7 (systematic) Mpc), and a total (Galactic plus internal) mean color excess E(B-V) = 0.08 +/- 0.05 (internal error) mag., Comment: 58 pages, 14 figures. Accepted for publication in the Astrophysical Journal. Figures 1, 4, 5, 8, 9 and 11 have been omitted in the astro-ph version

Using the Hubble Space Telescope (HST) 37 long-period Cepheid variables have been discovered in the Fornax Cluster spiral galaxy NGC 1365. The resulting V and I period-luminosity relations yield a true distance modulus of 31.35 +/- 0.07 mag, which corresponds to a distance of 18.6 +/- 0.6 Mpc. This measurement provides several routes for estimating the Hubble Constant. (1) Assuming this distance for the Fornax Cluster as a whole yields a local Hubble Constant of 70 +/-18_{random} [+/-7]_{systematic} km/s/Mpc. (2) Nine Cepheid-based distances to groups of galaxies out to and including the Fornax and Virgo clusters yield Ho = 73 (+/-16)_r [+/-7]_s km/s/Mpc. (3) Recalibrating the I-band Tully-Fisher relation using NGC 1365 and six nearby spiral galaxies, and applying it to 15 galaxy clusters out to 100 Mpc gives Ho = 76 (+/-3)_r [+/-8]_s km/s/Mpc. (4) Using a broad-based set of differential cluster distance moduli ranging from Fornax to Abell 2147 gives Ho = 72 (+/-)_r [+/-6]_s km/s/Mpc. And finally, (5) Assuming the NGC 1365 distance for the two additional Type Ia supernovae in Fornax and adding them to the SnIa calibration (correcting for light curve shape) gives Ho = 67 (+/-6)_r [+/-7]_s km/s/Mpc out to a distance in excess of 500 Mpc. All five of these Ho determinations agree to within their statistical errors. The resulting estimate of the Hubble Constant combining all these determinations is Ho = 72 (+/-5)_r [+/-12]_s km/s/Mpc., Comment: Accepted for publication in the Astrophysical Journal, Apr. 10 issue 28 pages, 3 tables, 12 figures (Correct figures and abstract)

Using WFPC2 on the Hubble Space Telescope, we have isolated a sample of 258 white dwarfs (WDs) in the Galactic globular cluster M4. Fields at three radial distances from the cluster center were observed and sizeable WD populations were found in all three. The location of these WDs in the color-magnitude diagram, their mean mass of 0.51($ \pm 0.03$)M$_{\odot}$, and their luminosity function confirm basic tenets of stellar evolution theory and support the results from current WD cooling theory. The WDs are used to extend the cluster main-sequence mass function upward to stars that have already completed their nuclear evolution. The WD/red dwarf binary frequency in M4 is investigated and found to be at most a few percent of all the main-sequence stars. The most ancient WDs found are about 9 Gyr old, a level which is set solely by the photometric limits of our data. Even though this is less than the age of M4, we discuss how these cooling WDs can eventually be used to check the turnoff ages of globular clusters and hence constrain the age of the Universe., Comment: 46 pages, latex, no figures included, figures available at ftp://ftp.astro.ubc.ca/pub/richer/wdfig.uu size 2.7Mb. To be published in the Astrophysical Journal