Search

Mass estimates of black holes (BHs) in the centers of Active Galactic Nuclei (AGN) often rely on the radius-luminosity relation. However, this relation, usually probed by reverberation mapping (RM), is poorly constrained in the high-luminosity and high-redshift ends due to the very long expected lag times. Multiply imaged AGN may thus offer a unique opportunity to explore the radius-luminosity relation at these ends. In addition to comprising several magnified images which enable a more efficient light-curve sampling, the time delay between multiple images of strongly lensed quasars can also aid in making such RM measurements feasible on reasonable timescales: If the time delay is, for example, of the order of the expected time lag, changes in the emission lines in the leading image can be observed around the same time as the changes in the continuum in the trailing image. In this work we probe the typical time-delay distribution in galaxy-cluster lenses and estimate the number of both high-mass ($\sim10^9-10^{10}$ M$_{\odot}$), and high-redshift ($z\gtrsim4-12$) quasars that are expected to be strongly lensed by clusters. We find that less than one, very massive (and luminous, L$_{UV}>10^{46.5}$ erg s$^{-1}$) multiply imaged quasar is expected across the sky down to 30 AB mag. Nonetheless, up to several dozen thousand M$_{BH}\sim10^{6}$-$10^{8}$ M$_{\odot}$ broad-line AGN at $z>4$ should be multiply imaged by galaxy clusters and detectable with JWST, hundreds with $\textit{Euclid}$ and several thousands with the $\textit{Roman}$ Space Telescope, across the whole sky. These could supply an important calibration for the BH mass scaling in the early Universe., Comment: To be submitted. Comments welcome. 14 pages, 3 figures, 1 table

We present new free-form and hybrid mass reconstructions of the galaxy cluster lens MACS J0416.1-2403 at $z=0.396$ using the lens inversion method GRALE. The reconstructions use 237 spectroscopically confirmed multiple images from Bergamini et. al. 2023 as the main input. Our primary model reconstructs images to a positional accuracy of 0.191", thus representing one of the most precise reconstructions of this lens to date. Our models find broad agreement with previous reconstructions, and identify two $\sim 10^{11} M_{\odot}$ light-unaffiliated substructures. We focus on two highly magnified arcs: Spock and Mothra. Our model features a unique critical curve structure around the Spock arc with 2 crossings. This structure enables sufficient magnification across this arc to potentially explain the large number of transients as microlensing events of supergiant stars. Additionally, we develop a model of the millilens substructure expected to be magnifying Mothra, which may be a binary pair of supergiants with $\mu \sim 6000$. This model accounts for flexibility in the millilens position while preserving the observed flux and minimizing image position displacements along the Mothra arc. We constrain the millilens mass and core radius to $\lesssim 10^6 M_{\odot}$ and $\lesssim 12$ pc, respectively, which would render it one of the smallest and most compact substructures constrained by lensing. If the millilens is dominated by wave dark matter, the axion mass is constrained to be $< 3.0 \times 10^{-21}$ eV. Further monitoring of this lens with JWST will uncover more transients, permitting tighter constraints on the structure surrounding these two arcs., Comment: Submitted. Comments Welcome

Strong gravitational magnification by massive galaxy clusters enable us to detect faint background sources, resolve their detailed internal structures, and in the most extreme cases identify and study individual stars in distant galaxies. Highly magnified individual stars allow for a wide range of applications, including studies of stellar populations in distant galaxies and constraining small-scale dark matter structures. However, these applications have been hampered by the small number of events observed, as typically one or a few stars are identified from each distant galaxy. Here, we report the discovery of 46 significant microlensed stars in a single strongly-lensed high-redshift galaxy behind the Abell 370 cluster at redshift of 0.725 when the Universe was half of its current age (dubbed the ``Dragon arc''), based on two observations separated by one year with the James Webb Space Telescope ({\it JWST}). These events are mostly found near the expected lensing critical curves, suggesting that these are magnified individual stars that appear as transients from intracluster stellar microlenses. Through multi-wavelength photometry and colors, we constrain stellar types and find that many of them are consistent with red giants/supergiants magnified by factors of thousands. This finding reveals an unprecedented high occurrence of microlensing events in the Dragon arc, and proves that {\it JWST}'s time-domain observations open up the possibility of conducting statistical studies of high-redshift stars and subgalactic scale perturbations in the lensing dark matter field., Comment: 15 pages, 4 figures, 1 table submitted to Nature Astronomy

The discovery of extremely luminous galaxies at ultra-high redshifts ($z\gtrsim 8$) has posed a challenge for galaxy formation models. Most statistical analyses of this tension to date have not properly accounted for the variance due to field-to-field clustering, which causes the number counts of galaxies to vary from field to field, greatly in excess of Poisson noise. This super-Poissonian variance is often referred to as cosmic variance. Since cosmic variance increases rapidly as a function of mass, redshift, and for small observing areas, the most massive objects in deep \textit{JWST} surveys are severely impacted by cosmic variance. In this paper, we introduce a simple model to predict the distribution of the mass of the most massive galaxy found for different survey designs, which includes cosmic variance. The distributions differ significantly from previous predictions using the Extreme Value Statistics formalism, changing both the position and shape of the distribution of most massive galaxies in a counter-intuitive way. We test our model using the \texttt{UniverseMachine} simulations, where the predicted effects of including cosmic variance are clearly identifiable. Moreover, we find that the highly significant skew in the distributions of galaxy number counts for typical deep \textit{JWST} surveys lead to a high "variance on the variance", which greatly impacts the calculation of the cosmic variance itself. We conclude that it is crucial to accurately account for the impact of cosmic variance in any future analysis of tension between extreme galaxies in the early universe and galaxy formation models., Comment: Submitted to ApJ. The main text contains 15 pages and 9 figures. The appendix contains 4 pages and 3 figures

