Open Access
Issue
A&A
Volume 666, October 2022
Article Number L9
Number of page(s) 10
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202244517
Published online 11 October 2022

© G. B. Caminha et al. 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model.

This Open access funding provided by Max Planck Society.

1. Introduction

Galaxy clusters are the most massive structures in the Universe and are therefore powerful cosmic telescopes for observing faint and distant sources via the gravitational lensing effect (e.g., Atek et al. 2015). When a background source galaxy lies behind a galaxy cluster, the light rays that we see from the source are deflected by the gravitational potential of the foreground galaxy cluster, that is, the source is gravitationally lensed by the cluster. In the regime of strong lensing, we see multiple magnified and distorted images of the background source. The lensing magnifications produced by galaxy clusters can reach factors of ∼100 − 1000, allowing us to detect the faintest sources, which would otherwise be impossible (see e.g., Caminha et al. 2016b; Vanzella et al. 2020, 2021). The positions and morphologies of the multiple lensed images also allow us to probe the distribution of matter, particularly dark matter, in the galaxy cluster, which is crucial for understanding the nature of dark matter (e.g., Bradač et al. 2006; Clowe et al. 2006). Furthermore, an individual galaxy cluster often lenses several background sources into several corresponding ‘families’ of multiple images that straddle different locations around the galaxy cluster. Families of multiple images formed by sources of different redshifts provide measurements of angular diameter distance ratios between the cluster and sources, allowing us to probe the geometry of the Universe (Jullo et al. 2010; Acebron et al. 2017; Caminha et al. 2016a, 2022). In cases where a background source is time varying, such as a quasar or supernova, the time delay(s) between the multiple images provide measurements of the ‘time-delay distance’ and the Hubble constant that sets the expansion rate of the Universe (Refsdal 1964; Suyu et al. 2010; Kelly et al. 2015; Grillo et al. 2018, 2020). Strong lensing galaxy clusters are therefore excellent laboratories for astrophysical and cosmological studies.

The James Webb Space Telescope (JWST; Rigby et al. 2022) operating in the infrared will provide unprecedented observations of high-redshift galaxies, in terms of both sensitivity and angular resolution. Studying these distant galaxies is crucial for understanding how the first galaxies formed and evolved into the structures that we see today. These first galaxies are inherently faint, and the combination of JWST and galaxy cluster lensing is therefore the best way to detect and study these faintest galaxies. As part of the Early Release Observations (ERO), the JWST team publicly released the observations of its first cosmic targets, including the galaxy cluster SMACS J0723.3−7327 (hereafter SMACS J0723, discovered by Ebeling et al. 2001 and Repp & Ebeling 2018). In this Letter we demonstrate the power of strong gravitational lensing and JWST to unveil and study faint distant galaxies. In particular, we identify new sources in the JWST images that were previously undetected in Hubble Space Telescope (HST) imaging, more than doubling the number of families of multiple images. Even though these new families of images do not yet have spectroscopic redshifts, we demonstrate the utility of our cluster lens mass model to constrain the redshifts of these new source galaxies, including one at z > 6.

Our paper is organised as follows. In Sect. 2 we present the pre-JWST archival data of the galaxy cluster SMACS J0723. The JWST observations used in this work are presented in Sect. 3. In Sect. 4 we present our lens models based on archival data only (pre-JWST), the updated version that includes new multiple images identified with JWST, and a comparison to other models in the literature. Finally, we provide a summary in Sect. 5.

2. Pre-JWST archival data

In this section we describe the data used in this work, consisting of HST photometry and Multi Unit Spectroscopic Explorer (MUSE) spectroscopy, to build the pre-JWST lens model.

2.1. HST

Observations from HST of SMACS J0723 were obtained under the treasury programme Reionization Lensing Cluster Survey (RELICS; Coe et al. 2019). Imagings were obtained in the filters F435W, F606W, and F814W using the Advanced Camera for Surveys, and in F105W, F125W, F140W, and F160W using Wide Field Camera 3. All reduced images and catalogues are publicly available at the Mikulski Archive for Space Telescope website1. The imaging and colour information were used to identify multiple images and cluster members with no spectroscopic information. More details are presented in Sect. 4.1.

2.2. VLT/MUSE

