Search

The Euclid mission is generating a vast amount of imaging data in four broadband filters at high angular resolution. This will allow the detailed study of mass, metallicity, and stellar populations across galaxies, which will constrain their formation and evolutionary pathways. Transforming the Euclid imaging for large samples of galaxies into maps of physical parameters in an efficient and reliable manner is an outstanding challenge. We investigate the power and reliability of machine learning techniques to extract the distribution of physical parameters within well-resolved galaxies. We focus on estimating stellar mass surface density, mass-averaged stellar metallicity and age. We generate noise-free, synthetic high-resolution imaging data in the Euclid photometric bands for a set of 1154 galaxies from the TNG50 cosmological simulation. The images are generated with the SKIRT radiative transfer code, taking into account the complex 3D distribution of stellar populations and interstellar dust attenuation. We use a machine learning framework to map the idealised mock observational data to the physical parameters on a pixel-by-pixel basis. We find that stellar mass surface density can be accurately recovered with a $\leq 0.130 {\rm \,dex}$ scatter. Conversely, stellar metallicity and age estimates are, as expected, less robust, but still contain significant information which originates from underlying correlations at a sub-kpc scale between stellar mass surface density and stellar population properties.

We introduce the Bias-free Extragalactic Analysis for Cosmic Origins with NIRCam (BEACON) survey, a JWST Cycle2 program allocated up to 600 pure-parallel hours of observations. BEACON explores high-latitude areas of the sky with JWST/NIRCam over $\sim100$ independent sightlines, totaling $\sim0.3$deg$^2$, reaching a median F444W depth of $\approx28.2$AB mag (5$\sigma$). Based on existing JWST observations in legacy fields, we estimate that BEACON will photometrically identify 25--150 galaxies at $z>10$ and 500--1000 at $z\sim7$--10 uniquely enabled by an efficient multiple filter configuration spanning $0.9$--5.0$\mu$m. The expected sample size of $z>10$ galaxies will allow us to obtain robust number density estimates and to discriminate between different models of early star formation. In this paper, we present an overview of the survey design and initial results using the first 19 fields. We present 129 galaxy candidates at $z>7$ identified in those fields, including 11 galaxies at $z>10$ and several UV-luminous ($M_{\rm UV}<-21$mag) galaxies at $z\sim8$. The number densities of $z<13$ galaxies inferred from the initial fields are overall consistent with those in the literature. Despite reaching a considerably large volume ($\sim10^5$Mpc$^3$), however, we find no galaxy candidates at $z>13$, providing us with a complimentary insight into early galaxy evolution with minimal cosmic variance. We publish imaging and catalog data products for these initial fields. Upon survey completion, all BEACON data will be coherently processed and distributed to the community along with catalogs for redshift and other physical quantities., Comment: Submitted to ApJ; DR1 data release will be made on the team website (https://beacon-jwst.github.io); Fig. 8 has been updated

The VLT Survey Telescope Survey of Mass Assembly and Structural Hierarchy (VST-SMASH) aims to detect tidal features and remnants around very nearby galaxies, a unique and essential diagnostic of the hierarchical nature of galaxy formation. Leveraging optimal sky conditions at ESO's Paranal Observatory, combined with the VST's multi-band optical filters, VST-SMASH aims to be the definitive survey of stellar streams and tidal remnants in the Local Volume, targeting a low surface-brightness limit of $\mu \sim$ 30 mag arcsec$^{-2}$ in the g and r bands, and $\mu \sim$ 28 mag arcsec$^{-2}$ in the i band, in a volume-limited sample of local galaxies within 11 Mpc and the Euclid footprint., Comment: 4 pages, 2 figures, 1 table, published in the ESO Messenger 193

Powerful, galactic outflows driven by Active Galactic Nuclei (AGNs) are commonly considered as a main mechanism to regulate star formation in massive galaxies. Ultra- and hyper-luminous IR galaxies (U/HyLIRGs) are thought to represent a transition phase of galaxies from a rapidly growing period to a quiescent status as gas swept out by outflows, providing a laboratory to investigate outflows and their feedback effects on the hosts. In this paper we report recent Gemini and ALMA observations of a HyLIRG, J1126 at $z=0.46842$, which has been identified with a puzzling co-existence of a fast ionized outflow ($>2000$ km s$^{-1}$) and an intense starburst (star formation rate of 800 $M_{\odot}$ yr$^{-1}$). The Gemini observation shows the fast ionized outflow is extended to several kpc with a mass-loss rate of 180 $M_{\odot}$ yr$^{-1}$. A massive molecular outflow with a high mass-loss rate (2500 $M_{\odot}$ yr$^{-1}$) is revealed by ALMA. The multi-phase outflows show large factors of momentum boost and loading of kinetic power, indicating a driving by thermal pressure of a nuclear hot wind and/or radiation pressure of a highly obscured AGN. In addition to ejection of kinetic energy, it is also found that the powerful outflow can induce an ionizing shock in the galaxy disk and enhance the excitation and dissociation of molecular gas. The powerful outflow probably results in an instantaneous negative feedback and shows potential to regulate the host growth in a long term., Comment: 35 pages, 26 figures, accepted to be published in ApJ

The advent of large integral field spectroscopic surveys has found that elliptical galaxies (EGs) can be classified into two classes: the fast rotators (whose kinematics are dominated by rotation) and the slow rotators (which exhibit slow or no rotation pattern). It is often suggested that while the slow rotators typically have boxy isophotal shapes, have a high $\alpha$-to-iron abundance ratio, and are quite massive, the fast rotators often exhibit the opposite properties (that is, having disky isophotes, lower $\alpha$-to-iron ratio, and of typical masses). Whether the EGs consist of two distinct populations (i.e., a dichotomy exists), remains an unsolved issue. To examine the existence of the dichotomy, we used a sample of 1,895 EGs from the SDSS-IV MaNGA survey, and measured robustly the stellar kinematics, isophotal shapes, and [Mg/Fe] ratio. We confirmed the previous finding that the bulk of the EGs are disky (65%) and fast rotators (67%), but found no evidence supporting a dichotomy, based on a principal component analysis. The different classes (boxy/disky and slow/fast rotators) of EGs occupy slightly different loci in the principal component space. This may explain the observed trends that led to the premature support of a dichotomy based on small samples of galaxies., Comment: 27 pages, 19 figures, 5 tables. Submitted to ApJ. Comments are welcome!

Investigating the metal enrichment in the early universe helps us constrain theories about the first stars and study the ages of galaxies. The lensed galaxy MACS0647$-$JD at $z=10.17$ is the brightest galaxy known at $z > 10$. Previous work analyzing JWST NIRSpec and MIRI data yielded a direct metallicity $\rm{12+log(O/H)}=7.79\pm0.09$ ($\sim$ 0.13 $Z_\odot$) and electron density $\rm{log}(n_e / \rm{cm^{-3}}) = 2.9 \pm 0.5$, the most distant such measurements to date. Here we estimate the direct C/O abundance for the first time at $z > 10$, finding a sub-solar ${\rm log(C/O)}=-0.44^{+0.06}_{-0.07}$. This is higher than other $z>6$ galaxies with direct C/O measurements, likely due to higher metallicity. It is also slightly higher than galaxies in the local universe with similar metallicity. This may suggest a very efficient and rapid burst of star formation, a low effective oxygen abundance yield, or the presence of unusual stellar populations including supermassive stars. Alternatively, the strong CIII]${\rm {\lambda}{\lambda}}$1907,1909 emission ($14\pm 3\,{\r{A}}$ rest-frame EW) may originate from just one of the two component star clusters JDB ($r \sim 20$ pc). Future NIRSpec IFU spectroscopic observations of MACS0647$-$JD will be promising for disentangling C/O in the two components to constrain the chemistry of individual star clusters just 460 Myr after the Big Bang., Comment: 9 pages, 2 figures, submitted

Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid telescope, using both simulated pre-launch images and the first early-release observations of the Fornax galaxy cluster. The Euclid Wide Survey will provide high-spatial resolution VIS imaging in the broad IE band as well as near-infrared photometry (YE, JE, and HE). We estimate that the galaxies within 100 Mpc in the footprint of the Euclid survey host around 830 000 EGCs of which about 350 000 are within the survey's detection limits. For about half of these EGCs, three infrared colours will be available as well. For any galaxy within 50Mpc the brighter half of its GC luminosity function will be detectable by the Euclid Wide Survey. The detectability of EGCs is mainly driven by the residual surface brightness of their host galaxy. We find that an automated machine-learning EGC-classification method based on real Euclid data of the Fornax galaxy cluster provides an efficient method to generate high purity and high completeness GC candidate catalogues. We confirm that EGCs are spatially resolved compared to pure point sources in VIS images of Fornax. Our analysis of both simulated and first on-sky data show that Euclid will increase the number of GCs accessible with high-resolution imaging substantially compared to previous surveys, and will permit the study of GCs in the outskirts of their hosts. Euclid is unique in enabling systematic studies of EGCs in a spatially unbiased and homogeneous manner and is primed to improve our understanding of many understudied aspects of GC astrophysics., Comment: Accepted in A&A