A new model is proposed in which typical galaxies form most of their stellar mass in a phase with an intrinsically red stellar population. In the standard picture, galaxies with intrinsically red stellar populations are believed to have old stellar populations, so that only galaxies with blue stellar populations have significant star formation, and subsequent changes to the stellar population come from predominantly from aging and merging populations which have already formed. However, several observational puzzles have developed which are difficult to reconcile with this standard scenario. The most massive blue star-forming galaxies, presumed to be at the end of their stellar mass growth, are $\sim 1$ dex less massive, have a $\sim 1$ dex lower $M_*/M_{BH}$ ratio, and have a bottom-lighter IMF than local quiescent galaxies. Here, a new solution is proposed: at low temperature and high metallicity, galaxies can continue to form stars efficiently without being able to form O and B stars. These red star-forming galaxies would have many of the same properties of the population currently described as post-starburst galaxies, allowing a new interpretation of their origin. Finally, additional falsifiable observational predictions of this model are also discussed., Comment: 10 pages, 2 figures

To study the role of environment in galaxy evolution, we reconstruct the underlying density field of galaxies based on COSMOS2020 (The Farmer catalog) and provide the density catalog for a magnitude limited ($K_{s}<24.5$) sample of $\sim 210 \, k$ galaxies at $0.42$., Comment: 22 pages, 11 figures, 2 tables, Submitted to ApJ

Several large JWST blank field observing programs have not yet discovered the first galaxies expected to form at $15 \leq z \leq 20$. This has motivated the search for more effective survey strategies that will be able to effectively probe this redshift range. Here, we explore the use of gravitationally lensed cluster fields, that have historically been the most effective discovery tool with HST. In this paper, we analyze the effectiveness of the most massive galaxy clusters that provide the highest median magnification factor within a single JWST NIRCam module in uncovering this population. The results of exploiting these lensing clusters to break the $z > 15$ barrier are compared against the results from large area, blank field surveys such as JADES and CEERS in order to determine the most effective survey strategy for JWST. We report that the fields containing massive foreground galaxy clusters specifically chosen to occupy the largest fraction of a single NIRCam module with high magnification factors in the source plane, whilst containing all multiple images in the image plane within a single module provide the highest probability of both probing the $15 \leq z \leq 20$ regime, as well as discovering the highest redshift galaxy possible with JWST. We also find that using multiple massive clusters in exchange for shallower survey depths is a more time efficient method of probing the $z > 15$ regime., Comment: 12 pages, 5 figures

This manuscript presents new astronomical source catalogs using data from the BUFFALO Survey. These catalogs contain detailed information for over 100,000 astronomical sources in the 6 BUFFALO clusters: Abell 370, Abell 2744, Abell S1063, MACS 0416, MACS 0717, and MACS 1149 spanning a total 240 arcmin^2. The catalogs include positions and forced photometry measurements of these objects in the F275W, F336W, F435W, F606W, F814W, F105W, F125W, F140W, and F160W HST -bands, Keck-NIRC2/VLT-HAWKI Ks band, and IRAC Channel 1 and 2 bands. Additionally, we include photometry measurements in the F475W, F625W, and F110W bands for Abell 370. This catalog also includes photometric redshift estimates computed via template fitting using LePhare. When comparing to spectroscopic reference, we obtain an outlier fraction of 9.2% and scatter, normalized median absolute deviation (NMAD), of 0.062. The catalogs are publicly available for their use by the community., Comment: 16 pages, to be submitted to ApJS

This work presents and releases a catalog of new photometrically-derived physical properties for the $\sim 10^5$ most well-measured galaxies in the COSMOS field on the sky. Using a recently developed technique, spectral energy distributions are modeled assuming a stellar initial mass function (IMF) that depends on the temperature of gas in star-forming regions. The method is applied to the largest current sample of high-quality panchromatic photometry, the COSMOS2020 catalog, that allows for testing this assumption. It is found that the galaxies exhibit a continuum of IMF, and gas temperatures, most of which are bottom-lighter than measured in the Milky Way. As a consequence, the stellar masses and star formation rates of most galaxies here are found to be lower than those measured by traditional techniques in the COSMOS2020 catalog by factors of $\sim 1.6-3.5$ and $2.5-70.0$, respectively, with the change being the strongest for the most active galaxies. The resulting physical properties provide new insights into variation of the IMF-derived gas temperature along the star-forming main sequence and at quiescence, produce a sharp and coherent picture of downsizing, as seen from the stellar mass functions, and hint at a possible high-temperature and high-density stage of early galactic evolution., Comment: ApJS in press, 25 pages, 10 figures. Catalog to be published online shortly (available on request now)