MUSE data of this field are available at the ESO Science Archive Facility website2. The field was observed as part of programme ID 0102.A-0718 (PI: A. Edge) under fair conditions, with seeing ≈0.7″ and a moon illumination of 68%, and with an exposure time of 2910 s on target. We processed the raw data following all steps in the standard reduction pipeline (version 2.8.3; Weilbacher et al. 2020). To improve the sky-subtraction, we used the auto-calibration mode and applied the Zurich Atmosphere Purge (Soto et al. 2016) to the data. The final datacube has a 1 arcmin2 field of view and covers the wavelength range 4750 Å–9350 Å.

We extracted the 1D spectra from all HST-detected sources to measure their redshifts. In a second step, we performed a blind search to detect faint emission lines with no clear photometric counterpart. We refer the reader to our previous works (Caminha et al. 2019, 2017a,b) for more details on this procedure. Our final redshift catalogue contains 78 secure and precise redshift measurements; it is presented in Table A.1 and is available at the CDS.

3. JWST Early Release observations

SMACS J0723 is the first gravitational lens observed by JWST. It was one of the targets of the ERO carried out during the telescope commissioning. The observations obtained photometric data from the Near Infrared Camera (NIRCam; e.g., Rieke et al. 2005) and the Mid-Infrared Instrument (MIRI; Rieke et al. 2015, Bouchet et al. 2015) and spectroscopy from the Near Infrared Spectrograph (NIRSpec; Jakobsen et al. 2022) and the Near Infrared Imager and Slitless Spectrograph (NIRISS; Doyon et al. 2012). Since MIRI imaging has a very large point spread function (> 1″ for wavelengths > 7.5 μm) compared to NIRCam, we do not consider mid-infrared data in this work. NIRSpec targeted multiple images: 4.1, 4.2, and 7.1 (see Table 1). No secure spectroscopic redshift has been obtained for images 4.1 or 4.2. For image 7.1, a secure redshift of z = 5.1727 is reported in Mahler et al. (2022) from O III and Hα detections. We incorporated this measurement into our model. We did not consider NIRISS data since slitless spectroscopy is subject to strong spectral contamination, especially in crowded fields such as the core of SMACS J0723. In the next sections we focus on NIRCam imaging to identify new multiple image families not detected with HST.

Table 1.

Multiple image list.

NIRCam ERO observations of SMACS J0723 were carried out in the filters F090W, F150W, F200W, F277W, F356W, and F444W, covering the wavelength range from 0.8 μm, to ≈5.0 μm, under programme ID 2736 (PI: K. Pontoppidan). The total exposure time in each filter was 7537 seconds divided over nine dither patterns, ensuring enough depth to detect and resolve faint sources. We used the standard NIRCam reduction pipeline with the calibration files available with the first data release. We produced mosaics with resolutions of 20 mas pixel−1 for the F090W, F150W, and F200W filters, and 40 mas pixel−1 for F277W, F356W, and F444W.

The final mosaics cover the core of SMACS J0723 and the entire region of multiple image formation. In Fig. 1 we show the colour composite image created using the software trilogy (Coe et al. 2012). We aligned the final images to the RELICS HST data to maintain the same coordinate reference as previous observations. This introduces a small rigid offset with respect to the nominal coordinates of the ERO release. The astrometric precision between the JWST and HST data we obtained has an rms of 0 . 015 $ 0{{\overset{\prime\prime}{.}}}015 $. Moreover, because of the nature of these early JWST observations, the calibration files are not ideal and the absolute flux calibration might have some issues. Nonetheless, these caveats do not influence the colours or morphology of multiple images in the imaging, with no significant effect on our model.

thumbnail Fig. 1.

Galaxy cluster SMACS J0723. Colour composed image using JWST/NIRCam imaging with the filters F090W, F150W in blue; F200W, F277W in green; and F356W and F444W in red. Cyan circles show the positions of multiple images identified in the pre-JWST lens model and the green circles are the newly identified multiple images with JWST. The red line is the critical line for a source at redshift z = 10.

4. Strong lens model