The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution., Comment: Submitted to A&A, 44 pages, 35 figures, Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

We present an analysis of Euclid observations of a 0.5 deg$^2$ field in the central region of the Fornax galaxy cluster that were acquired during the performance verification phase. With these data, we investigate the potential of Euclid for identifying GCs at 20 Mpc, and validate the search methods using artificial GCs and known GCs within the field from the literature. Our analysis of artificial GCs injected into the data shows that Euclid's data in $I_{\rm E}$ band is 80% complete at about $I_{\rm E} \sim 26.0$ mag ($M_{V\rm } \sim -5.0$ mag), and resolves GCs as small as $r_{\rm h} = 2.5$ pc. In the $I_{\rm E}$ band, we detect more than 95% of the known GCs from previous spectroscopic surveys and GC candidates of the ACS Fornax Cluster Survey, of which more than 80% are resolved. We identify more than 5000 new GC candidates within the field of view down to $I_{\rm E}$ mag, about 1.5 mag fainter than the typical GC luminosity function turn-over magnitude, and investigate their spatial distribution within the intracluster field. We then focus on the GC candidates around dwarf galaxies and investigate their numbers, stacked luminosity distribution and stacked radial distribution. While the overall GC properties are consistent with those in the literature, an interesting over-representation of relatively bright candidates is found within a small number of relatively GC-rich dwarf galaxies. Our work confirms the capabilities of Euclid data in detecting GCs and separating them from foreground and background contaminants at a distance of 20 Mpc, particularly for low-GC count systems such as dwarf galaxies., Comment: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level., Comment: Submitted to A&A, 44 pages, 36 figures - Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