Recent studies have reported tension between the presence of luminous, high-redshift galaxies and the halo mass functions predicted by standard cosmology. Here, an improved test is proposed using the presence of high-redshift Balmer breaks to probe the formation of early $10^4 - 10^5 M_\odot$ baryonic minihalos. Unlike previous tests, this does not depend upon the mass-to-light ratio, stellar initial mass function, or star-formation history, which are all weakly constrained at high redshift. We show that the strongest Balmer breaks allowed at $z = 9$ using the simplest $\Lambda$CDM cosmological model have $D_{4000} \leq 1.26$ under idealized circumstances and $D_{4000} \leq 1.14$ including realistic feedback models. Since current photometric template fitting to JWST sources infers the existence of stronger Balmer breaks out to $z \gtrsim 11$, upcoming spectroscopic followup will either demonstrate those templates are invalid at high redshift or imply new physics beyond `vanilla' $\Lambda$CDM., Comment: 10 pages, 5 figures

We select and characterise a sample of massive (log(M$_{*}/$M$_{\odot})>10.6$) quiescent galaxies (QGs) at $33$. We compute median rest-frame SEDs for our sample and find the median quiescent galaxy at $3

The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power law plus 'external shear', which notionally accounts for neighbouring galaxies and cosmic shear. We show that it does not. Except in a handful of rare systems, the best-fit values of external shear do not correlate with independent measurements of shear: from weak lensing in 45 Hubble Space Telescope images, or in 50 mock images of lenses with complex distributions of mass. Instead, the best-fit shear is aligned with the major or minor axis of 88% of lens galaxies; and the amplitude of the external shear increases if that galaxy is disky. We conclude that 'external shear' attached to a power law model is not physically meaningful, but a fudge to compensate for lack of model complexity. Since it biases other model parameters that are interpreted as physically meaningful in several science analyses (e.g. measuring galaxy evolution, dark matter physics or cosmological parameters), we recommend that future studies of galaxy-scale strong lensing should employ more flexible mass models., Comment: 13 pages, 11 figures, submitted to MNRAS and Review Response Submitted

Using a recently-developed technique to estimate gas temperatures ($T_\textrm{SF}$) in star-forming regions from large photometric surveys, we propose a diagram, analogous to the Hertzsprung-Russell diagram for individual stars, to probe the evolution of individual galaxies. On this $T_\textrm{SF}$-sSFR (specific star formation rate) diagram, a small fraction of star-forming galaxies appear to be dominated by different feedback mechanisms than typical star-forming galaxies. These galaxies generically have younger stellar populations, lower stellar masses and increase in relative abundance towards higher redshifts, so we argue that these objects are in an earlier stage of galactic star formation. Further, Hubble observations find that these "core-forming" galaxies also exhibit distinct morphology, and that tracks on the $T_\textrm{SF}$-sSFR diagram are also a morphological sequence. Thus, unlike starburst phases which can be triggered environmentally, these earliest, core-forming galaxies, appear to be a stage that typical galaxies go through early in their star formation history. We therefore argue that most galaxies first go through a core formation stage, then subsequently disk formation, and finally become quiescent., Comment: ApJL 949, L38

Although it is generally assumed that there are two dominant classes of gamma-ray bursts (GRB) with different typical durations, it has been difficult to unambiguously classify GRBs as short or long from summary properties such as duration, spectral hardness, and spectral lag. Recent work used t-distributed stochastic neighborhood embedding (t-SNE), a machine learning algorithm for dimensionality reduction, to classify all Swift gamma-ray bursts as short or long. Here, the method is expanded, using two algorithms, t-SNE and UMAP, to produce embeddings that are used to provide a classification for the 1911 BATSE bursts, 1321 Swift bursts, and 2294 Fermi bursts for which both spectra and metadata are available. Although the embeddings appear to produce a clear separation of each catalog into short and long bursts, a resampling-based approach is used to show that a small fraction of bursts cannot be robustly classified. Further, 3 of the 304 bursts observed by both Swift and Fermi have robust but conflicting classifications. A likely interpretation is that in addition to the two predominant classes of GRBs, there are additional, uncommon types of bursts which may require multi-wavelength observations in order to separate from more typical short and long GRBs., Comment: ApJ, in press

Recent data from the James Webb Space Telescope allow a first glimpse of galaxies at $z \gtrsim 11$. The most successful tool for identifying ultra-high-redshift candidates and inferring their properties is photometric template fitting. However, current methods rely on templates derived from much lower-redshift conditions, including stellar populations older than the age of the Universe at $z > 12$, a stellar initial mass function which is physically disallowed at $z > 6$, and weaker emission lines than currently observed at $z > 7.5$. Here, two sets of synthetic templates, optimized for the expected astrophysics of galaxies at $8 < z < 12$ and $z > 12$, are developed and used to fit three galaxies at $z > 12$ from the SMACS0723 field. Using these improved templates, quantitative estimates are produced of the bias in inferred properties from JWST observations at $z>8$ due to these effects. The best-fit redshifts are similar to those found with previous template sets, but the inferred stellar masses drop by as much as 1--1.6 dex, so that stellar masses are no longer seemingly inconsistent with $\Lambda$CDM. The two new template sets are released in formats compatible with EAZY and LePhare., Comment: Additional template files available via github

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, suggest most galaxies obey an IMF top-heavier than the Galactic IMF. The implications of these revised fits on mass functions, quiescence and turnoff are discussed. At all redshifts the highest mass galaxies become quiescent first with the turnoff mass decreasing towards the present. The synchronous turnoff mass across galaxies suggests quiescence is driven by universal mechanisms rather than by stochastic or environmental processes., Comment: 11 pages, ApJ, in press

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, allows galaxies to be fit with a range of IMFs. Three surprising new features appear: (1) most star-forming galaxies are best fit with a bottom-lighter IMF than the Milky Way; (2) most star-forming galaxies at fixed redshift are fit with a very similar IMF; and (3) the most massive star-forming galaxies at fixed redshift instead exhibit a less bottom-light IMF, similar to that measured in quiescent galaxies. Additionally, since stellar masses and star formation rates both depend on the IMF, these results slightly modify the resulting relationship, while yielding similar qualitative characteristics to previous studies., Comment: 11 pages, ApJ 931, 58

A universal stellar initial mass function (IMF) should not be expected from theoretical models of star formation, but little conclusive observational evidence for a variable IMF has been uncovered. In this paper, a parameterization of the IMF is introduced into photometric template fitting of the COSMOS2015 catalog. The resulting best-fit templates suggest systematic variations in the IMF, with most galaxies exhibiting top-heavier stellar populations than in the Milky Way. At fixed redshift, only a small range of IMFs are found, with the typical IMF becoming progressively top-heavier with increasing redshift. Additionally, subpopulations of ULIRGs, quiescent- and star-forming galaxies are compared with predictions of stellar population feedback and show clear qualitative similarities to the evolution of dust temperatures., Comment: Published in the 'The Astrophysical Journal'

One of the primary goals for the upcoming James Webb Space Telescope (JWST) is to observe the first galaxies. Predictions for planned and proposed surveys have typically focused on average galaxy counts, assuming a random distribution of galaxies across the observed field. The first and most massive galaxies, however, are expected to be tightly clustered, an effect known as cosmic variance. We show that cosmic variance is likely to be the dominant contribution to uncertainty for high-redshift mass and luminosity functions, and that median high-redshift and high-mass galaxy counts for planned observations lie significantly below average counts. Several different strategies are considered for improving our understanding of the first galaxies, including adding depth, area, and independent pointings. Adding independent pointings is shown to be the most efficient both for discovering the single highest-redshift galaxy and also for constraining mass and luminosity functions., Comment: 16 pages, ApJ, in press

We present a Spitzer/IRAC survey of H-faint ($H_{160} \gtrsim 26.4$, $<5\sigma$) sources in 101 lensing cluster fields. Across a CANDELS/Wide-like survey area of $\sim$648 arcmin$^2$ (effectively $\sim$221 arcmin$^2$ in the source plane), we have securely discovered 53 sources in the IRAC Channel-2 band (CH2, 4.5 $\mathrm{\mu m}$; median CH2$=22.46\pm0.11$ AB mag) that lack robust HST/WFC3-IR F160W counterparts. The most remarkable source in our sample, namely ES-009 in the field of Abell 2813, is the brightest H-faint galaxy at 4.5 $\mathrm{\mu m}$ known so far ($\mathrm{CH2}=20.48\pm0.03$ AB mag). We show that the H-faint sources in our sample are massive (median $M_\mathrm{star} = 10^{10.3\pm 0.3}$ $M_\odot$), star-forming (median star formation rate $=100_{-40}^{+60}$ $M_\odot$yr$^{-1}$) and dust-obscured ($A_V=2.6\pm0.3$) galaxies around a median photometric redshift of $z=3.9\pm0.4$. The stellar continua of 14 H-faint galaxies can be resolved in the CH2 band, suggesting a median circularized effective radius ($R_\mathrm{e,circ}$; lensing corrected) of $1.9\pm0.2$ kpc and $<1.5$ kpc for the resolved and whole samples, respectively. This is consistent with the sizes of massive unobscured galaxies at $z\sim4$, indicating that H-faint galaxies represent the dusty tail of the distribution of a wider galaxy population. Comparing with the ALMA dust continuum sizes of similar galaxies reported previously, we conclude that the heavy dust obscuration in H-faint galaxies is related to the compactness of both stellar and dust continua ($R_\mathrm{e,circ}\sim 1$ kpc). These H-faint galaxies make up $16_{-7}^{+13}$% of the galaxies in the stellar mass range of $10^{10}-10^{11.2}$ $M_\odot$ at $z=3\sim5$, contributing to $8_{-4}^{+8}$% of the cosmic star formation rate density in this epoch and likely tracing the early phase of massive galaxy formation., Comment: 26 pages, 10 figures, 3 tables. Accepted for publication in ApJ

Recent work reported the discovery of a gamma-ray burst (GRB) associated with the galaxy GN-z11 at $z\sim 11$. The extreme improbability of the transient source being a GRB in the very early Universe requires robust elimination of all plausible alternative hypotheses. We identify numerous examples of similar transient signals in separate archival MOSFIRE observations and argue that Solar system objects -- natural or artificial -- are a far more probable explanation for these phenomena. An appendix has been added in response to additional points raised in Jiang et al. (2021), which do not change the conclusion.

Our understanding of galaxy evolution is derived from large surveys designed to maximize efficiency by only observing the minimum amount needed to infer properties for a typical galaxy. However, for a few percent of galaxies in every survey, these observations are insufficient and derived properties can be catastrophically wrong. Further, it is currently difficult or impossible to determine which objects have failed, so that these contaminate every study of galaxy properties. We develop a novel method to identify these objects by combining the astronomical codes which infer galaxy properties with the dimensionality reduction algorithm t-SNE, which groups similar objects to determine which inferred properties are out of place. This method provides an improvement for the COSMOS catalog, which already uses existing techniques for catastrophic error removal, and therefore should improve the quality of large catalogs and any studies which are sensitive to large redshift errors., Comment: 11 pages, 7 figures. Accepted for publication in ApJ

The duration of a gamma-ray burst (GRB) is a key indicator of its physics origin, with long bursts perhaps associated with the collapse of massive stars and short bursts with mergers of neutron stars.However, there is substantial overlap in the properties of both short and long GRBs and neither duration nor any other parameter so far considered completely separates the two groups. Here we unambiguously classify every GRB using a machine-learning, dimensionality-reduction algorithm, t-distributed stochastic neighborhood embedding (t-SNE), providing a catalog separating all Swift GRBs into two groups. Although the classification takes place only using prompt emission light curves,every burst with an associated supernova is found in the longer group and bursts with kilonovae in the short, suggesting along with the duration distributions that these two groups are truly long and short GRBs. Two bursts with a clear absence of a supernova belong to the longer class, indicating that these might have been direct-collapse black holes, a proposed phenomenon that may occur in the deaths of more massive stars., Comment: Accepted to Astrophysical Journal Letters, 8 pages, 4 figures

Redshifts used in current cosmological supernova samples are measured using two primary techniques, one based on well-measured host galaxy spectral lines and the other based on supernova-dominated spectra. Here, we construct an updated Pantheon catalog with revised redshifts, redshift sources and estimated uncertainties for the entire sample to investigate whether these two techniques yield consistent results. The best-fit cosmological parameters using these two measurement techniques disagree, with a supernova-only sample producing $\Omega_m$ 3.2$\sigma$ higher and $H_0$ 2.5$\sigma$ lower than a hostz-only sample, and we explore several possible sources of bias which could result from using the lower-precision supernova-dominated redshifts. In a pilot study, we show that using a host redshift-only subsample will generically produce lower $\Omega_m$ and matter density $\Omega_m h^2$ and slightly higher $H_0$ than previous analysis which, for the Pantheon dataset, could result in supernova and CMB measurements agreeing on $\Omega_m h^2$ despite tension in $H_0$. To obtain rigorous results, though, the Pantheon catalog should be improved by obtaining host spectra for supernova that have faded and future surveys should be designed to use host galaxy redshifts rather than lower-precision methods., Comment: 13 pages

Large photometric surveys provide a rich source of observations of quiescent galaxies, including a surprisingly large population at z>1. However, identifying large, but clean, samples of quiescent galaxies has proven difficult because of their near-degeneracy with interlopers such as dusty, star-forming galaxies. We describe a new technique for selecting quiescent galaxies based upon t-distributed stochastic neighbor embedding (t-SNE), an unsupervised machine learning algorithm for dimensionality reduction. This t-SNE selection provides an improvement both over UVJ, removing interlopers which otherwise would pass color selection, and over photometric template fitting, more strongly towards high redshift. Due to the similarity between the colors of high- and low-redshift quiescent galaxies, under our assumptions t-SNE outperforms template fitting in 63% of trials at redshifts where a large training sample already exists. It also may be able to select quiescent galaxies more efficiently at higher redshifts than the training sample., Comment: ApJ, in press

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/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 multi-wavelength datasets, 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 which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a 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 followup. 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., Comment: Accepted ApJS; MAST archive will be live concurrent with publication

In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter distribution in both the cluster and parallel fields. The surface mass distribution derived for the parallel field shows clumpy overdensities connected by filament-like structures elongated in the direction of the cluster core. We also characterise the X-ray emission of the cluster, and compare it with the lensing mass distribution. We identify five substructures at the $>5\sigma$ level over the two fields, four of them being in the cluster one. Furthermore, three of them are located close to the edges of the field of view, and border issues can significantly hamper the determination of their physical parameters. Finally, we compare our results with the predicted subhalo distribution of one of the Hydrangea/C-EAGLE simulated cluster. Significant differences are obtained suggesting the simulated cluster is at a more advanced evolutionary state than MACS\,J0416.1-2403. Our results anticipate the upcoming BUFFALO observations that will link the two HFF fields, extending further the \emph{HST} coverage, and thus allowing a better characterisation of the reported substructures., Comment: 15 pages, 8 figures, 3 tables, submitted to MNRAS

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the 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 have 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 four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years., Comment: 30 pages, 10 figures, accepted in ApJ

We argue that the interplay between cosmic rays, the initial mass function, and star formation plays a crucial role in regulating the star-forming "main sequence". To explore these phenomena we develop a toy model for galaxy evolution in which star formation is regulated by a combination of a temperature-dependent initial mass function and heating due to starlight, cosmic rays and, at very high redshift, the cosmic microwave background. This produces an attractor, near-equilibrium solution which is consistent with observations of the star-forming main sequence over a broad redshift range. Additional solutions to the same equations may correspond to other observed phases of galaxy evolution including quiescent galaxies. This model makes several falsifiable predictions, including higher metallicities and dust masses than anticipated at high redshift and isotopic abundances in the Milky Way. It also predicts that stellar mass-to-light ratios are lower than produced using a Milky Way-derived IMF, so that inferences of stellar masses and star formation rates for high redshift galaxies are overestimated. In some cases, this may also transform inferred dark matter profiles from core-like to cusp-like., Comment: 12 pages

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

We argue that an increased temperature in star-forming clouds alters the stellar initial mass function to be more bottom-light than in the Milky Way. At redshifts $z \gtrsim 6$, heating from the cosmic microwave background radiation produces this effect in all galaxies, and it is also present at lower redshifts in galaxies with very high star formation rates (SFRs). A failure to account for it means that at present, photometric template fitting likely overestimates stellar masses and star formation rates for the highest-redshift and highest-SFR galaxies. In addition this may resolve several outstanding problems in the chemical evolution of galactic halos., Comment: 9 pages, 5 figures. Published in MNRAS. Added further references

Telescopes are much more expensive than astronomers, so it is essential to minimize required sample sizes by using the most data-efficient statistical methods possible. However, the most commonly used model-independent techniques for finding the relationship between two variables in astronomy are flawed. In the worst case they can lead without warning to subtly yet catastrophically wrong results, and even in the best case they require more data than necessary. Unfortunately, there is no single best technique for nonparametric regression. Instead, we provide a guide for how astronomers can choose the best method for their specific problem and provide a python library with both wrappers for the most useful existing algorithms and implementations of two new algorithms developed here., Comment: 19 pages, PASP

At redshift z = 2, when the Universe was just three billion years old, half of the most massive galaxies were extremely compact and had already exhausted their fuel for star formation(1-4). It is believed that they were formed in intense nuclear starbursts and that they ultimately grew into the most massive local elliptical galaxies seen today, through mergers with minor companions(5,6), but validating this picture requires higher-resolution observations of their centres than is currently possible. Magnification from gravitational lensing offers an opportunity to resolve the inner regions of galaxies(7). Here we report an analysis of the stellar populations and kinematics of a lensed z = 2.1478 compact galaxy, which surprisingly turns out to be a fast-spinning, rotationally supported disk galaxy. Its stars must have formed in a disk, rather than in a merger-driven nuclear starburst(8). The galaxy was probably fed by streams of cold gas, which were able to penetrate the hot halo gas until they were cut off by shock heating from the dark matter halo(9). This result confirms previous indirect indications(10-13) that the first galaxies to cease star formation must have gone through major changes not just in their structure, but also in their kinematics, to evolve into present-day elliptical galaxies., Comment: Nature (June 22, 2017)

We combine star formation along the `main sequence', quiescence, and clustering and merging to produce an empirical model for the evolution of individual galaxies. Main sequence star formation alone would significantly steepen the stellar mass function towards low redshift, in sharp conflict with observation. However, a combination of star formation and merging produces a consistent result for correct choice of the merger rate function. As a result, we are motivated to propose a model in which hierarchical merging is disconnected from environmentally-independent star formation. This model can be tested via correlation functions and would produce new constraints on clustering and merging., Comment: MNRAS, in press

