大栗 真宗

The Cosmic Gems arc is among the brightest and highly magnified galaxies observed at redshift z ≈ 10.2 (ref. 1). However, it is an intrinsically ultraviolet faint galaxy, in the range of those now thought to drive the reionization of the Universe2–4. Hitherto the smallest features resolved in a galaxy at a comparable redshift are between a few hundreds and a few tens of parsecs (pc)5,6. Here we report JWST observations of the Cosmic Gems. The light of the galaxy is resolved into five star clusters located in a region smaller than 70 pc. They exhibit minimal dust attenuation and low metallicity, ages younger than 50 Myr and intrinsic masses of about 106M⊙. Their lensing-corrected sizes are approximately 1 pc, resulting in stellar surface densities near 105M⊙ pc−2, three orders of magnitude higher than typical young star clusters in the local Universe7. Despite the uncertainties inherent to the lensing model, they are consistent with being gravitationally bound stellar systems, that is, proto-globular clusters. We conclude that star cluster formation and feedback likely contributed to shaping the properties of galaxies during the epoch of reionization.

Microlensing near macro-caustics is a complex phenomenon in which swarms of micro-images produced by micro-caustics form on both sides of a macro-critical curve. Recent discoveries of highly magnified images of individual stars in massive galaxy cluster lenses, predicted to be formed by these micro-image swarms, have stimulated studies on this topic. In this article, we explore microlensing near macro-caustics using both simulations and analytic calculations. We show that the mean total magnification of the micro-image swarms follows that of an extended source in the absence of microlensing. Micro-caustics join into a connected network in a region around the macro-critical line of a width proportional to the surface density of microlenses; within this region, the increase of the mean magnification toward the macro-caustic is driven by the increase of the number of micro-images rather than individual magnifications of micro-images. The maximum achievable magnification in micro-caustic crossings decreases with the mass fraction in microlenses. We conclude with a review of applications of this microlensing phenomenon, including limits to the fraction of dark matter in compact objects, and searches of Population III stars and dark matter subhalos. We argue that the discovered highly magnified stars at cosmological distances already imply that less than ∼ 10% of the dark matter may be in the form of compact objects with mass above ∼10−6M⊙.