Recent studies suggested that the ejecta velocity of Type Ia supernova (SN Ia) is a promising indicator in distinguishing the progenitor systems and explosion mechanisms. By classifying the SNe Ia based on their ejecta velocities, studies found SNe Ia with high Si II $\lambda$6355 velocities (HV SNe Ia; v>12000 km/s) tend to be physically different from their normal-velocity counterparts (NV SNe Ia). In this work, we revisit the low-$z$ sample studied in previous work and closely look into the spatially resolved environment local to the site of SN explosion. Our results reveal a possible trend (at $2.4\sigma$ significance) that HV SNe Ia are likely associated with older stellar populations than NV SNe Ia. While the trend is inconclusive, the local host-galaxy sample studied in this work is likely skewed toward massive galaxies, limiting the parameter space that we would like to investigate from the original parent sample. Nevertheless, our results do not rule out the possibility that parameters other than the host-galaxy age (such as metallicity) could be the underlying factors driving the differences between HV and NV SNe Ia due to the limitation of our dataset., Comment: 13 pages, 6 figures, accepted for publication in MNRAS

We study the impact of black hole nuclear activity on both the global and radial star formation rate (SFR) profiles in X-ray-selected active galactic nuclei (AGN) in the field of miniJPAS, the precursor of the much wider J-PAS project. Our sample includes 32 AGN with z < 0.3 detected via the XMM-Newton and Chandra surveys. For comparison, we assembled a control sample of 71 star-forming (SF) galaxies with similar magnitudes, sizes, and redshifts. To derive the global properties of both the AGN and the control SF sample, we used CIGALE to fit the spectral energy distributions derived from the 56 narrowband and 4 broadband filters from miniJPAS. We find that AGN tend to reside in more massive galaxies than their SF counterparts. After matching samples based on stellar mass and comparing their SFRs and specific SFRs (sSFRs), no significant differences appear. This suggests that the presence of AGN does not strongly influence overall star formation. However, when we used miniJPAS as an integral field unit (IFU) to dissect galaxies along their position angle, a different picture emerges. We find that AGN tend to be more centrally concentrated in mass with respect to SF galaxies. Moreover, we find a suppression of the sSFR up to 1Re and then an enhancement beyond 1Re , strongly contrasting with the decreasing radial profile of sSFRs in SF galaxies. This could point to an inside-out quenching of AGN host galaxies. These findings suggest that the reason we do not see differences on a global scale is because star formation is suppressed in the central regions and enhanced in the outer regions of AGN host galaxies. While limited in terms of sample size, this work highlights the potential of the upcoming J-PAS as a wide-field low-resolution IFU for thousands of nearby galaxies and AGN., Comment: 21 pages, 20 figures, Accepted for Publication in A&A