Using a grid of $\sim 2$ million elements ($\Delta z = 0.005$) adapted from COSMOS photometric redshift (photo-z) searches, we investigate the general properties of template-based photo-z likelihood surfaces. We find these surfaces are filled with numerous local minima and large degeneracies that generally confound rapid but "greedy" optimization schemes, even with additional stochastic sampling methods. In order to robustly and efficiently explore these surfaces, we develop BAD-Z [Brisk Annealing-Driven Redshifts (Z)], which combines ensemble Markov Chain Monte Carlo (MCMC) sampling with simulated annealing to sample arbitrarily large, pre-generated grids in approximately constant time. Using a mock catalog of 384,662 objects, we show BAD-Z samples $\sim 40$ times more efficiently compared to a brute-force counterpart while maintaining similar levels of accuracy. Our results represent first steps toward designing template-fitting photo-z approaches limited mainly by memory constraints rather than computation time., Comment: 14 pages, 8 figures; submitted to MNRAS; comments welcome

The current hierarchical merging paradigm and $\Lambda$CDM predict that the $z \sim 4-8$ universe should be a time in which the most massive galaxies are transitioning from their initial halo assembly to the later baryonic evolution seen in star-forming galaxies and quasars. However, no evidence of this transition has been found in many high redshift galaxy surveys including CFHTLS, CANDELS and SPLASH, the first studies to probe the high-mass end at these redshifts. Indeed, if halo mass to stellar mass ratios estimated at lower-redshift continue to $z \sim 6-8$, CANDELS and SPLASH report several orders of magnitude more $M \sim 10^{12-13} M_\odot$ halos than are possible to have formed by those redshifts, implying these massive galaxies formed impossibly early. We consider various systematics in the stellar synthesis models used to estimate physical parameters and possible galaxy formation scenarios in an effort to reconcile observation with theory. Although known uncertainties can greatly reduce the disparity between recent observations and cold dark matter merger simulations, even taking the most conservative view of the observations, there remains considerable tension with current theory., Comment: 10 pages, ApJ, in press

