Search

We study the stellar mass function (SMF) and the co-evolution with dark matter halos via abundance matching in the largest redshift range to date $0.25$, we find increased abundances of massive (log$\, M_{\star}/M_{\odot}>10.5$) implying integrated star formation efficiencies (SFE) $\epsilon_{\star}\equiv M_{\star}\, f_{\rm b}^{-1} M_{\rm halo}^{-1} \gtrsim 0.5$. We find a flattening of the SMF at the high-mass end that is better described by a double power law at $z>5.5$. At $z \lesssim 5.5$ it transitions to a Schechter law which coincides with the emergence of the first massive quiescent galaxies in the Universe. We trace the cosmic stellar mass density (SMD) and infer the star formation rate density (SFRD), which at $z>7.5$ agrees remarkably with recent \JWST{} UV luminosity function-derived estimates. However, at $z \lesssim 3.5$, we find significant tension ($\sim 0.3$ dex) with the cosmic star formation (SF) history from instantaneous SF measures, the causes of which remain poorly understood. We infer the stellar-to-halo mass relation (SHMR) and the SFE from abundance matching out to $z=12$, finding a non-monotonic evolution. The SFE has the characteristic strong dependence with mass in the range of $0.02 - 0.2$, and mildly decreases at the low mass end out to $z\sim3.5$. At $z\sim3.5$ the SFE increases sharply from $\sim 0.1$ to approach high SFE of $0.8-1$ by $z\sim 10$ for log$(M_{\rm h}/M_{\odot})\approx11.5$, albeit with large uncertainties. Finally, we use the SHMR to track the SFE and stellar mass growth throughout the halo history and find that they do not grow at the same rate -- from the earliest times up until $z\sim3.5$ the halo growth rate outpaces galaxy assembly, but at $z>3.5$ halo growth stagnates and accumulated gas reservoirs keep the SF going and galaxies outpace halos.

The stellar mass-star formation rate (M$_\star$ - SFR) plane is essential for distinguishing galaxy populations, but how galaxies move within this plane over cosmic time remains unclear. This study aims to describe galaxy migrations in the M$_\star$ - SFR plane by reconstructing star formation histories (SFHs) for a sample of galaxies out to redshift $z < 4$. This provides insights into the physical processes driving star formation. We use data from the COSMOS field, selecting 299131 galaxies at $z < 4$ with COSMOS-Web NIRCam data (m$_\mathrm{F444W} < 27$) over 0.54 deg$^2$. Using the SED modeling code CIGALE with non-parametric SFHs, we derive physical properties and migration vectors for these galaxies. These vectors describe the direction and velocity of evolutionary paths across the M$_\star$ - SFR plane. To assess the accuracy of these vectors, we compare them to results from the Horizon-AGN simulation. Galaxies within the main sequence show low migration amplitudes and dispersed movement directions, indicating oscillation within the main sequence. Most progenitors were already on the main sequence a billion years earlier. Starburst galaxies assembled half their mass in the last 350 Myr and originated from the main sequence. Passive galaxies show uniformly declining SFHs and include massive galaxies already in the passive region at $z = 3.5-4$, with increasing density over time. Using reconstructed SFHs up to $z < 4$, we propose a coherent picture of how galaxies migrate over cosmic time in the M$_\star$ - SFR plane, highlighting the connection between major phases in the star-formation history.

We use NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) to study the ionising properties of a sample of 15721 galaxies at $3 \leq z_{\rm{phot}} \leq 9$, 90\% complete in stellar mass down to log(M$_{\star}$/[M$_{\odot}$])$\approx 7.5$. Out of the full sample, 1620 of the galaxies have spectroscopic redshift measurements from the literature. We use the spectral energy distribution fitting code \texttt{Prospector} to fit all available photometry and infer galaxy properties. We find a significantly milder evolution of the ionising photon production efficiency (\xion\/) with redshift and UV magnitude than previously reported. Interestingly, we observe two distinct populations in \xion\/, distinguished by their burstiness (given by SFR$_{10}$/SFR$_{100}$). Both populations show the same evolution with $z$ and M$_{\rm{UV}}$, but have a different \xion\/ normalisation. We convolve the more representative $\log(\xi_{\rm{ion}} (z,\text{M}_{\rm{UV}}))$ relations (accounting for $\sim96$\% of the sample), with luminosity functions from literature, to place constraints on the cosmic ionising photon budget. By combining our results, we find that one of our models can match the observational constraints from the \lya\/ forest at $z\lesssim6$. We conclude that galaxies with M$_{\rm{UV}}$ between $-16$ and $-20$, adopting a reasonable escape fraction, can produce enough ionising photons to ionise the Universe, without exceeding the required ionising photon budget., Comment: Accepted in MNRAS. 23 pages, 21 figures