We present JWST/NIRSpec high-resolution spectroscopy G395H/F290LP of MACS0647-JD, a gravitationally lensed galaxy merger at $z=10.167$. The new spectroscopy, which is acquired for the two lensed images (JD1 and JD2), detects and resolves emission lines in the rest-frame ultraviolet (UV) and blue optical, including the resolved [OII]3726,3729 doublet, [NeIII]3870, [HeI]3890, H$\delta$, H$\gamma$, and [OIII]4363. This is the first observation of the resolved [OII]3726,3729 doublet for a galaxy at $z>8$. We measure a line flux ratio [OII]3729/3726 $= 0.9 \pm 0.3$, which corresponds to an estimated electron density of $\log(n_{e} / \rm{cm}^{-3}) = 2.9 \pm 0.5$. This is significantly higher than the electron densities of local galaxies reported in the literature. We compile the measurements from the literature and further analyze the redshift evolution of $n_{e}$. We find that the redshift evolution follows the power-law form of $n_{e} = A\times (1+z)^{p}$ with $A=54^{+31}_{-23}$ cm$^{-3}$ and $p=1.2^{+0.4}_{-0.4}$. This power-law form may be explained by a combination of metallicity and morphological evolution of galaxies, which become, on average, more metal-poor and more compact with increasing redshift., Comment: 20 pages, 11 figures. Accepted for publication in ApJ

JWST spectroscopy has revolutionized our understanding of galaxies in the early universe. Covering wavelengths up to $5.3\,{\rm \mu m}$, NIRSpec can detect rest-frame optical emission lines H$\alpha$ out to $z = 7$ and [O III] to $z = 9.5$. Observing these lines in more distant galaxies requires longer wavelength spectroscopy with MIRI. Here we present MIRI MRS IFU observations of the lensed galaxy merger MACS0647$-$JD at $z = 10.165$. With exposure times of 4.2 hours in each of two bands, we detect H$\alpha$ at $9\sigma$, [O III]$\,\lambda5008$ at $11\sigma$, and [O III]$\,\lambda4960$ at $3\sigma$. Combined with previously reported NIRSpec spectroscopy that yields seven emission lines including the auroral line [O III]$\,\lambda4363$, we present the first direct metallicity measurement of a $z > 10$ galaxy: $12+{\rm log(O/H)}= 7.79\pm0.09$, or $0.13^{+0.02}_{-0.03}\,Z_{\odot}$. This is similar to galaxies at $z \sim 4 - 9$ with direct metallicity measurements, though higher than expected given the high specific star formation rate ${\rm log(sSFR / yr^{-1})} = -7.4 \pm 0.3$. We further constrain the ionization parameter ${\rm log}(U)$ = $-1.9 \pm 0.1$, ionizing photon production efficiency ${\rm log}(\xi_{\rm ion})$ = $25.3\pm0.1$, and star formation rate $5.0\pm0.6\,M_{\odot}/{\rm yr}$ within the past $10\,{\rm Myr}$. These observations demonstrate the combined power of JWST NIRSpec and MIRI for studying galaxies in the first $500$ million years., Comment: 16 pages, 8 figures, accepted for publication in ApJ

We present recent JWST NIRCam imaging observations of SPT0615-JD (also known as the Cosmic Gems Arc), lensed by the galaxy cluster SPT-CL J0615-5746. The 5-arcsec-long arc is the most highly magnified $z>10$ galaxy known, straddling the lensing critical curve and revealing five star clusters with radii $\sim 1$ pc or less. We measure the full arc to have F200W 24.5 AB mag, consisting of two mirror images, each 25.3 AB mag with a magnification $\mu \sim 60$ (delensed 29.7 AB mag, $M_{UV} = -17.8$). The galaxy has an extremely strong Lyman break F115W$-$F200W $>3.2$ mag ($2\sigma$ lower limit), is undetected in all bluer filters ($< 2\sigma$), and has a very blue continuum slope redward of the break ($\beta = -2.7 \pm 0.1$), resulting in a photometric redshift $z_{phot} = 10.2 \pm 0.2$ (95% confidence) with no significant likelihood below $z < 9.8$. Based on SED fitting to the total photometry, we estimate an intrinsic stellar mass of $M_{*} \sim 2.4 - 5.6 \times 10^{7} M_{\odot}$, young mass-weighted age of $\sim 21 - 79$ Myr, low dust content ($A_V < 0.15$), and a low metallicity of $\lesssim 1\%~Z_{\odot}$. We identify a fainter third counterimage candidate within 2.2 arcsec of the predicted position, lensed to AB mag 28.4 and magnified by $\mu \sim 2$, suggesting the fold arc may only show $\sim60$% of the galaxy. SPT0615-JD is a unique laboratory to study star clusters observed within a galaxy just 460 Myr after the Big Bang., Comment: 22 pages, 8 figures, 4 tables, submitted to ApJ