We study how AGN activity changes across environments from galaxy pairs to clusters using $143\, 843$ galaxies with $z<0.2$ from the Sloan Digital Sky Survey (SDSS). Using a refined technique, we apply a continuous measure of AGN activity, characteristic of the ionization state of the narrow-line emitting gas. Changes in key emission-line ratios ([NII]$\lambda6548$/H$\alpha$, [OIII]$\lambda5007$/H$\beta$) between different samples allow us to disentangle different environmental effects while removing contamination. We confirm that galaxy interactions enhance AGN activity. However, conditions in the central regions of clusters are inhospitable for AGN activity even if galaxies are in pairs. These results can be explained through models of gas dynamics in which pair interactions stimulate the transfer of gas to the nucleus and clusters suppress gas availability for accretion onto the central black hole., Comment: Accepted for publication in ApJ

Recent observations indicate a remarkable similarity in the properties of evolving galaxies at fixed mass and redshift, prompting us to consider the possibility that most galaxies may evolve with a common history encompassing star formation, quasar accretion, and eventual quiescence. We quantify this by defining a "synchronization timescale" for galaxies as a function of mass and redshift that characterizes the extent to which different galaxies of a common mass are evolving in the same manner at various cosmic epochs. We measure this synchronization timescale using 9 different star-forming galaxy observations from the literature and SDSS quasar observations spanning $0 < z \lesssim 6$. Surprisingly, this synchronization timescale is a constant, approximately 1.5 Gyr for all combinations of mass and time. We also find that the ratio between the stellar mass of galaxies turning off star formation and black hole mass of turnoff quasars is approximately 30:1, much lower than the 500:1 for quiescent galaxies at low redshift. As a result, we propose a model in which the star-forming "main sequence", analogous quasar behavior, and other observations form a galactic evolution "main sequence", in which star formation occurs earliest, followed by supermassive black hole accretion, and feedback between the two are dominated by deterministic rather than stochastic processes., Comment: ApJ, in press