We used the software lenstool (Kneib et al. 1996; Jullo et al. 2007; Jullo & Kneib 2009) to model the mass distribution of SMACS J0723. The fiducial mass model consists of an elliptical cluster-scale dark matter halo (six free parameters), a truncated spherical isothermal mass profile for each cluster galaxy member (two free parameters with a constant M/L scaling relation for the members), and an external shear (two free parameters). The model constraints are the positions of multiple images and (when available) spectroscopic redshifts. For families of multiple images with no spectroscopic measurements, we leave the redshift values free to vary in the model. We built two versions of the lens model: The first uses information from HST and MUSE only (i.e., data available prior to the JWST ERO release). The second model was built upon JWST NIRCam imaging, which allowed us to confirm 30 additional, secure multiple images from 11 individual sources. Naturally, these new identifications are extremely faint or not detected in the HST optical and near-infrared imaging. We carefully checked the MUSE datacube around these positions, but as expected, no optical emission was detected. In the following sections we describe the two models we present in this work.

4.1. Pre-JWST lens model

Before the release of JWST/NIRCam imaging, the best dataset available for SMACS J0723 was photometry from the RELICS program and spectroscopy from MUSE (see Sects. 2.1 and 2.2). From the spectroscopic data, we confirm a total of 23 cluster members in the redshift range of 0.367 − 0.408, corresponding to a rest-frame velocity of ±3000 km s−1 from the cluster mean redshift of z = 0.387. However, the spectroscopically confirmed members are limited to the MUSE field of view of 1 arcmin2. We computed the 68% confidence intervals from the colour and photometric redshift distributions of spectroscopically confirmed members and used these intervals to select additional cluster members with no spectroscopic information. In Fig. 2 we show the colour–magnitude diagram for the HST filters F606W and F160W, where the red sequence is clearly defined by the spectroscopically confirmed members. In total, we selected 38 photometric members down to a magnitude of 24 in F160W.

thumbnail Fig. 2.

Colour–magnitude diagram of galaxies in the field of view of SMACS J0723. The spectroscopically confirmed cluster members (red points) define the red sequence. The colour and photometric redshift distribution of spectroscopically confirmed members are used to select 38 additional cluster members with no spectroscopic information (orange points) down to a magnitude of 24 in F160W.

From the MUSE spectroscopy, we confirmed the redshifts of four families of multiple images. In Table 1 we list all the positions and redshifts of multiple images used as constraints in the lens model. In this first version of the model, two families have no secure spectroscopic confirmations, but are clearly detected in HST and have colours and parities as expected in strong lensing. Our pre-JWST model thus has 19 multiple images as input from six families, and only two of the redshifts are free to vary in the model.

The mass profiles we adopted are the pseudo-isothermal elliptical parameterisation (PIEMD; Kassiola & Kovner 1993), to describe the extended cluster mass (dark matter), and the dual pseudo-isothermal mass profile with axial symmetry (Elíasdóttir et al. 2007; Suyu & Halkola 2010), to describe the 61 cluster members. As commonly done in the literature, we assumed a constant mass-to-light ratio for the cluster members. With this, the number of free parameters describing these components is reduced to two normalisations, namely σ v , gals norm . $ \sigma_{\mathrm{v,gals}}^{\mathrm{norm.}} $ and r cut , gals norm . $ r_{\mathrm{cut, gals}}^{\mathrm{norm.}} $. Finally, we also considered an external shear component to account for external massive perturbers that can affect the light deflection of background sources.

Our fiducial cluster mass model was obtained following the procedure in Caminha et al. (2019), in order to find the combination of mass profiles that can reproduce the position of multiple images well while not overfitting the data. We tested a range of combinations of mass profiles, from the simplest one, composed of a single PIEMD halo profile plus cluster members, to models with additional cluster halo components, external shear, and individual cluster members with mass parameters varying freely outside the mass-to-light scaling relation. We find that the parameterisation with the smallest values of the Bayesian information criterion (Schwarz 1978) and Akaike information criterion (Akaike 1974) is the one with ten parameters for the mass distribution (as described at the beginning of Sect. 4) and the two unknown redshifts of families 1 and 4 (see Table 1). In this model, the best-fit offset between the multiple observed images and the model predicted positions is δ rms = 0 . 39 $ \delta_{\mathrm{rms}} = 0{{\overset{\prime\prime}{.}}}39 $. We note that our best-fit δrms is significantly lower than that previously obtained by Golubchik et al. (2022), who quoted a value of δ rms = 2 . 3 $ \delta_{\mathrm{rms}}=2{{\overset{\prime\prime}{.}}}3 $.

4.2. JWST lens model