One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$\alpha$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$\alpha$ absorption (DLA) wings from high column densities of neutral atomic hydrogen (HI), signifying major gas accretion events in the formation of these galaxies. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at $z=5.5-13.4$. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly$\alpha$ damping parameter $D_{\rm Ly\alpha}$ to measure and quantify the Ly$\alpha$ emission strength, HI fraction in the IGM, or local HI column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories. The fraction of strong galaxy DLAs are found to be in the range $65-95\%$ at $z>5.5$. The fraction of strong Ly$\alpha$ emitters (LAEs) is found to increase with decreasing redshift, in qualitative agreement with previous observational results, and are predominantly associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine HI gas accretion. [abridged], Comment: 20 pages, 18 figures, 3 tables. Submitted to A&A. Comments welcome! All data and catalogs are available through the DAWN JWST Archive (DJA): https://dawn-cph.github.io/dja/ and https://github.com/keheintz/jwst-primal

Galaxy sizes correlate with many other important properties of galaxies, and the cosmic evolution of galaxy sizes is an important observational diagnostic for constraining galaxy evolution models. The effective radius is probably the most widely used indicator of galaxy size. We used the TNG50-SKIRT Atlas to investigate the wavelength dependence of the effective radius of galaxies at optical and near-infrared (NIR) wavelengths. We find that, on average, the effective radius in every band exceeds the stellar mass effective radius, and that this excess systematically decreases with increasing wavelength. The optical g-band (NIR Ks-band) effective radius is on average 58% (13%) larger than the stellar mass effective radius. Effective radii measured from dust-obscured images are systematically larger than those measured from dust-free images, although the effect is limited (8.7% in the g-band, 2.1% in the Ks-band). We find that stellar population gradients are the dominant factor (about 80%) in driving the wavelength dependence of the effective radius, and that differential dust attenuation is a secondary factor (20%). Comparing our results to recent observational data, we find offsets in the absolute values of the median effective radii, up to 50% for the population of blue galaxies. We find better agreement in the slope of the wavelength dependence of the effective radius, with red galaxies having a slightly steeper slope than green-blue galaxies. Comparing our effective radii with those of galaxies from the Siena Galaxy Atlas in separate bins in z-band absolute magnitude and g-z colour, we find excellent agreement for the reddest galaxies, but again significant offsets for the blue populations: up to 70% for galaxies around Mz=-21.5. This difference in median effective radius for the bluer galaxies is most probably due to (abridged...)., Comment: 11 pages, 10 figures, accepted for publication in A&A

Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a synthetic UV to near-infrared broadband image atlas for a complete stellar-mass selected sample of 1154 galaxies extracted from the TNG50 cosmological simulation at $z=0$. The images have a high spatial resolution (100 pc) and a wide field of view (160 kpc). In addition to the dust-obscured images, we also release dust-free images and physical parameter property maps with matching characteristics. As a sanity check and preview application we discuss the UVJ diagram of the galaxy sample. We investigate the effect of dust attenuation on the UVJ diagram and find that it affects both the star-forming and the quiescent galaxy populations. The quiescent galaxy region is polluted by younger and star-forming highly inclined galaxies, while dust attenuation induces a separation in inclination of the star-forming galaxy population, with low-inclination galaxies remaining at the blue side of the diagram and high-inclination galaxies systematically moving towards the red side. This image atlas can be used for a variety of other applications, including galaxy morphology studies and the investigation of local scaling relations. We publicly release the images and parameter maps, and we invite the community to use them., Comment: 13 pages, 7 figures, accepted for publication in A&A

The Cosmic Gems arc is among the brightest and highly magnified galaxies observed at redshift $z\sim10.2$. However, it is an intrinsically UV faint galaxy, in the range of those now thought to drive the reionization of the Universe. Hitherto the smallest features resolved in a galaxy at a comparable redshift are between a few hundreds and a few tens of parsecs. 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 parsec. They exhibit minimal dust attenuation and low metallicity, ages younger than 50 Myr and intrinsic masses of $\sim10^6$ M$_{\odot}$. Their lensing-corrected sizes are approximately 1 pc, resulting in stellar surface densities near $10^5$~M$_{\odot}$/pc$^2$, three orders of magnitude higher than typical young star clusters in the local universe. Despite the uncertainties inherent to the lensing model, they are consistent with being gravitationally bound stellar systems, i.e., proto-globular clusters. We conclude that star cluster formation and feedback likely contributed to shape the properties of galaxies during the epoch of reionization. [Abridged], Comment: Accepted for publication

