Search

We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5

New JWST/NIRCam wide-field slitless spectroscopy provides redshifts for four z>8 galaxies located behind the lensing cluster MACS J0416.1-2403. Two of them, "Y1" and "JD", have previously reported spectroscopic redshifts based on ALMA measurements of [OIII] 88 $\mu$m and/or [CII] 157.7 $\mu$m lines. Y1 is a merging system of three components, and the existing redshift z=8.31 is confirmed. However, JD is at z=8.34 instead of the previously claimed z=9.28. JD's close companion, "JD-N", which was a previously discovered z>8 candidate, is now identified at the same redshift as JD. JD and JD-N form an interacting pair. A new candidate at z>8, "f090d_018", is also confirmed and is at z=8.49. These four objects are likely part of an overdensity that signposts a large structure extending ~165 kpc in projected distance and ~48.7 Mpc in radial distance. They are magnified by less than one magnitude and have intrinsic $M_{UV}$ ranging from -19.57 to -20.83 mag. Their spectral energy distributions show that the galaxies are all very young with ages ~ 4-18 Myr and stellar masses about $10^{7-8}$ ${\rm M_\odot}$. These infant galaxies have very different star formation rates ranging from a few to over a hundred $\rm{M_\odot}$ yr$^{-1}$, but only two of them (JD and f090d_018) have blue rest-frame UV slopes $\beta<-2.0$ indicative of a high Lyman-continuum photon escape fraction that could contribute significantly to the cosmic hydrogen-reionizing background. Interestingly, these two galaxies are the least massive and least active ones among the four. The other two systems have much flatter UV slopes largely because of their high dust extinction ($A_{\rm V}$=0.9-1.0 mag). Their much lower indicated escape fractions show that even very young, actively star-forming galaxies can have negligible contribution to reionization when they quickly form dust throughout their bodies., Comment: 18 pages, 6 figures, after addressing the referee report

We present an analysis of rest-frame UV continuum slopes, $\beta$, using a sample of 1011 galaxies at $6.5

The first public 0.9-4.4{\mu}m NIRCam images of the North Ecliptic Pole (NEP) Time Domain Field (TDF) uncovered galaxies displaying point-source features in their cores as seen in the longer wavelength filters. We visually identified a sample of 66 galaxies (~1 galaxy per arcmin2) with point-like cores and have modeled their two-dimensional light profiles with GalFit, identifying 16 galactic nuclei with measurable point-source components. GalFit suggests the visual sample is a mix of both compact stellar bulge and point-source galaxy cores. This core classification is complemented by spectral energy distribution (SED) modeling to infer the sample's active galactic nucleus (AGN) and host-galaxy parameters. For galaxies with measurable point-source components, the median fractional AGN contribution to their 0.1-30.0{\mu}m flux is 0.44, and 14/16 are color-classified AGN. We conclude that near-infrared point-source galaxy cores are signatures of AGN. In addition, we define an automated sample-selection criterion to identify these point-source features. These criteria can be used in other extant and future NIRCam images to streamline the search for galaxies with unresolved IR-luminous AGN. The James Webb Space Telescope's superb angular resolution and sensitivity at infrared wavelengths is resurrecting the morphological identification of AGN., Comment: 17 pages, 9 figures, 2 tables

