Search

Abstract Extensive air showers induced from high-energy cosmic rays provide a window into understanding the most energetic phenomena in the universe. We present a new method for observing these showers using the silicon imaging detector Subaru Hyper Suprime-Cam (HSC). This method has the advantage of being able to measure individual secondary particles. When paired with a surface detector array, silicon imaging detectors like Subaru HSC will be useful for studying the properties of extensive air showers in detail. The following report outlines the first results of observing extensive air showers with Subaru HSC. The potential for reconstructing the incident direction of primary cosmic rays is demonstrated and possible interdisciplinary applications are discussed.

Cosmic hydrogen reionization and cosmic production of the first metals are major phase transitions of the Universe occurring during the first billion years after the Big Bang; however, these are still underexplored observationally. Using the JWST/NIRSpec prism spectroscopy, we report the discovery of a sub- L _* galaxy at z _spec = 8.1623 ± 0.0007, dubbed RX J2129–z8He II , via the detection of a series of strong rest-frame UV/optical nebular emission lines and the clear Lyman break. RX J2129–z8He II shows a pronounced UV continuum with an extremely steep (i.e., blue) spectral slope of $\beta =-{2.53}_{-0.07}^{+0.06}$ , the steepest among all spectroscopically confirmed galaxies at z _spec ≳ 7, in support of its very hard ionizing spectrum that could lead to a significant leakage of its ionizing flux. Therefore, RX J2129–z8He II is representative of the key galaxy population driving the cosmic reionization. More importantly, we detect a strong He II λ 1640 emission line in its spectrum, one of the highest redshifts at which such a line is robustly detected. Its high rest-frame equivalent width (EW = 21 ± 4 Å) and extreme flux ratios with respect to UV metal and Balmer lines raise the possibility that part of RX J2129–z8He II ’s stellar population could be Pop III (Pop III)-like. Through careful photoionization modeling, we show that the physically calibrated phenomenological models of the ionizing spectra of Pop III stars with strong mass loss can successfully reproduce the emission line flux ratios observed in RX J2129–z8He II . Assuming the Eddington limit, the total mass of the Pop III stars within this system is estimated to be 7.8 ± 1.4 × 10 ^5 M _⊙ . To date, this galaxy presents the most compelling case in the early Universe where trace Pop III stars might coexist with metal-enriched populations.

Stellar-mass primordial black holes (PBHs) from the early Universe can directly contribute to the gravitational wave (GW) events observed by LIGO, but can only comprise a subdominant component of the dark matter (DM). The primary DM constituent will generically form massive halos around seeding stellar-mass PBHs. We demonstrate that gravitational lensing of sources at cosmological (≳Gpc) distances can directly explore DM halo dresses engulfing PBHs, challenging for lensing of local sources in the vicinity of Milky Way. Strong lensing analysis of fast radio bursts detected by CHIME survey already starts to probe parameter space of dressed stellar-mass PBHs, and upcoming searches can efficiently explore dressed PBHs over ∼10−105M⊙ mass-range and provide a stringent test of the PBH scenario for the GW events. Our findings establish a general test for a broad class of DM models with stellar-mass PBHs, including those where QCD axions or WIMP-like particles comprise predominant DM. The results open a new route for exploring dressed PBHs with various types of lensing events at cosmological distances, such as supernovae and caustic crossings.

The luminosities ( L ) and velocity dispersions ( σ ) of the extinction-corrected Balmer emission lines of giant H ii regions in nearby galaxies exhibit a tight correlation (∼0.35 dex scatter). There are few constraints, however, on whether giant H ii regions at significant look-back times follow an L – σ relation, given the angular resolution and sensitivity required to study them individually. We measure the luminosities and velocity dispersions of H α and H β emission from 11 H ii regions in Sp1149, a spiral galaxy at redshift ( z ) z = 1.49 multiply imaged by the MACS J1149 galaxy cluster. Sp1149 is also the host galaxy of the first-known strongly lensed supernova with resolved images, SN Refsdal. We employ archival Keck I OSIRIS observations, and newly acquired Keck I MOSFIRE and Large Binocular Telescope LUCI long-slit spectra of Sp1149. When we use the GLAFIC simply parameterized lens model, we find that the H α luminosities of the H ii regions at z = 1.49 are a factor of ${6.4}_{-2.0}^{+2.9}$ brighter than predicted by the low-redshift L – σ relation we measure from Very Large Telescope MUSE spectroscopy. If the lens model is accurate, then the H ii regions in Sp1149 differ from their low-redshift counterparts. We identify an H ii region in Sp1149 that is dramatically brighter (by 2.03 ± 0.44 dex) than our low-redshift L – σ relation predicts given its low velocity dispersion. Finally, the H ii regions in Sp1149 are consistent, perhaps surprisingly, with the z ≈ 0 star-forming locus on the Baldwin–Phillips–Terlevich diagram.