We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of $z_{\mathrm{phot}}\simeq4.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 $z_{\mathrm{spec}}=4.758\pm0.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 $T_{\mathrm{eff,B}}\simeq15000$ K with either a redder F-type companion ($T_{\mathrm{eff,F}}\simeq6250$K) or significant dust attenuation ($A_V\simeq0.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\lesssim0.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 $\sim1.2$ Gyr old showcases the ability of JWST to study early stars through extreme lensing., Comment: Accepted for publication in MNRAS letters. v2 updated to match the published version

The onset of galaxy formation is thought to be initiated by the infall of neutral, pristine gas onto the first protogalactic halos. However, direct constraints on the abundance of neutral atomic hydrogen (HI) in galaxies have been difficult to obtain at early cosmic times. Here we present spectroscopic observations with JWST of three galaxies at redshifts $z=8.8 - 11.4$, about $400-600$ Myr after the Big Bang, that show strong damped Lyman-$\alpha$ absorption ($N_{\rm HI} > 10^{22}$ cm$^{-2}$) from HI in their local surroundings, an order of magnitude in excess of the Lyman-$\alpha$ absorption caused by the neutral intergalactic medium at these redshifts. Consequently, these early galaxies cannot be contributing significantly to reionization, at least at their current evolutionary stages. Simulations of galaxy formation show that such massive gas reservoirs surrounding young galaxies so early in the history of the universe is a signature of galaxy formation in progress., Comment: Submitted

We present JWST/NIRSpec prism spectroscopy of MACS0647-JD, the triply-lensed $z \sim 11$ candidate discovered in HST imaging and spatially resolved by JWST imaging into two components A and B. Spectroscopy of component A yields a spectroscopic redshift $z=10.17$ based on 7 detected emission lines: CIII] $\lambda\lambda$1907,1909, [OII] $\lambda$3727, [NeIII] $\lambda$3869, [NeIII] $\lambda$3968, H$\delta$ $\lambda$4101, H$\gamma$ $\lambda$4340, and [OIII] $\lambda$4363. These are the second-most distant detections of these emission lines to date, in a galaxy observed just 460 million years after the Big Bang. Based on observed and extrapolated line flux ratios we derive a gas-phase metallicity $Z =$ log(O/H) = $7.5 - 8.0$, or $(0.06 - 0.2)$ $Z_\odot$, ionization parameter log($U$) $\sim -1.9\pm0.2$, and an ionizing photon production efficiency ${\rm log}(\xi_{\rm ion})=25.2\pm0.2\,$erg$^{-1}$ Hz. The spectrum has a softened Lyman-$\alpha$ break, evidence for a strong Ly$\alpha$ damping wing, suggesting that MACS0647-JD was unable to ionize its surroundings beyond its immediate vicinity ($R_{\text{HII}} \ll 1$ pMpc). The Ly$\alpha$ damping wing also suppresses the F150W photometry, explaining the slightly overestimated photometric redshift $z = 10.6 \pm 0.3$. MACS0647-JD has a stellar mass log($M/M_\odot$) = $8.1 \pm 0.3$, including $\sim$ 6$\times 10^7 M_\odot$ in component A, most of which formed recently (within $\sim$ 20 Myr) with a star formation rate $2\pm1 M_\odot$ / yr, all within an effective radius $70\pm24\,$pc. The smaller component B ($r \sim 20$) pc is likely older ($\sim$100 Myr) with more dust ($A_V \sim 0.1$ mag), as found previously. Spectroscopy of a fainter companion galaxy C separated by a distance of \about\ 3$\,$kpc reveals a Lyman break consistent with $z = 10.17$. MACS0647-JD is likely the most distant galaxy merger known., Comment: 22 pages, 6 figures, accepted for publication in ApJ

We present new observations of MRG-M2129, a quiescent galaxy at z = 2.15 with the Atacama Large Millimeter/submillimeter Array (ALMA). With the combination of the gravitational lensing effect by the foreground cluster and the angular resolution provided by ALMA, our data reveal 1.2 mm continuum emission at $\sim130$ pc angular resolution. Compact dust continuum is detected at 7.9 $\sigma$ in the target but displaced from its stellar peak position by $62 \pm 38$ mas, or $\sim169 \pm 105$ pc in the source plane. We find considerably high dust-to-stellar mass ratio, $4 \times 10^{-4}$. From non-detection of the [C i] 3P2 -> 3P1 line, we then derive $3 \sigma$ upper limits on the molecular gas-to-dust mass ratio $\delta_\mathrm{GDR} < 60$ and the molecular gas-to-stellar mass ratio fH2 < 2.3%. The derived $\delta_\mathrm{GDR}$ is $>2\times$ smaller than the typical value assumed for quiescent galaxies in the literature. Our study supports that there exists a broad range of $\delta_\mathrm{GDR}$ and urges submillimeter follow-up observations of quenching/recently quenched galaxies at similar redshifts. Based on the derived low $\delta_\mathrm{GDR}$ and other observed dust properties, we argue that the central black hole is still active and regulates star formation in the system. Our study exhibits a rare case of a gravitationally lensed type 2 QSO harbored by a quiescent galaxy., Comment: Accepted for publication in ApJ