With the new JWST/NIRCam imaging, at wavelengths not reachable with HST and the high spatial resolution, we were able to significantly improve our lens model constraints. We carefully inspected the near-infrared imaging in order to find new multiple images that were not previously identified. Thanks to the colour information and resolution, we obtained a sample of 30 new secure multiple images from 11 background galaxies, increasing by a factor of two the number of model constraints compared to the pre-JWST (HST and MUSE only) lens model.

Given that these sources are very faint in observed optical wavelengths, none have spectroscopic confirmation from MUSE and, in most cases, they have no clear HST counterpart. In Table 1 we list all new multiple images as well as their model redshifts. Moreover, in Fig. A.3 we show the JWST NIRCam cutouts around all multiple images used in our lens model.

The lens model with the new multiple images has a best-fit δrms of 0 . 51 $ 0{{\overset{\prime\prime}{.}}}51 $. Considering that the number of model constraints is around twice that of the pre-JWST model, this small increase in δrms is expected. Nonetheless, it is significantly lower than what was reported in previous published works on SMACS J0723. In Fig. 1 we show the critical line for a source at z = 10 overlaid on the JWST/NIRCam imaging. The lens model, including magnification maps and lenstool configuration files, are available at the CDS and upon request to the authors. We encourage the community to further explore the lens model presented in this work.

In Table A.2 we list the constrained total mass parameters, and in Fig. A.1 we show the total mass and density profiles. We find a mass of 8 . 7 0.2 + 0.2 × 10 13 M $ 8.7_{-0.2}^{+0.2} \times 10^{13}\,M_{\odot} $ at a radius of 128 kpc, approximately the Einstein radius for a source at z → ∞. At this radius, this mass estimate is 51% more accurate when compared to our pre-JWST lens model.

One special case is family 8. Images 8.1 and 8.2, shown in Fig. A.3, are secure multiple images thanks to their colours and morphology. They are composed of two main clumps and some smaller possible substructures. Image #8.3 is very likely a counter-image, showing similar colour and the double-clump morphology. In order to avoid any bias in our lens model, we did not include #8.3 as an observational constraint. The total lensing magnifications are 6 . 1 0.6 + 0.7 $ 6.1_{-0.6}^{+0.7} $, 9 . 5 0.8 + 0.9 $ 9.5_{-0.8}^{+0.9} $, and 4 . 5 0.4 + 0.5 $ 4.5_{-0.4}^{+0.5} $ for 8.1, 8.2, and #8.3, respectively. The multiple images are located in the easternmost region of the cluster (see Fig. 1). This makes it difficult to precisely constrain its redshift from the model. However, the lens model indicates it might be at z > 7.5 from the posterior probability distribution.

4.3. Comparison to models in the literature

Mahler et al. (2022) and Pascale et al. (2022) also used the new JWST observations, particularly the imaging, to model SMACS J0723. The model by Pascale et al. (2022) was built upon the model of Golubchik et al. (2022), which was based on HST observations and created before the JWST ERO. We briefly compare our model ingredients and results to these studies.

All studies identified new families of multiple images in the NIRCam imaging: 16 in Mahler et al. (2022), 14 in Pascale et al. (2022), and 11 in this work. Since no families have spectroscopic redshift confirmations (except one identified by Mahler et al. 2022 with NIRSpec), we deliberately used only highly secure families of multiple images (in Table 1) to constrain our mass model in order to avoid misidentifications of families that could bias our mass model parameters. This explains our lower number of reported families.

All three models are parametric in having a cluster halo component with cluster members, and we have a further external shear component that is not present in the other two studies. Our robust selection of 61 cluster members is based on both MUSE spectroscopy and HST photometry, whereas Mahler et al. (2022) and Pascale et al. (2022) relied on HST photometry to identify ≈130 members (which will be refined in future work by including spectroscopy).

Pascale et al. (2022) have posted three versions of their model with measured photometric redshifts based on the NIRCam photometry and also predicted the redshifts based on their mass model. The first version from Pascale et al. (2022)3 produced ten redshift estimates for the new JWST lensed sources, as did our model, although only six of the ten are the same. While the model redshifts from both models for five of these sources are all in the range ≈1 − 3, only one of the five agrees within the 1σ uncertainty. The sixth source is our family 7 at z > 7.5, whereas its redshift is unconstrained by Pascale et al. (2022). The model redshifts are therefore highly model dependent and rely on the goodness of the model fits. Future spectroscopic redshift measurements of some of these sources would be valuable. The third version of Pascale et al. (2022) was posted after the release of our model and now has nine redshift estimates in common with ours, all agreeing within the 1σ uncertainty.