We present a new parametric lens model for the G165.7+67.0 galaxy cluster, which was discovered with $Planck$ through its bright submillimeter flux, originating from a pair of extraordinary dusty star-forming galaxies (DSFGs) at $z\approx 2.2$. Using JWST and interferometric mm/radio observations, we characterize the intrinsic physical properties of the DSFGs, which are separated by only $\sim 1^{\prime\prime}$ (8 kpc) and a velocity difference $\Delta V \lesssim 600~{\rm km}~{\rm s}^{-1}$ in the source plane, and thus likely undergoing a major merger. Boasting intrinsic star formation rates ${\rm SFR}_{\rm IR} = 320 \pm 70$ and $400 \pm 80~ M_\odot~{\rm yr}^{-1}$, stellar masses ${\rm log}[M_\star/M_\odot] = 10.2 \pm 0.1$ and $10.3 \pm 0.1$, and dust attenuations $A_V = 1.5 \pm 0.3$ and $1.2 \pm 0.3$, they are remarkably similar objects. We perform spatially-resolved pixel-by-pixel SED fitting using rest-frame near-UV to near-IR imaging from JWST/NIRCam for both galaxies, resolving some stellar structures down to 100 pc scales. Based on their resolved specific SFRs and $UVJ$ colors, both DSFGs are experiencing significant galaxy-scale star formation events. If they are indeed interacting gravitationally, this strong starburst could be the hallmark of gas that has been disrupted by an initial close passage. In contrast, the host galaxy of the recently discovered triply-imaged SN H0pe has a much lower SFR than the DSFGs, and we present evidence for the onset of inside-out quenching and large column densities of dust even in regions of low specific SFR. Based on the intrinsic SFRs of the DSFGs inferred from UV through FIR SED modeling, this pair of objects alone is predicted to yield an observable $1.1 \pm 0.2~{\rm CCSNe~yr}^{-1}$, making this cluster field ripe for continued monitoring., Comment: 47 pages, 21 figures, 5 tables. Submitted to ApJ, comments welcome!

We present a new lens model for the $z=0.396$ galaxy cluster MACS J0416.1$-$2403 based on a previously known set of 77 spectroscopically confirmed, multiply imaged galaxies plus an additional set of 42 candidate multiply imaged galaxies from past HST and new JWST data. The new galaxies lack spectroscopic redshifts but have geometric and/or photometric redshift estimates that are presented here. The new model predicts magnifications and time delays for all multiple images. The full set of constraints totals 343, constituting the largest sample of multiple images lensed by a single cluster to date. Caustic-crossing galaxies lensed by this cluster are especially interesting. Some of these galaxies show transient events, most of which are interpreted as micro-lensing of stars at cosmological distances. These caustic-crossing arcs are expected to show similar events in future, deeper JWST observations. We provide time delay and magnification models for all these arcs. The time delays and the magnifications for different arcs are generally anti-correlated, as expected from $N$-body simulations. In the major sub-halos of the cluster, the dark-matter mass from our lens model correlates well with the observed number of globular clusters. This confirms earlier results, derived at lower redshifts, which suggest that globular clusters can be used as powerful mass proxies for the halo masses when lensing constraints are scarce or not available., Comment: 21 pages and 11 figures. Accepted for publication in A&A

The Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS), a JWST GTO program, obtained a set of unique NIRCam observations that have enabled us to significantly improve the default photometric calibration across both NIRCam modules. The observations consisted of three epochs of 4-band (F150W, F200W, F356W, and F444W) NIRCam imaging in the Spitzer IRAC Dark Field (IDF). The three epochs were six months apart and spanned the full duration of Cycle 1. As the IDF is in the JWST continuous viewing zone, we were able to design the observations such that the two modules of NIRCam, modules A and B, were flipped by 180 degrees and completely overlapped each other's footprints in alternate epochs. We were therefore able to directly compare the photometry of the same objects observed with different modules and detectors, and we found significant photometric residuals up to ~ 0.05 mag in some detectors and filters, for the default version of the calibration files that we used (jwst_1039.pmap). Moreover, there are multiplicative gradients present in the data obtained in the two long-wavelength bands. The problem is less severe in the data reduced using the latest pmap (jwst_1130.pmap as of September 2023), but it is still present, and is non-negligible. We provide a recipe to correct for this systematic effect to bring the two modules onto a more consistent calibration, to a photometric precision better than ~ 0.02 mag., Comment: 12 pages, 8 figures, 3 tables. Accepted to PASP