Using the first 50% of data collected for the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH) observations on the 1.8 deg$^2$ Cosmological Evolution Survey (COSMOS) we estimate the masses and star formation rates of 3398 $M_*>10^{10}M_\odot $ star-forming galaxies at $4 < z < 6$ with a substantial population up to $M_* \gtrsim 10^{11.5} M_\odot$. We find that the strong correlation between stellar mass and star formation rate seen at lower redshift (the "main sequence" of star-forming galaxies) extends to $z\sim6$. The observed relation and scatter is consistent with a continued increase in star formation rate at fixed mass in line with extrapolations from lower-redshift observations. It is difficult to explain this continued correlation, especially for the most massive systems, unless the most massive galaxies are forming stars near their Eddington-limited rate from their first collapse. Furthermore, we find no evidence for moderate quenching at higher masses, indicating quenching either has not occurred prior to $z \sim 6$ or else occurs rapidly, so that few galaxies are visible in transition between star-forming and quenched., Comment: ApJL, accepted

Using a compilation of 25 studies from the literature, we investigate the evolution of the star-forming galaxy (SFG) Main Sequence (MS) in stellar mass and star formation rate (SFR) out to $z \sim 6$. After converting all observations to a common set of calibrations, we find a remarkable consensus among MS observations ($\sim 0.1$ dex 1$\sigma$ interpublication scatter). By fitting for time evolution of the MS in bins of constant mass, we deconvolve the observed scatter about the MS within each observed redshift bins. After accounting for observed scatter between different SFR indicators, we find the width of the MS distribution is $\sim 0.2$ dex and remains constant over cosmic time. Our best fits indicate the slope of the MS is likely time-dependent, with our best fit $\log\textrm{SFR}(M_*,t) = \left(0.84 \pm 0.02 - 0.026 \pm 0.003 \times t\right) \log M_* - \left(6.51 \pm 0.24 - 0.11 \pm 0.03 \times t\right)$, with $t$ the age of the Universe in Gyr. We use our fits to create empirical evolutionary tracks in order to constrain MS galaxy star formation histories (SFHs), finding that (1) the most accurate representations of MS SFHs are given by delayed-$\tau$ models, (2) the decline in fractional stellar mass growth for a "typical" MS galaxy today is approximately linear for most of its lifetime, and (3) scatter about the MS can be generated by galaxies evolving along identical evolutionary tracks assuming an initial $1\sigma$ spread in formation times of $\sim 1.4$ Gyr., Comment: 59 pages, 10 tables, 12 figures, accepted to ApJS; v2, slight changes to text, added new figure and fits

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2., Comment: Submitted to ApJ, 9 pages, 7 figures