Even though a direct comparison of the different models is not straightforward given the different families of multiple images used, the δrms of the models provides an indication of the goodness of the model fit. Our δrms of 0 . 51 $ 0{{\overset{\prime\prime}{.}}}51 $ is a factor of ≈2 lower than the 1 . 08 $ 1{{\overset{\prime\prime}{.}}}08 $ from Mahler et al. (2022) and the 0 . 93 $ 0{{\overset{\prime\prime}{.}}}93 $ from the first version of Pascale et al. (2022). We attribute our better fitting model to (1) our use of only the secure families of multiple images and (2) our robust cluster member selection based on both spectroscopy and photometry throughout our work since the beginning. The third version of Pascale et al. (2022) quotes δ rms = 0 . 48 $ \delta_{\mathrm{rms}}=0{{\overset{\prime\prime}{.}}}48 $, similar to the value obtained by our model.

5. Conclusions

We present our strong lens model of SMACS J0723, the first gravitational lens ever observed by JWST, using newly available data from NIRCam. With the new data, we have doubled the number of model constraints, increasing the number of multiple images from the 19 identified with HST and MUSE data to 49. This demonstrates the unique capabilities of JWST in discovering faint and/or high-redshift galaxies that were previously undetected in HST observations. Our lens model is capable of reproducing the positions of all multiple images with a high accuracy of δ rms = 0 . 51 $ \delta_{\mathrm{rms}} = 0{{\overset{\prime\prime}{.}}}51 $. We report the lens model redshifts for 12 lensed galaxies, most of which are very faint; it remains challenging to obtain spectroscopic information for them. While most of the lensed galaxies have modelled redshifts ≲3, one object stands out as a high-redshift galaxy at z > 7.5, with image 8.2 having a magnification of | μ | 9 . 5 0.8 + 0.9 $ |\mu| \approx 9.5_{-0.8}^{+0.9} $. This might be the first galaxy with spatially resolved substructures at the epoch of re-ionisation of the Universe (see Fig. A.3) and an interesting target for follow-up deep spectroscopy. We report spectroscopic redshifts from MUSE archival data and their lensing magnifications. Finally, our lens models, including magnification maps and configuration files for the software lenstool, are made publicly available.

These first observations by JWST of SMACS J0723 have proven to be a treasure trove, allowing us to reduce the uncertainties in our cluster mass measurements by 51% even without incorporating any spectroscopic information on the newly identified lensed galaxies. We anticipate future JWST observations of other strong lensing galaxy clusters to be transformative in both constraining cluster mass distributions and studying high-redshift galaxies.


Acknowledgments

The authors thank the anonymous referee for the useful comments on the manuscript. We thank C. Grillo, P. Rosati and D. Sluse for constructive discussions during the preparation of this work. GBC and SHS acknowledge the Max Planck Society for financial support through the Max Planck Research Group for SHS and the academic support from the German Centre for Cosmological Lensing. This work is partially funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2094 – 390783311. AM, PB, AA and EV acknowledge financial contributions by PRIN-MIUR 2017WSCC32 “Zooming into dark matter and proto-galaxies with massive lensing clusters” (PI: P.Rosati), INAF “main-stream” 1.05.01.86.20: “Deep and wide view of galaxy clusters” (PI: M. Nonino) and INAF “main-stream” 1.05.01.86.31 “The deepest view of high-redshift galaxies and globular cluster precursors in the early Universe” (PI: E. Vanzella). AA has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 101024195 – ROSEAU.

