Search

The Euclid, Rubin/LSST and Roman (WFIRST) projects will undertake flagship optical/near-infrared surveys in the next decade. By mapping thousands of square degrees of sky and covering the electromagnetic spectrum between 0.3 and 2 microns with sub-arcsec resolution, these projects will detect several tens of billions of sources, enable a wide range of astrophysical investigations by the astronomical community and provide unprecedented constraints on the nature of dark energy and dark matter. The ultimate cosmological, astrophysical and time-domain science yield from these missions will require joint survey processing (JSP) functionality at the pixel level that is outside the scope of the individual survey projects. The JSP effort scoped here serves two high-level objectives: 1) provide precise concordance multi-wavelength images and catalogs over the entire sky area where these surveys overlap, which accounts for source confusion and mismatched isophotes, and 2) provide a science platform to analyze concordance images and catalogs to enable a wide range of astrophysical science goals to be formulated and addressed by the research community. For the cost of about 200WY, JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also allow for the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

SPHEREx, the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and ices Explorer, is a NASA MIDEX mission planned for launch in 2024. SPHEREx will carry out the first all-sky spectral survey at wavelengths between 0.75µm and 5µm with spectral resolving power ~40 between 0.75 and 3.8µm and ~120 between 3.8 and 5µm At the end of its two-year mission, SPHEREx will provide 0.75-to-5µm spectra of each 6.”2 x 6.”2 pixel on the sky - 14 billion spectra in all. This paper updates an earlier description of SPHEREx presenting changes made during the mission's Preliminary Design Phase, including a discussion of instrument integration and test ow and a summary of the data processing, analysis, and distribution plans.

The infrared (IR) spectral energy distributions (SEDs) of main-sequence galaxies in the early Universe (z > 4) is currently unconstrained as IR continuum observations are time-consuming and not feasible for large samples. We present Atacama Large Millimetre Array Band 8 observations of four main-sequence galaxies at z similar to 5.5 to study their IR SED shape in detail. Our continuum data (rest-frame 110 mu m, close to the peak of IR emission) allows us to constrain luminosity-weighted dust temperatures and total IR luminosities. With data at longer wavelengths, we measure for the first time the emissivity index at these redshifts to provide more robust estimates of molecular gas masses based on dust continuum. The Band 8 observations of three out of four galaxies can only be reconciled with optically thin emission redward of rest-frame 100 mu m. The derived dust peak temperatures at z similar to 5.5 (30-43 K) are elevated compared to average local galaxies, however, similar to 10K below what would be predicted from an extrapolation of the trend at z <4. This behaviour can be explained by decreasing dust abundance (or density) towards high redshifts, which would cause the IR SED at the peak to be more optically thin, making hot dust more visible to the external observer. From the 850-mu m dust continuum, we derive molecular gas masses between 10(10) and 10(11)M(circle dot) and gas fractions (gas over total mass) of 30-80 per cent (gas depletion times of 100-220 Myr). All in all, our results provide a first measured benchmark SED to interpret future millimetre observations of normal, main-sequence galaxies in the early Universe.

We present the physical extent of [C II] 158 mu m line-emitting gas from 46 star-forming galaxies at z.=.4-6 from the ALMA Large Program to INvestigate C II at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [C II] line (r(e),([C II])) for individual galaxies and compare them with the rest-frame ultraviolet (UV) continuum (re, UV) from Hubble Space Telescope images. The effective radius r(e),([C II]) exceeds r(e, UV) by factors of similar to 2-3, and the ratio of r(e),([C II])/r(e, UV) increases as a function of M-star. We do not find strong evidence that the [C II] line, rest-frame UV, and far-infrared (FIR) continuum are always displaced over similar or equal to 1 kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [C II] component over 10 kpc scales detected at 4.1s-10.9s beyond the size of rest-frame UV and FIR continuum. One object has tentative rotating features up to similar to 10 kpc, where the 3D model fit shows the rotating [C II]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [C II] line structure have high star formation rate, high stellar mass (M-star), low Ly alpha equivalent width, and more blueshifted (redshifted) rest-frame UV metal absorption (Ly alpha line), as compared to galaxies without such extended [C II] structures. Although we cannot rule out the possibility that a selection bias toward luminous objects may be responsible for such trends, the star-formation-driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [C II] line structures are ubiquitous to high-z star-forming galaxies.