We report on the discovery of a Type 1 quasar, SDSS 0956+5128, with a surprising combination of extreme velocity offsets. SDSS 0956+5128 is a broad-lined quasar exhibiting emission lines at three substantially different redshifts: a systemic redshift of z ~ 0.714 based on narrow emission lines, a broad MgII emission line centered 1200 km/s bluer than the systemic velocity, at z ~ 0.707, and broad H\alpha and H\beta emission lines centered at z ~ 0.690. The Balmer line peaks are 4100 km/s bluer than the systemic redshift. There are no previously known objects with such an extreme difference between broad MgII and broad Balmer emission. The two most promising explanations are either an extreme disk emitter or a high-velocity black hole recoil. However, neither explanation appears able to explain all of the observed features of SDSS 0956+5128, so the object may provide a challenge to our general understanding of quasar physics., Comment: ApJ, accepted

Anomalous Narrow-Line Quasars (ANLs) are a population of quasars with narrow H\beta, and sometimes [O III] broader than ~1000 km/s, in total comprising \sim 10-30% (most likely ~25%) of Type I quasars at 0.2 < z < 0.8. We find that virial masses using the H\beta and Mg II lines systematically differ for ANLs by an average of as much as 0.5 dex. Because the broad H\beta component width increases in ANLs but Mg II does not, we might suspect H\beta-based virial masses for ANLs are wrong but Mg II masses are correct. If this is due to an outflow reaching the lower-ionization potential H\beta line, C IV masses will be similarly flawed. However, we cannot be certain of this explanation without followup work, and may be unable to identify which quasars are ANLs at z > 0.8. Therefore, it is essential that ANLs be well-understood and well-modeled in order to allow the use of virial mass estimators on large optical spectroscopic catalogs, particularly at z < 0.4 and z > 2.0 where only one broad line is available for use in mass estimation., Comment: 5 pages, submitted