Galaxy-cluster gravitational lenses enable the study of faint galaxies even at large lookback times, and, recently, time-delay constraints on the Hubble constant. There have been few tests, however, of lens model predictions adjacent to the critical curve (≲8″) where the magnification is greatest. In a companion paper, we use the GLAFIC lens model to constrain the Balmer L – σ relation for H ii regions in a galaxy at redshift z = 1.49 strongly lensed by the MACS J1149 galaxy cluster. Here we perform a detailed comparison between the predictions of 10 cluster lens models that employ multiple modeling assumptions with our measurements of 11 magnified, giant H ii regions. We find that that the models predict magnifications an average factor of 6.2 smaller, a ∼2 σ tension, than that inferred from the H ii regions under the assumption that they follow the low-redshift L – σ relation. To evaluate the possibility that the lens model magnifications are strongly biased, we next consider the flux ratios among knots in three images of Sp1149, and find that these are consistent with model predictions. Moreover, while the mass-sheet degeneracy could in principle account for a factor of ∼6 discrepancy in magnification, the value of H _0 inferred from SN Refsdal’s time delay would become implausibly small. We conclude that the lens models are not likely to be highly biased, and that instead the H ii regions in Sp1149 are substantially more luminous than the low-redshift Balmer L – σ relation predicts.

Herein, we present the statistical properties of active galactic nuclei (AGNs) for approximately 1 million member galaxies of galaxy groups and clusters with 0.1 < cluster redshift ( z _cl ) < 1.4 selected using the Subaru Hyper Suprime-Cam, the so-called CAMIRA clusters. In this research, we focused on the AGN power fraction ( f _AGN ), which is defined as the proportion of the contribution of AGNs to the total infrared (IR) luminosity, L _IR (AGN)/ L _IR , and examined how f _AGN depends on (i) z _cl and (ii) the distance from the cluster center. We compiled multiwavelength data using the ultraviolet–mid-IR range. Moreover, we performed spectral energy distribution fits to determine f _AGN using the CIGALE code with the SKIRTOR AGN model. We found that (i) the value of f _AGN in the CAMIRA clusters is positively correlated with z _cl , with the correlation slope being steeper than that for field galaxies, and (ii) f _AGN exhibits a high value at the cluster outskirts. These results indicate that the emergence of the AGN population depends on the redshift and environment and that galaxy groups and clusters at high redshifts are important in AGN evolution. Additionally, we demonstrated that cluster–cluster mergers may enhance AGN activity at the outskirts of particularly massive galaxy clusters. Our findings are consistent with a related study on the CAMIRA clusters that was based on the AGN number fraction.

Sub/millimeter galaxies are a key population for the study of galaxy evolution because the majority of star formation at high redshifts occurred in galaxies deeply embedded in dust. To search for this population, we have performed an extensive survey with Atacama Large Millimeter/submillimeter Array (ALMA), called the ALMA Lensing Cluster Survey (ALCS). This survey covers 133 arcmin ^2 area and securely detects 180 sources at z ∼ 0.5–6 with a flux limit of ∼0.2 mJy at 1.2 mm. Here, we report the results of multiwavelength spectral energy distribution analysis of the whole ALCS sample, utilizing the observed-frame UV to millimeter photometry. We find that the majority of the ALCS sources lie on the star-forming main sequence, with a smaller fraction showing intense starburst activities. The ALCS sample contains high infrared-excess sources ( $\mathrm{IRX}=\mathrm{log}({L}_{\mathrm{dust}}/{L}_{\mathrm{UV}})\gt 1$ ), including two extremely dust-obscured galaxies (IRX > 5). We also confirm that the ALCS sample probes a broader range in lower dust mass than conventional submillimeter galaxy samples in the same redshift range. We identify six heavily obscured active galactic nucleus (AGN) candidates that are not detected in the archival Chandra data in addition to the three X-ray AGNs reported by Uematsu et al. (2023). The inferred AGN luminosity density shows a possible excess at z = 2–3 compared with that determined from X-ray surveys below 10 keV.