We consider the effect of including an active galactic nuclei (AGN) component when fitting spectral energy distributions of 109 spectroscopically confirmed $z\approx 3.5-12.5$ galaxies with JWST. Remarkably, we find that the resulting cosmic star formation history is $\approx 0.4$ dex lower at $z\gtrsim 9.5$ when an AGN component is included in the fitting. This alleviates previously reported excess star formation at $z\gtrsim 9.5$ compared to models based on typical baryon conversion efficiencies inside dark matter halos. We find that the individual stellar masses and star formation rates can be as much as $\approx 4$ dex lower when fitting with an AGN component. These results highlight the importance of considering both stellar mass assembly and supermassive black hole growth when interpreting the light distributions of among the first galaxies to ever exist., Comment: 15 pages, 5 figures. Accepted for publication in ApJL after addressing reviewer comments

A wealth of observations have long suggested that the vast majority of isolated classical dwarf galaxies ($M_*=10^7$-$10^9$ M$_\odot$) are currently star-forming. However, recent observations of the large abundance of "Ultra-Diffuse Galaxies" beyond the reach of previous large spectroscopic surveys suggest that our understanding of the dwarf galaxy population may be incomplete. Here we report the serendipitous discovery of an isolated quiescent dwarf galaxy in the nearby Universe, which was imaged as part of the PEARLS GTO program. Remarkably, individual red-giant branch stars are visible in this near-IR imaging, suggesting a distance of $30\pm4$ Mpc, and a wealth of archival photometry point to an sSFR of $2\times10^{-11}$ yr$^{-1}$ and SFR of $4\times10^{-4}$ M$_\odot$ yr$^{-1}$. Spectra obtained with the Lowell Discovery Telescope find a recessional velocity consistent with the Hubble Flow and ${>}1500$ km/s separated from the nearest massive galaxy in SDSS, suggesting that this galaxy was either quenched from internal mechanisms or had a very high-velocity ($>1000$ km/s) interaction with a nearby massive galaxy in the past. This analysis highlights the possibility that many nearby quiescent dwarf galaxies are waiting to be discovered and that JWST has the potential to resolve them., Comment: Accepted to ApJ Letters

A Type Ia supernova (SN) at $z=1.78$ was discovered in James Webb Space Telescope Near Infrared Camera imaging of the galaxy cluster PLCK G165.7+67.0 (G165; $z = 0.35$). The SN is situated 1.5-2 kpc from the host-galaxy nucleus and appears in three different locations as a result of gravitational lensing by G165. These data can yield a value for Hubble's constant using time delays from this multiply-imaged SN Ia that we call "SN H0pe." Over the cluster, we identified 21 image multiplicities, confirmed five of them using the Near-Infrared Spectrograph, and constructed a new lens model that gives a total mass within 600 kpc of ($2.6 \pm 0.3) \times 10^{14}$ $M_{\odot}$. The photometry uncovered a galaxy overdensity coincident with the SN host galaxy. NIRSpec confirmed six member galaxies, four of which surround the SN host galaxy with relative velocity $\lesssim$900 km s$^{-1}$ and projected physical extent $\lesssim$33 kpc. This compact galaxy group is dominated by the SN host galaxy, which has a stellar mass of $(5.0 \pm 0.1) \times 10^{11}$ $M_{\odot}$. The group members have specific star-formation rates of 2-260 Gyr$^{-1}$ derived from the H$\alpha$-line fluxes corrected for stellar absorption, dust extinction, and slit losses. Another group centered on a strongly-lensed dusty star forming galaxy is at $z=2.24$. The total (unobscured and obscured) SFR of this second galaxy group is estimated to be ($\gtrsim$100 $M_{\odot}$ yr$^{-1}$), which translates to a supernova rate of $\sim$1 SNe yr$^{-1}$, suggesting that regular monitoring of this cluster may yield additional SNe., Comment: 29 pages, Accepted to ApJ on November 24, 2023