We present the JWST discovery of SN 2023adsy, a transient object located in a host galaxy JADES-GS$+53.13485$$-$$27.82088$ with a host spectroscopic redshift of $2.903\pm0.007$. The transient was identified in deep James Webb Space Telescope (JWST)/NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) program. Photometric and spectroscopic followup with NIRCam and NIRSpec, respectively, confirm the redshift and yield UV-NIR light-curve, NIR color, and spectroscopic information all consistent with a Type Ia classification. Despite its classification as a likely SN Ia, SN 2023adsy is both fairly red (E(B-V)$\sim0.9$) despite a host galaxy with low-extinction and has a high Ca II velocity ($19,000\pm2,000$km/s) compared to the general population of SNe Ia. While these characteristics are consistent with some Ca-rich SNe Ia, particularly SN 2016hnk, SN 2023adsy is intrinsically brighter than the low-z Ca-rich population. Although such an object is too red for any low-z cosmological sample, we apply a fiducial standardization approach to SN 2023adsy and find that the SN 2023adsy luminosity distance measurement is in excellent agreement ($\lesssim1\sigma$) with $\Lambda$CDM. Therefore unlike low-z Ca-rich SNe Ia, SN 2023adsy is standardizable and gives no indication that SN Ia standardized luminosities change significantly with redshift. A larger sample of distant SNe Ia is required to determine if SN Ia population characteristics at high-z truly diverge from their low-z counterparts, and to confirm that standardized luminosities nevertheless remain constant with redshift., Comment: Submitted to ApJL

We present JWST NIRCam and NIRSpec observations of a Type Ic supernova (SN Ic) and its host galaxy (JADES-GS+53.13533-27.81457) at $z = 2.83$. This SN (named SN 2023adta) was identified in deep James Webb Space Telescope (JWST)/NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) Program. Follow-up observations with JWST/NIRSpec provided a spectroscopic redshift of $z = 2.83$ and the classification as a SN Ic-BL. The light curve of SN 2023adta matches well with other stripped envelope supernovae and we find a high peak luminosity, $M_V = -19.0 \pm 0.2$ mag, based on the distribution of best-fit SNe. The broad absorption features in its spectrum are consistent with other SNe Ic-BL 1-3 weeks after peak brightness. We measure a Ca II NIR triplet expansion velocity of $29{,}000 \pm 2{,}000$ km s$^{-1}$. The host galaxy of SN 2023adta is irregular, and modeling of its spectral energy distribution (SED) indicates a metallicity of $Z = 0.35^{+0.16}_{-0.08} Z_{\odot}$. This environment is consistent with the population of low-$z$ SNe Ic-BL which prefer lower metallicities relative to other stripped envelope supernovae, and track long duration $\gamma$-ray burst (LGRB) environments. We do not identify any GRBs that are coincident with SN 2023adta. Given the rarity of SNe Ic-BL in the local universe, the detection of a SN Ic-BL at $z = 2.83$ could indicate that their rates are enhanced at high redshift., Comment: 16 pages, 7 figures, Published in ApJL

We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find an increased fraction of interacting galaxies in the starburst regime at all redshifts out to z=4. This increase correlates with star formation efficiency (SFE), but not with gas fraction. However, the correlation is weak (and only significant out to z=2), which could be explained by the short duration of SFE increase during interaction. In addition, we find that isolated disk galaxies make up a significant fraction of the starburst population. The fraction of such galaxies with star-forming clumps ("clumpy disks") is significantly increased compared to the main-sequence disk population. Furthermore, this fraction directly correlates with SFE. This is direct observational evidence for a long-term increase of SFE maintained due to disk instabilities, contributing to the majority of starburst galaxies in our sample and hence to substantial mass growth in these systems. This result could also be of importance for explaining the growth of the most massive galaxies at z>6., Comment: 14 pages, 5 figures. Submitted to ApJ. Comments welcome