We present Atacama Large Millimeter/submillimeter Array (ALMA) deep spectroscopy for a lensed galaxy at z _spec = 8.496 with $\mathrm{log}({M}_{\mathrm{star}}/{M}_{\odot })\sim 7.8$ whose optical nebular lines and stellar continuum are detected by JWST/NIRSpec and NIRCam Early Release Observations in the field of SMACS J0723.3–7327. Our ALMA spectrum shows [O iii ] 88 μ m and [C ii ] 158 μ m line detections at 4.0 σ and 4.5 σ , respectively. The redshift and position of the [O iii ] line coincide with those of the JWST source, while the [C ii ] line is blueshifted by 90 km s ^−1 with a spatial offset of 0.″5 (≈0.5 kpc in the source plane) from the centroid of the JWST source. The NIRCam F444W image, including [O iii ] λ 5007 and H β line emission, spatially extends beyond the stellar components by a factor of >8. This indicates that the z = 8.5 galaxy has already experienced strong outflows as traced by extended [O iii ] λ 5007 and offset [C ii ] emission, which would promote ionizing photon escape and facilitate reionization. With careful slit-loss corrections and the removal of emission spatially outside the galaxy, we evaluate the [O iii ] 88 μ m/ λ 5007 line ratio, and derive the electron density n _e by photoionization modeling to be ${220}_{-130}^{+230}$ cm ^−3 , which is comparable with those of z ∼ 2–3 galaxies. We estimate an [O iii ] 88 μ m/[C ii ] 158 μ m line ratio in the galaxy of >4, as high as those of known z ∼ 6–9 galaxies. This high [O iii ] 88 μ m/[C ii ] 158 μ m line ratio is generally explained by the high n _e as well as the low metallicity ( ${Z}_{\mathrm{gas}}/{Z}_{\odot }={0.04}_{-0.02}^{+0.02}$ ), high ionization parameter ( $\mathrm{log}U\gt -2.27$ ), and low carbon-to-oxygen abundance ratio (log(C/O) = [−0.52: −0.24]) obtained from the JWST/NIRSpec data; further [C ii ] follow-up observations will constrain the covering fraction of photodissociation regions.

We report Hubble Space Telescope (HST) Wide Field Camera 3 deep IR (F160W) imaging of SDSS J1608+2716. This system, located at a redshift of z = 2.575, was recently reported as a triple-quasar candidate with subarcsecond separations (∼0.″25) based on selection from Gaia astrometry and follow-up Keck adaptive optics–assisted integral field unit spectroscopy. Our new HST deep IR imaging reveals the presence of a fourth point-like component located ∼0.″9 away from the triple system. Additionally, we detect an edge-on disk galaxy located in between the four point sources. The entire system exhibits a characteristic cusp structure in the context of strong gravitational lensing, and the observed image configuration can be successfully reproduced using a lens model based on a singular isothermal ellipsoid mass profile. These findings indicate that this system is a quadruply lensed quasar. Our results highlight the challenges associated with identifying dual/multiple quasars on ∼kiloparsec scales at high redshifts and emphasize the crucial role of deep, high-resolution IR imaging in robustly confirming such systems.