We present a detailed analysis of a large sample of spectroscopically confirmed massive quiescent galaxies (MQGs; log(M */M ⊙) ∼ 11.5) at z ≳ 2. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with Very Large Telescope (VLT) X-shooter spectroscopy and Hubble Space Telescope (HST)/WFC3 H F160W imaging. These observations allow us to unambiguously confirm their redshifts, ascertain their quiescent nature and stellar ages, and reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high-mass end of the stellar mass–size relation for quiescent galaxies at z = 2. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey, showing that evolution in these galaxies is dominated by changes in size. The HST images reveal, as surprisingly high, that 40% of the sample has tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from z = 2 to 0, coupled with a doubling of the stellar mass, with a factor of 4 size increase and the observed disturbed stellar morphologies, supports dry minor mergers as the primary drivers of the evolution of the MQGs over the last 10 billion yr.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope (HST) Treasury program taking data from 2018 to 2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in Wide Field Camera 3/IR F105W, F125W, and F160W and Advanced Camera for Surveys/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multiwavelength data sets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies that would otherwise lie below HST detection limits and model foreground clusters to study the properties of dark matter and galaxy assembly. The expanded area will provide the first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient follow-up. BUFFALO data including mosaics, value-added catalogs, and cluster mass distribution models will be released via MAST on a regular basis as the observations and analysis are completed for the six individual clusters.

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves., Comment: 20 pages (+appendix), 11 figures, accepted for publication in the Astrophysical Journal (ApJ). This work comes with a companion letter by Tanaka, Valentino, Toft et al. 2019, ApJL, 885, L34

We calibrated a technique to measure dust attenuation in star-forming galaxies. The technique utilizes the stellar-wind lines in Wolf-Rayet stars, which are widely observed in galaxy spectra. The He II 1640 and 4686 features are recombination lines whose ratio is largely determined by atomic physics. Therefore they can serve as a stellar dust probe in the same way as the Balmer lines are used as a nebular probe. We measured the strength of the He II 1640 line in 97 Wolf-Rayet stars in the Galaxy and the Large Magellanic Cloud. The reddening corrected fluxes follow a tight correlation with a fixed ratio of 7.76 for the He II 1640 to 4686 line ratio. Dust attenuation decreases this ratio. We provide a relation between the stellar E(B-V) and the observed line ratio for several attenuation laws. Combining this technique with the use of the nebular Balmer decrement allows the determination of the stellar and nebular dust attenuation in galaxies and can probe its effects at different stellar age and mass regimes, independently of the initial mass function and the star-formation history. We derived the dust reddening from the He II line fluxes and compared it to the reddening from the Balmer decrement and from the slope of the ultraviolet continuum in two star-forming galaxies. The three methods result in dust attenuations which agree to within the errors. Future application of this technique permits studies of the stellar dust attenuation compared to the nebular attenuation in a representative galaxy sample., Comment: Accepted by The Astronomical Journal

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers., Comment: Published in the Astrophysical Journal letters. Fixed an error in dynamical mass

Machine-learning (ML) algorithms will play a crucial role in studying the large datasets delivered by new facilities over the next decade and beyond. Here, we investigate the capabilities and limits of such methods in finding galaxies with brightness-variable active galactic nuclei (AGN). Specifically, we focus on an unsupervised method based on self-organizing maps (SOM) that we apply to a set of nonparametric variability estimators. This technique allows us to maintain domain knowledge and systematics control while using all the advantages of ML. Using simulated light curves that match the noise properties of observations, we verify the potential of this algorithm in identifying variable light curves. We then apply our method to a sample of ~8300 WISE color-selected AGN candidates in Stripe 82, in which we have identified variable light curves by visual inspection. We find that with ML we can identify these variable classified AGN with a purity of 86% and a completeness of 66%, a performance that is comparable to that of more commonly used supervised deep-learning neural networks. The advantage of the SOM framework is that it enables not only a robust identification of variable light curves in a given dataset, but it is also a tool to investigate correlations between physical parameters in multi-dimensional space - such as the link between AGN variability and the properties of their host galaxies. Finally, we note that our method can be applied to any time-sampled light curve (e.g., supernovae, exoplanets, pulsars, and other transient events)., Comment: 9 pages, 4 figures, accepted for publication in ApJL

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1

We perform a comprehensive analysis of the redshift evolution of the rest-UV spectra of star-forming galaxies out to z~5. We combine new z~5 measurements of HI Ly$\alpha$ and low- and high-ionization interstellar metal absorption features with comparable measurements at z~2-4. We measure the equivalent widths of interstellar absorption features using stacked spectra in bins of Ly$\alpha$ equivalent width, performing corrections to Ly$\alpha$ strengths based on a model for the transmission of the intergalactic medium. We find a strong correlation between decreasing low-ionization absorption strength and increasing Ly$\alpha$ emission strength over the redshift range z~2-5, suggesting that both of these quantities are fundamentally linked to neutral gas covering fraction. At the highest Ly$\alpha$ equivalent widths, we observe evolution at $z\sim5$ towards greater Ly$\alpha$ emission strength at fixed low-ionization absorption strength. If we interpret the non-evolving relationship of Ly$\alpha$ emission strength and low-ionization line strength at z~2-4 as primarily reflecting the radiative transfer of Ly$\alpha$ photons, this evolution at z~5 suggests a higher intrinsic production rate of Ly$\alpha$ photons than at lower redshift. Our conclusion is supported by the joint evolution of the relationships among Ly$\alpha$ emission strength, interstellar absorption strength, and dust reddening. We perform additional analysis in bins of stellar mass, star-formation rate, UV luminosity, and age, examining how the relationships between galaxy properties and Ly$\alpha$ emission evolve towards higher redshift. We conclude that increasing intrinsic Ly$\alpha$ photon production and strong detection of nebular CIV emission (signaling lower metallicity) at z~5 indicate an elevated ionized photon production efficiency ($\xi_{\rm ion}$)., Comment: 20 pages, 17 figures, accepted for publication in MNRAS