References

  1. Acebron, A., Jullo, E., Limousin, M., et al. 2017, MNRAS, 470, 1809 [Google Scholar]
  2. Akaike, H. 1974, IEEE Trans. Auto. Control, 19, 716 [Google Scholar]
  3. Atek, H., Richard, J., Kneib, J.-P., et al. 2015, ApJ, 800, 18 [Google Scholar]
  4. Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bouchet, P., García-Marín, M., Lagage, P. O., et al. 2015, PASP, 127, 612 [CrossRef] [Google Scholar]
  6. Bradač, M., Clowe, D., Gonzalez, A. H., et al. 2006, ApJ, 652, 937 [CrossRef] [Google Scholar]
  7. Caminha, G. B., Grillo, C., Rosati, P., et al. 2016a, A&A, 587, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Caminha, G. B., Karman, W., Rosati, P., et al. 2016b, A&A, 595, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Caminha, G. B., Grillo, C., Rosati, P., et al. 2017a, A&A, 600, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Caminha, G. B., Grillo, C., Rosati, P., et al. 2017b, A&A, 607, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Caminha, G. B., Rosati, P., Grillo, C., et al. 2019, A&A, 632, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Caminha, G. B., Suyu, S. H., Grillo, C., & Rosati, P. 2022, A&A, 657, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Clowe, D., Bradač, M., Gonzalez, A. H., et al. 2006, ApJ, 648, L109 [NASA ADS] [CrossRef] [Google Scholar]
  14. Coe, D., Umetsu, K., Zitrin, A., et al. 2012, ApJ, 757, 22 [NASA ADS] [CrossRef] [Google Scholar]
  15. Coe, D., Salmon, B., Bradač, M., et al. 2019, ApJ, 884, 85 [Google Scholar]
  16. Doyon, R., Hutchings, J. B., Beaulieu, M., et al. 2012, in Space Telescopes and Instrumentation 2012: Optical, Infrared, and Millimeter Wave, eds. M. C. Clampin, G. G. Fazio, H. A. MacEwen, J. Oschmann, & M. Jacobus, SPIE Conf. Ser., 8442, 84422R [NASA ADS] [CrossRef] [Google Scholar]
  17. Ebeling, H., Edge, A. C., & Henry, J. P. 2001, ApJ, 553, 668 [Google Scholar]
  18. Elíasdóttir, Á., Limousin, M., Richard, J., et al. 2007, ArXiv e-prints [arXiv:0710.5636] [Google Scholar]
  19. Golubchik, M., Furtak, L. J., Meena, A. K., & Zitrin, A. 2022, ArXiv e-prints [arXiv:2207.05007] [Google Scholar]
  20. Grillo, C., Rosati, P., Suyu, S. H., et al. 2018, ApJ, 860, 94 [Google Scholar]
  21. Grillo, C., Rosati, P., Suyu, S. H., et al. 2020, ApJ, 898, 87 [Google Scholar]
  22. Jakobsen, P., Ferruit, P., Alves de Oliveira, C., et al. 2022, A&A, 661, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Jullo, E., & Kneib, J. P. 2009, MNRAS, 395, 1319 [NASA ADS] [CrossRef] [Google Scholar]
  24. Jullo, E., Kneib, J. P., Limousin, M., et al. 2007, New J. Phys., 9, 447 [Google Scholar]
  25. Jullo, E., Natarajan, P., Kneib, J. P., et al. 2010, Science, 329, 924 [Google Scholar]
  26. Kassiola, A., & Kovner, I. 1993, ApJ, 417, 450 [Google Scholar]
  27. Kelly, P. L., Rodney, S. A., Treu, T., et al. 2015, Science, 347, 1123 [Google Scholar]
  28. Kneib, J. P., Ellis, R. S., Smail, I., Couch, W. J., & Sharples, R. M. 1996, ApJ, 471, 643 [Google Scholar]
  29. Mahler, G., Jauzac, M., Richard, J., et al. 2022, ApJ, submitted [arXiv:2207.07101] [Google Scholar]
  30. Pascale, M., Frye, B., Diego, J., et al. 2022, ApJ, submitted [arXiv:2207.07102] [Google Scholar]
  31. Refsdal, S. 1964, MNRAS, 128, 307 [NASA ADS] [CrossRef] [Google Scholar]
  32. Repp, A., & Ebeling, H. 2018, MNRAS, 479, 844 [NASA ADS] [Google Scholar]
  33. Rieke, M. J., Kelly, D., & Horner, S. 2005, in Cryogenic Optical Systems and Instruments XI, eds. J. B. Heaney, & L. G. Burriesci, SPIE Conf. Ser., 5904, 1 [NASA ADS] [Google Scholar]
  34. Rieke, G. H., Wright, G. S., Böker, T., et al. 2015, PASP, 127, 584 [NASA ADS] [CrossRef] [Google Scholar]
  35. Rigby, J., Perrin, M., McElwain, M., et al. 2022, ArXiv e-prints [arXiv:2207.05632] [Google Scholar]
  36. Schwarz, G. 1978, Ann. Statist., 6, 461 [Google Scholar]
  37. Soto, K. T., Lilly, S. J., Bacon, R., Richard, J., & Conseil, S. 2016, MNRAS, 458, 3210 [Google Scholar]
  38. Suyu, S. H., & Halkola, A. 2010, A&A, 524, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Suyu, S. H., Marshall, P. J., Auger, M. W., et al. 2010, ApJ, 711, 201 [Google Scholar]
  40. Vanzella, E., Meneghetti, M., Caminha, G. B., et al. 2020, MNRAS, 494, L81 [Google Scholar]
  41. Vanzella, E., Caminha, G. B., Rosati, P., et al. 2021, A&A, 646, A57 [EDP Sciences] [Google Scholar]
  42. Weilbacher, P. M., Palsa, R., Streicher, O., et al. 2020, A&A, 641, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Supporting material

