Search

We present results from the JWST First Reionization Epoch Spectroscopically Complete Observations survey on the star-forming sequence (SFS) of galaxies at 1.0 < z < 1.7, around the peak of the cosmic star formation history. Star formation rates (SFRs) are measured from the redshifted, relatively dust-insensitive Paschen- α emission line, and stellar mass measurements include the F444W (4.4 μ m; rest-frame H) band. We find SFRs of galaxies with log(M */ M ⊙) > 9.5 that are lower than found in many earlier studies by up to 0.6 dex, but in good agreement with recent results obtained with the Prospector fitting framework. The difference (log(SFR(Pa α)-SFR(Prospector)) is −0.09 ± 0.04 dex at 1010−11 M ⊙. We also measure the empirical relation between Paschen- α luminosity and rest-frame H- band magnitude and find that the scatter is only 0.04 dex lower than that of the SFR– M * relation and is much lower than the systematic differences among relations in the literature due to various methods of converting observed measurements to physical properties. We additionally identify examples of sources—that, with standard cutoffs via the UVJ diagram, would be deemed quiescent—with significant (log(sSFR)> −11 yr−1), typically extended, Paschen- α emission. Our results may be indicative of the potential unification of methods used to derive the SFS with careful selection of star-forming galaxies and independent SFR and stellar mass indicators. [ABSTRACT FROM AUTHOR]

This theoretical work investigates different models to predict the redshift of fast radio bursts (FRBs) from their observed dispersion measure (DM) and other reported properties. We performed an extended revision of the FRBs with confirmed galaxy hosts in the literature and built the most updated catalog to date. With this sample of FRBs, we propose four models that relate the DM and z: a linear trend (inspired by the Macquart relation), a log-parabolic function, a power law, and an interpolation for DM that includes z and the position of the host galaxy of the transient. The latter model has the highest success rate according to the metrics implemented: likelihood, median of the z difference, and the Akaike and Bayesian information criteria. Although the performance of model D is closely followed by the power law and linear models, the former has the advantage of accounting for anisotropies in the dispersion measure due to the angular coordinates. Conversely, the log-parabolic formula performs poorly in this task but provides a good prediction for FRBs with low DM at a low redshift. Additionally, we use the reported galaxy properties of the hosts to establish a connection between the observed DM with the star formation rate (SFR) and stellar mass (Ms) of the galaxies where the FRBs reside. In both cases, we find a weak correlation. Although the studied correlations are well-motivated, the sample of FRBs is not statistically significant enough to draw solid conclusions in this second part of our work. With the advent of new facilities devoted to studying the localization and nature of these transients, we will get access to new data that will enrich the proposed models and give us hints on the astrophysical origin and evolution of FRBs. [ABSTRACT FROM AUTHOR]

This theoretical work investigates different models to predict the redshift of fast radio bursts (FRBs) from their observed dispersion measure (DM) and other reported properties. We performed an extended revision of the FRBs with confirmed galaxy hosts in the literature and built the most updated catalog to date. With this sample of FRBs, we propose four models that relate the DM and z: a linear trend (inspired by the Macquart relation), a log-parabolic function, a power law, and an interpolation for DM that includes z and the position of the host galaxy of the transient. The latter model has the highest success rate according to the metrics implemented: likelihood, median of the z difference, and the Akaike and Bayesian information criteria. Although the performance of model D is closely followed by the power law and linear models, the former has the advantage of accounting for anisotropies in the dispersion measure due to the angular coordinates. Conversely, the log-parabolic formula performs poorly in this task but provides a good prediction for FRBs with low DM at a low redshift. Additionally, we use the reported galaxy properties of the hosts to establish a connection between the observed DM with the star formation rate (SFR) and stellar mass (Ms) of the galaxies where the FRBs reside. In both cases, we find a weak correlation. Although the studied correlations are well-motivated, the sample of FRBs is not statistically significant enough to draw solid conclusions in this second part of our work. With the advent of new facilities devoted to studying the localization and nature of these transients, we will get access to new data that will enrich the proposed models and give us hints on the astrophysical origin and evolution of FRBs.

Understanding galaxy properties may be the key to unlocking some of the most intriguing mysteries of modern cosmology. Recent work relied on machine learning to extract cosmological constraints on Ω _m using only one galaxy. But if this is true, how should we select the galaxy to use for cosmology inference? In this Letter, we consider selecting a galaxy that lies in cosmic voids, the underdense regions of the cosmic web, and compare the constraints obtained with the ones obtained when randomly selecting a galaxy in the whole sample. We use the IllustrisTNG galaxy catalog from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and the VIDE void finder to identify galaxies inside voids. We show that void galaxies provide stronger constraints on Ω _m compared to randomly selected galaxies. This result suggests that the distinctive characteristics of void galaxies may provide a cleaner and more effective environment for extracting cosmological information.

Photometrically derived stellar masses are known to suffer from systematic uncertainties, particularly due to nebular emission contributions to the spectral energy distribution (SED). Using James Webb Space Telescope NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey, we introduce a comparison study of photometrically derived redshifts and stellar masses based on two photometric catalogs of the same field spanning ∼0.4–4.5 μ m: one consisting solely of wide band photometry, and another employing a combination of wide and medium band photometry. We find that ∼70% of likely line emitters have consistent photometric redshifts between both catalogs, with a median stellar mass difference between the two catalogs of 5. These systematic differences caused by the poor spectral resolution of wide bands have implications for both ongoing and future planned observing programs that determine stellar mass and other physical properties of high-redshift galaxies solely via wide band photometry.

We explore the morphological features and star formation activities of [O ii ] emitters in the COSMOS UltraDeep field at z ∼ 1.5 using JWST NIRCam data from the COSMOS-Web survey and Subaru Hyper Suprime-Cam. We also report the discovery of large filamentary structures traced by [O ii ] emitters surrounding an extremely overdense core with a galaxy number density ∼11× higher than the field average. These structures span over 50 cMpc, underscoring their large scale in the cosmic web at this epoch. After matching the stellar-mass distributions, the core galaxies show a higher frequency of disturbances (50% ± 9%) than those in the outskirts (41% ± 9%) and the field (21% ± 5%), indicative of more frequent mergers and interactions in the innermost ≲1.′5 region. Additionally, we observe that specific star formation rates are elevated in denser environments. A Kolmogorov–Smirnov test comparing the distribution of specific star formation rates of core and field galaxies yields a p -value of 0.02, suggesting an enhancement of star formation activity driven by the dense environment. Our findings underscore the environmental impact on galaxy evolution during a pivotal cosmic epoch and set the stage for further investigation with the increasing larger data from upcoming surveys.

We study a sample of 14 spectroscopically confirmed Ly α emitters (LAEs) in the late era of reionization (at redshift z ≈ 6) based on the JWST/NIRCam imaging data set. These LAEs with high Ly α luminosity of L (Ly α ) ∼10 ^42.4 –10 ^43.4 erg s ^−1 have been covered by the (ongoing) COSMOS-Web survey over 0.28 deg ^2 in four NIRCam bands (F115W, F150W, F277W, and F444W). With JWST/NIRCam imaging, we determine the UV continua with M _UV ranging from −20.5 to −18.5 mag. The UV slopes have a median value of β ≈ −2.35, and the steepest slopes can reach β < −3. Under the excellent spatial resolution of JWST, we identify three objects in the sample as potential merging/interacting systems. The 14 LAEs (and their components) are compact in morphology, residing substantially below the mass–size relation of high- z galaxies. We further investigate their physical properties, including the stellar mass ( M _* ) and star formation rates (SFRs). Most of the LAEs lie on the SFR– M _* main-sequence relation, while two of them, featured as “little red dots,” likely host active galactic nuclei (AGNs), implying a ∼10% AGN fraction. Moreover, we reveal that a new correlation may exist between Ly α equivalent width and the offset between Ly α and UV emission (Δ d _Ly _α ), with a median Δ d _Ly _α ∼ 1 kpc. This could be explained by the Ly α radiative transfer process in both the interstellar medium and circumgalactic medium. The results usher in a new era of detailed analysis on high- z LAEs with the JWST capability.

We present a study of Lyman continuum (LyC) emission in a sample of ∼150 Ly α emitters (LAEs) at z ≈ 3.1 in the Subaru-XMM Deep Survey field. These LAEs were previously selected using the narrowband technique and spectroscopically confirmed with Ly α equivalent widths (EWs) ≥ 45 Å. We obtain deep UV images using a custom intermediate-band filter U J that covers a wavelength range of 3330 ∼ 3650 Å, corresponding to 810 ∼ 890 Å in the rest frame. We detect five individual LyC galaxy candidates in the U J band, and their escape fractions (f esc) of LyC photons are roughly between 40% and 80%. This supports a previous finding that a small fraction of galaxies may have very high f esc. We find that the f esc values of the five LyC galaxies are not apparently correlated with other galaxy properties such as Ly α luminosity and EW, UV luminosity and slope, and star formation rate (SFR). This is partly due to the fact that these galaxies only represent a small fraction (∼3%) of our LAE sample. For the remaining LAEs that are not detected in U J, we stack their U J-band images and constrain their average f esc. The upper limit of the average f esc value is about 16%, consistent with the results in the literature. Compared with the non-LyC LAEs, the LyC LAEs tend to have higher Ly α luminosities, Ly α EWs, and SFRs, but their UV continuum slopes are similar to those of other galaxies. [ABSTRACT FROM AUTHOR]

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, the basic demographics of clumps and their host galaxies remain uncertain. Using the Hubble Space Telescope/Wide Field Camera 3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, we detect and analyze giant star-forming clumps in galaxies at 0.5 ≤ z ≤ 1, connecting two epochs when clumps are common (at cosmic high noon, z ∼ 2) and rare (in the local Universe). We construct a clump sample whose rest-frame 1600 Å luminosity is 3 times higher than the most luminous local H ii regions (M UV ≤ −16 AB). In our sample, 35% ± 3% of low-mass galaxies (log[ M ∗/ M ⊙] < 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22% ± 3% and 22% ± 4% for intermediate (10 ≤ log[ M ∗/ M ⊙] ≤ 10.5) and high-mass (log[ M ∗/ M ⊙] > 10.5) galaxies, in agreement with previous studies. When compared to similar-mass nonclumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher star formation rates (SFRs) and bluer rest-frame U − V colors, while high-mass clumpy SFGs tend to be larger than nonclumpy SFGs. However, clumpy and nonclumpy SFGs have similar Sérsic index, indicating a similar underlying density profile. Furthermore, we investigate how the UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher specific SFR galaxies and are found closer to their hosts' galactic centers. [ABSTRACT FROM AUTHOR]

Recent observations have reignited interest in a population of dwarf galaxies that are large and diffuse for their mass, often called ultra-diffuse galaxies (UDGs). However, the origin and evolution of these mass–size outliers and the role of the environment are still unclear. Using the exquisitely deep and wide Hyper Suprime-Cam Strategic Survey images, we search for ultra-puffy galaxies (UPGs), defined as being 1.5 σ larger than the average size for their mass, around Milky Way–like galaxies. We present the sizes and radial distributions of mass–size outliers and derive their quenched fraction to explore the impact of the environment. Surprisingly, despite being outliers in size, the UPGs have a similar quenched fraction as normal-sized satellites of Milky Way analogs in both observations and simulations, suggesting that quenching is not tied to being a mass–size outlier. The quenched fraction is higher for the UPGs associated with redder hosts, as well as those that are closer to the host in projection. In contrast, the UDGs are overall redder and more quiescent compared with normal satellites. We show that the classic definition of UDGs is heavily weighted toward quenched galaxies and thus cannot be used for a study of quenching of mass–size outliers. [ABSTRACT FROM AUTHOR]

Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru's Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and modeling them. As a case study, we surveyed 922 Milky Way analogs in the nearby Universe (0.01 < z < 0.04) and built a large sample of satellite galaxies that are outliers in the mass–size relation. These "ultra-puffy" galaxies (UPGs), defined to be 1.5 σ above the average mass–size relation, represent the tail of the satellite size distribution. We find that each MW analog hosts N UPG = 0.31 ± 0.05 UPGs on average, which is consistent with but slightly lower than the observed abundance at this halo mass in the Local Volume. We also construct a sample of ultra-diffuse galaxies (UDGs) in MW analogs and find an abundance of N UDG = 0.44 ± 0.05 per host. With literature results, we confirm that the UDG abundance scales with the host halo mass following a sublinear power law. We argue that our definition of UPGs, which is based on the mass–size relation, is more physically motivated than the common definition of UDGs, which depends on the surface brightness and size cuts and thus yields different surface mass density cuts for quenched and star-forming galaxies. [ABSTRACT FROM AUTHOR]

We measure the galaxy-ellipticity (GI) correlations for the Sloan Digital Sky Survey Data Release 12 LOWZ and CMASS samples with the shape measurements from the DESI Legacy Imaging Surveys. We model the GI correlations in an N -body simulation with our recent accurate stellar–halo mass relation from the Photometric object Around Cosmic webs (PAC) method. The large data set and our accurate modeling turns out an accurate measurement of the alignment angle between central galaxies and their host halos. We find that the alignment of central elliptical galaxies with their host halos increases monotonically with galaxy stellar mass or host halo mass, which can be well described by a power law for the massive galaxies. We also find that central elliptical galaxies are more aligned with their host halos in LOWZ than in CMASS, which might indicate an evolution of galaxy–halo alignment, though future studies are needed to verify this is not induced by the sample selections. In contrast, central disk galaxies are aligned with their host halos about 10 times more weakly in the GI correlation. These results have important implications for intrinsic alignment (IA) correction in weak lensing studies, IA cosmology, and theory of massive galaxy formation. [ABSTRACT FROM AUTHOR]

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at 1.3 ≤ z ≤ 2.6 from the MOSFIRE Deep Evolution Field survey. We divide galaxies with a detected H α emission line and coverage of H β into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), and we then stack their spectra. From each stack, we measure the Balmer decrement and gas-phase metallicity, and then we compute the median A V and UV continuum spectral slope (β). First, we find that none of the dust properties (Balmer decrement, A V, or β) varies with the axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model in which attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases the SFR and lowers the metallicity, but leaves the dust column density mostly unchanged. We quantify this idea using the Kennicutt–Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation. [ABSTRACT FROM AUTHOR]

We combine an unprecedented MaNGA sample of over 3000 passive galaxies in the stellar mass range 109–1012 M ⊙ with the Sloan Digital Sky Survey group catalog by Tinker to quantify how central and satellite formation, quantified by radial profiles in stellar age, [Fe/H], and [Mg/Fe], depends on the stellar mass of the galaxy (M *) and the mass of the host halo (M h ). After controlling for M * and M h , the stacked spectra of centrals and satellites beyond the effective radius (r e ) show small, yet significant differences in multiple spectral features at the 1% level. According to spectral fitting with the code alf, a primary driver of these differences appears to be [Mg/Fe] variations, suggesting that stellar populations in the outskirts of satellites formed more rapidly than the outer populations of centrals. To probe the physical mechanisms that may be responsible for this signal, we examined how satellite stellar populations depend on M h . We find that satellites in high- M h halos show older stellar ages, lower [Fe/H], and higher [Mg/Fe] compared to satellites in low- M h halos, especially for M * = 109.5–1010.5 M ⊙. These signals lend support to environmentally driven processes that quench satellite galaxies, although variations in the merger histories of central and satellite galaxies also emerge as a viable explanation. [ABSTRACT FROM AUTHOR]

Modern spectroscopic surveys can only target a small fraction of the vast amount of photometrically cataloged sources in wide-field surveys. Here, we report the development of a generative artificial intelligence (AI) method capable of predicting optical galaxy spectra from photometric broadband images alone. This method draws from the latest advances in diffusion models in combination with contrastive networks. We pass multiband galaxy images into the architecture to obtain optical spectra. From these, robust values for galaxy properties can be derived with any methods in the spectroscopic toolbox, such as standard population synthesis techniques and Lick indices. When trained and tested on 64 × 64 pixel images from the Sloan Digital Sky Survey, the global bimodality of star-forming and quiescent galaxies in photometric space is recovered, as well as a mass–metallicity relation of star-forming galaxies. The comparison between the observed and the artificially created spectra shows good agreement in overall metallicity, age, Dn4000, stellar velocity dispersion, and E ( B − V ) values. Photometric redshift estimates of our generative algorithm can compete with other current, specialized deep learning techniques. Moreover, this work is the first attempt in the literature to infer velocity dispersion from photometric images. Additionally, we can predict the presence of an active galactic nucleus up to an accuracy of 82%. With our method, scientifically interesting galaxy properties, normally requiring spectroscopic inputs, can be obtained in future data sets from large-scale photometric surveys alone. The spectra prediction via AI can further assist in creating realistic mock catalogs.

Investigating the impact of galaxy properties on emergent Ly α emission is crucial for reionization studies, given the sensitivity of Ly α to neutral hydrogen. This study presents an analysis of the physical characteristics of 155 star-forming galaxies, 29 with Ly α detected, and 126 with Ly α not detected with Ly α EW < 20 Å, at z = 1.9–3.5, drawn from the MOSFIRE Deep Evolution Field survey, that have overlapping observations from the Hobby–Eberly Telescope Dark Energy Experiment survey. To unravel the interstellar medium (ISM) conditions in our sample, we developed a custom nebular line modeling algorithm based on the MAPPINGS V photoionization model grid and the emcee framework. Combining nebular-based ISM properties with photometry-based global properties, constrained via Bagpipes , we explore distinctions in the stellar and gas properties between Ly α -detected and Ly α -nondetected galaxies. Our analysis reveals statistically significant differences between the two samples in terms of stellar mass and dust attenuation ( A _V ) at >2 σ significance, as determined via a Kolmogorov–Smirnov test. Moreover, there are weaker (≲1 σ significance) indications that the ionization parameter and metallicity differ between the two samples. Our results demonstrate that the escape fraction of Ly α ( ${f}_{\mathrm{esc}}^{\mathrm{Ly}\alpha }$ ) is inversely correlated with stellar mass, star formation rate, and dust attenuation, while it is positively correlated with the ionization parameter, with significance levels exceeding 2 σ . Our findings suggest that the interstellar environments of Ly α -detected galaxies, characterized by low mass, low dust, low gas-phase metallicity, and high ionization parameters, play a pivotal role in promoting the escape of Ly α radiation.

Lyman continuum (LyC) emitters at z > 3 provide critical samples for studying the contribution of galaxies to the ionizing background in the epoch of reionization. We collect a sample of z > 3 LyC emitters, a dominant fraction (∼60%–70%) of which show spatial offsets between LyC emission and the nonionizing continuum. From this sample, especially, we find a case of an offset LyC emitter, CDFS-6664 ( z = 3.797), which shows two components in the high-resolution Hubble Space Telescope and James Webb Space Telescope images. The exceptionally rich data set of CDFS-6664 enables us to extract the two components across multiple wavelengths and estimate their physical properties. We show that CDFS-6664 is consistent with a major merger system with boosted star formation in both components and that the offset LyC emission is most likely associated with the bluer and younger component in this merging system. Our result offers an example in which the offset can be caused by a merger. Future observations of more offset LyC emitters would elucidate the role that mergers play in the escape of LyC photons.

We present results from the JWST First Reionization Epoch Spectroscopically Complete Observations survey on the star-forming sequence (SFS) of galaxies at 1.0 < z < 1.7, around the peak of the cosmic star formation history. Star formation rates (SFRs) are measured from the redshifted, relatively dust-insensitive Paschen- α emission line, and stellar mass measurements include the F444W (4.4 μ m; rest-frame H ) band. We find SFRs of galaxies with log( M _* / M _⊙ ) > 9.5 that are lower than found in many earlier studies by up to 0.6 dex, but in good agreement with recent results obtained with the Prospector fitting framework. The difference (log(SFR(Pa α )-SFR( Prospector )) is −0.09 ± 0.04 dex at 10 ^10−11 M _⊙ . We also measure the empirical relation between Paschen- α luminosity and rest-frame H- band magnitude and find that the scatter is only 0.04 dex lower than that of the SFR– M * relation and is much lower than the systematic differences among relations in the literature due to various methods of converting observed measurements to physical properties. We additionally identify examples of sources—that, with standard cutoffs via the UVJ diagram, would be deemed quiescent—with significant (log(sSFR)> −11 yr ^−1 ), typically extended, Paschen- α emission. Our results may be indicative of the potential unification of methods used to derive the SFS with careful selection of star-forming galaxies and independent SFR and stellar mass indicators.

The detections of Ly α emission in galaxies with redshifts above 5 are of utmost importance for constraining the cosmic reionization timeline, yet such detections are usually based on slit spectroscopy. Here we investigate the significant bias induced by slit placement on the estimate of Ly α escape fraction ( ${f}_{\mathrm{esc}}^{\mathrm{Ly}\alpha }$ ) by presenting a galaxy (dubbed A2744-z6Lya) at z = 5.66, where its deep JWST/NIRSpec prism spectroscopy completely misses the strong Ly α emission detected in the MUSE data. A2744-z6Lya exhibits a pronounced UV continuum with an extremely steep spectral slope of $\beta =-{2.640}_{-0.008}^{+0.008}$ , and it has a stellar mass of ∼10 ^8.75 M _⊙ , a star formation rate of ∼5.95 M _⊙ yr ^−1 , and a gas-phase metallicity of $12+\mathrm{log}({\rm{O}}/{\rm{H}})\sim 7.90$ . The observed flux and rest-frame equivalent width of its Ly α from MUSE spectroscopy are 1.2 × 10 ^−16 erg s ^−1 cm ^−2 and 75 Å, equivalent to ${f}_{\mathrm{esc}}^{\mathrm{Ly}\alpha }=42 \% \pm 1 \% $ . However, its Ly α nondetection from JWST NIRSpec gives a 5 σ upper limit of

The universal rotation curve (URC) of disk galaxies was originally proposed to predict the shape and amplitude of any rotation curve (RC) based solely on photometric data. Here, the URC is investigated with an extensive set of spatially resolved RCs drawn from the PROBES-I, PROBES-II, and MaNGA databases with matching multiband surface brightness profiles from the DESI-LIS and Wide-Field Infrared Survey Explorer surveys for 3846 disk galaxies. Common URC formulations fail to achieve an adequate level of accuracy to qualify as truly universal over fully sampled RCs. We develop neural network (NN) equivalents for the proposed URCs that predict RCs with higher accuracy, showing that URC inaccuracies are not due to insufficient data but rather nonoptimal formulations or sampling effects. This conclusion remains even if the total RC sample is pruned for symmetry. The latest URC prescriptions and their NN equivalents trained on our subsample of 579 disk galaxies with symmetric RCs perform similarly to the URC/NN trained on the complete data sample. We conclude that a URC with an acceptable level of accuracy (Δ V _circ ≲ 15%) at all radii would require a detailed modeling of a galaxy’s central regions and outskirts (e.g., for baryonic effects leading to contraction or expansion of any dark-matter-only halo).

Empirical correlations connecting starlight to galaxy dynamics (e.g., the fundamental plane (FP) of elliptical/quiescent (Q) galaxies and the Tully–Fisher relation of spiral/star-forming (SF) galaxies) provide cosmology-independent distance estimation and are central to local Universe cosmology. In this work, we introduce the mass hyperplane (MH), which is the stellar-to-dynamical mass relation ( M _⋆ / M _dyn ) recast as a linear distance indicator. Building on recent FP studies, we show that both SF and Q galaxies follow the same empirical MH, then use this to measure the peculiar velocities (PVs) for a sample of 2496 galaxies at z < 0.12 from GAMA. The limiting precision of MH-derived distance/PV estimates is set by the intrinsic scatter in size, which we find to be ≈0.1 dex for both Q and SF galaxies (when modeled independently) and ≈0.11 dex when all galaxies are modeled together, showing that the MH is as good as the FP. To empirically validate our framework and distance/PV estimates, we compare the inferred distances to groups as derived using either Q or SF galaxies. A good agreement is obtained with no discernible bias or offset, having a scatter of ≈0.05 dex ≈12% in distance. Further, we compare our PV measurements for the Q galaxies to the previous PV measurements of the galaxies in common between GAMA and the Sloan Digital Sky Survey, which shows similarly good agreement. Finally, we provide comparisons of PV measurements made with the FP and the MH, then discuss possible improvements in the context of upcoming surveys such as the 4MOST Hemisphere Survey.

We present the Galaxy Assembly and Interaction Neural Networks ( Gainn ), a series of artificial neural networks for predicting the redshift, stellar mass, halo mass, and mass-weighted age of simulated galaxies based on James Webb Space Telescope (JWST) photometry. Our goal is to determine the best neural network for predicting these variables at 11 < z < 15. The parameters of the optimal neural network can then be used to estimate these variables for real, observed galaxies. The inputs of the neural networks are JWST filter magnitudes of a subset of five broadband filters (F150W, F200W, F277W, F356W, and F444W) and two medium-band filters (F162M and F182M). We compare the performance of the neural networks using different combinations of these filters, as well as different activation functions and numbers of layers. The best neural network predicted redshift with a normalized rms error of ${0.010}_{-0.001}^{+0.003}$ , stellar mass with rms = ${0.089}_{-0.022}^{+0.044}$ , halo mass with a mean-squared error of ${0.022}_{-0.008}^{+0.014}$ , and mass-weighted age with rms = ${12.466}_{-2.408}^{+5.065}$ . We also test the performance of Gainn on real data from MACS0647JD, an object observed by JWST. Predictions from Gainn for the first projection of the object (JD1) have normalized bias 〈Δ z 〉 < 0.00228, which is significantly smaller than found with template-fitting methods. We find that the optimal filter combination is F277W, F356W, F162M, and F200W when considering both theoretical accuracy and observational resources from JWST.

Recent works have discovered a relatively tight correlation between Ω _m and the properties of individual simulated galaxies. Because of this, it has been shown that constraints on Ω _m can be placed using the properties of individual galaxies while accounting for uncertainties in astrophysical processes such as feedback from supernovae and active galactic nuclei. In this work, we quantify whether using the properties of multiple galaxies simultaneously can tighten those constraints. For this, we train neural networks to perform likelihood-free inference on the value of two cosmological parameters (Ω _m and σ _8 ) and four astrophysical parameters using the properties of several galaxies from thousands of hydrodynamic simulations of the CAMELS project. We find that using properties of more than one galaxy increases the precision of the Ω _m inference. Furthermore, using multiple galaxies enables the inference of other parameters that were poorly constrained with one single galaxy. We show that the same subset of galaxy properties are responsible for the constraints on Ω _m from one and multiple galaxies. Finally, we quantify the robustness of the model and find that without identifying the model range of validity, the model does not perform well when tested on galaxies from other galaxy formation models.

Using the Systematically Measuring Ultra-Diffuse Galaxies and Sloan Digital Sky Survey catalogs and our own reprocessing of the Legacy Survey imaging, we investigate the properties of nuclear star clusters (NSCs) in galaxies having central surface brightnesses as low as 27 mag arcsec ^−2 . We identify 273 (123 with known redshift) and 32 NSC-bearing galaxies in the two samples, respectively, where we require candidate NSCs to have a separation of less than 0.10 r _e from the galaxy center. We find that galaxies with low central surface brightness ( μ _0, _g > 24 mag arcsec ^−2 ) are more likely to contain an NSC if they (1) have a higher stellar mass, (2) have a higher stellar-to-total mass ratio, (3) have a brighter central surface brightness, (4) have a larger axis ratio, or (5) lie in a denser environment. Because of the correlations among these various quantities, it is likely that only one or two are true physical drivers. We also find scaling relations for the NSC mass with stellar mass ( M _NSC / ${M}_{\odot }={10}^{6.02\pm 0.03}{({M}_{* ,\mathrm{gal}}/{10}^{8}{M}_{\odot })}^{0.77\pm 0.04}$ ) and halo mass ( M _NSC / ${M}_{\odot }={10}^{6.11\pm 0.05}{({M}_{h,\mathrm{gal}}/{10}^{10}{M}_{\odot })}^{0.92\pm 0.05}$ ), although it is the scaling with halo mass that is consistent with a direct proportionality. In galaxies with an NSC, M _NSC ≈ 10 ^−4 M _h _,gal . This proportionality echoes the finding of a direct proportionality between the mass (or number) of globular clusters (GCs) in galaxies and the galaxy’s total mass. These findings favor a related origin for GCs and NSCs.

The current and upcoming large data volume galaxy surveys require the use of machine-learning techniques to maximize their scientific return. This study explores the use of Self-Organizing Maps (SOMs) to estimate galaxy parameters with a focus on handling cases of missing data and providing realistic probability distribution functions for the parameters. We train an SOM with a simulated mass-limited lightcone assuming a ugrizY JHK _s +IRAC data set, mimicking the Hyper Suprime-Cam Deep joint data set. For parameter estimation, we derive SOM likelihood surfaces considering photometric errors to derive total (statistical and systematic) uncertainties. We explore the effects of missing data, including which bands are particularly critical to the accuracy of the derived parameters. We demonstrate that the parameter recovery is significantly better when the missing bands are “filled in” rather than if they are completely omitted. We propose a practical method for such recovery of missing data.

To enhance our understanding of the impact of galaxy mergers on the kinematics of early-type galaxies (ETGs), we examine differences in specific stellar angular momentum within the half-light radius ( ${\lambda }_{{R}_{e}}$ ) among ETGs with different types of tidal features and those without such features. This is accomplished by categorizing tidal features, which serve as direct evidence of recent mergers, into shells, streams, and tails, through deep images from the DESI Legacy Survey, and by using MaNGA data for the analysis of the kinematics of 1244 ETGs at z < 0.055. We find that ETGs with tidal features typically have reduced ${\lambda }_{{R}_{e}}$ values that are lower by 0.12 dex than ETGs without tidal features. ETGs with shells contribute most to the reduction in ${\lambda }_{{R}_{e}}$ . Consequently, nearly half of ETGs with shells are classified as slow rotators, a fraction that is more than twice as high as that of ETGs with tails or streams, and over three times higher than that of ETGs without tidal features. These trends generally remain valid even when ETGs are divided into several mass bins. Our findings support the idea that radial mergers, which are more effective at reducing ${\lambda }_{{R}_{e}}$ than circular mergers, are more closely associated with the formation of shells rather than streams or tails. The detection of shells in slightly more massive ETGs compared to streams and tails may be attributed to the fact that massive satellite galaxies are more likely to be accreted through radial orbits, due to the nature of dynamical friction.

We examine the central stellar velocity dispersion of subhalos based on IllustrisTNG cosmological hydrodynamic simulations. The central velocity dispersion is a fundamental observable that links galaxies with their dark matter subhalos. We carefully explore simulated stellar velocity dispersions derived with different definitions to assess possible systematics. We explore the impact of variation in the identification of member stellar particles, the viewing axes, the velocity dispersion computation technique, and simulation resolution. None of these issues impact the velocity dispersion significantly; any systematic uncertainties are smaller than the random error. We examine the stellar mass–velocity dispersion relation as an observational test of the simulations. At fixed stellar mass, the observed velocity dispersions significantly exceed the simulation results. This discrepancy is an interesting benchmark for the IllustrisTNG simulations because the simulations are not explicitly tuned to match this relation. We demonstrate that the stellar velocity dispersion provides measures of the dark matter velocity dispersion and the dark matter subhalo mass.

We present a study of a sample of 45 spectroscopically confirmed, UV luminous galaxies at z ∼ 6. They were selected as bright Lyman-break galaxies (LBGs) using deep multiband optical images in more than 2 deg ^2 of the sky, and subsequently identified via their strong Ly α emission. The majority of these LBGs span an absolute UV magnitude range from −22.0 to −20.5 mag with Ly α equivalent width (EW) between ∼10 and ∼200 Å, representing the most luminous galaxies at z ∼ 6 in terms of both UV continuum emission and Ly α line emission. We model the spectral energy distributions of 10 LBGs that have deep infrared observations from Hubble Space Telescope, JWST, and/or Spitzer, and find that they have a wide range of stellar masses and ages. They also have high star formation rates ranging from a few tens to a few hundreds of solar mass per year. Five of the LBGs have JWST or HST images, and four of them show compact morphology in these images, including one that is roughly consistent with a point source, suggesting that UV luminous galaxies at this redshift are generally compact. The fraction of our photometrically selected LBGs with strong Ly α emission (EW > 25 Å) is about 0.2, which is consistent with previous results and supports a moderate evolution of the intergalactic medium opacity at the end of cosmic reionization. Using deep X-ray images, we do not find evidence of strong active galactic nucleus (AGN) activity in these galaxies, but our constraint is loose, and we are not able to rule out the possibility of any weak AGN activity.

Robust estimation of star formation rates (SFRs) at higher redshifts ( z ≳ 1) using UV–optical–near-infrared (NIR) photometry is contingent on the ability of spectral energy distribution (SED) fitting to constrain the dust attenuation, stellar metallicity, and star formation history (SFH) simultaneously. IR-derived dust luminosities can help break the degeneracy between these parameters, but IR data are often not available. Here, we explore strategies for SED fitting at z ≳ 1 in the absence of IR data using a sample of log M _* > 10.2 star-forming galaxies from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) for which 24 μ m data are available. We adopt the total IR luminosity ( L _TIR ) obtained from 24 μ m as the “ground truth,” which allows us to assess how well it can be recovered (as L _dust ) from UV–optical–NIR SED fitting. We test a variety of dust attenuation models, stellar population synthesis models, metallicity assumptions, and SFHs separately to identify which assumptions maximize the agreement (correlation and linearity) between L _TIR and L _dust . We find that a flexible dust attenuation law performs best. For stellar populations, we find that Bruzual & Charlot models are favored over those of Eldridge et al. Fixing the stellar metallicity at solar value is preferred to other fixed values or leaving it as a free parameter. For SFHs, we find that minimizing the variability in the recent (

We present our photometric search for potential nuclear star clusters (NSCs) in ultra-diffuse galaxies (UDGs) as an extension of the SMUDGes catalog. We identify 325 SMUDGes galaxies with NSCs and, from the 144 with existing distance estimates, identify 33 NSC hosts as UDGs ( μ _0, _g ≥ 24 mag arcsec ^−2 , r _e ≥ 1.5 kpc). The SMUDGes with NSCs lie on the galaxy red sequence, satisfy the relationship between NSC and host galaxy stellar masses, have a mean NSC stellar mass fraction of 0.02 but reach as high as 0.1, have NSCs that are displaced from the host center with a standard deviation of 0.10 r _e , and weakly favor higher-density environments. All of these properties are consistent with previous results from higher surface brightness galaxy samples, allowing for at most a relatively weak dependence of NSC behavior on host galaxy surface brightness.

The PRObabilistic Value-added Bright Galaxy Survey (PROVABGS) catalog will provide measurements of galaxy properties, such as stellar mass (M *), star formation rate (SFR), stellar metallicity (Z), and stellar age (t age), for >10 million galaxies of the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey. Full posterior distributions of the galaxy properties will be inferred using state-of-the-art Bayesian spectral energy distribution (SED) modeling of DESI spectroscopy and Legacy Surveys photometry. In this work, we present the SED model, the neural emulator for the model, and the Bayesian inference framework of PROVABGS. Furthermore, we apply the PROVABGS SED modeling on realistic synthetic DESI spectra and photometry, constructed using the L-Galaxies semi-analytic model. We compare the inferred galaxy properties to the true values of the simulation using a hierarchical Bayesian framework to quantify accuracy and precision. Overall, we accurately infer the true M *, SFR, Z, and t age of the simulated galaxies. However, the priors on galaxy properties induced by the SED model have a significant impact on the posteriors, which we characterize in detail. This work also demonstrates that a joint analysis of spectra and photometry significantly improves the constraints on galaxy properties over photometry alone and is necessary to mitigate the impact of the priors. With the methodology presented and validated in this work, PROVABGS will maximize information extracted from DESI observations and extend current galaxy studies to new regimes and unlock cutting-edge probabilistic analyses. https://github.com/changhoonhahn/provabgs/ [ABSTRACT FROM AUTHOR]

We present results from the Cosmic Evolution Early Release Survey on the stellar population parameters for 28 galaxies with redshifts 4 < z < 9 using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 μ m data extend the coverage of the rest-frame spectral energy distribution to nearly 1 μ m for galaxies in this redshift range. By modeling the galaxies’ SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å)— I -band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower ratios of stellar mass to light. The MIRI data reduce the stellar masses by $\langle {\rm{\Delta }}\mathrm{log}{M}_{* }\rangle =0.25$ dex at 4 < z < 6 and 0.37 dex at 6 < z < 9. This also reduces the star formation rates (SFRs) by 〈ΔlogSFR〉 = 0.14 dex at 4 < z < 6 and 0.27 dex at 6 < z < 9. The MIRI data also improve constraints on the allowable stellar mass formed in early star formation. We model this using a star formation history that includes both a “burst” at z _f = 100 and a slowly varying (“delayed- τ ”) model. The MIRI data reduce the allowable stellar mass by 0.6 dex at 4 < z < 6 and by ≈1 dex at 6 < z < 9. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early Universe ( z ≈ 9). Therefore, observations of rest-frame ≳1 μ m are paramount for constraining the stellar-mass buildup in galaxies at very high redshifts.

We present a study of seven spectroscopically confirmed (Ly α -emitting) galaxies at redshift z ≃ 6 using the JWST/NIRCam imaging data. These galaxies, with a wide range of Ly α luminosities, were recently observed in a series of NIRCam broad and medium bands. We constrain the rest-frame UV/optical continua and measure the H α line emission of the galaxies using the combination of JWST/NIRCam and archival HST/WFC3 infrared photometry. We further estimate the escape fractions of their Ly α photons ( ${f}_{\mathrm{esc}}^{\mathrm{Ly}\alpha }$ ) and the production efficiency of ionizing photons ( ξ _ion ). Among the sample, six out of seven galaxies have Ly α escape fractions of ≲10%, which might be the status for most of the star-forming galaxies at z ≃ 6. One UV-faint Ly α galaxy with an extremely blue UV slope has a large value of ${f}_{\mathrm{esc}}^{\mathrm{Ly}\alpha }$ reaching ≃50%. These galaxies have a broad range of ξ _ion over log _10 ξ _ion,0 (Hz erg ^−1 ) ∼ 25.0–26.5. We find that UV-fainter galaxies with bluer UV-continuum slopes likely have higher escape fractions of Ly α photons. We also find that galaxies with higher Ly α line emission tend to produce ionizing photons more efficiently. The most Ly α -luminous galaxy in the sample has a very high ξ _ion,0 of log _10 ξ _ion,0 (Hz erg ^−1 ) > 26. Our results support the scenario that Ly α galaxies may have served as an important contributor to cosmic reionization. Blue and bright Ly α galaxies are excellent targets for JWST follow-up spectroscopic observations.

We present a spectroscopic survey of Ly α emitters (LAEs) at z ∼ 3.7 and z ∼ 4.8. The LAEs are selected using the narrowband technique based on the combination of deep narrowband and broadband imaging data in two deep fields, and then spectroscopically confirmed with the MMT multifiber spectrograph Hectospec. The sample consists of 71 LAEs at z ∼ 3.7 and 69 LAEs at z ∼ 4.8 over ∼1.5 deg ^2 , making it one of the largest spectroscopically confirmed samples of LAEs at the two redshifts. Their Ly α luminosities are measured using the secure redshifts and deep photometric data, and span a range of ∼10 ^42.5 –10 ^43.6 erg s ^−1 , so these LAEs represent the most luminous galaxies at the redshifts in terms of Ly α luminosity. We estimate and correct sample incompletenesses and derive reliable Ly α luminosity function (LFs) at z ∼ 3.7 and 4.8 based on the two spectroscopic samples. We find that our Ly α LFs are roughly consistent (within a factor of 2−3) with previous measurements at similar redshifts that were derived from either photometric samples or spectroscopic samples. By comparing with previous studies in different redshifts, we find that the Ly α LFs decrease mildly from z ∼ 3.1 to z ∼ 5.7, supporting the previous claim of the slow LF evolution between z ∼ 2 and z ∼ 6. At z > 5.7, the LF declines rapidly toward higher redshift, partly due to the effect of cosmic reionization.

By utilizing the spatially resolved photometry of galaxies at 0.2 < z < 3.0 in the CEERS field, we estimate the resolved and unresolved stellar mass via fitting of the spectral energy distribution (SED) to study the discrepancy between them. We first compare M _* derived from photometry with and without the JWST wavelength coverage and find that M _* can be overestimated by up to 0.2 dex when lacking rest-frame near-infrared (NIR) data. The SED fitting process tends to overestimate both stellar age and dust attenuation in the absence of rest-frame NIR data, consequently leading to a larger observed mass-to-light ratio and hence an elevated M _* . With the inclusion of the JWST NIR photometry, we find no significant disparity between the estimates of resolved and unresolved stellar mass, providing a plausible solution to the conflict between them out to z ∼ 3. Further investigation demonstrates that reliable M _* estimates can be obtained, regardless of whether they are derived from spatially resolved or spatially unresolved photometry, so long as the reddest filter included in the SED fitting has a rest-frame wavelength longer than 10000 Å.

We present a study of Lyman continuum (LyC) emission in a sample of ∼150 Ly α emitters (LAEs) at z ≈ 3.1 in the Subaru-XMM Deep Survey field. These LAEs were previously selected using the narrowband technique and spectroscopically confirmed with Ly α equivalent widths (EWs) ≥ 45 Å. We obtain deep UV images using a custom intermediate-band filter U _J that covers a wavelength range of 3330 ∼ 3650 Å, corresponding to 810 ∼ 890 Å in the rest frame. We detect five individual LyC galaxy candidates in the U _J band, and their escape fractions ( f _esc ) of LyC photons are roughly between 40% and 80%. This supports a previous finding that a small fraction of galaxies may have very high f _esc . We find that the f _esc values of the five LyC galaxies are not apparently correlated with other galaxy properties such as Ly α luminosity and EW, UV luminosity and slope, and star formation rate (SFR). This is partly due to the fact that these galaxies only represent a small fraction (∼3%) of our LAE sample. For the remaining LAEs that are not detected in U _J , we stack their U _J -band images and constrain their average f _esc . The upper limit of the average f _esc value is about 16%, consistent with the results in the literature. Compared with the non-LyC LAEs, the LyC LAEs tend to have higher Ly α luminosities, Ly α EWs, and SFRs, but their UV continuum slopes are similar to those of other galaxies.

We investigate the alignment of galaxy and halo orientations using the TNG300-1 hydrodynamical simulation. Our analysis reveals that the distribution of the 2D misalignment angle θ _2D can be well described by a truncated shifted exponential distribution with only one free parameter across different redshifts and galaxy/halo properties. We demonstrate that the galaxy–ellipticity (GI) correlations of galaxies can be reproduced by perturbing halo orientations with the obtained θ _2D distribution, with only a small bias (

To better understand the formation of large, low-surface-brightness galaxies, we measure the correlation function between ultradiffuse galaxy (UDG) candidates and Milky Way analogs (MWAs). We find that: (1) the projected radial distribution of UDG satellites (projected surface density ∝ r ^−0.84±0.06 ) is consistent with that of normal satellite galaxies; (2) the number of UDG satellites per MWA ( S _UDG ) is ∼0.5 ± 0.1 over projected radii from 20 to 250 kpc and −17 < M _r < −13.5; (3) S _UDG is consistent with a linear extrapolation of the relationship between the number of UDGs per halo versus halo mass obtained over galaxy group and cluster scales; (4) red UDG satellites dominate the population of UDG satellites (∼80%); (5) over the range of satellite magnitudes studied, UDG satellites comprise ∼10% of the satellite galaxy population of MWAs; and (6) a significant fraction of these (∼13%) have estimated total masses >10 ^10.9 M _⊙ or, equivalently, at least half the halo mass of the LMC, and populate a large fraction (∼18%) of the expected subhalos down to these masses. All of these results suggest a close association between the overall low-mass galaxy population and UDGs, which we interpret as favoring models where UDG formation principally occurs within the general context of low-mass galaxy formation over models invoking more exotic physical processes specifically invoked to form UDGs.

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, the basic demographics of clumps and their host galaxies remain uncertain. Using the Hubble Space Telescope/Wide Field Camera 3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, we detect and analyze giant star-forming clumps in galaxies at 0.5 ≤ z ≤ 1, connecting two epochs when clumps are common (at cosmic high noon, z ∼ 2) and rare (in the local Universe). We construct a clump sample whose rest-frame 1600 Å luminosity is 3 times higher than the most luminous local H ii regions ( M _UV ≤ −16 AB). In our sample, 35% ± 3% of low-mass galaxies (log[ M _∗ / M _⊙ ] < 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22% ± 3% and 22% ± 4% for intermediate (10 ≤ log[ M _∗ / M _⊙ ] ≤ 10.5) and high-mass (log[ M _∗ / M _⊙ ] > 10.5) galaxies, in agreement with previous studies. When compared to similar-mass nonclumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher star formation rates (SFRs) and bluer rest-frame U − V colors, while high-mass clumpy SFGs tend to be larger than nonclumpy SFGs. However, clumpy and nonclumpy SFGs have similar Sérsic index, indicating a similar underlying density profile. Furthermore, we investigate how the UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher specific SFR galaxies and are found closer to their hosts’ galactic centers.

Large diffuse galaxies are hard to find, but understanding the environments where they live, their numbers, and ultimately their origins, is of intense interest and importance for galaxy formation and evolution. Using Subaru’s Hyper Suprime-Cam Strategic Survey Program, we perform a systematic search for low surface brightness galaxies and present novel and effective methods for detecting and modeling them. As a case study, we surveyed 922 Milky Way analogs in the nearby Universe (0.01 < z < 0.04) and built a large sample of satellite galaxies that are outliers in the mass–size relation. These “ultra-puffy” galaxies (UPGs), defined to be 1.5 σ above the average mass–size relation, represent the tail of the satellite size distribution. We find that each MW analog hosts N _UPG = 0.31 ± 0.05 UPGs on average, which is consistent with but slightly lower than the observed abundance at this halo mass in the Local Volume. We also construct a sample of ultra-diffuse galaxies (UDGs) in MW analogs and find an abundance of N _UDG = 0.44 ± 0.05 per host. With literature results, we confirm that the UDG abundance scales with the host halo mass following a sublinear power law. We argue that our definition of UPGs, which is based on the mass–size relation, is more physically motivated than the common definition of UDGs, which depends on the surface brightness and size cuts and thus yields different surface mass density cuts for quenched and star-forming galaxies.

Recent observations have reignited interest in a population of dwarf galaxies that are large and diffuse for their mass, often called ultra-diffuse galaxies (UDGs). However, the origin and evolution of these mass–size outliers and the role of the environment are still unclear. Using the exquisitely deep and wide Hyper Suprime-Cam Strategic Survey images, we search for ultra-puffy galaxies (UPGs), defined as being 1.5 σ larger than the average size for their mass, around Milky Way–like galaxies. We present the sizes and radial distributions of mass–size outliers and derive their quenched fraction to explore the impact of the environment. Surprisingly, despite being outliers in size, the UPGs have a similar quenched fraction as normal-sized satellites of Milky Way analogs in both observations and simulations, suggesting that quenching is not tied to being a mass–size outlier. The quenched fraction is higher for the UPGs associated with redder hosts, as well as those that are closer to the host in projection. In contrast, the UDGs are overall redder and more quiescent compared with normal satellites. We show that the classic definition of UDGs is heavily weighted toward quenched galaxies and thus cannot be used for a study of quenching of mass–size outliers.

This work presents and releases a catalog of new photometrically derived physical properties for the ∼10 ^5 most well-measured galaxies in the COSMOS field on the sky. Using a recently developed technique, spectral energy distributions are modeled assuming a stellar initial mass function (IMF) that depends on the temperature of gas in star-forming regions. The method is applied to the largest current sample of high-quality panchromatic photometry, the COSMOS2020 catalog, that allows for testing this assumption. It is found that the galaxies exhibit a continuum of IMF and gas temperatures, most of which are bottom-lighter than measured in the Milky Way. As a consequence, the stellar masses and star formation rates of most galaxies here are found to be lower than those measured by traditional techniques in the COSMOS2020 catalog by factors of ∼1.6–3.5 and 2.5–70.0, respectively, with the change being the strongest for the most active galaxies. The resulting physical properties provide new insights into variation of the IMF-derived gas temperature along the star-forming main sequence and at quiescence, produce a sharp and coherent picture of downsizing, as seen from the stellar mass functions, and hint at a possible high-temperature and high-density stage of early galactic evolution.

We measure the galaxy-ellipticity (GI) correlations for the Sloan Digital Sky Survey Data Release 12 LOWZ and CMASS samples with the shape measurements from the DESI Legacy Imaging Surveys. We model the GI correlations in an N -body simulation with our recent accurate stellar–halo mass relation from the Photometric object Around Cosmic webs (PAC) method. The large data set and our accurate modeling turns out an accurate measurement of the alignment angle between central galaxies and their host halos. We find that the alignment of central elliptical galaxies with their host halos increases monotonically with galaxy stellar mass or host halo mass, which can be well described by a power law for the massive galaxies. We also find that central elliptical galaxies are more aligned with their host halos in LOWZ than in CMASS, which might indicate an evolution of galaxy–halo alignment, though future studies are needed to verify this is not induced by the sample selections. In contrast, central disk galaxies are aligned with their host halos about 10 times more weakly in the GI correlation. These results have important implications for intrinsic alignment (IA) correction in weak lensing studies, IA cosmology, and theory of massive galaxy formation.

In this empirical work, we aim to quantify the systematic uncertainties in stellar-mass ( M _⋆ ) estimates made from spectral energy distribution (SED) fitting through stellar population synthesis (SPS) for galaxies in the local Universe by using the dynamical mass ( M _dyn ) estimator as an SED-independent check on stellar mass. We first construct a statistical model of the high-dimensional space of galaxy properties; including size ( R _e ), velocity dispersion ( σ _e ), surface brightness ( I _e ), mass-to-light ratio ( M _⋆ / L ), rest-frame color, Sérsic index ( n ), and dynamical mass ( M _dyn ), and accounting for selection effects and covariant errors. We disentangle the correlations among galaxy properties and find that the variation in M _⋆ / M _dyn is driven by σ _e , Sérsic index and color. We use these parameters to calibrate an SED-independent M _⋆ estimator, ${\hat{M}}_{\star }$ . We find the random scatter of the relation ${M}_{\star }-{\hat{M}}_{\star }$ to be 0.108 dex and 0.147 dex for quiescent and star-forming galaxies, respectively. Finally, we inspect the residuals as a function of SPS parameters (dust, age, metallicity, and star formation rate) and spectral indices (H α , H δ , and D _n 4000). For quiescent galaxies, ∼65% of the scatter can be explained by the uncertainty in SPS parameters, with dust and age being the largest sources of uncertainty. For star-forming galaxies, while age and metallicity are the leading factors, SPS parameters account for only ∼13% of the scatter. These results leave us with remaining unmodelled scatters of 0.055 dex and 0.122 dex for quiescent and star-forming galaxies, respectively. This can be interpreted as a conservative limit on the precision in M _⋆ that can be achieved via simple SPS modeling.

We present the completed catalog of ultradiffuse galaxy (UDG) candidates (7070 objects) from our search of the DR9 Legacy Survey images, including distance and total mass estimates for 1529 and 1436 galaxies, respectively, that we provide and describe in detail. From the sample with estimated distances, we obtain a sample of 585 UDGs ( μ _0, _g ≥ 24 mag arcsec ^−2 and r _e ≥ 1.5 kpc) over 20,000 square degrees of sky in various environments. We conclude that UDGs in our sample are limited to 10 ^10 ≲ M _h / M _⊙ ≲ 10 ^11.5 and are on average a factor of 1.5–7 deficient in stars relative to the general population of galaxies of the same total mass. That factor increases with increasing galaxy size and mass up to a factor of ∼10 when the total mass of the UDG increases beyond M _h = 10 ^11 M _⊙ . We do not find evidence that this factor has a dependence on the UDGs large-scale environment.

We present a new selection of 358 blue compact dwarf galaxies (BCDs) from 5000 square degrees in the Dark Energy Survey, and the spectroscopic follow-up of a subsample of 68 objects. For the subsample of 34 objects with deep spectra, we measure the metallicity via the direct T _e method using the auroral [O iii ] λ 4363 emission line. These BCDs have an average oxygen abundance of 12+log(O/H) = 7.8, with stellar masses between 10 ^7 and 10 ^8 M _⊙ and specific star-formation rates between ∼10 ^−9 and 10 ^−7 yr ^−1 . We compare the position of our BCDs with the mass–metallicity ( M – Z ) and luminosity–metallicity ( L – Z ) relation derived from the Local Volume Legacy sample. We find the scatter about the M – Z relation is smaller than the scatter about the L – Z relation. We identify a correlation between the offsets from the M – Z and L – Z relation that we suggest is due to the contribution of metal-poor inflows. Finally, we explore the validity of the mass–metallicity–SFR fundamental plane in the mass range probed by our galaxies. We find that BCDs with stellar masses smaller than 10 ^8 M _⊙ do not follow the extrapolation of the fundamental plane. This result suggests that mechanisms other than the balance between inflows and outflows may be at play in regulating the position of low-mass galaxies in the M – Z –SFR space.

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at 1.3 ≤ z ≤ 2.6 from the MOSFIRE Deep Evolution Field survey. We divide galaxies with a detected H α emission line and coverage of H β into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), and we then stack their spectra. From each stack, we measure the Balmer decrement and gas-phase metallicity, and then we compute the median A _V and UV continuum spectral slope ( β ). First, we find that none of the dust properties (Balmer decrement, A _V , or β ) varies with the axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model in which attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases the SFR and lowers the metallicity, but leaves the dust column density mostly unchanged. We quantify this idea using the Kennicutt–Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation.

We present an updated version of Lightning , a galaxy spectral energy distribution (SED) fitting code that can model X-ray to submillimeter observations. The models in Lightning include the options to contain contributions from stellar populations, dust attenuation and emission, and active galactic nuclei (AGNs). X-ray emission, when utilized, can be modeled as originating from stellar compact binary populations with the option to include emission from AGNs. We have also included a variety of algorithms to fit the models to observations and sample parameter posteriors; these include adaptive Markov Chain Monte Carlo (MCMC), affine-invariant MCMC, and Levenberg–Marquardt gradient-descent ( MPFIT ) algorithms. To demonstrate some of the capabilities of Lightning , we present several examples using a variety of observational data. These examples include (1) deriving the spatially resolved stellar properties of the nearby galaxy M81, (2) demonstrating how X-ray emission can provide constraints on the properties of the supermassive black hole of a distant AGN, (3) exploring how to rectify the attenuation effects of inclination on the derived the star formation rate of the edge-on galaxy NGC 4631, (4) comparing the performance of Lightning to similar Bayesian SED-fitting codes when deriving physical properties of the star-forming galaxy NGC 628, and (5) comparing the derived X-ray and UV-to-IR AGN properties from Lightning and CIGALE for a distant AGN. Lightning is an open-source application developed in IDL and is available at https://github.com/rafaeleufrasio/lightning .

We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ( z < 0.4) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as q _PAH , which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous mid-infrared photometric data points covering 7–24 μ m and far-infrared flux densities, q _PAH values are derived through spectral energy distribution fitting. The correlation between q _PAH and other physical properties of galaxies, i.e., gas-phase metallicity ( $12+\mathrm{log}({\rm{O}}/{\rm{H}})$ ), stellar mass, and specific star-formation rate (sSFR) are explored. As in previous studies, q _PAH values of galaxies with high metallicity are found to be higher than those with low metallicity. The strength of PAH emission is also positively correlated with the stellar mass and negatively correlated with the sSFR. The correlation between q _PAH and each parameter still exists even after the other two parameters are fixed. In addition to the PAH strength, the application of metallicity-dependent gas-to-dust mass ratio appears to work well to estimate gas mass that matches the observed relationship between molecular gas and physical parameters. The result obtained will be used to calibrate the observed PAH luminosity-total infrared luminosity relation, based on the variation of MIR-FIR SED, which is used in the estimation of hidden star formation.

Taking advantage of the Photometric objects Around Cosmic webs method developed in Paper I, we measure the excess surface density ${\bar{n}}_{2}{w}_{{\rm{p}}}$ of the photometric objects around spectroscopic objects down to stellar masses 10 ^8.0 M _⊙ , 10 ^9.2 M _⊙ , and 10 ^9.8 M _⊙ in the redshift ranges of z _s < 0.2, 0.2 < z _s < 0.4, and 0.5 < z _s < 0.7, respectively, using data from the DESI Legacy Imaging Surveys and spectroscopic samples from the Sloan Digital Sky Survey (i.e., the Main, LOWZ, and CMASS samples). We model the measured ${\bar{n}}_{2}{w}_{{\rm{p}}}$ in an N -body simulation, using the abundance matching method, and we constrain the stellar–halo mass relations (SHMRs) in the three redshift ranges to percent levels. With the accurate modeling, we demonstrate that the stellar mass scatter for a given halo mass is nearly a constant, and that the empirical form of Behroozi et al. describes the SHMR better than the double-power-law form at low mass. Our SHMR accurately captures the downsizing of massive galaxies from z _s = 0.7, while it also indicates that small galaxies are still growing faster than their host halos. The galaxy stellar mass functions (GSMFs) from our modeling are in perfect agreement with the model-independent measurements in Paper III, although the current work extends the GSMFs to a much smaller stellar mass. Based on the GSMFs and the SHMRs, we derive the stellar mass completeness and halo occupation distributions for the LOWZ and CMASS samples, which are useful for correctly interpreting their cosmological measurements, such as galaxy–galaxy lensing and redshift space distortion.