The presence of evolved stars in high-redshift galaxies can place valuable indirect constraints on the onset of star formation in the Universe. Thus we use PEARLS GTO and public NIRCam photometric data to search for Balmer-break candidate galaxies at $7 < z < 12$. We find that our Balmer-break candidates at $z \sim 10.5$ tend to be older (115 Myr), have lower inferred [O III] + H$\beta$ emission line equivalent widths (120 \r{A}), have lower specific star formation rates (6 Gyr$^{-1}$) and redder UV slopes ($\beta = -1.8$) than our control sample of galaxies. However, these trends all become less strong at $z \sim 8$, where the F444W filter now probes the strong rest-frame optical emission lines, thus providing additional constraints on the current star formation activity of these galaxies. Indeed, the bursty nature of Epoch of Reionisation galaxies can lead to a disconnect between their current SED profiles and their more extended star-formation histories. We discuss how strong emission lines, the cumulative effect of weak emission lines, dusty continua and AGN can all contribute to the photometric excess seen in the rest-frame optical, thus mimicking the signature of a Balmer break. Additional medium-band imaging will thus be essential to more robustly identify Balmer-break galaxies. However, the Balmer break alone cannot serve as a definitive proxy for the stellar age of galaxies, being complexly dependent on the star-formation history. Ultimately, deep NIRSpec continuum spectroscopy and MIRI imaging will provide the strongest indirect constraints on the formation era of the first galaxies in the Universe, thereby revealing when cosmic dawn breaks., Comment: 25 pages, 16 figures, 2 tables. Updated to published version in MNRAS

Direct-collapse black holes (DCBHs) of mass $\sim 10^4$-$10^5 {M}_\odot$ that form in HI-cooling halos in the early Universe are promising progenitors of the $\gtrsim 10^9 {M}_\odot$ supermassive black holes that fuel observed $z \gtrsim 7$ quasars. Efficient accretion of the surrounding gas onto such DCBH seeds may render them sufficiently bright for detection with the JWST up to $z\approx 20$. Additionally, the very steep and red spectral slope predicted across the $\approx 1$-5 $\mu$m wavelength range of the JWST/NIRSpec instrument during their initial growth phase should make them photometrically identifiable up to very high redshifts. In this work, we present a search for such DCBH candidates across the 34 arcmin$^{2}$ in the first two spokes of the JWST cycle-1 PEARLS survey of the north ecliptic pole time-domain field covering eight NIRCam filters down to a maximum depth of $\sim$ 29 AB mag. We identify two objects with spectral energy distributions consistent with the Pacucci et al. (2016) DCBH models. However, we also note that even with data in eight NIRCam filters, objects of this type remain degenerate with dusty galaxies and obscured active galactic nuclei over a wide range of redshifts. Follow-up spectroscopy would be required to pin down the nature of these objects. Based on our sample of DCBH candidates and assumptions on the typical duration of the DCBH steep-slope state, we set a conservative upper limit of $\lesssim 5\times 10^{-4}$ comoving Mpc$^{-3}$ (cMpc$^{-3}$) on the comoving density of host halos capable of hosting DCBHs with spectral energy distributions similar to the Pacucci et al. (2016) models at $z\approx 6$-14., Comment: 10 pages, 6 figures, 1 table, accepted for publication in A&A

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 towards 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 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 twelve could not have been found by HST because they are 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 de-magnified supernova surface density, when monitored at a time cadence of a few months to a ~3--4 micron survey limit of AB ~ 28.5 mag, of ~0.5 per sq. arcmin integrated to z ~ 2. This survey depth is beyond the capability of HST but can be easily reached by JWST., Comment: ApJS accepted

