We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST “Prime Extragalactic Areas for Reionization and Lensing Science” (PEARLS) project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift protoclusters; and the iconic backlit VV 191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, active galactic nucleus (AGN) growth, and First Light. Five fields—the JWST NEP Time-Domain Field (TDF), IRAC Dark Field, and three lensing clusters—will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9–4.5 μm galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9–4.5 μm. PEARLS is designed to be of lasting benefit to the community.

We present an analysis of the ultraviolet luminosity function (UV LF) and star formation rate density of distant galaxies (7.5 < z < 13.5) in the “blank” fields of the Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) survey combined with Early Release Science data from the CEERS, GLASS, and NGDEEP surveys/fields and the first data release of JADES. We use strict quality cuts on EAZY photometric redshifts to obtain a reliable selection and characterization of high-redshift (z > 6.5) galaxies from a consistently processed set of deep, near-infrared imaging. Within an area of 180 arcmin², we identify 1046 candidate galaxies at redshifts z > 6.5 and we use this sample to study the UV LF in four redshift bins between 7.5 < z < 13.5. The measured number density of galaxies at z = 8 and z = 9 matches those of past observations undertaken by the Hubble Space Telescope (HST). Our z = 10.5 measurements lie between early James Webb Space Telescope (JWST) results and past HST results, indicating cosmic variance may be the cause of previous high density measurements. However, the number densities of UV-luminous galaxies at z = 12.5 are high compared to predictions from simulations. When examining the star formation rate density of galaxies at this period, our observations are still largely consistent with a constant star formation efficiency, are slightly lower than previous early estimations using JWST, and support galaxy driven reionization at z ≤ 8.

On 2022 July 13, NASA released to the whole world the data obtained by the James Webb Space Telescope (JWST) Early Release Observations (ERO). These are the first set of science-grade data from this long-awaited facility, marking the beginning of a new era in astronomy. In the study of the early universe, JWST will allow us to push far beyond z ≈ 11, the redshift boundary previously imposed by the 1.7 μm red cutoff of the Hubble Space Telescope (HST). In contrast, JWST’s NIRCam reaches ∼5 μm. Among the JWST ERO targets there is a nearby galaxy cluster SMACS 0723-73, which is a massive cluster and has been long recognized as a potential “cosmic telescope” in amplifying background galaxies. The ERO six-band NIRCam observations on this target have covered an additional flanking field not boosted by gravitational lensing, which also sees far beyond HST. Here we report the result from our search of candidate objects at z > 11 using these ERO data. In total, there are 87 such objects identified by using the standard “dropout” technique. These objects are all detected in multiple bands and therefore cannot be spurious. For most of them, their multiband colors are inconsistent with known types of contaminants. If the detected dropout signature is interpreted as the expected Lyman break, it implies that these objects are at z ≈ 11–20. The large number of such candidate objects at such high redshifts is not expected from the previously favored predictions and demands further investigations. JWST spectroscopy on such objects will be critical.

We utilize deep JWST Near Infrared Camera (NIRCam) observations for the first direct constraints on the Galaxy Stellar Mass Function (GSMF) at z > 10. Our EPOCHS v1 sample includes 1120 galaxy candidates at 6.5 < z < 13.5 taken from a consistent reduction and analysis of publicly available deep JWST NIRCam data covering the Prime Extragalactic Areas for Reionization Science, CEERS, GLASS, JADES GOOD-S, NGDEEP, and SMACS0723 surveys, totaling 187 arcmin². We investigate the impact of spectral energy distribution fitting methods, assumed star formation histories (SFHs), dust laws, and priors on galaxy masses and the resultant GSMF. While our fiducial GSMF agrees with the literature at z < 13.5, we find that the assumed SFH model has a large impact on the GSMF and stellar mass density (SMD), finding a 0.75 dex increase in the SMD at z = 10.5 between a flexible nonparametric and standard parametric SFH. Overall, we find a flatter SMD evolution at z ≥ 9 than some studies predict, suggesting a rapid buildup of stellar mass in the early Universe. We find no incompatibility between our results and those of standard cosmological models, as suggested previously, although the most massive galaxies may require a high star formation efficiency. We find that the “little red dot” galaxies dominate the z = 7 GSMF at high masses, necessitating a better understanding of the relative contributions of active galactic nucleus and stellar emission. We show that assuming a theoretically motivated top-heavy initial mass function (IMF) reduces stellar mass by 0.5 dex without affecting fit quality, but our results remain consistent with existing cosmological models with a standard IMF.

We present a multiwavelength analysis using the Submillimeter Array (SMA), James Clerk Maxwell Telescope, NOEMA, JWST, the Hubble Space Telescope (HST), and the Spitzer Space Telescope of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A 1489. These SMA-located sources both lie at z = 4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at ≲ 2 μm, even with deep JWST/NIRCam observations. We investigate their stellar, interstellar medium, and dynamical properties, including a pixel-level spectral energy distribution analysis to derive subkiloparsec-resolution stellar-mass and A _V maps. We find that 850.1 is one of the most massive and highly obscured, A _V ∼ 5, galaxies known at z > 4 with M _* ∼10^11.8 M _⊙ (likely forming at z > 6), and 850.2 is one of the least massive and least obscured, A _V ∼ 1, members of the z > 4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z ≳ 3–4 based on imaging at ≲ 2 μm, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar-mass distribution, R e mass ∼1 kpc, but its expected evolution means that it matches both the properties of massive, quiescent galaxies at z ∼ 1.5 and ultramassive early-type galaxies at z ∼ 0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature that we link to the active dynamical environment it resides within.

With its unprecedented sensitivity and spatial resolution, the James Webb Space Telescope (JWST) has opened a new window for time-domain discoveries in the infrared. Here we report observations in the only field that has received four epochs (spanning 126 days) of JWST NIRCam observations in Cycle 1. This field is toward MACS J0416.1−2403, which is a rich galaxy cluster at redshift z = 0.4 and is one of the Hubble Frontier Fields. We have discovered 14 transients from these data. Twelve of these transients happened in three galaxies (with z = 0.94, 1.01, and 2.091) crossing a lensing caustic of the cluster, and these transients are highly magnified by gravitational lensing. These 12 transients are likely of a similar nature to those previously reported based on the Hubble Space Telescope (HST) data in this field, i.e., individual stars in the highly magnified arcs. However, these 12 could not have been found by HST because they were too red and too faint. The other two transients are associated with background galaxies (z = 2.205 and 0.7093) that are only moderately magnified, and they are likely supernovae. They indicate a demagnified supernova surface density, when monitored at a time cadence of a few months to a ∼3–4 μm survey limit of AB ∼28.5 mag, of ∼0.5 arcmin⁻² integrated to z ≈ 2. This survey depth is beyond the capability of HST but can be easily reached by JWST.