In this appendix we present the supporting material of this work. The full redshift catalogue from the MUSE data and the lensing magnification estimates from the JWST model are presented in Table A.1. In Table A.2 we show the recovered mass model parameters for the JWST lens model, and the total mass and density profiles in Fig. A.1. In Fig. A.2 we present the normalised probability density function of the lens model-z for all multiple image families, except for family 8 because its redshift value is above 7.5 and otherwise unconstrained by our model. Figure A.3 shows the cutouts of all multiple images used in our lens models, and Fig. A.4 the spectra of confirmed multiple images.

thumbnail Fig. A.1.

Total mass and density profiles of SMACS J0723 from our pre-JWST (red) and JWST (green) lens models. The regions correspond to the 95% confidence, and the positions of multiple images are indicated by vertical lines.

thumbnail Fig. A.2.

Probability distribution functions (PDF) for lens model redshifts for multiple image families. Median and 68% confidence levels are listed in Table 1. The PDF for family 8 is not shown in this figure for better visualisation of the other families’ PDFs.

thumbnail Fig. A.3.

Multiple image cutouts. The white scale bar has a length of 1″.

thumbnail Fig. A.3.

continued.

thumbnail Fig. A.3.

continued.

thumbnail Fig. A.4.

MUSE spectrum of confirmed multiple images. Blues lines are the spectra extracted from the data cube, and the shaded regions are the re-scaled variance. Vertical lines indicate the O II emission at the source redshift.

Table A.1.

MUSE spectroscopic catalogue.

Table A.2.

Recovered mass parameters from the JWST lens model.

All Tables

Table 1.

Multiple image list.

Table A.1.

MUSE spectroscopic catalogue.

Table A.2.

Recovered mass parameters from the JWST lens model.

All Figures

thumbnail Fig. 1.

Galaxy cluster SMACS J0723. Colour composed image using JWST/NIRCam imaging with the filters F090W, F150W in blue; F200W, F277W in green; and F356W and F444W in red. Cyan circles show the positions of multiple images identified in the pre-JWST lens model and the green circles are the newly identified multiple images with JWST. The red line is the critical line for a source at redshift z = 10.

In the text
thumbnail Fig. 2.

Colour–magnitude diagram of galaxies in the field of view of SMACS J0723. The spectroscopically confirmed cluster members (red points) define the red sequence. The colour and photometric redshift distribution of spectroscopically confirmed members are used to select 38 additional cluster members with no spectroscopic information (orange points) down to a magnitude of 24 in F160W.

In the text
thumbnail Fig. A.1.

Total mass and density profiles of SMACS J0723 from our pre-JWST (red) and JWST (green) lens models. The regions correspond to the 95% confidence, and the positions of multiple images are indicated by vertical lines.

In the text
thumbnail Fig. A.2.

Probability distribution functions (PDF) for lens model redshifts for multiple image families. Median and 68% confidence levels are listed in Table 1. The PDF for family 8 is not shown in this figure for better visualisation of the other families’ PDFs.

In the text
thumbnail Fig. A.3.

Multiple image cutouts. The white scale bar has a length of 1″.

In the text
thumbnail Fig. A.4.

MUSE spectrum of confirmed multiple images. Blues lines are the spectra extracted from the data cube, and the shaded regions are the re-scaled variance. Vertical lines indicate the O II emission at the source redshift.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.