Search

We combine the deepest X-ray survey from the Chandra Deep Field-South (CDF-S) `7-Ms' survey with the deepest mid-infrared (5.6$ \mu m$) image from the JWST/MIRI Deep Imaging Survey (MIDIS) in the Hubble Ultra-Deep Field (HUDF) to study the infrared counterparts and point-source emission of 31 X-ray sources with a median, intrinsic, rest-frame X-ray luminosity of $\log_{10}(L_{\rm Xc}^{\rm 0.5-7keV})$=42.04$\pm$0.22 erg $\rm s^{-1}$. The sample includes 24 AGN with a redshift range, as set by the X-ray detectability, of $z \simeq 0.5-3$. Through a multi-wavelength morphological decomposition, employing three separate classifications (visual, parametric and non-parametric) we separate (where present) the luminosity of the point-like AGN component from the remainder of the host-galaxy emission. The unprecedented mid-infrared sensitivity and imaging resolution of MIRI allows, in many cases, the direct characterisation of point-like (i.e. unresolved) components in the galaxies' emission. We establish a broad agreement between the three morphological classifications. At least 70% of the X-ray sources, including some classified as galaxies, show unresolved emission in the MIRI images, with the unresolved-to-total flux fraction at rest-frame 2$\mu m$ ranging from $\sim$0.2 to $\sim$0.9. At high X-ray luminosities ($\log_{10}(L_{\rm Xc}$)>43 erg $\rm s^{-1}$) we derive a consistent rest-frame near-infrared 2$ \mu m$ point-source luminosity to that derived for local AGN, whilst at lower X-ray luminosity we identify an excess in the 2$ \mu m$ emission compared to pre-JWST studies. We speculate this offset may be driven by a combination of Compton-thick AGN components and nuclear starburst, merger driven activity. Our observations highlight the complex nature of X-ray sources in the distant Universe and demonstrate the power of JWST/MIRI in quantifying their nuclear infrared emission. (Abridged), Comment: 21 pages, 7 figures. Submitted to A&A

The recently launched James Webb Space Telescope (JWST) is opening new observing windows on the distant universe. Among JWST's instruments, the Mid Infrared Instrument (MIRI) offers the unique capability of imaging observations at wavelengths $\lambda > 5\mu$m. This enables unique access to the rest frame near infra-red (NIR, $\lambda \ge 1$\mum) emission from galaxies at redshifts $z>4$ and the visual ($\lambda \gtrsim 5000$\AA) rest frame for $z>9$. We here report on the guaranteed time observations (GTO) from the MIRI European Consortium, of the Hubble Ultra Deep Field (HUDF), forming the MIRI Deep Imaging Survey (MIDIS), consisting of an on source integration time of $\sim41$ hours in the MIRI/F560W (5.6 $\mu$m) filter. To our knowledge, this constitutes the longest single filter exposure obtained with JWST of an extragalactic field as yet., Comment: submitted to A&A on July 30, 2024

The circumgalactic medium (CGM) plays a vital role in galaxy evolution, however, studying the emission from CGM is challenging due to its low surface brightness and the complexities involved in interpreting resonant lines such as Ly$\alpha$. The near-infrared coverage, unprecedented sensitivity, and high spatial resolution of JWST enable us to study the optical strong lines associated with the extended Ly$\alpha$ "nebulae" at redshifts of 2--3. These lines serve as diagnostic tools to infer the physical conditions in the CGM gas reservoir of these systems. In deep medium-band images taken by the JWST, we serendipitously discovered the [O III] emission from the CGM around a massive interacting galaxy system at a redshift z~2.8, known to be embedded in a bright extended (100 kpc) Ly$\alpha$ "nebula." This is the first time that the [O III] lines have been detected from a Ly$\alpha$ "nebula." The JWST images reveal that the CGM gas actually resides in narrow (~ 2.5 kpc) filamentary structures with strong [O III] emission, tracing the same extent as the Ly$\alpha$ emission. An analysis of the [O III] suggests that the emitting CGM is fully ionized and is energetically dominated by mechanical heating. We also find that the density and pressure are higher than those commonly predicted by simulations of the CGM. We conclude that the observed CGM emission originates from the gas expelled by the episodic feedback processes, cooling down and enriching the CGM, while traveling a distance of at least 60 kpc. These observations demonstrate how intensive feedback processes shape gas distribution and properties in the CGM around massive halos. While access to such deep, high-resolution imaging opens up a new discovery space for investigating the CGM, it also challenges numerical simulations with respect to explaining and reproducing the exquisitely complex structures revealed by the observations., Comment: 15 pages, 7 figures, 1 table, accepted by A&A Letter

Making use of JWST NIRSpec and NIRCam data, we conduct a detailed analysis of the spectroscopic and photometric properties of GN-z8-LAE, a strong Lya emitter at z=8.279. Our goal is to investigate the interstellar medium (ISM) physical conditions that enable the Lya detection in this source at the Epoch of Reionization (EoR) and scrutinize GN-z8-LAE as an early reionizer. In broad agreement with previous results, we find that GN-z8-LAE is a young galaxy (age ~ 10 Myr) with a low stellar mass (M* ~ 10^7.66 Msun), significantly lower than those of most Lya emitters known at similarly high redshifts. The derived stellar mass and star formation rate surface densities are 355 Msun/pc^2 and 88 Msun/yr/kpc^2, respectively. Our spectral analysis indicates that: the Lya line peak has a small velocity offset 133+-72 km/s with respect to the galaxy systemic redshift; CIV] / CIII] ~ 3.3; the ISM is characterized by a hard ionization field, although no signature of AGN is present. Moreover, we report the presence of NIII]1750 emission with super-solar N abundance, which makes GN-z8-LAE one of the first known cases of a simultaneous strong Lya and nitrogen emitter at the EoR. Based on all these properties, we apply a wide range of methods to constrain the absolute Lyman continuum escape fraction of GN-z8-LAE, and find that it is >14% in all cases. Therefore, we conclude that GN-z8-LAE is a robust candidate for a Lyman continuum (LyC) leaker at the EoR which is being caught at the moment of efficiently reionizing its surrounding medium., Comment: Submitted to ApJ. 22 pages, 8 Figures

We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$\alpha$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observations reveal an unexpected extremely red color at rest-frame near-infrared wavelengths, $\rm F444W - F1000W = 2.33 \pm 0.06$. Such steep rise in the near-infrared, completely missed without MIRI imaging, is poorly reproduced by models including only stellar populations and hints towards the presence of an Active Galactic Nucleus (AGN). According to the shape of its overall SED, Virgil belongs to the recently discovered population of Little Red Dots (LRDs) but displays an extended rest-frame UV-optical wavelengths morphology following a 2D-S\'ersic profile with average index $n = 0.93^{+0.85}_{-0.31}$ and $r_e = 0.49^{+0.05}_{-0.11}$ pkpc. Only at MIRI wavelengths Virgil is unresolved due to the coarser PSF. We also estimate a bolometric luminosity $L_{\rm bol, AGN} = (8.9-11)\times 10^{44}\rm~erg~s^{-1}$ and a supermassive black hole mass $M_{\rm BH} = (7-9)\times 10^6\rm ~ M_\odot$ (if $\lambda_{\rm Edd} = 1$) in agreement with recently reported values for LRDs. This discovery demonstrates the crucial importance of deep MIRI surveys to find AGN amongst high-$z$ galaxies that otherwise would be completely missed and raises the question of how common Virgil-like objects could be in the early Universe., Comment: 19 pages, 11 figures, 3 tables. Submitted to ApJ

We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5\sigma$ detection limits in individual bands, as well as in the MIDIS F560W ultra-deep data ($\sim$29 mag, $5\sigma$). Analyzing the spectral energy distribution built with low-$S/N$ ($<5$) measurements in individual optical-to-mid-infrared filters and higher $S/N$ ($\gtrsim5$) in stacked NIRCam data, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a Solar System body based on the $<0.016$" proper motion in the 1-year apart JADES and MIDIS observations. A sub-stellar Galactic nature is deemed unlikely, given that the Cerberus' relatively flat NIRCam-to-NIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers 3 possibilities: (1) A $z\sim0.4$ galaxy with strong emission from polycyclic aromatic hydrocarbons; the very low inferred stellar mass, $\mathrm{M}_\star=10^{5-6}$ M$_\odot$, makes this possibility highly improbable. (2) A dusty galaxy at $z\sim4$ with an inferred stellar mass $\mathrm{M}_\star\sim10^{8}$ M$_\odot$. (3) A galaxy with observational properties similar to those of the reddest little red dots discovered around $z\sim7$, but Cerberus lying at $z\sim15$, with the rest-frame optical dominated by emission from a dusty torus or a dusty starburst., Comment: Accepted for publication in ApJL, results unchanged, manuscript shortened to publish as a letter

Observations with the James Webb Space Telescope (JWST) have uncovered numerous faint active galactic nuclei (AGN) at $z\sim5$ and beyond. These objects are key to our understanding of the formation of supermassive black holes (SMBHs), their co-evolution with host galaxies, as well as the role of AGN in cosmic reionization. Using photometric colors and size measurements, we perform a search for compact red objects in an array of blank deep JWST/NIRCam fields totaling $\sim640$ arcmin$^{2}$. Our careful selection yields 260 reddened AGN candidates at $4

We investigate the properties of strong (Hb+[OIII]) emitters before and after the end of the Epoch of Reionization from z=8 to z=5.5. We make use of ultra-deep JWST/NIRCam imaging in the Parallel Field of the MIRI Deep Imaging Survey (MIDIS) in the Hubble eXtreme Deep Field (P2-XDF), in order to select prominent (Hb+[OIII]) emitters (with rest EW_0 > 100 Angstroms) at z=5.5-7, based on their flux density enhancement in the F356W band with respect to the spectral energy distribution continuum. We complement our selection with other (Hb+[OIII]) emitters from the literature at similar and higher (z=7-8) redshifts. We find (non-independent) anti-correlations between EW_0(Hb+[OIII]) and both galaxy stellar mass and age, in agreement with previous studies, and a positive correlation with specific star formation rate (sSFR). On the SFR-M* plane, the (Hb+[OIII]) emitters populate both the star-formation main sequence and the starburst region, which become indistinguishable at low stellar masses (log10(M*) < 7.5). We find tentative evidence for a non-monotonic relation between EW_0(Hb+[OIII]) and SFR, such that both parameters correlate with each other at SFR > 1 Msun/yr, while the correlation flattens out at lower SFRs. This suggests that low metallicities producing high EW_0(Hb+[OIII]) could be important at low SFR values. Interestingly, the properties of the strong emitters and other galaxies (33% and 67% of our z=5.5-7 sample, respectively) are similar, including, in many cases, high sSFR. Therefore, it is crucial to consider both emitters and non-emitters to obtain a complete picture of the cosmic star formation activity around the Epoch of Reionization., Comment: 14 pages, 8 figures. Accepted for publication at the ApJ. Version v3 includes minor update to match journal version

We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$\alpha$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MIRI playing a crucial role in robustly assessing the LAE's stellar mass and ages. Our LAEs are low-mass objects (log$_{10}$(M$_\star$[M$_\odot$]) ~ 7.5), with little or no dust extinction (E(B - V) ~ 0.1) and a blue UV continuum slope ($\beta$ ~ -2.2). While 75% of our LAEs are young (< 100 Myr), the remaining 25% have significantly older stellar populations (> 100 Myr). These old LAEs are statistically more massive, less extinct and have lower specific star formation rate (sSFR) compared to young LAEs. Besides, they populate the M$_\star$ - SFR plane along the main-sequence (MS) of star-forming galaxies, while young LAEs populate the starburst region. The comparison between the LAEs properties to those of a stellar-mass matched sample of LBGs shows no statistical difference between these objects, except for the LBGs redder UV continuum slope and marginally larger E(B - V) values. Interestingly, 48% of the LBGs have ages < 10 Myr and are classified as starbursts, but lack detectable Ly$\alpha$ emission. This is likely due to HI resonant scattering and/or selective dust extinction. Overall, we find that JWST observations are crucial in determining the properties of LAEs and shedding light on the properties and similarities between LAEs and LBGs., Comment: 22 pages, 15 figures, 5 tables. Accepted for publication in APJ

We present deep James Webb Space Telescope (JWST)/MIRI F560W observations of a flux-limited, ALMA-selected sample of 28 galaxies at z=0.5-3.6 in the Hubble Ultra Deep Field (HUDF). The data from the MIRI Deep Imaging Survey (MIDIS) reveal the stellar structure of the HUDF galaxies at rest-wavelengths of >1 micron for the first time. We revise the stellar mass estimates using new JWST photometry and find good agreement with pre-JWST analysis; the few discrepancies can be explained by blending issues in the earlier lower-resolution Spitzer data. At z~2.5, the resolved rest-frame near-infrared (1.6 micron) structure of the galaxies is significantly more smooth and centrally concentrated than seen by HST at rest-frame 450 nm (F160W), with effective radii of Re(F560W)=1-5 kpc and S\'ersic indices mostly close to an exponential (disk-like) profile (n~1), up to n~5 (excluding AGN). We find an average size ratio of Re(F560W)/Re(F160W)~0.7 that decreases with stellar mass. The stellar structure of the ALMA-selected galaxies is indistinguishable from a HUDF reference sample of galaxies with comparable MIRI flux density. We supplement our analysis with custom-made, position-dependent, empirical PSF models for the F560W observations. The results imply that a smoother stellar structure is in place in massive gas-rich, star-forming galaxies at Cosmic Noon, despite a more clumpy rest-frame optical appearance, placing additional constraints on galaxy formation simulations. As a next step, matched-resolution, resolved ALMA observations will be crucial to further link the mass- and light-weighted galaxy structures to the dusty interstellar medium., Comment: 19 pages, 12 figures, 1 table, Accepted for publication in ApJ

The Cosmic Evolution Early Release Science Survey (CEERS), targeting the Extended Groth Strip extragalactic field, is one of the JWST Director's Discretionary Early Release Science programs. To date, all observations have been executed and include NIRCam/MIRI imaging and NIRSpec/NIRCam spectroscopic exposures. Here, we discuss the MIRI imaging, which includes eight pointings, four of which provide deep imaging with the bluer bands (F560W, F770W) and four with contiguous wavelength coverage in F1000W, F1280W, F1500W, and F1800W, where two of these also include coverage in F770W and F2100W. We present a summary of the data, the data quality, and data reduction. The data reduction is based on the JWST Calibration Pipeline combined with custom modifications and additional steps designed to enhance the output quality, including improvements in astrometry and the removal of detector artifacts. We estimate the image depth of the reduced mosaics, and show that these generally agree with expectations from the Exposure Time Calculator. We compare the MIRI F560W and F770W flux densities for bright sources to measurements from Spitzer/IRAC Ch3 (5.8 $\mu$m) and Ch4 (8.0 $\mu$m), and we find that they agree with systematic differences of $<0.1$ mag. For the redder MIRI bands, we assess their quality by studying the spectral energy distributions (SEDs) of Galactic stars. The SEDs are consistent with the expected Rayleigh-Jeans law with a deviation $\sim 0.03$ mag, indicating that the MIRI colors are reliable. We also discuss all publicly released data products (images and source catalogs), which are available on the CEERS website (https://ceers.github.io/)., Comment: 17 pages, 11 figures, and 4 tables. ApJL in press

We present a JWST/MRS spectrum of the quasar J1120+0641 at z=7.0848, the first spectroscopic observation of a reionisation-era quasar in the rest-frame infrared ($0.6<\lambda<3.4\mu$m). In the context of the mysterious fast assembly of the first supermassive black holes at z>7, our observations enable for the first time the detection of hot torus dust, the H$\alpha$ emission line, and the Paschen-series broad emission lines in a quasar at z>7. Hot torus dust is clearly detected as an upturn in the continuum emission at $\lambda_{\text{rest}}\simeq1.3\mu$m, leading to a black-body temperature of $T=1413.5^{+5.7}_{-7.4}$K. Compared to similarly-luminous quasars at 06 ($2\sigma$ significance). We measure the black hole mass of J1120+0641 based on the H$\alpha$ Balmer line, $M_{\text{BH}}=1.52\pm0.17\cdot 10^9 M_\odot$, which is in good agreement with the previous rest-UV MgII black hole mass measurement. The black hole mass based on the Paschen-series lines is also consistent, indicating no significant extinction in the rest-frame UV measurement. The broad H$\alpha$, Pa-$\alpha$ and Pa-$\beta$ emission lines are consistent with an origin in a common broad-line region (BLR) with density log$N_H/\text{cm}^{-3}\geq 12$, ionisation parameter $-7<$log$U<-4$, and extinction E(B-V)$\lesssim 0.1$mag. These BLR parameters are consistent with similarly-bright quasars at 0

We analyze a sample of 3300 galaxies between redshifts z~3.5 and z~8.5 selected from JWST images in the Hubble Ultra Deep Field (HUDF) and UKIDSS Ultra Deep Survey field, including objects with stellar masses as low as ~ 10^8 Msun up to z~8. The depth and wavelength coverage of the JWST data allow us, for the first time, to derive robust stellar masses for such high-z, low stellar-mass galaxies on an individual basis. We compute the galaxy stellar mass function (GSMF), after complementing our sample with ancillary data from CANDELS to constrain the GMSF at high stellar masses (M > M*). Our results show a steepening of the low stellar-mass end slope (a) with redshift, with a = -1.61 (+/-0.05) at z~4 and a = -1.98 (+/-0.14) at z~7. We also observe an evolution of the normalization phi* from z~7 to z~4, with phi*(z~4)/phi*(z~7)= 130 (+210/-50). Our study incorporates a novel method for the estimation of the Eddington bias that takes into account its possible dependence both on stellar mass and redshift, while allowing for skewness in the error distribution. We finally compute the resulting cosmic stellar mass density and find a flatter evolution with redshift than previous studies., Comment: 21 pages, including 11 figures and 6 tables. The manuscript has been accepted for publication in the ApJ

We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $88$ galaxy candidates based on their dropout nature in the F115W and/or F150W filters, a high probability for their photometric redshifts, estimated with three different codes, being at $z>8$, good fits based on $\chi^2$ calculations, and predominant solutions compared to $z<8$ alternatives. We find mild evolution in the luminosity function from $z\sim13$ to $z\sim8$, i.e., only a small increase in the average number density of $\sim$0.2 dex, while the faint-end slope and absolute magnitude of the knee remain approximately constant, with values $\alpha=-2.2\pm0.1$ and $M^*=-20.8\pm0.2$ mag. Comparing our results with the predictions of state-of-the-art galaxy evolution models, we find two main results: (1) a slower increase with time in the cosmic star formation rate density compared to a steeper rise predicted by models; (2) nearly a factor of 10 higher star formation activity concentrated in scales around 2 kpc in galaxies with stellar masses $\sim10^8$ M$_\odot$ during the first 350 Myr of the Universe, $z\sim12$, with models matching better the luminosity density observational estimations $\sim$150 Myr later, by $z\sim9$., Comment: Accepted for publication in ApJL, in memoriam Hans Ulrik Noorgaard-Nielsen and Olivier Le F\`evre

Using the novel $JWST$/NIRCam observations in the Abell 2744 field, we present a first spatially resolved overview of an $HST$-dark galaxy, spectroscopically confirmed at $z=2.58$ with magnification $\mu\approx1.9$. While being largely invisible at $\sim$1 $\mu$m with NIRCam, except for sparse clumpy sub-structures, the object is well-detected and resolved in the long-wavelength bands with a spiral shape clearly visible in F277W. By combining ancillary ALMA and $Herschel$ data, we infer that this object is an edge-on dusty spiral with an intrinsic stellar mass log$(M_*/M_\odot)\sim11.3$ and a dust-obscured SFR$\sim 300~M_\odot$~yr$^{-1}$. A massive quiescent galaxy (log$(M_*/M_\odot)\sim10.8$) with tidal features lies 2\farcs{0} away ($r$$\sim$9 kpc), at a consistent redshift as inferred by $JWST$ photometry, indicating a potential major merger. The dusty spiral lies on the main-sequence of star formation, and shows high dust attenuation in the optical ($34$ would be only detectable in F356W and F444W in UNCOVER-like survey, and become totally $JWST$-dark at $z\sim6$. This suggests that detecting highly attenuated galaxies in the Epoch of Reionization might be a challenging task for $JWST$., Comment: 15 pages, 5 figures, 1 table. Accepted to ApJL

Until now, our knowledge of the extragalactic Universe at mid-IR wavelengths (> 5 microns) was limited to rare active galactic nuclei (AGN) and the brightest normal galaxies up to z~3. The advent of the JWST with its Mid-Infrared Instrument (MIRI) will revolutionise the ability of the mid-IR regime as a key wavelength domain to probe the high-z Universe. In this work we present a first study of JWST MIRI 7.7 micron sources selected with > 3 sigma significance from the lensing cluster field SMACS J0723.3-7327. We model their spectral energy distribution fitting with 13 JWST and HST broad bands, in order to obtain photometric redshifts and derived physical parameters for all these sources. We find that this 7.7 micron galaxy sample is mainly composed of normal galaxies up to z=4 and has a tail of about 2% of sources at higher redshifts to z~ 9-10. The vast majority of our galaxies have [3.6]-[7.7]<0 colours and very few of them need high dust extinction values (A_V=3 -6 mag) for their SED fitting. The resulting lensing-corrected stellar masses span the range 10^7-10^{11} Msun. Overall, our results clearly show that the first MIRI 7.7 micron observations of deep fields are already useful to probe the high-redshift Universe and suggest that the deeper 7.7 micron observations to be available very soon will open up, for the first time, the epoch of reionisation at mid-IR wavelengths., Comment: 13 pages, 9 figures, 2 tables. Accepted for publication at the ApJ. In the current version we make the 7.7 micron source catalogue publicly available

A wealth of extragalactic populations completely missed at UV-optical wavelengths has been identified in the last decade, combining the deepest HST and Spitzer observations. These dark sources are thought to be very dusty and star-forming systems at 3

We study the relation between stellar mass (M*) and star formation rate (SFR) for star-forming galaxies over approximately five decades in stellar mass (5.5 <~ log10(M*/Msun) <~ 10.5) at z ~ 3-6.5. This unprecedented coverage has been possible thanks to the joint analysis of blank non-lensed fields (COSMOS/SMUVS) and cluster lensing fields (Hubble Frontier Fields) which allow us to reach very low stellar masses. Previous works have revealed the existence of a clear bimodality in the SFR-M* plane with a star-formation Main Sequence and a starburst cloud at z ~ 4-5. Here we show that this bimodality extends to all star-forming galaxies and is valid in the whole redshift range z ~ 3-6.5. We find that starbursts constitute at least 20% of all star-forming galaxies with M* >~ 10^9 Msun at these redshifts and reach a peak of 40% at z=4-5. More importantly, 60% to 90% of the total SFR budget at these redshifts is contained in starburst galaxies, indicating that the starburst mode of star-formation is dominant at high redshifts. Almost all the low stellar-mass starbursts with log10(M*/Msun) <~ 8.5 have ages comparable to the typical timescales of a starburst event, suggesting that these galaxies are being caught in the process of formation. Interestingly, galaxy formation models fail to predict the starburst/main-sequence bimodality and starbursts overall, suggesting that the starburst phenomenon may be driven by physical processes occurring at smaller scales than those probed by these models., Comment: 24 pages, including 15 figures (17 files in total) and 4 tables. The manuscript has been accepted for publication in the ApJ