We present the quasar luminosity function (LF) at z = 7, measured with 35 spectroscopically confirmed quasars at 6.55 < z < 7.15. The sample of 22 quasars from the Subaru High- z Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over −28 < M _1450 < −23. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude ${M}_{1450}^{* }=-{25.60}_{-0.30}^{+0.40}$ , a characteristic density ${{\rm{\Phi }}}^{* }={1.35}_{-0.30}^{+0.47}$ Gpc ^−3 mag ^−1 , and a bright-end slope $\beta =-{3.34}_{-0.57}^{+0.49}$ , when the faint-end slope is fixed to α = −1.2 as observed at z ≤ 6. The overall LF shape remains remarkably similar from z = 4 to 7, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at z ≥ 6. The estimated ionizing photon density, 10 ^48.2±0.1 s ^−1 Mpc ^−3 , is less than 1% of the critical rate to keep the intergalactic medium ionized at z = 7, and thus indicates that quasars are not a major contributor to cosmic reionization.

MACS0647–JD is a triply lensed z ∼ 11 galaxy originally discovered with the Hubble Space Telescope. The three lensed images are magnified by factors of ∼8, 5, and 2 to AB mag 25.1, 25.6, and 26.6 at 3.5 μ m. The brightest is over a magnitude brighter than other galaxies recently discovered at similar redshifts z > 10 with JWST. Here, we report new JWST imaging that clearly resolves MACS0647–JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. The brighter larger component “A” is intrinsically very blue ( β ∼ −2.6 ± 0.1), likely due to very recent star formation and no dust, and is spatially extended with an effective radius ∼70 ± 24 pc. The smaller component “B” ( r ∼ 20 ${}_{-5}^{+8}\,$ pc) appears redder ( β ∼ −2 ± 0.2), likely because it is older (100–200 Myr) with mild dust extinction ( A _V ∼ 0.1 mag). With an estimated stellar mass ratio of roughly 2:1 and physical projected separation ∼400 pc, we may be witnessing a galaxy merger 430 million years after the Big Bang. We identify galaxies with similar colors in a high-redshift simulation, finding their star formation histories to be dissimilar, which is also suggested by the spectral energy distribution fitting, suggesting they formed further apart. We also identify a candidate companion galaxy “C” ∼3 kpc away, likely destined to merge with A and B. Upcoming JWST Near Infrared Spectrograph observations planned for 2023 January will deliver spectroscopic redshifts and more physical properties for these tiny magnified distant galaxies observed in the early universe.

We report the discovery of two extremely magnified lensed star candidates behind the galaxy cluster MACS J0647.7+015 using recent multiband James Webb Space Telescope (JWST) NIRCam observations. The star candidates are seen in a previously known, z _phot ≃ 4.8 dropout giant arc that straddles the critical curve. The candidates lie near the expected critical curve position, but lack clear counter-images on the other side of it, suggesting these are possibly stars undergoing caustic crossings. We present revised lensing models for the cluster, including multiply imaged galaxies newly identified in the JWST data, and use them to estimate background macro-magnifications of at least ≳90 and ≳50 at the positions of the two candidates, respectively. With these values, we expect effective, caustic-crossing magnifications of ∼[10 ^3 –10 ^5 ] for the two star candidates. The spectral energy distributions of the two candidates match well the spectra of B-type stars with best-fit surface temperatures of ∼10,000 K, and ∼12,000 K, respectively, and we show that such stars with masses ≳20 M _⊙ and ≳50 M _⊙ , respectively, can become sufficiently magnified to be observable. We briefly discuss other alternative explanations and conclude that these objects are likely lensed stars, but also acknowledge that the less-magnified candidate may alternatively reside in a star cluster. These star candidates constitute the second highest-redshift examples to date after Earendel at z _phot ≃ 6.2, establishing further the potential of studying extremely magnified stars at high redshifts with JWST. Planned future observations, including with NIRSpec, will enable a more detailed view of these candidates in the near future.