The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$\mu$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 $\mu > 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 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_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, 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., Comment: Accepted to ApJL. Data products, lens models, and analysis code will be available online at https://cosmic-spring.github.io

We present a multiwavelength study of IC 860, a nearby post-starburst galaxy at the early stage of transitioning from blue and star-forming to red and quiescent. Optical images reveal a galaxy-wide, dusty outflow originating from a compact core. We find evidence for a multiphase outflow in the molecular and neutral gas phase from the CO position-velocity diagram and NaD absorption features. We constrain the neutral mass outflow rate to be ~0.5 M$_{\odot}/$yr, and the total hydrogen mass outflow rate to be ~12 M$_{\odot}$/yr. Neither outflow component seems able to escape the galaxy. We also find evidence for a recent merger in the optical images, CO spatial distribution, and kinematics, and evidence for a buried AGN in the optical emission line ratios, mid-IR properties, and radio spectral shape. The depletion time of the molecular gas reservoir under the current star formation rate is ~7 Gyr, indicating that the galaxy could stay at the intermediate stage between the blue and red sequence for a long time. Thus the timescales for a significant decline in star formation rate ("quenching") and gas depletion are not necessarily the same. Our analysis supports the quenching picture where outflows help suppress star formation by disturbing rather than expelling the gas and shed light on possible ongoing activities in similar quenching galaxies., Comment: 31 pages, 15 figures, accepted for publication in The Astrophysical Journal. Online material will be available from the publisher

We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro'uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($\Sigma_{\rm H_{2}}$--$\Sigma_{\rm SFR}$--$\Sigma_{*}$) and dust scaling relations. While the relations with all sub-galactic regions of the galaxies are reasonably tight ($\sigma \lesssim 0.3$ dex), we find that most of the scaling relations exhibit galaxy-to-galaxy variations in normalization and shape. Only two relations of $\Sigma_{\rm dust}$--$\Sigma_{\rm gas}$ and $\Sigma_{\rm dust}$--$\Sigma_{\rm H_{2}}$ do not show noticeable galaxy-to-galaxy variations among our sample galaxies. We further investigate correlations among the scaling relations. We find significant correlations among the normalization of the $\Sigma_{\rm H_{2}}$--$\Sigma_{\rm SFR}$--$\Sigma_{*}$ relations, which suggest that galaxies with higher levels of resolved $\text{H}_{2}$ fraction ($f_{\rm H_{2}}$) tend to have higher levels of resolved star formation efficiency (SFE) and specific star formation rate (sSFR). We also observe that galaxies with higher levels of resolved dust-to-stellar mass ratios tend to have higher levels of resolved sSFR, SFE, and $f_{\rm H_{2}}$. Moreover, we find that galaxies with higher global sSFR and less compact morphology tend to have higher levels of the resolved sSFR, SFE, and $f_{\rm H_{2}}$, which can explain the variations in the normalization of the $\Sigma_{\rm H_{2}}$--$\Sigma_{\rm SFR}$--$\Sigma_{*}$ relationships. Overall, we observe indications of the contributions of both global and local factors in governing the star formation process in galaxies., Comment: 28 pages, 13 figures, accepted for publication in ApJ. piXedfit v1.0 is publicly available at https://github.com/aabdurrouf/piXedfit. Documentation is available at https://pixedfit.readthedocs.io/en/latest/. Some tutorials of practical usages in jupyter notebooks can be found at https://github.com/aabdurrouf/piXedfit/tree/main/examples

This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys., Comment: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