We study galaxies in JADES Deep to study the evolution of the ionising photon production efficiency, $\xi_{\rm{ion}}$, observed to increase with redshift. We estimate $\xi_{\rm{ion}}$ for a sample of 677 galaxies at $z \sim 4 - 9$ using NIRCam photometry. Specifically, combinations of the medium and wide bands F335M-F356W and F410M-F444W to constrain emission lines that trace $\xi_{\rm{ion}}$: H$\alpha$ and [OIII]. Additionally, we use the spectral energy distribution fitting code \texttt{Prospector} to fit all available photometry and infer galaxy properties. The flux measurements obtained via photometry are consistent with FRESCO and NIRSpec-derived fluxes. Moreover, the emission-line-inferred measurements are in tight agreement with the \texttt{Prospector} estimates. We also confirm the observed $\xi_{\rm{ion}}$ trend with redshift and M$_{\rm{UV}}$, and find: $\log \xi_{\rm{ion}} (z,\text{M}_{\rm{UV}}) = (0.05 \pm 0.02)z + (0.11 \pm 0.02) \text{M}_{\rm{UV}} + (27.33 \pm 0.37)$. We use \texttt{Prospector} to investigate correlations of $\xi_{\rm{ion}}$ with other galaxy properties. We see a clear correlation between $\xi_{\rm{ion}}$ and burstiness in the star formation history of galaxies, given by the ratio of recent to older star formation, where burstiness is more prevalent at lower stellar masses. We also convolve our $\xi_{\rm{ion}}$ relations with luminosity functions from the literature, and constant escape fractions of 10 and 20\%, to place constraints on the cosmic ionising photon budget. By combining our results, we find that if our sample is representative of the faint low-mass galaxy population, galaxies with bursty star formation are efficient enough in producing ionising photons and could be responsible for the reionisation of the Universe., Comment: Submitted to MNRAS. 19 pages, 14 figures

Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED fitting codes. Using JWST/NIRCam images, we also produce two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. The lens is a very massive and quiescent (sSFR < 10^(-13) yr-1) elliptical galaxy at z = 2.02 \pm 0.02 with a total mass Mtot(

Interactions between galaxies leave distinguishable imprints in the form of tidal features which hold important clues about their mass assembly. Unfortunately, these structures are difficult to detect because they are low surface brightness features so deep observations are needed. Upcoming surveys promise several orders of magnitude increase in depth and sky coverage, for which automated methods for tidal feature detection will become mandatory. We test the ability of a convolutional neural network to reproduce human visual classifications for tidal detections. We use as training $\sim$6000 simulated images classified by professional astronomers. The mock Hyper Suprime Cam Subaru (HSC) images include variations with redshift, projection angle and surface brightness ($\mu_{lim}$ =26-35 mag arcsec$^{-2}$). We obtain satisfactory results with accuracy, precision and recall values of Acc=0.84, P=0.72 and R=0.85, respectively, for the test sample. While the accuracy and precision values are roughly constant for all surface brightness, the recall (completeness) is significantly affected by image depth. The recovery rate shows strong dependence on the type of tidal features: we recover all the images showing shell features and 87% of the tidal streams; these fractions are below 75% for mergers, tidal tails and bridges. When applied to real HSC images, the performance of the model worsens significantly. We speculate that this is due to the lack of realism of the simulations and take it as a warning on applying deep learning models to different data domains without prior testing on the actual data., Comment: 13 pages, 10 figures, accepted for publication in MNRAS

Surveys with James Webb Space Telescope (JWST) have discovered candidate galaxies in the first 400 Myr of cosmic time. Preliminary indications have suggested these candidate galaxies may be more massive and abundant than previously thought. However, without confirmed distances, their inferred properties remain uncertain. Here we identify four galaxies located in the JWST Advanced Deep Extragalactic Survey (JADES) Near-Infrared Camera (NIRCam) imaging with photometric redshifts z~10-13. These galaxies include the first redshift z>12 systems discovered with distances spectroscopically confirmed by JWST in a companion paper. Using stellar population modelling, we find the galaxies typically contain a hundred million solar masses in stars, in stellar populations that are less than one hundred million years old. The moderate star formation rates and compact sizes suggest elevated star formation rate surface densities, a key indicator of their formation pathways. Taken together, these measurements show that the first galaxies contributing to cosmic reionisation formed rapidly and with intense internal radiation fields., Comment: Author version of manuscript, please visit Nature Astronomy for the version published 04 April 2023

The detection of variations of fundamental constants of the Standard Model would provide us with compelling evidence of new physics, and could lift the veil on the nature of dark matter and dark energy. In this work, we discuss how a network of atomic and molecular clocks can be used to look for such variations with unprecedented sensitivity over a wide range of time scales. This is precisely the goal of the recently launched QSNET project: A network of clocks for measuring the stability of fundamental constants. QSNET will include state-of-the-art atomic clocks, but will also develop next-generation molecular and highly charged ion clocks with enhanced sensitivity to variations of fundamental constants. We describe the technological and scientific aims of QSNET and evaluate its expected performance. We show that in the range of parameters probed by QSNET, either we will discover new physics, or we will impose new constraints on violations of fundamental symmetries and a range of theories beyond the Standard Model, including dark matter and dark energy models.

The SNO+ experiment is located 2 km underground at SNOLAB in Sudbury, Canada. A low background search for neutrinoless double beta ($0\nu\beta\beta$) decay will be conducted using 780 tonnes of liquid scintillator loaded with 3.9 tonnes of natural tellurium, corresponding to 1.3 tonnes of $^{130}$Te. This paper provides a general overview of the SNO+ experiment, including detector design, construction of process plants, commissioning efforts, electronics upgrades, data acquisition systems, and calibration techniques. The SNO+ collaboration is reusing the acrylic vessel, PMT array, and electronics of the SNO detector, having made a number of experimental upgrades and essential adaptations for use with the liquid scintillator. With low backgrounds and a low energy threshold, the SNO+ collaboration will also pursue a rich physics program beyond the search for $0\nu\beta\beta$ decay, including studies of geo- and reactor antineutrinos, supernova and solar neutrinos, and exotic physics such as the search for invisible nucleon decay. The SNO+ approach to the search for $0\nu\beta\beta$ decay is scalable: a future phase with high $^{130}$Te-loading is envisioned to probe an effective Majorana mass in the inverted mass ordering region., Comment: 61 pages, 23 figures, 4 tables

We discuss the spectral energy distributions and physical properties of six galaxies whose photometric redshifts suggest they lie beyond a redshift $z\simeq$9. Each was selected on account of a prominent excess seen in the Spitzer/IRAC 4.5$\mu$m band which, for a redshift above $z=9.0$, likely indicates the presence of a rest-frame Balmer break and a stellar component that formed earlier than a redshift $z\simeq10$. In addition to constraining the earlier star formation activity on the basis of fits using stellar population models with BAGPIPES, we have undertaken the necessary, but challenging, follow-up spectroscopy for each candidate using various combinations of Keck/MOSFIRE, VLT/X-shooter, Gemini/FLAMINGOS2 and ALMA. Based on either Lyman-$\alpha$ or [OIII] 88 $\mu$m emission, we determine a convincing redshift of $z$=8.78 for GN-z-10-3 and a likely redshift of $z$=9.28 for the lensed galaxy MACS0416-JD. For GN-z9-1, we conclude the case remains promising for a source beyond $z\simeq$9. Together with earlier spectroscopic data for MACS1149-JD1, our analysis of this enlarged sample provides further support for a cosmic star formation history extending beyond redshifts $z\simeq$10. We use our best-fit stellar population models to reconstruct the past rest-frame UV luminosities of our sources and discuss the implications for tracing earlier progenitors of such systems with the James Webb Space Telescope., Comment: updated version with correction of a numerical error on MACS1149-JD1 age estimate. The predicted earlier luminosity for this galaxy has been corrected on Figure 11

A search has been performed for neutrinos from two sources, the $hep$ reaction in the solar $pp$ fusion chain and the $\nu_e$ component of the diffuse supernova neutrino background (DSNB), using the full dataset of the Sudbury Neutrino Observatory with a total exposure of 2.47 kton-years after fiducialization. The $hep$ search is performed using both a single-bin counting analysis and a likelihood fit. We find a best-fit flux that is compatible with solar model predictions while remaining consistent with zero flux, and set a one-sided upper limit of $\Phi_{hep} < 30\times10^{3}~\mathrm{cm}^{-2}~\mathrm{s}^{-1}$ [90% credible interval (CI)]. No events are observed in the DSNB search region, and we set an improved upper bound on the $\nu_e$ component of the DSNB flux of $\Phi^\mathrm{DSNB}_{\nu_e} < 19~\textrm{cm}^{-2}~\textrm{s}^{-1}$ (90% CI) in the energy range $22.9 < E_\nu < 36.9$~MeV., Comment: 11 pages, 6 figures

Neutrons produced in nuclear interactions initiated by cosmic-ray muons present an irreducible background to many rare-event searches, even in detectors located deep underground. Models for the production of these neutrons have been tested against previous experimental data, but the extrapolation to deeper sites is not well understood. Here we report results from an analysis of cosmogenically produced neutrons at the Sudbury Neutrino Observatory. A specific set of observables are presented, which can be used to benchmark the validity of GEANT4 physics models. In addition, the cosmogenic neutron yield, in units of $10^{-4}\;\text{cm}^{2}/\left(\text{g}\cdot\mu\right)$, is measured to be $7.28 \pm 0.09\;\text{stat.} ^{+1.59}_{-1.12}\;\text{syst.}$ in pure heavy water and $7.30 \pm 0.07\;\text{stat.} ^{+1.40}_{-1.02}\;\text{syst.}$ in NaCl-loaded heavy water. These results provide unique insights into this potential background source for experiments at SNOLAB.

We propose an extension of the LSST survey to cover the northern sky to DEC < +30 (accessible at airmass <1.8). This survey will increase the LSST sky coverage by ~9,600 square degrees from 18,900 to 28,500 square degrees (a 50% increase) but use only 0.6-2.5% of the time depending on the synergies with other surveys. This increased area addresses a wide range of science cases that enhance all of the primary LSST science goals by significant amounts. The science enabled includes: increasing the area of the sky accessible for follow-up of multi-messenger transients including gravitational waves, mapping the milky way halo and halo dwarfs including discovery of RR Lyrae stars in the outer galactic halo, discovery of z>7 quasars in combination Euclid, enabling a second generation DESI and other spectroscopic surveys, and enhancing all areas of science by improving synergies with Euclid, WFIRST, and unique northern survey facilities. This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and presents a unified plan for northern coverage. The range of time estimates reflects synergies with other surveys. If the modified DESC WFD survey, the ecliptic plane mini survey, and the north galactic spur mini survey are executed this plan would only need 0.6% of the LSST time, however if none of these are included the overall request is 2.5% of the 10 year survey life. In other words, the majority of these observations are already suggested as part of these other surveys and the intent of this white paper is to propose a unified baseline plan to carry out a broad range of objectives to facilitate a combination of multiple science objectives. A companion white paper gives Euclid specific science goals, and we support the white papers for southern extensions of the LSST survey., Comment: LSST Survey White Paper Submitted in Dec 2018

This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey., Comment: LSST Survey White Paper Submitted in Dec 2018

Neutron production in GeV-scale neutrino interactions is a poorly studied process. We have measured the neutron multiplicities in atmospheric neutrino interactions in the Sudbury Neutrino Observatory experiment and compared them to the prediction of a Monte Carlo simulation using GENIE and a minimally modified version of GEANT4. We analyzed 837 days of exposure corresponding to Phase I, using pure heavy water, and Phase II, using a mixture of Cl in heavy water. Neutrons produced in atmospheric neutrino interactions were identified with an efficiency of $15.3\%$ and $44.3\%$, for Phase I and II respectively. The neutron production is measured as a function of the visible energy of the neutrino interaction and, for charged current quasi-elastic interaction candidates, also as a function of the neutrino energy. This study is also performed classifying the complete sample into two pairs of event categories: charged current quasi-elastic and non charged current quasi-elastic, and $\nu_{\mu}$ and $\nu_e$. Results show good overall agreement between data and Monte Carlo for both phases, with some small tension with a statistical significance below $2\sigma$ for some intermediate energies.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

The long baseline between the Earth and the Sun makes solar neutrinos an excellent test beam for exploring possible neutrino decay. The signature of such decay would be an energy-dependent distortion of the traditional survival probability which can be fit for using well-developed and high precision analysis methods. Here a model including neutrino decay is fit to all three phases of $^8$B solar neutrino data taken by the Sudbury Neutrino Observatory. This fit constrains the lifetime of neutrino mass state $\nu_2$ to be ${>8.08\times10^{-5}}$ s/eV at $90\%$ confidence. An analysis combining this SNO result with those from other solar neutrino experiments results in a combined limit for the lifetime of mass state $\nu_2$ of ${>1.04\times10^{-3}}$ s/eV at $99\%$ confidence.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science., Comment: White Paper in response to LSST Call for Observing Strategy Input

Experimental tests of Lorentz symmetry in systems of all types are critical for ensuring that the basic assumptions of physics are well-founded. Data from all phases of the Sudbury Neutrino Observatory, a kiloton-scale heavy water Cherenkov detector, are analyzed for possible violations of Lorentz symmetry in the neutrino sector. Such violations would appear as one of eight possible signal types in the detector: six seasonal variations in the solar electron neutrino survival probability differing in energy and time dependence, and two shape changes to the oscillated solar neutrino energy spectrum. No evidence for such signals is observed, and limits on the size of such effects are established in the framework of the Standard Model Extension, including 40 limits on perviously unconstrained operators and improved limits on 15 additional operators. This makes limits on all minimal, Dirac-type Lorentz violating operators in the neutrino sector available for the first time.

Tests on $B-L$ symmetry breaking models are important probes to search for new physics. One proposed model with $\Delta(B-L)=2$ involves the oscillations of a neutron to an antineutron. In this paper a new limit on this process is derived for the data acquired from all three operational phases of the Sudbury Neutrino Observatory experiment. The search was concentrated in oscillations occurring within the deuteron, and 23 events are observed against a background expectation of 30.5 events. These translate to a lower limit on the nuclear lifetime of $1.48\times 10^{31}$ years at 90% confidence level (CL) when no restriction is placed on the signal likelihood space (unbounded). Alternatively, a lower limit on the nuclear lifetime was found to be $1.18\times 10^{31}$ years at 90% CL when the signal was forced into a positive likelihood space (bounded). Values for the free oscillation time derived from various models are also provided in this article. This is the first search for neutron-antineutron oscillation with the deuteron as a target., Comment: 14 pages, 8 figures

We review the uncertainties in high-z star-formation rate (SFR) measures and the constraints that one obtains from high-z gamma-ray burst (GRB) rates on them. We show that at the present time, the GRB rates per unit star-formation at z>3 are higher than at lower redshift. We also compare metallicity predictions made using a hierarchical model of cosmic chemical evolution based on two recently proposed SFRs, one based on the observed galaxy luminosity function at high redshift and one based on the GRB rate and find that within the considerable scatter in metal abundance measures, they both are consistent with the data. Analyzing the ensemble of different measurements together, we conclude that despite metallicity biases, GRBs may be a less biased probe of star-formation at z>3 than at z<2. There is likely to be a common origin to the high GRB rate per unit star-formation and the high observed Lyman-continuum production rate in high redshift galaxies and that this may be due to a relative overabundance of stars with mass $>$25Msun which are likely GRB progenitors. We also find that to reconcile these measurements with the Thomson scattering cross section of cosmic microwave background (CMB) photons measured by Planck, the escape fraction of Lyman-continuum photons from galaxies must be low, about 15 percent or less and that the clumping factor of the IGM is likely to be small, about 3. Finally, we demonstrate that GRBs are unique probes of metallicity evolution in low-mass galaxy samples and that GRB hosts likely lost a significant fraction of metals to the intergalactic medium (IGM) due to feedback processes such as stellar winds and supernovae., Comment: 17 pages, 3 figures, accepted to Space Science Reviews

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge., Comment: 11 pages, 4 figures, accepted for publication in ApJ

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts ( = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4., Comment: 26 pages, 15 figures, updated to match version accepted for publication in MNRAS

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events., Comment: submitted to ApJ; 18 pages; see http://www.stsci.edu/hst/campaigns/frontier-fields/ and http://ssc.spitzer.caltech.edu/warmmission/scheduling/approvedprograms/ddt/frontier/ for data and more information

Long-duration Gamma-Ray Bursts (GRBs) allow us to pinpoint and study star-forming galaxies in the early universe, thanks to their orders of magnitude brighter peak luminosities compared to other astrophysical sources, and their association with deaths of massive stars. We present Hubble Space Telescope Wide Field Camera 3 detections of three Swift GRB host galaxies lying at redshifts $z = 5.913$ (GRB 130606A), $z = 6.295$ (GRB 050904), and $z = 6.327$ (GRB 140515A) in the F140W (wide-$JH$ band, $\lambda_{\rm{obs}}\sim1.4\,\mu m$) filter. The hosts have magnitudes (corrected for Galactic extinction) of $m_{\rm{\lambda_{obs},AB}}= 26.34^{+0.14}_{-0.16}, 27.56^{+0.18}_{-0.22},$ and $28.30^{+0.25}_{-0.33}$ respectively. In all three cases the probability of chance coincidence of lower redshift galaxies is $\lesssim2\,\%$, indicating that the detected galaxies are most likely the GRB hosts. These are the first detections of high redshift ($z > 5$) GRB host galaxies in emission. The galaxies have luminosities in the range $0.1-0.6\,L^{*}_{z=6}$ (with $M_{1600}^{*}=-20.95\pm0.12$), and half-light radii in the range $0.6-0.9\,\rm{kpc}$. Both their half-light radii and luminosities are consistent with existing samples of Lyman-break galaxies at $z\sim6$. Spectroscopic analysis of the GRB afterglows indicate low metallicities ($[\rm{M/H}]\lesssim-1$) and low dust extinction ($A_{\rm{V}}\lesssim0.1$) along the line of sight. Using stellar population synthesis models, we explore the implications of each galaxy's luminosity for its possible star formation history, and consider the potential for emission-line metallicity determination with the upcoming James Webb Space Telescope., Comment: Published in ApJ

Decision trees are a popular technique in statistical data classification. They recursively partition the feature space into disjoint sub-regions until each sub-region becomes homogeneous with respect to a particular class. The basic Classification and Regression Tree (CART) algorithm partitions the feature space using axis parallel splits. When the true decision boundaries are not aligned with the feature axes, this approach can produce a complicated boundary structure. Oblique decision trees use oblique decision boundaries to potentially simplify the boundary structure. The major limitation of this approach is that the tree induction algorithm is computationally expensive. In this article we present a new decision tree algorithm, called HHCART. The method utilizes a series of Householder matrices to reflect the training data at each node during the tree construction. Each reflection is based on the directions of the eigenvectors from each classes' covariance matrix. Considering axis parallel splits in the reflected training data provides an efficient way of finding oblique splits in the unreflected training data. Experimental results show that the accuracy and size of the HHCART trees are comparable with some benchmark methods in the literature. The appealing feature of HHCART is that it can handle both qualitative and quantitative features in the same oblique split., Comment: 13 Pages, 1 Figure, 4 Tables, 1 Algorithm

We present the results of a study investigating the sizes and morphologies of redshift 4 < z < 8 galaxies in the CANDELS GOODS-S, HUDF and HUDF parallel fields. Based on non-parametric measurements and incorporating a careful treatment of measurement biases, we quantify the typical size of galaxies at each redshift as the peak of the log-normal size distribution, rather than the arithmetic mean size. Parameterizing the evolution of galaxy half-light radius as $r_{50} \propto (1+z)^n$, we find $n = -0.20 \pm 0.26$ at bright UV-luminosities ($0.3L_{*(z=3)} < L < L_*$) and $n = -0.47 \pm 0.62$ at faint luminosities ($0.12L_* < L < 0.3L_*$). Furthermore, simulations based on artificially redshifting our z~4 galaxy sample show that we cannot reject the null hypothesis of no size evolution. We show that this result is caused by a combination of the size-dependent completeness of high-redshift galaxy samples and the underestimation of the sizes of the largest galaxies at a given epoch. To explore the evolution of galaxy morphology we first compare asymmetry measurements to those from a large sample of simulated single S\'ersic profiles, in order to robustly categorise galaxies as either `smooth' or `disturbed'. Comparing the disturbed fraction amongst bright ($M_{UV} \leq -20$) galaxies at each redshift to that obtained by artificially redshifting our z~4 galaxy sample, while carefully matching the size and UV-luminosity distributions, we find no clear evidence for evolution in galaxy morphology over the redshift interval 4 < z < 8. Therefore, based on our results, a bright ($M_{UV} \leq -20$) galaxy at z~6 is no more likely to be measured as `disturbed' than a comparable galaxy at z~4, given the current observational constraints., Comment: 29 pages, 25 figures, 4 tables, published in MNRAS

The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astrophysical bursts of widely varying intensity. No correlations were found between neutrinos detected in SNO and such astrophysical sources.

We present a new determination of the UV galaxy luminosity function (LF) at redshift z ~ 7 and z ~ 8, and a first estimate at z ~ 9. An accurate determination of the form and evolution of the LF at high z is crucial for improving our knowledge of early galaxy evolution and cosmic reionization. Our analysis exploits fully the new, deepest WFC3/IR imaging from our HST UDF12 campaign, and includes a new, consistent analysis of all appropriate, shallower/wider-area HST data. Our new measurement of the evolving LF at z ~ 7-8 is based on a final catalogue of ~600 galaxies, and involves a step-wise maximum likelihood determination based on the redshift probability distribution for each object; this makes full use of the 11-band imaging now available in the HUDF, including the new UDF12 F140W data, and the deep Spitzer IRAC imaging. The final result is a determination of the z ~ 7 LF extending down to M_UV = -16.75, and the z ~ 8 LF down to M_UV = -17.00. Fitting a Schechter function, we find M* = -19.90 (+0.23/-0.28), log phi* = -2.96 (+0.18/-0.23), and a faint-end slope alpha=-1.90 (+0.14/-0.15) at z~7, and M* = -20.12 (+0.37/-0.48), log phi* = -3.35 (+0.28/-0.47), alpha=-2.02 (+0.22/-0.23) at z~8. These results strengthen suggestions that the evolution at z > 7 is more akin to `density evolution' than the apparent `luminosity evolution' seen at z ~ 5-7. We also provide the first meaningful information on the LF at z ~ 9, explore alternative extrapolations to higher z, and consider the implications for the evolution of UV luminosity density. Finally, we provide catalogues (including z_phot, M_UV and all photometry) for the 100 most robust z~6.5-11.9 galaxies in the HUDF used in this analysis. We discuss our results in the context of earlier work and the results of an independent analysis of the UDF12 data based on colour-colour selection (Schenker et al. 2013)., Comment: 25 pages, 7 figures, accepted for publication in MNRAS. Full data tables for all eight survey fields included in the appendix

We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a new sample of galaxies at z ~ 9. Utilising galaxy selection in the new F140W imaging to minimize colour bias, and applying both colour and power-law estimators of beta, we find beta = -2.1 (+/-0.2) at z ~ 7 for galaxies with M_UV ~ -18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z ~ 8 we find a similar value, beta = -1.9 (+/-0.3). At z ~ 9, we find beta = -1.8 (+/-0.6), essentially unchanged from z ~ 6 - 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in beta within our new, robust z ~ 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either ~ solar metallicity and zero dust, or moderately sub-solar (~ 10-20%) metallicity with modest dust obscuration (A_V ~ 0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with beta ~ -2.5 may not emerge until M_UV > -16, a regime likely to remain inaccessible until the James Webb Space Telescope., Comment: 15 pages, 7 figures, updated to match accepted version, now including an appendix tabulating the individual beta values for the objects used in this study

We report the results of a study exploring the stellar populations of 13 luminous (L>L*), spectroscopically confirmed, galaxies in the redshift interval 5.5= 300 Myr, the degeneracies introduced by dust extinction mean that only two of these objects actually require a >300 Myr old stellar population to reproduce the observed photometry. Moreover, when considering only smoothly-varying star-formation histories, we observe a clear tension between the data and models such that a galaxy SED template with an old age is often chosen in order to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break this tension we explore SED fitting with two-component models (burst plus on-going star-formation) and allow for nebular emission. On average, the inclusion of nebular emission leads to lower stellar-mass estimates (median offset 0.18 dex), moderately higher specific star-formation rates, and allows for a wider range of plausible stellar ages. However, based on our SED modelling, we find no strong evidence for extremely young ages in our sample (<50 Myr). Finally, considering all of the different star-formation histories explored, we find that the median best-fitting ages are of the order 200-300 Myr and that the objects with the tightest constraints indicate ages in the range 50-200 Myr (Abridged)., Comment: 23 pages, 12 figures, for submission to MNRAS

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive., Comment: Nature, in press

We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8B decays in the Sun of (5.25 \pm 0.16(stat.)+0.11-0.13(syst.))\times10^6 cm^{-2}s^{-1}. A two-flavor neutrino oscillation analysis yielded \Deltam^2_{21} = (5.6^{+1.9}_{-1.4})\times10^{-5} eV^2 and tan^2{\theta}_{12}= 0.427^{+0.033}_{-0.029}. A three-flavor neutrino oscillation analysis combining this result with results of all other solar neutrino experiments and the KamLAND experiment yielded \Deltam^2_{21} = (7.41^{+0.21}_{-0.19})\times10^{-5} eV^2, tan^2{\theta}_{12} = 0.446^{+0.030}_{-0.029}, and sin^2{\theta}_{13} = (2.5^{+1.8}_{-1.5})\times10^{-2}. This implied an upper bound of sin^2{\theta}_{13} < 0.053 at the 95% confidence level (C.L.).

This paper details the solar neutrino analysis of the 385.17-day Phase-III data set acquired by the Sudbury Neutrino Observatory (SNO). An array of $^3$He proportional counters was installed in the heavy-water target to measure precisely the rate of neutrino-deuteron neutral-current interactions. This technique to determine the total active $^8$B solar neutrino flux was largely independent of the methods employed in previous phases. The total flux of active neutrinos was measured to be $5.54^{+0.33}_{-0.31}(stat.)^{+0.36}_{-0.34}(syst.)\times 10^{6}$ cm$^{-2}$ s$^{-1}$, consistent with previous measurements and standard solar models. A global analysis of solar and reactor neutrino mixing parameters yielded the best-fit values of $\Delta m^2 = 7.59^{+0.19}_{-0.21}\times 10^{-5}{eV}^2$ and $\theta = 34.4^{+1.3}_{-1.2}$ degrees., Comment: 104 pages, 26 figures, 20 tables, submitted to Phys. Rev. C

Using deep Keck spectroscopy of Lyman break galaxies selected from infrared imaging data taken with WFC3/IR onboard the Hubble Space Telescope, we present new evidence for a reversal in the redshift-dependent fraction of star forming galaxies with detectable Lyman alpha emission in the redshift range 6.3 < z < 8.8. Our earlier surveys with the DEIMOS spectrograph demonstrated a significant increase with redshift in the fraction of line emitting galaxies over the interval 4 < z < 6, particularly for intrinsically faint systems which dominate the luminosity density. Using the longer wavelength sensitivities of LRIS and NIRSPEC, we have targeted 19 Lyman break galaxies selected using recent WFC3/IR data whose photometric redshifts are in the range 6.3 < z < 8.8 and which span a wide range of intrinsic luminosities. Our spectroscopic exposures typically reach a 5-sigma sensitivity of < 50 A for the rest-frame equivalent width (EW) of Lyman alpha emission. Despite the high fraction of emitters seen only a few hundred million years later, we find only 2 convincing and 1 possible line emitter in our more distant sample. Combining with published data on a further 7 sources obtained using FORS2 on the ESO VLT, and assuming continuity in the trends found at lower redshift, we discuss the significance of this apparent reversal in the redshift-dependent Lyman alpha fraction in the context of our range in continuum luminosity. Assuming all the targeted sources are at their photometric redshift and our assumptions about the Lyman alpha EW distribution are correct, we would expect to find so few emitters in less than 1% of the realizations drawn from our lower redshift samples. Our new results provide further support for the suggestion that, at the redshifts now being probed spectroscopically, we are entering the era where the intergalactic medium is partially neutral., Comment: 8 pages, 5 figures, Accepted to ApJ 10/1/11

It has recently been reported that the faintest galaxies at z~6-7 display extremely blue UV continuum slopes, with a UV power-law index beta ~ -3. Such slopes are bluer than previously reported for any other galaxy population, and imply extinction-free, young, very low-metallicity stellar populations with a high ionizing photon escape fraction. Here we undertake a critical study of the evidence for such extreme values of beta, combining three new WFC3/IR-selected samples of galaxies spanning ~2 decades in UV luminosity over the redshift range z~4.5-8. We explore the impact of inclusion/exclusion of less robust high-z candidates, and use the varying depths of the samples to explore the effects of noise and selection bias. Simple data-consistency arguments suggest that artificially blue average values of beta can result when the analysis is extended into the deepest ~ 0.5-mag bin of these WFC3/IR-selected samples, regardless of the actual luminosity or z range probed. By confining attention to robust, well-detected high-z galaxy candidates, we find that the average value of beta is consistent with -2.05 +/- 0.10 for z=5-7, and -22 < M_UV < -18. We create and analyse a set of simulations which demonstrate that a bias towards artifically low/blue average values of beta is indeed expected when the UV slope analysis is extended towards the source detection threshold, and conclude that there is as yet no clear evidence for UV slopes significantly bluer than beta ~ -2, the typical value displayed by the bluest star-forming galaxies at more modest z. A robust measurement of beta for the faintest galaxies at z~7-8 remains a key observational goal, as it provides a fundamental test for high escape fractions from a potentially abundant source of reionizing photons. This goal is achievable with HST, but requires still deeper WFC3/IR imaging in the HUDF., Comment: 12 pages, 10 figures, minor modifications made to make arXiv version match the accepted/published MNRAS version