We present the evolution of the mass–metallicity (MZ) relation at z = 4–10 derived with 135 galaxies identified in JWST/NIRSpec data taken from the three major public spectroscopy programs of ERO, GLASS, and CEERS. Because there are many discrepancies between the flux measurements reported by the early ERO studies, we first establish our NIRSpec data reduction procedure for reliable emission-line flux measurements and errors, successfully explaining Balmer decrements with no statistical tensions thorough comparisons with the early ERO studies. Applying the reduction procedure to the 135 galaxies, we obtain emission-line fluxes for physical property measurements. We confirm that 10 out of the 135 galaxies with [O iii ] λ 4363 lines have electron temperatures of ≃(1.1–2.3) × 10 ^4 K, similar to lower- z star-forming galaxies, which can be explained by heating by young massive stars. We derive the metallicities of the 10 galaxies by a direct method and the rest of the galaxies with strong lines using the metallicity calibrations of Nakajima et al. applicable for these low-mass metal-poor galaxies, anchoring the metallicities with the direct-method measurements. We thus obtain the MZ relation and star formation rate (SFR)–MZ relation over z = 4–10. We find that there is a small evolution of the MZ relation from z ∼ 2–3 to z = 4–10, while interestingly the SFR–MZ relation shows no evolution up to z ∼ 8 but a significant decrease at z > 8 beyond the errors This SFR–MZ relation decrease at z > 8 may suggest a break of the metallicity equilibrium state via star formation, inflow, and outflow, while further statistical and local-baseline studies are needed for a conclusion.