We study spatially resolved properties (on spatial scales of $\sim 1-2$ kpc out to at least $3$ effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX--$\beta$ relation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kpc scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center., Comment: 35 pages, 19 figures, 4 tables. Accepted by ApJ (results remain the same). Complementary figures to Fig. 4 (for 10 galaxies) and high resolution version of Fig. 8 are available at https://github.com/aabdurrouf/maps-nearby-galaxies

We present piXedfit, pixelized spectral energy distribution (SED) fitting, a Python package that provides tools for analyzing spatially resolved properties of galaxies using multiband imaging data alone or in combination with integral field spectroscopy (IFS) data. piXedfit has six modules that can handle all tasks in the spatially resolved SED fitting. The SED fitting module uses the Bayesian inference technique with two kinds of posteriors sampling methods: Markov Chain Monte Carlo (MCMC) and random densely-sampling of parameter space (RDSPS). We test the performance of the SED fitting module using mock SEDs of simulated galaxies from IllustrisTNG. The SED fitting with both posteriors sampling methods can recover physical properties and star formation histories of the IllustrisTNG galaxies well. We further test the performance of piXedfit modules by analyzing 20 galaxies observed by the CALIFA and MaNGA surveys. The data comprises of 12-band imaging data from GALEX, SDSS, 2MASS, and WISE, and the IFS data from CALIFA or MaNGA. piXedfit can spatially match (in resolution and sampling) of the imaging and IFS data. By fitting only the photometric SEDs, piXedfit can predict the spectral continuum, $\text{D}_{\rm n}4000$, $H_{\alpha}$, and $H_{\beta}$ well. The star formation rate (SFR) derived by piXedfit is consistent with that derived from $H_{\alpha}$ emission. The RDSPS method gives equally good fitting results as the MCMC and it is much faster than the MCMC. piXedfit is a versatile tool equipped with a parallel computing module for efficient analysis of large datasets, and will be made publicly available (https://github.com/aabdurrouf/piXedfit)., Comment: 40 pages, 21 figures, 1 table, Accepted for publication in The Astrophysical Journal Supplement Series (ApJS). piXedfit will be available at https://github.com/aabdurrouf/piXedfit and documented at https://pixedfit.readthedocs.io/en/latest/index.html

We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 5 observations of a galaxy at $z=1.91$, GDS24569, in search of molecular gas in its vicinity via the [C I] $^3$P$_1$-$^3$P$_0$ line. GDS24569 is a massive ($\log M_*/M_\odot=11$) passively evolving galaxy, and characterized by compact morphology with an effective radius of $\sim0.5$ kpc. We apply two blind detection algorithms to the spectral data cubes, and find no promising detection in or around GDS24569 out to projected distance of $\sim320$ kpc, while a narrow tentative line ($4.1 \sigma$) is identified at $+1200$ km/s by one of the algorithms. From the non-detection of [C I], we place a $3\sigma$ upper limit on molecular hydrogen mass, $\sim 7.1 \times 10^9 M_\odot$, which converts to an extremely low gas-to-stellar mass fraction, $< 5 \%$. We conduct a spectral energy distribution modeling by including optical-to-far-infrared data, and find a considerably high ($\sim0.1\%$) dust-to-stellar mass ratio, $\sim10$-$100\times$ higher than those of local early-type galaxies. In combination with a previous result of an insufficient number of surrounding satellite galaxies, it is suggested that GDS24569 is unlikely to experience significant size evolution via satellite mergers. We discuss possible physical mechanisms that quenched GDS24569., Comment: Accepted in ApJ

We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs., Comment: 28 pages, 25 figures, 7 tables; accepted by ApJ

Despite decreasing cosmic star formation rate density over the last 10 Gyr, the stellar mass ($M_{*}$) buildups in galaxies were still progressing during this epoch. About 50\% of the current $M_{*}$ density in the universe was built over the last $\sim 8.7$ Gyr. In this research, we investigated the stellar mass buildup and quenching of spatially resolved regions within massive disk galaxies over the last 10 Gyr. We apply the spectral energy distribution (SED) fitting method to SEDs of sub-galactic regions in galaxies to derive the spatially resolved distributions of SFR and $M_{*}$ in the galaxies. This namely \textit{pixel-to-pixel SED fitting} method is applied to massive disk galaxies at $0.01

We investigate a relation between surface densities of star formation rate (SFR) and stellar mass ($M_{*}$) at a $\sim 1$ kpc scale namely spatially resolved star formation main sequence (SFMS) in massive ($\log(M_{*}/M_{\odot})>10.5$) face-on disc galaxies at $0.01

We investigate the relation between star formation rate (SFR) and stellar mass ($M_*$) in sub-galactic ($\sim 1$kpc) scale of 93 local ($0.0110^{10.5}M_{\odot}$) spiral galaxies. To derive spatially-resolved SFR and stellar mass, we perform so-called pixel-to-pixel SED fitting, which fits an observed spatially-resolved multiband SED with a library of model SEDs using Bayesian statistics approach. We use 2 bands (FUV and NUV) and 5 bands ($u$, $g$, $r$, $i$, and $z$) imaging data from Galaxy Evolution Explorer (GALEX) and Sloan Digital Sky Survey (SDSS), respectively. We find a tight nearly linear relation between the local surface density of SFR ($\Sigma_{\rm{SFR}}$) and stellar mass ($\Sigma_{*}$) which has flattening in high $\Sigma_{*}$. The near linear relation between $\Sigma_{*}$ and $\Sigma_{\rm SFR}$ suggests constant sSFR throughout the galaxies, and the scatter of the relation is directly related to that of sSFR. Therefore, we analyse the variation of sSFR in various scales. More massive galaxies on average have lower sSFR throughout them than less massive galaxies. We also find that barred galaxies have lower sSFR in a core region than non-barred galaxies. However, in the outside region, sSFR of barred and non-barred galaxies are similar and lead to the similar total sSFR., Comment: 16 pages, 19 figures. Accepted for publication in MNRAS