We present an analysis of a sample of robust high redshift galaxies selected photometrically from the `blank' fields of the Prime Extragalactic Areas for Reionization Science (PEARLS) survey and Early Release Observations (ERO) data of the James Webb Space Telescope (JWST) with the aim of selecting candidate high redshift active galactic nuclei (AGN). Sources were identified from the parent sample using a threefold selection procedure, which includes spectral energy distribution (SED) fitting to identify sources that are best fitted by AGN SED templates, a further selection based on the relative performance of AGN and non-AGN models, and finally morphological fitting to identify compact sources of emission, resulting in a purity-oriented procedure. Using this procedure, we identify a sample of nine AGN candidates at $6.5 < z < 12$, from which we constrain their physical properties as well as measure a lower bound on the AGN fraction in this redshift range of $5 \pm 1$\%. As this is an extreme lower limit due to our focus on purity and our SEDs being calibrated for unobscured Type 1 AGN, this demonstrates that AGN are perhaps quite common at this early epoch. The rest-frame UV colors of our candidate objects suggest that these systems are potentially candidate obese black hole galaxies (OBG), or AGN with very little galaxy component. We also investigate emission from our sample sources from fields overlapping with Chandra and VLA surveys, allowing us to place X-ray and 3 GHz radio detection limits on our candidates. Of note is a $z = 11.9$ candidate source exhibiting an abrupt morphological shift in the reddest band as compared to the bluer bands, indicating a potential merger or an unusually strong outflow., Comment: Accepted by MNRAS, 12 pages, 11 figures, updated to the final accepted version

We present a multi-wavelength analysis using SMA, JCMT, NOEMA, JWST, HST, and SST of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A1489. 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 <2um, even with deep JWST/NIRCam observations. We investigate their stellar, ISM, and dynamical properties, including a pixel-level SED analysis to derive sub-kpc-resolution stellar-mass and Av maps. We find that 850.1 is one of the most massive and highly obscured, Av~5, galaxies known at z>4 with M*~10^11.8 Mo (likely forming at z>6), and 850.2 is one of the least massive and least obscured, Av~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 <2um, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar mass distribution, Re(mass)~1kpc, but its expected evolution to z~1.5 and then z~0 matches both the properties of massive, quiescent galaxies at z~1.5 and ultra-massive 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 which we link to the active dynamical environment it resides within., Comment: Submitted to ApJ, comments welcome!

We use the Galaxy and Mass Assembly (GAMA) and the Deep Extragalactic Visible Legacy Survey (DEVILS) observational data sets to calculate the cosmic star formation rate (SFR) and active galactic nuclei (AGN) bolometric luminosity history (CSFH/CAGNH) over the last 12.5 billion years. SFRs and AGN bolometric luminosities were derived using the spectral energy distribution fitting code ProSpect, which includes an AGN prescription to self consistently model the contribution from both AGN and stellar emission to the observed rest-frame ultra-violet to far-infrared photometry. We find that both the CSFH and CAGNH evolve similarly, rising in the early Universe up to a peak at look-back time $\approx 10$~Gyr ($z \approx 2$), before declining toward the present day. The key result of this work is that we find the ratio of CAGNH to CSFH has been flat ($\approx 10^{42.5}\mathrm{erg \, s^{-1}M_{\odot}^{-1}yr}$) for $11$~Gyr up to the present day, indicating that star formation and AGN activity have been coeval over this time period. We find that the stellar masses of the galaxies that contribute most to the CSFH and CAGNH are similar, implying a common cause, which is likely gas inflow. The depletion of the gas supply suppresses cosmic star formation and AGN activity equivalently to ensure that they have experienced similar declines over the last 10 Gyr. These results are an important milestone for reconciling the role of star formation and AGN activity in the life cycle of galaxies., Comment: 16 pages, 10 figures. Figures 9 and 10 are the main results. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We identify 71 distant stars in JWST/NIRCam ERO images of the field of galaxy cluster SMACS J0723.3-7327 (SMACS 0723). Given the relatively small ($\sim$$10^{\circ}$) angular separation between SMACS 0723 and the Large Magellanic Cloud, it is likely that these stars are associated with the LMC outskirts or Leading Arm. This is further bolstered by a spectral energy distribution analysis, which suggests an excess of stars at a physical distance of $40-100$ kpc, consistent with being associated with or located behind the Magellanic system. In particular, we find that the overall surface density of stars brighter than 27.0 mag in the field of SMACS 0723 is $\sim$2.3 times that of stars in a blank field with similar galactic latitude (the North Ecliptic Pole Time Domain Field), and that the density of stars in the SMACS 0723 field with SED-derived distances consistent with the Magellanic system is $\sim$6.1 times larger than that of the blank field. The candidate stars at these distances are consistent with a stellar population at the same distance modulus with [Fe/H] $= -1.0$ and an age of $\sim$$5.0$ Gyr. On the assumption that all of the 71 stars are associated with the LMC, then the stellar density of the LMC at the location of the SMACS 0723 field is $\sim$$740$ stars kpc$^{-3}$, which helps trace the density of stars in the LMC outskirts., Comment: Published in ApJ

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 Science (PEARLS) survey combined with Early Release Science (ERS) data from the CEERS, GLASS, 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 characterisation of high-redshift ($z>6.5$) galaxies from a consistently processed set of deep, near-infrared imaging. Within an area of 180 arcmin$^{2}$, we identify 1046 candidate galaxies at redshifts $z>6.5$ and we use this sample to study the ultraviolet luminosity function (UV LF) in four redshift bins between $7.5

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 proto-clusters; 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, AGN growth, and First Light. Five fields, the JWST NEP Time-Domain Field (TDF), IRAC Dark Field (IDF), 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 $\mu$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 {\mu}m. PEARLS is designed to be of lasting benefit to the community., Comment: Accepted to AJ, comments welcome. We ask anyone who uses our public PEARLS (NEP TDF) data to refer to this overview paper

Using the first epoch of four-band NIRCam observations obtained by the James Webb Space Telescope (JWST) Prime Extragalactic Areas for Reionization and Lensing Science Program in the Spitzer IRAC Dark Field, we search for F150W and F200W dropouts. In 14.2 arcmin^2, we have found eight F150W dropouts and eight F200W dropouts, all brighter than 27.5 mag (the brightest being ~24 mag) in the band to the red side of the break. As they are detected in multiple bands, these must be real objects. Their nature, however, is unclear, and characterizing their properties is important for realizing the full potential of JWST. If the observed color decrements are due to the Lyman break, these objects should be at z >~ 11.7 and z >~ 15.4, respectively. The color diagnostics show that at least four F150W dropouts are far away from the usual contaminators encountered in dropout searches (red galaxies at much lower redshifts or brown dwarf stars). While the diagnostics of the F200W dropouts are less certain due to the limited number of passbands, at least one of them is likely not a known type of contaminant, and the rest are consistent with either high-redshift galaxies with evolved stellar populations or old galaxies at z ~ 3 to 8. If a significant fraction of our dropouts are indeed at z ~ 12, we have to face the severe problem of explaining their high luminosities and number densities. Spectroscopic identifications of such objects are urgently needed., Comment: Accepted for publication in ApJL

We use two independent, galaxy formation simulations, FLARES, a cosmological hydrodynamical simulation, and SHARK, a semi-analytic model, to explore how well the James Webb Space Telescope (JWST) will be able to uncover the existence and parameters of the star-forming main sequence (SFS) at $z=5\to10$, i.e. shape, scatter, normalisation. Using two independent simulations allows us to isolate predictions (e.g., stellar mass, star formation rate, SFR, luminosity functions) that are robust to or highly dependent on the implementation of the physics of galaxy formation. Both simulations predict that JWST can observe $\ge 70-90\%$ (for SHARK and FLARES respectively) of galaxies up to $z\sim10$ (down to stellar masses of $\approx 10^{8.3}\,\rm M_{\odot}$ and SFRs of $\approx 10^{0.5}\,\rm M_{\odot}\, yr^{-1}$) in modest integration times and given current proposed survey areas (e.g. the Web COSMOS $0.6\,\rm deg^2$) to accurately constrain the parameters of the SFS. Although both simulations predict qualitatively similar distributions of stellar mass and SFR, there are important quantitative differences, such as the abundance of massive, star-forming galaxies, with FLARES predicting a higher abundance than SHARK; the early onset of quenching as a result of black hole growth in FLARES (at $z\approx 8$), not seen in SHARK until much lower redshifts; and the implementation of synthetic photometry, with FLARES predicting more JWST-detected galaxies ($\sim 90\%$) than SHARK ($\sim 70\%$) at $z=10$. JWST observations will distinguish between these models, leading to a significant improvement upon our understanding of the formation of the very first galaxies., Comment: 21 pages, 13 figures including 1 in appendix, accepted to MNRAS