This study examines the impact of cluster environments on galaxy properties using data from the Hyper Suprime-Cam Subaru Strategic Program and an optically selected CAMIRA cluster sample. Specifically, the study analyzes the fractions of quiescent and green valley galaxies with stellar masses above 10 ^8.6 M _⊙ at z ∼ 0.2 and 10 ^9.8 M _⊙ at z ∼ 1.1, investigating their trends in radius and density. The results indicate that a slow quenching mechanism is at work, as evidenced by a radially independent specific star formation rate reduction of 0.1 dex for star-forming galaxies in a cluster environment. The study also finds that slow quenching dominates fast quenching only for low-mass galaxies (

A tight positive correlation between the stellar mass and the gas-phase metallicity of galaxies has been observed at low redshifts. The redshift evolution of this correlation can strongly constrain theories of galaxy evolution. The advent of JWST allows probing the mass–metallicity relation at redshifts far beyond what was previously accessible. Here we report the discovery of two emission line galaxies at redshifts 8.15 and 8.16 in JWST NIRCam imaging and NIRSpec spectroscopy of targets gravitationally lensed by the cluster RX J2129.4+0005. We measure their metallicities and stellar masses along with nine additional galaxies at 7.2 < z _spec < 9.5 to report the first quantitative statistical inference of the mass–metallicity relation at z ≈ 8. We measure ∼0.9 dex evolution in the normalization of the mass–metallicity relation from z ≈ 8 to the local universe; at a fixed stellar mass, galaxies are 8 times less metal enriched at z ≈ 8 compared to the present day. Our inferred normalization is in agreement with the predictions of FIRE simulations. Our inferred slope of the mass–metallicity relation is similar to or slightly shallower than that predicted by FIRE or observed at lower redshifts. We compare the z ≈ 8 galaxies to extremely low-metallicity analog candidates in the local universe, finding that they are generally distinct from extreme emission line galaxies or “green peas,” but are similar in strong emission line ratios and metallicities to “blueberry galaxies.” Despite this similarity, at a fixed stellar mass, the z ≈ 8 galaxies have systematically lower metallicities compared to blueberry galaxies.

Dual quasars at small physical separations are an important precursor phase of galaxy mergers, ultimately leading to the coalescence of the two supermassive black holes. Starting from a sample of dual and/or lensed quasar candidates discovered using astrometric jitter in Gaia data, we present a pilot case study of one of the most promising yet puzzling candidate dual quasars at cosmic noon ( z ∼ 1.8). Using multiwavelength imaging and spectroscopy from X-ray to radio, we test whether the SDSS J0823+2418 system is two individual quasars in a bound pair at separation ∼0.″64, or instead a single quasar being gravitationally lensed by a foreground galaxy. We find consistent flux ratios (∼1.25−1.45) between the two sources in optical, near-IR (NIR), UV, and radio, and thus similar spectral energy distributions, suggesting a strong-lensing scenario. However, differences in the radio spectral index, as well as changing X-ray fluxes, hint at either a dual quasar with otherwise nearly identical properties or perhaps lensing-based time lag of ∼3 days paired with intrinsic variability. We find with lens mass modeling that the relative NIR positions and magnitudes of the two quasars and a marginally detected central galaxy are consistent with strong lensing. Archival Sloan Digital Sky Survey spectra likewise suggest a foreground absorber via Mg ii absorption lines. We conclude that SDSS J0823+2418 is likely a lensed quasar, and therefore that the VODKA sample contains a population of these lensed systems (perhaps as high as 50%) as well as dual quasars.

We report the discovery of four galaxy candidates observed 450–600 Myr after the Big Bang with photometric redshifts between z ∼ 8.3 and 10.2 measured using James Webb Space Telescope (JWST) NIRCam imaging of the galaxy cluster WHL0137−08 observed in eight filters spanning 0.8–5.0 μ m, plus nine Hubble Space Telescope filters spanning 0.4–1.7 μ m. One candidate is gravitationally lensed with a magnification of μ ∼ 8, while the other three are located in a nearby NIRCam module with expected magnifications of μ ≲ 1.1. Using SED fitting, we estimate the stellar masses of these galaxies are typically in the range $\mathrm{log}{M}_{\star }/{M}_{\odot }$ = 8.3–8.7. All appear young, with mass-weighted ages

In late 2014, four images of supernova (SN) “Refsdal,” the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images’ discovery, the SN was predicted to reappear within hundreds of days at a new position ∼8″ away in the field. The observed reappearance in late 2015 makes it possible to carry out Refsdal’s original proposal to use a multiply imaged SN to measure the Hubble constant H _0 , since the time delay between appearances should vary inversely with H _0 . Moreover, the position, brightness, and timing of the reappearance enable a novel test of the blind predictions of galaxy-cluster models, which are typically constrained only by the positions of multiply imaged galaxies. We have developed a new photometry pipeline that uses DOLPHOT to measure the fluxes of the five images of SN Refsdal from difference images. We apply four separate techniques to perform a blind measurement of the relative time delays and magnification ratios between the last image SX and the earlier images S1–S4. We measure the relative time delay of SX–S1 to be $\displaystyle {376.0}_{-5.5}^{+5.6}$ days and the relative magnification to be $\displaystyle {0.30}_{-0.3}^{+0.5}$ . This corresponds to a 1.5% precision on the time delay and 17% precision for the magnification ratios and includes uncertainties due to millilensing and microlensing. In an accompanying paper, we place initial and blind constraints on the value of the Hubble constant.

The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but they have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of ∼6.2 ± 0.5 × 10 ^−4 integrated over ∼0.″3–3″ separations (projected physical separations of ∼3–30 kpc at z ∼ 2) in luminous ( L _bol > 10 ^45.8 erg s ^−1 ) unobscured quasars at 1.5 < z < 3.5 using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The measurement was based on a sample of 60 Gaia-resolved double quasars (out of 487 Gaia pairs dominated by quasar+star superpositions) at these separations, corrected for pair completeness in Gaia, which we quantify as functions of pair separation, magnitude of the primary, and magnitude contrast. The double quasar fraction increases toward smaller separations by a factor of ∼5 over these scales. The division between physical quasar pairs and lensed quasars in our sample is currently unknown, requiring dedicated follow-up observations (in particular, deep, subarcsecond-resolution IR imaging for the closest pairs). Intriguingly, at this point, the observed pair statistics are in rough agreement with theoretical predictions both for the lensed quasar population in mock catalogs and for dual quasars in cosmological hydrodynamic simulations. Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST, and Roman, combined with targeted follow-up observations, will conclusively measure the abundances and host galaxy properties of galactic-scale quasar pairs, offset AGNs, and subarcsecond lensed quasars across cosmic time.

Star cluster formation in the early universe and its contribution to reionization remains largely unconstrained to date. Here we present JWST/NIRCam imaging of the most highly magnified galaxy known at z ∼ 6, the Sunrise arc. We identify six young massive star clusters (YMCs) with measured radii spanning from ∼20 down to ∼1 pc (corrected for lensing magnification), estimated stellar masses of ∼10 ^6–7 M _⊙ , and ages of 1–30 Myr based on SED fitting to photometry measured in eight filters extending to rest frame 7000 Å. The resulting stellar mass surface densities are higher than 1000 M _⊙ pc ^−2 (up to a few 10 ^5 M _⊙ pc ^−2 ), and their inferred dynamical ages qualify the majority of these systems as gravitationally bound stellar clusters. The star cluster ages map the progression of star formation along the arc, with two evolved systems (≳10 Myr old) followed by very young clusters. The youngest stellar clusters (1000 Å rest frame and are hosted in a 200 pc sized star-forming complex. Such a region dominates the ionizing photon production with a high efficiency $\mathrm{log}({\xi }_{\mathrm{ion}}[\mathrm{Hz}\,{\mathrm{erg}}^{-1}])\sim 25.7$ . A significant fraction of the recently formed stellar mass of the galaxy (10%–30%) occurred in these YMCs. We speculate that such sources of ionizing radiation boost the ionizing photon production efficiency, which eventually carves ionized channels that might favor the escape of Lyman continuum radiation. The survival of some of the clusters would make them the progenitors of massive and relatively metal-poor globular clusters in the local universe.

We study the spatially resolved stellar populations of 444 galaxies at 0.3 < z < 6.0 in two clusters (WHL 0137–08 and MACS 0647+70) and a blank field, combining imaging data from the Hubble Space Telescope and JWST to perform spatially resolved spectral energy distribution (SED) modeling using piXedfit . The high spatial resolution of the imaging data combined with magnification from gravitational lensing in the cluster fields allows us to resolve a large fraction of our galaxies (109) to subkiloparsec scales. At redshifts around cosmic noon and higher (2.5 ≲ z ≲ 6.0), we find mass-doubling times to be independent of radius, inferred from flat specific star formation rate (sSFR) radial profiles and similarities between the half-mass and half-SFR radii. At lower redshifts (1.5 ≲ z ≲ 2.5), a significant fraction of our star-forming galaxies shows evidence for nuclear starbursts, inferred from a centrally elevated sSFR and a much smaller half-SFR radius compared to the half-mass radius. At later epochs, we find more galaxies suppress star formation in their centers but are still actively forming stars in the disk. Overall, these trends point toward a picture of inside-out galaxy growth consistent with theoretical models and simulations. We also observe a tight relationship between the central mass surface density and global stellar mass with ∼0.38 dex scatter. Our analysis demonstrates the potential of spatially resolved SED analysis with JWST data. Future analysis with larger samples will be able to further explore the assembly of galaxy mass and the growth of their structures.

We conduct a comprehensive study on dropout galaxy candidates at z ∼ 9–16 using the first 90 arcmin ^2 James Webb Space Telescope (JWST) Near Infrared Camera images taken by the early release observations (ERO) and early release science programs. With the JWST simulation images, we find that a number of foreground interlopers are selected with a weak photo- z determination (Δ χ ^2 > 4). We thus carefully apply a secure photo -z selection criterion (Δ χ ^2 > 9) and conventional color criteria with confirmations of the ERO Near Infrared Spectrograph spectroscopic redshifts, and obtain a total of 23 dropout galaxies at z ∼ 9–16, including two candidates at ${z}_{\mathrm{phot}}={16.25}_{-0.46}^{+0.24}$ and ${16.41}_{-0.55}^{+0.66}$ . We perform thorough comparisons of dropout galaxies found in our work with recent JWST studies, and conclude that our galaxy sample is reliable enough for statistical analyses. We derive the UV luminosity functions at z ∼ 9–16, and confirm that our UV luminosity functions at z ∼ 9 and 12 agree with those determined by other Hubble Space Telescope and JWST studies. The cosmic star formation rate (SFR) density decreases from z ∼ 9 to 12, and perhaps to 16, but the densities at z ∼ 12–16 are higher than the constant star formation efficiency model. Interestingly, there are six bright galaxy candidates at z ∼ 10–16 with M _UV < −19.5 mag and M _* ∼ 10 ^8−9 M _⊙ . Because a majority (∼80%) of these galaxies show no signatures of active galactic nuclei in their morphologies, the high cosmic SFR densities and the existence of these UV-luminous galaxies are explained by the lack of suppression of star formation by the UV background radiation at the pre-reionization epoch and/or an efficient UV radiation production by a top-heavy initial mass function with Population III–like star formation.