The H\beta\ emission line in a typical Type I quasar is composed of a broad base and a narrow core, with the core velocity typical of narrow-line region emission, and line-fitting routines typically assume this picture. We test the effects of removing this constraint, and find a substantial group of Type I quasars in the Sloan Digital Sky Survey catalog with H\beta\ emission line cores broader than 1200 km/s, above the velocity believed possible for gas in the quasar narrow-line region. We identify this group of "anomalous H\beta\ quasars" (AHQs) as a distinct population because of a variety of spectral and photometric signatures common to these AHQs but atypical of other quasars. These features are similar to some aspects of narrow-line Seyfert 1s and correlations identified by Eigenvector 1, but also contain distinct features that make AHQs difficult to classify. We demonstrate that AHQs comprise at least 11% and most likely approximately one quarter of the SDSS Type I quasar population at 0.2 < z < 0.8. For AHQs, the [O III]\lambda 4959,5007 profile is often better fit by de-linking it from the H\beta\ core, while a more standard linked fit produces a tight correlation between narrow- and broad-line velocities. We find that [O III] in AHQs sometimes has a standard narrow-line profile and other times matches the H\beta core, but is rarely in between the two, implying that the broadened core emission arises from a distinct physical region. Another feature of AHQs is a diminished [O II] line, which might indicate a connection between AHQs and the interstellar mediums of their host galaxies, through reduced photoionization or star formation. We find that it is difficult to produce AHQs using the current quasar standard model., Comment: 12 pages, accepted to PASJ

Recent work using virial mass estimates and the quasar mass-luminosity plane has yielded several new puzzles regarding quasar accretion, including a sub-Eddington boundary on most quasar accretion, near-independence of the accretion rate from properties of the host galaxy, and a cosmic synchronization of accretion among black holes of a common mass. We consider how these puzzles might change if virial mass estimation turns out to have a systematic bias. As examples, we consider two recent claims of mass-dependent biases in MgII masses. Under any such correction, the surprising cosmic synchronization of quasar accretion rates and independence from the host galaxy remain. The slope and location of the sub-Eddington boundary are very sensitive to biases in virial mass estimation, and various mass calibrations appear to favor different possible physical explanations for feedback between the central black hole and its environment. The alternative mass estimators considered do not simply remove puzzling quasar behavior, but rather replace it with new puzzles that may be more difficult to solve than those using current virial mass estimators and the Shen et al. (2008) catalog., Comment: 8 pages, 4 figures

Previous work on the quasar mass-luminosity plane indicates the possibility that quasars of the same central black hole mass might follow a common evolutionary track, independent of the properties of the host galaxy. We consider two simple models for the evolution of individual quasars. Requiring these tracks to lie within the observed quasar locus at all redshifts strongly constrains the model parameters, but does allow some solutions. These solutions include a family of tracks with similar shape but different initial masses that might match the observed quasar distributions at all redshifts z < 2.0. This family of solutions is characterized by short (1-2 Gyr) lifetimes, a duty cycle in which the quasar is on at least 25% of the time, and a rapid decline in Eddington ratio, perhaps with L/L_Edd ~ t^-6 or steeper., Comment: Accepted by MNRAS