We report the discovery of a mid-infrared variable AGN which is hosted by an ultraluminous infrared galaxy (ULIRG) candidate in the Sloan Stripe 82 field. \textit{WISE} \textit{J030654.88+010833.6} is a red, extended galaxy, which we estimate to be at a photometric redshift of 0.28 $\leq$ z $\leq$ 0.31, based on its optical and near-infrared spectral energy distribution (SED). The factor of two variability over 8 years seen in the \textit{WISE} 3.4 and 4.6 $\mu$m wavelength channels is not clearly correlated with optical variability in archival data. Based on our estimation of the physical parameters of the host galaxy, \textit{J030654.88+010833.6} is possibly a composite AGN/starburst ULIRG in a phase where high star formation $\sim$ 70 M$_{\odot}$ year$^{-1}$ is occurring. Our estimate of the black hole mass to stellar mass ratio also appears to be consistent with that of broad-line AGN in the local universe. The long-term variability of \textit{J030654.88+010833.6} as seen in the \textit{WISE} \textit{W1} and \textit{W2} light curves is likely due to variations in the accretion rate, with the energy being reprocessed by dust in the vicinity of the AGN., Comment: Accepted for publication in the Astrophysical Journal, 11 pages

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$., Comment: 17 pages, 11 figures, Accepted for Publication in ApJ

Thanks to deep optical to near-IR imaging and spectroscopy, significant progress is made in characterizing the rest-frame UV to optical properties of galaxies in the early universe (z > 4). Surveys with Hubble, Spitzer, and ground-based facilities (Keck, Subaru, and VLT) provide spectroscopic and photometric redshifts, measurements of the spatial structure, stellar masses, and optical emission lines for large samples of galaxies. Recently, the Atacama Large (Sub) Millimeter Array (ALMA) has become a major player in pushing studies of high redshift galaxies to far-infrared wavelengths, hence making panchromatic surveys over many orders of frequencies possible. While past studies focused mostly on bright sub-millimeter galaxies, the sensitivity of ALMA now enables surveys like ALPINE, which focuses on measuring the gas and dust properties of a large sample of normal main-sequence galaxies at z > 4. Combining observations across different wavelengths into a single, panchromatic picture of galaxy formation and evolution is currently and in the future an important focus of the astronomical community., Comment: 4 pages, 2 figures. Submitted to Proceedings IAU Symposium No. 341, 2018

Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

The redshift range z = 4-6 marks a transition phase between primordial and mature galaxy formation in which galaxies considerably increase their stellar mass, metallicity, and dust content. The study of galaxies in this redshift range is therefore important to understanding early galaxy formation and the fate of galaxies at later times. Here, we investigate the burstiness of the recent star formation history (SFH) of 221z similar to 4.5 main-sequence galaxies at log(M/M-circle dot) > 9.7 by comparing their ultra-violet (UV) continuum, H alpha luminosity, and H alpha equivalent-width (EW). The H alpha properties are derived from the Spitzer [3.6 mu m]-[4.5 mu m] broadband color, thereby properly taking into account model and photometric uncertainties. We find a significant scatter between H alpha- and UV-derived luminosities and star formation rates (SFRs). About half of the galaxies show a significant excess in Ha.compared to expectations from a constant smooth SFH. We also find a tentative anticorrelation between Ha.EW and stellar mass, ranging from 1000 angstrom at log(M/M-circle dot) < 10 to below 100 angstrom at log(M/M-circle dot) > 11. Consulting models suggests that most z similar to 4.5 galaxies had a burst of star formation within the last 50 Myr, increasing their SFRs by a factor of >5. The most massive galaxies on the other hand might decrease their SFRs and may be transitioning to a quiescent stage by z = 4. We identify differential dust attenuation (f) between stars and nebular regions as the main contributor to the uncertainty. With local galaxies selected by increasing H alpha EW (reaching values similar to high-z galaxies), we predict that f approaches unity at z > 4, consistent with the extrapolation of measurements out to z = 2.