We present the size-mass relation of the brightest cluster galaxies (BCGs) at 0.1 < z < 1.4 using the imaging data obtained by the Hyper Suprime-Cam Subaru Strategic Program. Our sample consists of 471 photometrically selected BCGs with stellar mass logM∗/M⊙ = 11-12. We measure the size of the BCGs using i-band imaging and model fits based on a single SCrossed D sign.

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.

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 β = − 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 × 105 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.

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_{\rm cl}$) $<$ 1.4, selected using Subaru Hyper Suprime-Cam, the so-called CAMIRA clusters. In this research, we focused on the AGN power fraction ($f_{\rm AGN}$), which is defined as the proportion of the contribution of AGNs to the total infrared (IR) luminosity, $L_{\rm IR}$ (AGN)/$L_{\rm IR}$, and examined how $f_{\rm AGN}$ depends on (i) $z_{\rm 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_{\rm AGN}$ using the CIGALE code with the SKIRTOR AGN model. We found that (i) the value of $f_{\rm AGN}$ in the CAMIRA clusters is positively correlated with $z_{\rm cl}$, with the correlation slope being steeper than that for field galaxies, and (ii) $f_{\rm AGN}$ exhibits a high value at the cluster outskirts. These results indicate that the emergence of 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.

Astrometric perturbations of critical curves in strong lens systems are thought to be one of the most promising probes of substructures down to small-mass scales. While a smooth mass distribution creates a symmetric geometry of critical curves with radii of curvature about the Einstein radius, substructures introduce small-scale distortions on critical curves, which can break the symmetry of gravitational lensing events near critical curves, such as highly magnified individual stars. We derive a general formula that connects the fluctuation of critical curves with the fluctuation of the surface density caused by substructures, which is useful when constraining models of substructures from observed astrometric perturbations of critical curves. We numerically check that the formula is valid and accurate as long as substructures are not dominated by a small number of massive structures. As a demonstration of the formula, we also explore the possibility that an anomalous position of an extremely magnified star, recently reported as "Mothra,"can be explained by fluctuations in the critical curve due to substructures. We find that cold dark matter subhalos with masses ranging from 5×107M⊙/h to 109M⊙/h can well explain the anomalous position of Mothra, while in the fuzzy dark matter model, the very small mass of ∼10-24 eV is needed to explain it.

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 arcmin2 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 ( IRX = log ( L dust / L UV ) > 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 log ( M star / M ⊙ ) ∼ 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 gas / Z ⊙ = 0.04 − 0.02 + 0.02 ), high ionization parameter ( log U > − 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 perform the weak lensing mass mapping analysis to identify troughs, which are defined as local minima in the mass map. Since weak lensing probes the projected matter distribution along the line of sight, these troughs can be produced by single or multiple voids projected along the line of sight. To scrutinize the origins of the weak lensing troughs, we systematically investigate the line-of-sight structure of troughs selected from the latest Subaru Hyper Suprime-Cam (HSC) Year 3 weak lensing data covering 433.48 deg2. From a curved sky mass map constructed with the HSC data, we identify 15 troughs with the signal-to-noise ratio higher than 5.7 and address their line-of-sight density structure utilizing redshift distributions of two galaxy samples, photometric luminous red galaxies observed by HSC and spectroscopic galaxies detected by Baryon Oscillation Spectroscopic Survey. While most weak lensing signals due to the troughs are explained by multiple voids aligned along the line of sight, we find that two of the 15 troughs potentially originate from single voids at redshift ∼0.3. The single void interpretation appears to be consistent with our three-dimensional mass mapping analysis. We argue that single voids can indeed reproduce observed weak lensing signals at the troughs if these voids are not spherical but are highly elongated along the line-of-sight direction.

We report the discovery of new lensed quasar candidates in the imaging data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Data Release 4, covering of the sky with median seeing in the i band of ≈0.6 arcsec. In addition to two catalogues of Million Quasars Catalog v7.2 and AllWISE Catalog of Mid-Infra Red Active Galactic Nuclei, containing confirmed and candidate quasars, we preselect quasar sources using colour cuts from the HSC (grizy) and unWISE (W1 + W2) photometric data based on SDSS spectroscopic catalogues. We search for the presence of multiple point sources with similar colour through the convolution of the Laplacian of the preselected quasar image cutouts with the Laplacian of the point spread function, resulting in a reduction of lens candidates from 1 652 329 to 121 511 (7.4 per cent). After visual binary classification, we grade 6 199 (0.4 per cent) potential lenses on a scale of 0 to 3, with 3 indicating a lens and 0 indicating a non-lens. Finally we obtain 162 lens candidates with an average grade of ≥2, and among them, we successfully recover 18 known lenses. By fitting the light distribution and removing the known contaminants, we discover that 57 new systems contain at least two point sources and a galaxy in between, including 10 possible quadruply lensed quasars. This new sample exhibits a median separation of 1.26 arcsec and a magnitude limit of i ≈ 22. Spectroscopic or high-resolution imaging follow up on these newly discovered lensed quasar candidates will further allow their natures to be confirmed.

We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of zphot ≃ 4.8, previously detected in JWST/NIRCam imaging of the strong lensing (SL) cluster MACS J0647+7015 (z = 0.591). The spectroscopic observation allows us to precisely measure the redshift of the host arc at zspec = 4.758 ± 0.004, and the star’s spectrum displays clear Lyman- and Balmer-breaks commensurate with this redshift. A fit to the spectrum suggests a B-type super-giant star of surface temperature Teff,B ≃ 15 000 K with either a redder F-type companion (Teff,F ≃ 6 250 K) or significant dust attenuation (AV ≃ 0.82) along the line of sight. We also investigate the possibility that this object is a magnified young globular cluster rather than a single star. We show that the spectrum is in principle consistent with a star cluster, which could also accommodate the lack of flux variability between the two epochs. However, the lack of a counter image and the strong upper limit on the size of the object from lensing symmetry, r ≲ 0.5 pc, could indicate that this scenario is somewhat less likely – albeit not completely ruled out by the current data. The presented spectrum seen at a time when the Universe was only ∼1.2 Gyr old showcases the ability of JWST to study early stars through extreme lensing.

We present cosmology results from a blinded joint analysis of cosmic shear, ζ±, galaxy-galaxy weak lensing, Δς(R), and projected galaxy clustering, wp(R), measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog - a 3×2 pt cosmology analysis. We define luminosity-cut, and therefore nearly volume-limited, samples of SDSS galaxies to serve as the tracers of wp and as the lens samples for Δς in three spectroscopic redshift bins spanning the range 0.15<z<0.7. For the ζ± and Δς measurements, we use a single sample of about seven million source galaxies over 416 deg2, selected from HSC-Y3 based on having photometric redshifts (photo-z) greater than 0.75. The deep, high-quality HSC-Y3 data enable significant detections of the Δς signals, with integrated signal-to-noise ratio S/N∼24 in the range 3≤R/[h-1 Mpc]≤30 over the three lens samples. ζ± has S/N∼19 in the range 8′≤ ≤50′ and 30′≤ ≤150′ for ζ+ and ζ-, respectively. For cosmological parameter inference, we use the dark emulator package, combined with a halo occupation distribution prescription for the relation between galaxies and halos, to model wp and Δς down to quasinonlinear scales, and we estimate cosmological parameters after marginalizing over nuisance parameters. In our baseline analysis we employ an uninformative flat prior of the residual photo-z error, given by Π(Δzph)=U(-1,1), to model a residual bias in the mean redshift of HSC source galaxies. Comparing the relative lensing amplitudes for Δς in the three redshift bins and for ζ± with the single HSC source galaxy sample allows us to calibrate the photo-z parameter Δzph to the precision of σ(Δzph)≃0.09. With these methods, we obtain a robust constraint on the cosmological parameters for the flat ΛCDM model: S8=σ8(ωm/0.3)0.5=0.763-0.036+0.040, or the best-constrained parameter given by S8′=σ8(ωm/0.3)0.22=0.721±0.028, determined with about 4% fractional precision. Based on multidimensional tension metrics, HSC-Y3 data exhibits about 2.5σ tension with the cosmological constraint inferred by Planck for the ΛCDM model, and hints at a nonzero residual photo-z bias implying that the true mean redshift of the HSC galaxies at z≳0.75 is higher than that implied by the original photo-z estimates.

We measure weak lensing cosmic shear power spectra from the 3-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers 416 deg2 of the northern sky, with a mean i-band seeing of 0.59 arcsec and an effective galaxy number density of 15 arcmin-2 within our adopted redshift range. With an i-band magnitude limit of 24.5 mag, and four tomographic redshift bins spanning 0.3≤zph≤1.5 based on photometric redshifts, we obtain a high-significance measurement of the cosmic shear power spectra, with a signal-to-noise ratio of approximately 26.4 in the multipole range 300<ℓ<1800. The accuracy of our power spectrum measurement is tested against realistic mock shear catalogs, and we use these catalogs to get a reliable measurement of the covariance of the power spectrum measurements. We use a robust blinding procedure to avoid confirmation bias, and model various uncertainties and sources of bias in our analysis, including point spread function systematics, redshift distribution uncertainties, the intrinsic alignment of galaxies and the modeling of the matter power spectrum. For a flat ΛCDM model, we find S8σ8(ωm/0.3)0.5=0.776-0.033+0.032, which is in excellent agreement with the constraints from the other HSC Year 3 cosmology analyses, as well as those from a number of other cosmic shear experiments. This result implies a ∼2σ-level tension with the Planck 2018 cosmology. We study the effect that various systematic errors and modeling choices could have on this value, and find that they can shift the best-fit value of S8 by no more than ∼0.5σ, indicating that our result is robust to such systematics.

We perform a blinded cosmology analysis with cosmic shear two-point correlation functions measured from more than 25 million galaxies in the Hyper Suprime-Cam three-year shear catalog in four tomographic redshift bins ranging from 0.3 to 1.5. After conservative masking and galaxy selection, the survey covers 416 deg2 of the northern sky with an effective galaxy number density of 15 arcmin-2 over the four redshift bins. The 2PCFs adopted for cosmology analysis are measured in the angular range; 7.1<θ/arcmin<56.6 for ζ+ and 31.2<θ/arcmin<248 for ζ-, with a total signal-to-noise ratio of 26.6. We apply a conservative, wide, flat prior on the photometric redshift errors on the last two tomographic bins, and the relative magnitudes of the cosmic shear amplitude across four redshift bins allow us to calibrate the photometric redshift errors. With this flat prior on redshift errors, we find ωm=0.256-0.044+0.056 and S8σ8ωm/0.3=0.769-0.034+0.031 (both 68% C.I.) for a flat Λ cold dark matter cosmology. We find, after unblinding, that our constraint on S8 is consistent with the Fourier space cosmic shear and the 3×2 pt analyses on the same HSC dataset. We carefully study the potential systematics from astrophysical and systematic model uncertainties in our fiducial analysis using synthetic data, and report no biases (including projection bias in the posterior space) greater than 0.5σ in the estimation of S8. Our analysis hints that the mean redshifts of the two highest tomographic bins are higher than initially estimated. In addition, a number of consistency tests are conducted to assess the robustness of our analysis. Comparing our result with Planck-2018 cosmic microwave background observations, we find a ∼2σ tension for the ΛCDM model.

We present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering about 416 deg2, and the SDSS DR11 spectroscopic galaxies spanning the redshift range [0.15, 0.70]. We subdivide the SDSS galaxies into three luminosity-cut, and therefore nearly volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, wp(R). We also use the measurements of the galaxy-galaxy weak-lensing signal Δς(R) for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, ζ± measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat ΛCDM cosmology. We use conservative scale cuts, R>12 and 8h-1 Mpc for Δς and wp, respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. We present various validation tests of our model as well as analysis methods. Our baseline analysis yields S8=0.775-0.038+0.043 (68% C.I.) for the ΛCDM model, after marginalizing over uncertainties in other parameters. Our value of S8 is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. We show that various internal consistency tests based on different splits of the data are passed. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model, using Δς(R) and wp(R) down to R>3 and 2h-1 Mpc, respectively.

We utilize the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (SDSS-BOSS) galaxies and its overlap with approximately 416 sq degrees of deep grizy-band imaging from the Subaru Hyper Suprime-Cam Survey (HSC). We perform measurements of three two-point correlations which form the basis of the cosmological inference presented in our companion papers, Miyatake et al. and Sugiyama et al. We use three approximately volume limited subsamples of spectroscopic galaxies by their i-band magnitude from the SDSS-BOSS: LOWZ (0.1<z<0.35), CMASS1 (0.43<z<0.55) and CMASS2 (0.55<z<0.7), respectively. We present high signal-to-noise ratio measurements of the projected correlation functions of these galaxies, which is expected to be proportional to the projected matter correlation function on large scales with a proportionality constant dependent on the bias of galaxies. In order to help break the degeneracy between the amplitude of the matter correlation and the bias of these spectroscopic galaxies, we use the distortions of the shapes of fainter galaxies in HSC due to weak gravitational lensing, to measure the galaxy-galaxy lensing signal, which probes the projected galaxy-matter cross-correlation function of the SDSS-BOSS galaxies. We also measure the cosmic shear correlation functions from HSC galaxies which is related to the projected matter correlation function. We demonstrate the robustness of our measurements by subjecting each of them to a variety of systematic tests. Our use of a single sample of HSC source galaxies is crucial to calibrate any residual systematic biases in the inferred redshifts of our galaxies. We also describe the construction of a suite of mocks: (i) spectroscopic galaxy catalogs which obey the clustering and abundance of each of the three SDSS-BOSS subsamples, and (ii) galaxy shape catalogs which obey the footprint of the HSC survey and have been appropriately sheared by the large-scale structure expected in a Λ Cold Dark Matter model. We use these mock catalogs to compute the covariance of each of our observables.

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.

One of the key questions on active galactic nuclei (AGN) in galaxy clusters is how AGN could affect the formation and evolution of member galaxies and galaxy clusters in the history of the Universe. To address this issue, we investigate the dependence of AGN number fraction (fAGN) on cluster redshift (zcl) and distance from the cluster center (RR200). We focus on more than 27000 galaxy groups and clusters at 0.1 < zcl < 1.4 with more than 1 million member galaxies selected from the Subaru Hyper Suprime-Cam. By combining various AGN selection methods based on infrared (IR), radio, and X-ray data, we identify 2688 AGN. We find that (i) fAGN increases with zcl and (ii) fAGN decreases with RR200. The main contributors to the rapid increase of fAGN towards high-z and cluster center are IR- and radio-selected AGN, respectively. These results indicate that the emergence of the AGN population depends on the environment and redshift, and galaxy groups and clusters at high z play an important role in AGN evolution. We also find that cluster-cluster mergers may not drive AGN activity in at least the cluster center, while we have tentative evidence that cluster-cluster mergers could enhance AGN activity in the outskirts of (particularly massive) galaxy clusters.

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) × 104 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.

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.

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.

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 108.6 M ⊙ at z ∼ 0.2 and 109.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 (<109.2 M ⊙) near the cluster edge, based on their contributions to the quiescent fraction. After controlling for M *, z, and local overdensity, the study still finds a significant radial gradient in the quiescent fraction, indicating active ram pressure stripping in dense environments. That said, analyzing the density trend of the quiescent fraction with other fixed parameters suggests that radial and density-related quenching processes are equally crucial for low-mass cluster galaxies. The study concludes that ram pressure stripping is the primary environmental quenching mechanism for high stellar mass galaxies in clusters. By contrast, ram pressure stripping and density-related quenching processes act comparably for low-mass cluster galaxies around the center. Near the cluster boundary, starvation and harassment become the leading quenching processes for low stellar mass galaxies.

As a fundamental parameter for modern cosmology, the Hubble constant H0 is experiencing a serious crisis. In this paper, we explore an independent approach to measure H0 based on the time-delay cosmography with strong gravitational lensing of a quasar by a galaxy cluster. Specifically we focus on the strong lensing system SDSS J1004+4112 with the maximum image separation of 14.62′′, the first system of a quasar lensed by a galaxy cluster with five multiple images. Incorporating the latest time-delay measurements, we investigate the lens model dependence from the combination of 16 different lens mass models. We find that the lens model dependence is indeed large, with the combined measurement of the Hubble constant of H0=67.5-8.9+14.5 km s-1 Mpc-1 that is obtained by summing posteriors of the Hubble constant from the 16 models with equal weighting. Interestingly, our results show that the value of Hubble constant decreases as the complexity of the perturbation around the lens increases, although weighting based on positional errors of quasar images does not significantly improve the H0 constraint. We find that the 16 different mass models predict largely different shapes of the lensed quasar host galaxy as well as other lensed galaxies behind the cluster. By selecting two mass models that best reproduces those shapes, the constraint on the Hubble constant is significantly tightened to H0=59.1-3.5+3.6 km s-1 Mpc-1. While we caution that our analysis still does not fully explore all the possible mass model uncertainty, our results highlight the importance of including as many constraints as possible such as extended shapes of lensed galaxies for obtaining tight constraints on the Hubble constant from cluster-lensed quasar lens systems.

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 log M ⋆ / M ⊙ = 8.3-8.7. All appear young, with mass-weighted ages <240 Myr, low dust content A V < 0.15 mag, and specific star formation rates sSFR ∼0.25-10 Gyr−1 for most. One z ∼ 9 candidate is consistent with an age <5 Myr and an sSFR ∼10 Gyr−1, as inferred from a strong F444W excess, implying [O iii ]+H β rest-frame equivalent width ∼2000 Å, although an older z ∼ 10 object is also allowed. Another z ∼ 9 candidate is lensed into an arc 2.″4 long with a magnification of μ ∼ 8. This arc is the most spatially resolved galaxy at z ∼ 9 known to date, revealing structures ∼30 pc across. Follow-up spectroscopy of WHL0137−08 with JWST/NIRSpec will be useful to spectroscopically confirm these high-redshift galaxy candidates and to study their physical properties in more detail.

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.

In the version of this article initially published, there was a typographical error in “MARACHAS” in the Acknowledgements text now reading “G.M. received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement no. MARACHAS – DLV-896778”. The error has been corrected in the HTML and PDF versions of the article.

The gravitationally lensed supernova Refsdal appeared in multiple images produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted that an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measurement of the expansion rate of the Universe, quantified by the Hubble constant (H0). Using eight cluster lens models, we infer H0 = 64:8+4.4-4.3 kilometers per second per megaparsec. Using the two models most consistent with the observations, we find H0 = 66:6+4.1-3.3 kilometers per second per megaparsec. The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster.

Unveiling the true nature of dark matter, which manifests itself only through gravity, is one of the principal quests in physics. Leading candidates for dark matter are weakly interacting massive particles or ultralight bosons (axions), at opposite extremes in mass scales, that have been postulated by competing theories to solve deficiencies in the Standard Model of particle physics. Whereas dark matter weakly interacting massive particles behave like discrete particles (ϱDM), quantum interference between dark matter axions is manifested as waves (ψDM). Here, we show that gravitational lensing leaves signatures in multiply lensed images of background galaxies that reveal whether the foreground lensing galaxy inhabits a ϱDM or ψDM halo. Whereas ϱDM lens models leave well documented anomalies between the predicted and observed brightnesses and positions of multiply lensed images, ψDM lens models correctly predict the level of anomalies remaining with ϱDM lens models. More challengingly, when subjected to a battery of tests for reproducing the quadruply lensed triplet images in the system HS 0810+2554, ψDM is able to reproduce all aspects of this system whereas ϱDM often fails. The ability of ψDM to resolve lensing anomalies even in demanding cases such as HS 0810+2554, together with its success in reproducing other astrophysical observations, tilt the balance toward new physics invoking axions.

We report the discovery of a transient seen in a strongly lensed arc at redshift zs = 1.2567 in Hubble Space Telescope imaging of the Abell 370 galaxy cluster. The transient is detected at 29.51 ± 0.14 AB mag in a WFC3/UVIS F200LP difference image made using observations from two different epochs, obtained in the framework of the Flashlights programme, and is also visible in the F350LP band (mF350LP ≈ 30.53 ± 0.76 AB mag). The transient is observed on the negative-parity side of the critical curve at a distance of ∼0.6 arcsec from it, greater than previous examples of lensed stars. The large distance from the critical curve yields a significantly smaller macromagnification, but our simulations show that bright, O/B-type supergiants can reach sufficiently high magnifications to be seen at the observed position and magnitude. In addition, the observed transient image is a trailing image with an observer-frame time delay of ∼+0.8 d from its expected counterpart, so that any transient lasting for longer than that should have also been seen on the minima side and is thus excluded. This, together with the blue colour we measure for the transient (mF200LP - mF350LP ≈ [-0.3, -1.6] AB), rules out most other transient candidates such as (kilo)novae, for example, and makes a lensed star the prime candidate. Assuming that the transient is indeed a lensed star as suggested, many more such events should be detected in the near future in cluster surveys with the Hubble Space Telescope and JWST.

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 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 M1450∗=-25.60-0.30+0.40, a characteristic density φ∗=1.35-0.30+0.47 Gpc-3 mag-1, and a bright-end slope β=-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, 1048.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.

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 376.0 − 5.5 + 5.6 days and the relative magnification to be 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.

We report the spectroscopic follow-up of 175 lensed quasar candidates selected using Gaia Data Release 2 observations following Paper III of this series. Systems include 86 confirmed lensed quasars and a further 17 likely lensed quasars based on imaging and/or similar spectra. We also confirm 11 projected quasar pairs and 11 physical quasar pairs, while 25 systems are left as unclassified quasar pairs – pairs of quasars at the same redshift, which could be either distinct quasars or potential lensed quasars. Especially interesting objects include eight quadruply imaged quasars of which two have BAL sources, an apparent triple, and a doubly lensed LoBaL quasar. The source redshifts and image separations of these new lenses range between 0.65–3.59 and 0.78–6.23 arcsec, respectively. We compare the known population of lensed quasars to an updated mock catalogue at image separations between 1 and 4 arcsec, showing a very good match at z < 1.5. At z > 1.5, only 47 per cent of the predicted number are known, with 56 per cent of these missing lenses at image separations below 1.5 arcsec. The missing higher redshift, small-separation systems will have fainter lensing galaxies, and are partially explained by the unclassified quasar pairs and likely lenses presented in this work, which require deeper imaging. Of the 11 new reported projected quasar pairs, 5 have impact parameters below 10 kpc, almost tripling the number of such systems, which can probe the innermost regions of quasar host galaxies through absorption studies. We also report four new lensed galaxies discovered through our searches, with source redshifts ranging from 0.62 to 2.79.

We present high angular resolution measurements of the thermal Sunyaev–Zel’dovich effect (SZE) toward two galaxy clusters, RCS J2319+0038 at z = 0.9 and HSC J0947−0119 at z = 1.1, by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. They are supplemented with available Chandra X-ray data, optical data taken by Hyper Suprime-Cam on Subaru, and millimeter-wave SZE data from the Atacama Cosmology Telescope. Taking into account departures from spherical symmetry, we have reconstructed non-parametrically the inner pressure profile of two clusters as well as electron temperature and density profiles for RCS J2319+0038. This is one of the first such measurements for an individual cluster at z ≳ 0.9. We find that the inner pressure profile of both clusters is much shallower than that of local cool-core clusters. Our results consistently suggest that RCS J2319+0038 hosts a weak cool core, where radiative cooling is less significant than in local cool cores. On the other hand, HSC J0947−0119 exhibits an even shallower pressure profile than RCS J2319+0038 and is more likely to be a non-cool-core cluster. The SZE centroid position is offset by more than 140 h−170 kpc from the peaks of galaxy distribution in HSC J0947−0119, suggesting a stronger influence of mergers in this cluster. We conclude that these distant clusters are at a very early stage of developing the cool cores typically found in clusters at lower redshifts.

Galaxy mergers produce pairs of supermassive black holes (SMBHs), which may be witnessed as dual quasars if both SMBHs are rapidly accreting. The kiloparsec (kpc)-scale separation represents a physical regime sufficiently close for merger-induced effects to be important1 yet wide enough to be directly resolvable with the facilities currently available. Whereas many kpc-scale, dual active galactic nuclei—the low-luminosity counterparts of quasars—have been observed in low-redshift mergers2, no unambiguous dual quasar is known at cosmic noon (z ≈ 2), the peak of global star formation and quasar activity3,4. Here we report multiwavelength observations of Sloan Digital Sky Survey (SDSS) J0749 + 2255 as a kpc-scale, dual-quasar system hosted by a galaxy merger at cosmic noon (z = 2.17). We discover extended host galaxies associated with the much brighter compact quasar nuclei (separated by 0.46″ or 3.8 kpc) and low-surface-brightness tidal features as evidence for galactic interactions. Unlike its low-redshift and low-luminosity counterparts, SDSS J0749 + 2255 is hosted by massive compact disk-dominated galaxies. The apparent lack of stellar bulges and the fact that SDSS J0749 + 2255 already follows the local SMBH mass–host stellar mass relation, suggest that at least some SMBHs may have formed before their host stellar bulges. While still at kpc-scale separations where the host-galaxy gravitational potential dominates, the two SMBHs may evolve into a gravitationally bound binary system in around 0.22 Gyr.

Ultraviolet light from early galaxies is thought to have ionized gas in the intergalactic medium. However, there are few observational constraints on this epoch because of the faintness of those galaxies and the redshift of their optical light into the infrared. We report the observation, in JWST imaging, of a distant galaxy that is magnified by gravitational lensing. JWST spectroscopy of the galaxy, at rest-frame optical wavelengths, detects strong nebular emission lines that are attributable to oxygen and hydrogen. The measured redshift is z = 9.51 ± 0.01, corresponding to 510 million years after the Big Bang. The galaxy has a radius of 16.2 +4.6-7.2 parsecs, which is substantially more compact than galaxies with equivalent luminosity at z ~ 6 to 8, leading to a high star formation rate surface density.

We conduct a comprehensive study on dropout galaxy candidates at z ∼ 9-16 using the first 90 arcmin2 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 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 * ∼ 108−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.

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 ∼106-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 105 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 (<5 Myr) show evidence of prominent Hβ+[O iii] emission based on photometry with equivalent widths larger than >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 log ( ξ ion [ Hz erg − 1 ] ) ∼ 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 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 ∼[103-105] 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.

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 > 1045.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.

Context. We present the results of a systematic X-ray analysis of optically rich galaxy clusters detected by the Subaru Hyper Suprime-Cam (HSC) survey in the eROSITA Final Equatorial-Depth Survey (eFEDS) field. Aims. Through a joint analysis of the SRG (Spectrum Roentgen Gamma)/eROSITA and Subaru/HSC surveys, we aim to investigate the dynamical status of the optically selected clusters and to derive the cluster scaling relations. Methods. The sample consists of 43 optically selected galaxy clusters with a richness &gt;40 in the redshift range of 0.16–0.89. We systematically analyzed the X-ray images and emission spectra using the eROSITA data. We identified the brightest cluster galaxy (BCG) using the optical and far-infrared databases. We evaluated the cluster’s dynamical status by measuring three quantities: offset between the X-ray peak and BCG position, the gas concentration parameter, and the number of galaxy-density peaks. We investigated the luminosity–temperature and mass–luminosity relations based on eROSITA X-ray spectra and HSC weak-lensing data analyses. Results. Based on these three measurements, we estimated the fraction of relaxed clusters to be 2(&lt; 39)%, which is smaller than that of the X-ray-selected cluster samples. After correcting for a selection bias due to the richness cut, we obtained a shallow L−T slope of 2.1 ± 0.5, which is consistent with the predictions of the self-similar model and the baseline model incorporating a mass–concentration relation. The L−M slope of 1.5 ± 0.3 is in agreement with the above-cited theoretical models as well as the data on the shear-selected clusters in the eFEDs field. Conclusions. Our analysis of high-richness optical clusters yields a small fraction of relaxed clusters and a shallow slope for the luminosity–temperature relation. This suggests that the average X-ray properties of the optical clusters are likely to be different from those observed in the X-ray samples. Thus, the joint eROSITA and HSC observations are a powerful tool in extending the analysis to a larger sample and understanding the selection effect with a view to establish cluster scaling relations.

We present a set of multiwavelength mosaics and photometric catalogs in the Atacama Large Millimeter/ submillimeter Array (ALMA) lensing cluster survey fields. The catalogs were built by the reprocessing of archival data from the Complete Hubble Archive for Galaxy Evolution compilation, taken by the Hubble Space Telescope (HST) in the Reionization Lensing Cluster Survey, Cluster Lensing And Supernova survey with Hubble, and Hubble Frontier Fields. Additionally, we have reconstructed the Spitzer Infrared Array Camera 3.6 and 4.5 μm mosaics, by utilizing all the available archival IPAC Infrared Science Archive/Spitzer Heritage Archive exposures. To alleviate the effect of blending in such a crowded region, we have modeled the Spitzer photometry by convolving the HST detection image with the Spitzer point-spread function using the novel GOLFIR software. The final catalogs contain 218,000 sources, covering a combined area of 690 arcmin2, a factor of ∼2 improvement over the currently existing photometry. A large number of detected sources is a result of reprocessing of all available and sometimes deeper exposures, in conjunction with a combined optical–near-IR detection strategy. These data will serve as an important tool in aiding the search of the submillimeter galaxies in future ALMA surveys, as well as follow-ups of the HST dark and high-z sources with JWST. Coupled with the available HST photometry, the addition of the 3.6 and 4.5 μm bands will allow us to place a better constraint on the photometric redshifts and stellar masses of these objects, thus giving us an opportunity to identify high-redshift candidates for spectroscopic follow-ups and to answer the important questions regarding the Epoch of Reionization and formation of the first galaxies. The mosaics, photometric catalogs, and the best-fit physical properties are publicly available at https:// github.com/dawn-cph/alcs-clusters.

The gravitationally lensed star WHL 0137-LS, nicknamed Earendel, was identified with a photometric redshift z phot = 6.2 ± 0.1 based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera images of Earendel in eight filters spanning 0.8-5.0 μm. In these higher-resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to μ > 4000 and restricting the source plane radius further to r < 0.02 pc, or ∼4000 au. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of T eff ≃ 13,000-16,000 K, assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from log ( L ) = 5.8 to 6.6 L ⊙, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe.

The core-collapse supernova of a massive star rapidly brightens when a shock, produced following the collapse of its core, reaches the stellar surface. As the shock-heated star subsequently expands and cools, its early-time light curve should have a simple dependence on the size of the progenitor1 and therefore final evolutionary state. Measurements of the radius of the progenitor from early light curves exist for only a small sample of nearby supernovae2-14, and almost all lack constraining ultraviolet observations within a day of explosion. The several-day time delays and magnifying ability of galaxy-scale gravitational lenses, however, should provide a powerful tool for measuring the early light curves of distant supernovae, and thereby studying massive stellar populations at high redshift. Here we analyse individual rest-frame exposures in the ultraviolet to the optical taken with the Hubble Space Telescope, which simultaneously capture, in three separate gravitationally lensed images, the early phases of a supernova at redshift z ≈ 3 beginning within 5.8 ± 3.1 hours of explosion. The supernova, seen at a lookback time of approximately 11.5 billion years, is strongly lensed by an early-type galaxy in the Abell 370 cluster. We constrain the pre-explosion radius to be [Formula: see text] solar radii, consistent with a red supergiant. Highly confined and massive circumstellar material at the same radius can also reproduce the light curve, but because no similar low-redshift examples are known, this is unlikely.

We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing (Δς) and projected galaxy clustering (wp) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range 0.15<z<0.7. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of wp in three spectroscopic redshift bins, and as the lens samples for Δς. For the Δς measurements, we select a single sample of 4×106 source galaxies over 140 deg2 from HSC-Y1 with photometric redshifts (photo z) greater than 0.75, enabling a better handle of photo-z errors by comparing the Δς amplitudes for the three lens redshift bins. The deep, high-quality HSC-Y1 data enable significant detections of the Δς signals, with integrated signal-to-noise ratio S/N∼15 in the range 3≤R/[h-1 Mpc]≤30 for the three lens samples, despite the small area coverage. For cosmological parameter inference, we use an input galaxy-halo connection model built on the dark emulator package (which uses an ensemble set of high-resolution N-body simulations and enables fast, accurate computation of the clustering observables) with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the Δς and wp measurements on scales from R≃3 and 2 h-1 Mpc, respectively, up to 30 h-1 Mpc (therefore excluding the baryon acoustic oscillations information) assuming a flat ΛCDM cosmology, marginalizing over about 20 nuisance parameters and demonstrating the robustness of our results to them. With various tests using mock catalogs described in Miyatake et al. [preceding paper, Phys. Rev. D 106, 083519 (2022)10.1103/PhysRevD.106.083519], we show that any bias in the clustering amplitude S8σ8(ωm/0.3)0.5 due to uncertainties in the galaxy-halo connection is less than ∼50% of the statistical uncertainty of S8, unless the assembly biaseffect is unexpectedly large. Our best-fit models have S8=0.795-0.042+0.049 (mode and 68% credible interval) for the flat ΛCDM model; we find tighter constraints on the quantity S8(α=0.17)σ8(ωm/0.3)0.17=0.745-0.031+0.039.

Abstract We conduct a search for galaxy-scale strong gravitational lens systems in Data Release 4 of the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP), consisting of data taken up to the S21A semester. We select 103191 luminous red galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS) sample that have deep multiband imaging from the HSC SSP and use the YattaLens algorithm to identify lens candidates with blue arc-like features automatically. The candidates are visually inspected and graded based on their likelihood of being a lens. We find eight definite lenses, 28 probable lenses, and 138 possible lenses. The new lens candidates generally have lens redshifts in the range 0.3 ≲ zL ≲ 0.9, a key intermediate redshift range to study the evolution of galaxy structure. Follow-up spectroscopy will confirm these new lenses and measure source redshifts to enable detailed lens modeling.

We discuss how small-scale density perturbations on the Fresnel scale affect amplitudes and phases of gravitational waves that are magnified by gravitational lensing in geometric optics. We derive equations that connect the small-scale density perturbations with the amplitude and phase fluctuations to show that such perturbative wave optics effects are significantly boosted in the presence of macro model magnifications such that amplitude and phase fluctuations can easily be observed for highly magnified gravitational waves. We discuss expected signals due to microlensing by stars, dark matter substructure, fuzzy dark matter, and primordial black holes.