Search

Active galactic nuclei (AGN) are the signposts of black hole growth, and likely play an important role in galaxy evolution. An outstanding question is whether AGN of different spectral types indicate different evolutionary stages in the coevolution of black holes and galaxies. We present the angular correlation function between an AGN sample selected from Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) optical photometry and Wide-field Infrared Survey Explorer mid-IR photometry and a luminous red galaxy (LRG) sample from HSC-SSP. We investigate AGN clustering strength as a function of luminosity and spectral features across three independent HSC fields totaling ∼600 deg2, for z ∈ 0.6 −1.2 and AGN with L 6 μ m > 3 × 1044 erg s−1. There are ∼28,500 AGN and ∼1.5 million LRGs in our primary analysis. We determine the average halo mass for the full AGN sample (M h ≈ 1012.9 h −1 M ⊙), and note that it does not evolve significantly as a function of redshift (over this narrow range) or luminosity. We find that, on average, unobscured AGN (M h ≈ 1013.3 h −1 M ⊙) occupy ∼4.5× more massive halos than obscured AGN (M h ≈ 1012.6 h −1 M ⊙), at 5 σ statistical significance using 1D uncertainties, and at 3 σ using the full covariance matrix, suggesting a physical difference between unobscured and obscured AGN, beyond the line-of-sight viewing angle. Furthermore, we find evidence for a halo mass dependence on reddening level within the Type I AGN population, which could support the existence of a dust-obscured phase. However, we also find that quite small systematic shifts in the redshift distributions of the AGN sample could explain current and previously observed differences in M h . [ABSTRACT FROM AUTHOR]

We study the stellar properties of a sample of simulated ultradiffuse galaxies (UDGs) with stellar mass M ⋆ = 107.5–109 M ⊙, selected from the TNG50 simulation, where UDGs form mainly in high-spin dwarf-mass halos. We divide our sample into star-forming and quenched UDGs, finding good agreement with the stellar assembly history measured in observations. Star-forming UDGs and quenched UDGs with M ⋆ ≥ 108 M ⊙ in our sample are particularly inefficient at forming stars, having 2–10 times less stellar mass than non-UDGs for the same virial mass halo. These results are consistent with recent mass inferences in UDG samples and suggest that the most inefficient UDGs arise from a late assembly of the dark matter mass followed by a stellar growth that is comparatively slower (for star-forming UDGs) or that was interrupted due to environmental removal of the gas (for quenched UDGs). Regardless of efficiency, UDGs are 60% poorer in [Fe/H] than the population of non-UDGs at a fixed stellar mass, with the most extreme objects having metal content consistent with the simulated mass–metallicity relation at z ∼ 2. Quenched UDGs stop their star formation in shorter timescales than non-UDGs of similar mass and are, as a consequence, alpha enhanced with respect to non-UDGs. We identify metallicity profiles in UDGs as a potential avenue to distinguish between different formation paths for these galaxies, where gentle formation as a result of high-spin halos would present well-defined declining metallicity radial profiles while powerful-outflows or tidal stripping formation models would lead to flatter or constant metallicity as a function of radius due to the inherent mixing of stellar orbits. [ABSTRACT FROM AUTHOR]

We investigate size variation with rest-frame wavelength for star-forming galaxies based on the second James Webb Space Telescope Advanced Deep Extragalactic Survey data release. Star-forming galaxies are typically smaller at longer wavelengths from the UV to the near-IR at z < 3.5, especially for more massive galaxies, indicating the inside-out assembly with in situ star formation if ignoring dust attenuation. The size variation with wavelength shows a strong dependence on stellar mass and shows little or no dependence on redshift, specific star formation rate, and galaxy environment. This suggests that the size growth of star-forming galaxies is a self-regulated process primarily governed by stellar mass. We model size as a function of both mass and redshift simultaneously, obtaining R e ∝ M * 0.23 (1 + z) − 1.04 at a wavelength of 0.45 μ m, and R e ∝ M * 0.20 (1 + z) − 1.08 at 1.0 μ m. Based on this size evolution and the star formation main sequence from the literature, we obtain the locus of typical size growth for individual galaxies of different masses on the mass–size plane. The moving trend of galaxies on the mass–size plane, which indicates the slopes of their locus, strongly correlates with the size ratio between 0.45 μ m and 1.0 μ m, supporting the idea that the size variation with wavelength provides important information on the size growth of galaxies on short timescales. [ABSTRACT FROM AUTHOR]

We use a sample of 73 simulated satellite and central dwarf galaxies spanning a stellar mass range of 105.3–109.1 M ⊙ to investigate the origin of their stellar age gradients. We find that dwarf galaxies often form their stars "inside-out," i.e., the stars form at successively larger radii over time. However, the oldest stars get reshuffled beyond the star-forming radius by fluctuations in the gravitational potential well caused by stellar feedback (the same mechanisms that cause dwarfs to form dark matter cores). The result is that many dwarfs appear to have an "outside-in" age gradient at z = 0, with younger stellar populations more centrally concentrated. However, for the reshuffled galaxies with the most extended star formation, young stars can form out to the large radii to which the old stars have been reshuffled, erasing the age gradient. We find that major mergers do not play a significant role in setting the age gradients of dwarfs. We find similar age gradient trends in satellites and field dwarfs, suggesting that environment plays only a minor role, if any. Finally, we find that the age gradient trends are imprinted on the galaxies at later times, suggesting that the stellar reshuffling dominates after the galaxies have formed 50% of their stellar mass. The later reshuffling is at odds with results from the fire-2 simulations. Hence, age gradients offer a test of current star formation and feedback models that can be probed via observations of resolved stellar populations. [ABSTRACT FROM AUTHOR]

We describe new ultradeep James Webb Space Telescope (JWST) NIRSpec PRISM and grating spectra for the galaxies JADES-GS-z11-0 ( z spec = 11.122 − 0.003 + 0.005 ) and JADES-GS-z13-0 ( z spec = 13.20 − 0.04 + 0.03 ), the most distant spectroscopically confirmed galaxy discovered in the first year of JWST observations. The extraordinary depth of these observations (75 hr and 56 hr, respectively) provides a unique opportunity to explore the redshifts, stellar properties, UV magnitudes, and slopes for these two sources. For JADES-GS-z11-0, we find evidence for multiple emission lines, including [O ii ] λ λ 3726, 3729 and [Ne iii ] λ 3869, resulting in a spectroscopic redshift we determine with 94% confidence. We present stringent upper limits on the emission-line fluxes and line equivalent widths for JADES-GS-z13-0. At this spectroscopic redshift, the Ly α break in JADES-GS-z11-0 can be fit with a damped Ly α absorber with log (N HI / cm − 2) = 22.42 − 0.120 + 0.093 . These results demonstrate how neutral hydrogen fraction and Lyman-damping wings may impact the recovery of spectroscopic redshifts for sources like these, providing insight into the overprediction of the photometric redshifts seen for distant galaxies observed with JWST. In addition, we analyze updated NIRCam photometry to calculate the morphological properties of these resolved sources, and find a secondary source 0.″3 south of JADES-GS-z11-0 at a similar photometric redshift, hinting at how galaxies grow through interactions in the early Universe. [ABSTRACT FROM AUTHOR]

The scaling of galaxy properties with halo mass suggests that feedback loops regulate star formation, but there is no consensus yet about how those feedback loops work. To help clarify discussions of galaxy-scale feedback, Paper I presented a very simple model for supernova feedback that it called the minimalist regulator model. This follow-up paper interprets that model and discusses its implications. The model itself is an accounting system that tracks all of the mass and energy associated with a halo's circumgalactic baryons—the central galaxy's atmosphere. Algebraic solutions for the equilibrium states of that model reveal that star formation in low-mass halos self-regulates primarily by expanding the atmospheres of those halos, ultimately resulting in stellar masses that are insensitive to the mass-loading properties of galactic winds. What matters most is the proportion of supernova energy that couples with circumgalactic gas. However, supernova feedback alone fails to expand galactic atmospheres in higher-mass halos. According to the minimalist regulator model, an atmospheric contraction crisis ensues, which may be what triggers strong black hole feedback. The model also predicts that circumgalactic medium properties emerging from cosmological simulations should depend largely on the specific energy of the outflows they produce, and we interpret the qualitative properties of several numerical simulations in light of that prediction. [ABSTRACT FROM AUTHOR]

This paper presents a new framework for understanding the relationship between a galaxy and its circumgalactic medium (CGM). It focuses on how imbalances between heating and cooling cause either expansion or contraction of the CGM. It does this by tracking all of the mass and energy associated with a halo's baryons, including their gravitational potential energy, even if feedback has pushed some of those baryons beyond the halo's virial radius. We show how a star-forming galaxy's equilibrium state can be algebraically derived within the context of this framework, and we analyze how the equilibrium star formation rate depends on supernova feedback. We consider the consequences of varying the mass loading parameter η M ≡ M ̇ wind / M ̇ * relating a galaxy's gas mass outflow rate ( M ̇ wind ) to its star formation rate ( M ̇ * ) and obtain results that challenge common assumptions. In particular, we find that equilibrium star formation rates in low-mass galaxies are generally insensitive to mass loading, and when mass loading does matter, increasing it actually results in more star formation because more supernova energy is needed to resist atmospheric contraction. [ABSTRACT FROM AUTHOR]

Massive quiescent galaxies in the young Universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 14 massive quiescent galaxies ( log (M ⋆ / M ⊙) > 10 ) at z ∼ 2 selected from a representative sample of the Blue Jay survey. We reconstruct their star formation histories (SFHs) by fitting spectral energy distribution models to the JWST/NIRSpec R ∼ 1000 spectra. We find that massive quiescent galaxies can be split into three categories with roughly equal numbers of galaxies according to their SFHs: (1) relatively old galaxies quenched at early epochs; (2) galaxies that are rapidly and recently quenched after a flat or bursty formation history (depending on the assumed prior); and (3) galaxies that are rapidly and recently quenched after a major starburst. Most recently quenched galaxies show neutral gas outflows, probed by blueshifted Na i D absorption, and ionized gas emission, with line ratios consistent with active galactic nucleus (AGN) diagnostics. This suggests that AGN activity drives multiphase gas outflows, leading to rapid quenching. By tracing back the SFHs of the entire sample, we predict the number density of massive quiescent galaxies at z = 4–6: n = (1.5–6.0) × 10−5 Mpc−3. The two old massive quiescent galaxies in our sample appear to have extremely early formation and quenching (z ≳ 6) and are possibly descendants of early post-starbursts at z > 3. These galaxies still show neutral gas reservoirs and weak H α emission, perhaps because the ejective AGN feedback that caused rapid quenching has weakened over time. [ABSTRACT FROM AUTHOR]

Using a combination of Hubble Space Telescope (HST), JWST, and Atacama Large Millimeter/submillimeter Array (ALMA) data, we perform spatially resolved spectral energy distributions (SED) fitting of fourteen 4 < z < 6 ultraviolet (UV)-selected main-sequence galaxies targeted by the ALMA Large Program [C ii ] Resolved ISM in Star-forming Galaxies. We consistently model the emission from stars and dust in ∼0.5–1 kpc spatial bins to obtain maps of their physical properties. We find no offsets between the stellar masses (M *) and star formation rates (SFRs) derived from their global emission and those from adding up the values in our spatial bins, suggesting there is no bias of outshining by young stars on the derived global properties. We show that ALMA observations are important to derive robust parameter maps because they reduce the uncertainties in L dust (hence, AV and SFR). Using these maps, we explore the resolved star-forming main sequence for z ∼ 5 galaxies, finding that this relation persists in typical star-forming galaxies in the early Universe. We find less obscured star formation where the M * (and SFR) surface densities are highest, typically in the central regions, contrary to the global relation between these parameters. We speculate this could be caused by feedback driving gas and dust out of these regions. However, more observations of IR luminosities with ALMA are needed to verify this. Finally, we test empirical SFR prescriptions based on the UV+IR and [C ii ] line luminosity, finding they work well at the scales probed (approximately kiloparsec). Our work demonstrates the usefulness of joint HST-, JWST-, and ALMA-resolved SED modeling analyses at high redshift. [ABSTRACT FROM AUTHOR]

We investigate the connection between galaxies, dark matter halos, and their large-scale environments at z = 0 with Illustris TNG300 hydrodynamic simulation data. We predict stellar masses from subhalo properties to test two types of machine learning (ML) models: explainable boosting machines (EBMs) with simple galaxy environment features and E (3) -invariant graph neural networks (GNNs). The best-performing EBM models leverage spherically averaged overdensity features on 3 Mpc scales. Interpretations via SHapley Additive exPlanations also suggest that in the context of the TNG300 galaxy–halo connection, simple spherical overdensity on ∼3 Mpc scales is more important than cosmic web distance features measured using the DisPerSE algorithm. Meanwhile, a GNN with connectivity defined by a fixed linking length, L, outperforms the EBM models by a significant margin. As we increase the linking length scale, GNNs learn important environmental contributions up to the largest scales we probe (L = 10 Mpc). We conclude that 3 Mpc distance scales are most critical for describing the TNG galaxy–halo connection using the spherical overdensity parameterization, but that information on larger scales, which is not captured by simple environmental parameters or cosmic web features, can further augment these models. Our study highlights the benefits of using interpretable ML algorithms to explain models of astrophysical phenomena, and the power of using GNNs to flexibly learn complex relationships directly from data while imposing constraints from physical symmetries. [ABSTRACT FROM AUTHOR]

Poststarburst galaxies (PSBs) are young quiescent galaxies that have recently experienced a rapid decrease in star formation, allowing us to probe the fast-quenching period of galaxy evolution. In this work, we obtained Hubble Space Telescope (HST)/WFC3 F110W imaging to measure the sizes of 171 massive ( log (M * / M ⊙) ∼ 11) spectroscopically identified PSBs at 1 < z 1.3 selected from the DESI Survey Validation luminous red galaxy sample. This statistical sample constitutes an order of magnitude increase from the ∼20 PSBs with space-based imaging and deep spectroscopy. We perform structural fitting of the target galaxies with pysersic and compare them to quiescent and star-forming galaxies in the 3D-HST survey. We find that these PSBs are more compact than the general population of quiescent galaxies, lying systematically ∼0.1 dex below the established size–mass relation. However, their central surface mass densities are similar to those of their quiescent counterparts ( log (Σ 1 kpc / (M ⊙ kpc − 2)) ∼ 10.1 ). These findings are easily reconciled by later ex situ growth via minor mergers or a slight progenitor bias. These PSBs are round in projection (b / a median ∼ 0.8), suggesting that they are primarily spheroids, not disks, in 3D. We find no correlation between the time since quenching and light-weighted PSB sizes or central densities. This disfavors apparent structural growth due to the fading of centralized starbursts in this galaxy population. Instead, we posit that the fast quenching of massive galaxies at this epoch occurs preferentially in galaxies with preexisting compact structures. [ABSTRACT FROM AUTHOR]

We investigate the multiphase structure of gas flows in galaxies. We study 80 galaxies during the epoch of peak star formation (1.4 ≤ z ≤ 2.7) using data from the Keck/Low-Resolution Imaging Spectrometer (LRIS) and the Very Large Telescope/ K -Band Multi-Object Spectrograph (KMOS). Our analysis provides a simultaneous probe of outflows using UV emission and absorption features and H α emission. With this unprecedented data set, we examine the properties of gas flows estimated from LRIS and KMOS in relation to other galaxy properties, such as star formation rate (SFR), SFR surface density (ΣSFR), stellar mass (M *), and main-sequence offset (ΔMS). We find no strong correlations between outflow velocity measured from rest-UV line centroids and galaxy properties. However, we find that galaxies with detected outflows show higher averages in SFR, ΣSFR, and ΔMS than those lacking outflow detections, indicating a connection between outflow and galaxy properties. Furthermore, we find a lower average outflow velocity than previously reported, suggesting greater absorption at the systemic redshift of the galaxy. Finally, we detect outflows in 49% of our LRIS sample and 30% in the KMOS sample and find no significant correlation between outflow detection and inclination. These results may indicate that outflows are not collimated and that H α outflows have a lower covering fraction than low-ionization interstellar absorption lines. Additionally, these tracers may be sensitive to different physical scales of outflow activity. A larger sample size with a wider dynamic range in galaxy properties is needed to further test this picture. [ABSTRACT FROM AUTHOR]

Galaxies in the early Universe appear to have grown too big too fast, assembling into massive, monolithic objects more rapidly than anticipated in the hierarchical Lambda cold dark matter (ΛCDM) structure formation paradigm. The available photometric data are consistent with there being a population of massive galaxies that form early (z ≳ 10) and quench rapidly over a short (≲1 Gyr) timescale, consistent with the traditional picture for the evolution of giant elliptical galaxies. Similarly, kinematic observations as a function of redshift show that massive spirals and their scaling relations were in place at early times. Explaining the early emergence of massive galaxies requires either an extremely efficient conversion of baryons into stars at z > 10 or a more rapid assembly of baryons than anticipated in ΛCDM. The latter possibility was explicitly predicted in advance by modified Newtonian dynamics (MOND). We discuss some further predictions of MOND, such as the early emergence of clusters of galaxies and early reionization. [ABSTRACT FROM AUTHOR]

We use resolved spectroscopy from MaNGA to investigate the significance of both local and global properties of galaxies to the cessation of star formation at kpc scales. Quenched regions are identified from a sample of isolated disk galaxies by a single-parameter criterion D n (4000) - log EW(H α) > 1.6 − log 2 = 1.3 , and are divided into gas-rich quenched regions (GRQRs) and gas-poor quenched regions according to the surface density of cold gas (Σgas). Both types of quenched regions tend to be hosted by non-AGN galaxies with relatively high mass (M * ≳ 1010 M ⊙) and red colors (NUV − r ≳ 3), as well as low star formation rate and high central density at fixed mass. They span wide ranges in other properties including structural parameters that are similar to the parent sample, indicating that the conditions responsible for quenching in gas-rich regions are largely independent of the global properties of galaxies. We train random forest classifiers and regressors for predicting quenching in our sample with 15 local/global properties. Σ* is the most important property for quenching, especially for GRQRs. These results strongly indicate the important roles of low-mass hot evolved stars, which are numerous and long-lived in quenched regions and can provide substantial radiation pressure to support the surrounding gas against gravitational collapse. The different feature importance for quenching, as found previously by A. F. L. Bluck et al., is partly due to the different definitions of quenched regions, particularly the different requirements on EW(H α). [ABSTRACT FROM AUTHOR]

We introduce a GPU-accelerated hybrid hydro/ N -body code (Enzo-N) designed to address the challenges of concurrently simulating star clusters and their parent galaxies. This task has been exceedingly challenging, primarily due to the considerable computational time required, which stems from the substantial scale difference between galaxies (∼0.1 Mpc) and star clusters (∼parsecs). Yet, this significant scale separation means that particles within star clusters perceive those outside the star cluster in a semistationary state. By leveraging this aspect, we integrate a direct N -body code (Nbody6++GPU) into a cosmological (magneto)hydrodynamic code (Enzo) through utilization of the semistationary background acceleration approximation. We solve the dynamics of particles within star clusters using a direct N -body solver with regularization for few-body interactions, while evolving particles outside—dark matter, gas, and stars—using a particle-mesh gravity solver and hydrodynamic methods. We demonstrate that Enzo-N successfully simulates the coevolution of star clusters and their parent galaxies, capturing phenomena such as core collapse of the star cluster and tidal stripping due to galactic tides. This comprehensive framework opens up new possibilities for studying the evolution of star clusters within galaxies, offering insights that were previously inaccessible. [ABSTRACT FROM AUTHOR]

Mergers of and interactions between galaxies imprint a wide diversity of morphological, dynamical, and chemical characteristics in stellar halos and tidal streams. Measuring these characteristics elucidates aspects of the progenitors of the galaxies we observe today. The M81 group is the perfect galaxy group to understand the past, present, and future of a group of galaxies in the process of merging. Here, we measure the end of star formation (t 90) and metallicity ([M/H]) of the stellar halo of M82 and the eastern tidal stream of NGC 3077 to: (1) test the idea that M82 possesses a genuine stellar halo, formed before any interaction with M81; (2) determine if NGC 3077's tidal disruption is related to the star formation history in its tails; and (3) create a timeline of the assembly history of the central trio in the M81 group. We argue that M82 possesses a genuine, metal-poor ([M/H] ∼ −1.62 dex) stellar halo, formed from the merger of a small satellite galaxy roughly 6.6 Gyr ago. We also find that the stars present in NGC 3077's tails formed before tidal disruption with M81, and possess a roughly uniform metallicity as shown in S. Okamoto et al., implying that NGC 3077's progenitor had significant population gradients. Finally, we present a timeline of the central trio's merger/interaction history. [ABSTRACT FROM AUTHOR]

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]

Stars in the Galactic disk are born on cold, nearly circular orbits with small vertical excursions. After their birth, their orbits evolve, driven by small- or large-scale perturbations in the Galactic disk's gravitational potential. Here, we study the vertical motions of young stars over their first few orbital periods, using a sample of OBA stars from Gaia E/DR3, which includes radial velocities and ages τ from LAMOST. We construct a parametric model for the time evolution of the stellar orbits' mean vertical actions J z as a function of Galactocentric radius, R GC. Accounting for data uncertainties, we use Markov Chain Monte Carlo analysis in annuli of Galactocentric radius (R GC) to constrain the model parameters. Our best-fit model shows a remarkably linear increase of vertical actions with age across all Galactocentric radii examined. Orbital heating by random scattering could offer a straightforward interpretation for this trend. However, various other dynamical aspects of the Galactic disk, such as stars being born in a warped disk, might offer alternative explanations that could be tested in the future. [ABSTRACT FROM AUTHOR]

We present CloudFlex, an open-source tool for predicting absorption-line signatures of cool gas in galaxy halos with small-scale structure. Motivated by analyses of ∼104 K material in hydrodynamical simulations of turbulent, multiphase media, we model cool gas structures as complexes of cloudlets sampled from a power-law distribution of mass ∝ m cl − α with velocities drawn from a turbulent velocity field. The user may specify α, the lower limit of the cloudlet mass distribution ( m cl , min ), and several other parameters that set the mass, size, and velocity distribution of the complex. This permits investigation of the relation between these parameters and absorption-line observables. As a proof-of-concept, we calculate the Mg ii λ 2796 absorption induced by the cloudlets in background quasi-stellar object (QSO) spectra. We demonstrate that, at fixed metallicity, the covering fraction of sight lines with equivalent widths W 2796 < 0.3 Å increases significantly with decreasing m cl , min , cloudlet number density (n cl), and complex size. We then use this framework to predict the halo-scale W 2796 distribution around ∼ L * galaxies. We show that the observed incidences of W 2796 > 0.3 Å sight lines with impact parameters 10 kpc < R ⊥ < 50 kpc in projected QSO–galaxy studies are consistent with our model over much of parameter space. However, they are underpredicted by models with m cl , min ≥ 100 M ⊙ and n cl ≥ 0.03 cm−3, in keeping with a picture in which the inner cool circumgalactic medium (CGM) is dominated by numerous low-mass cloudlets (m cl ≲ 100 M ⊙) with a volume filling factor ≲1%. When used to model absorption-line data sets built from multi-sight line and/or spatially extended background probes, CloudFlex enables detailed constraints on the size and velocity distributions of structures comprising the photoionized CGM. [ABSTRACT FROM AUTHOR]

We present the galaxy stellar mass–size relation in the rest-frame near-IR (1.5 μ m) and its evolution with redshift up to z = 2.5. Sérsic profiles are measured for ∼26,000 galaxies with stellar masses M ⋆ > 109 M ⊙ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys using coordinates, redshifts, colors, and stellar mass estimates from the COSMOS2020 catalog. The new rest-frame near-IR effective radii are generally smaller than previously measured rest-frame optical sizes, on average by 0.14 dex, with no significant dependence on redshift. For quiescent galaxies, this size offset does not depend on stellar mass, but for star-forming galaxies, the offset increases from −0.1 dex at M ⋆ = 109.5 M ⊙ to −0.25 dex at M ⋆ > 1011 M ⊙. That is, we find that the near-IR stellar mass–size relation for star-forming galaxies is flatter in the rest-frame near-IR than in the rest-frame optical at all redshifts 0.5 < z < 2.5. The general pace of size evolution is the same in the near-IR as previously demonstrated in the optical, with slower evolution (R e ∝ (1 + z)−0.7) for L * star-forming galaxies and faster evolution (R e ∝ (1 + z)−1.3) for L * quiescent galaxies. Massive (M ⋆ > 1011 M ⊙) star-forming galaxies evolve in size almost as fast as quiescent galaxies. Low-mass (M ⋆ < 1010 M ⊙) quiescent galaxies evolve as slow as star-forming galaxies. Our main conclusion is that the size evolution narrative as it has emerged over the past two decades does not radically change when accessing the rest-frame near-IR with JWST, a better proxy of the underlying stellar mass distribution. [ABSTRACT FROM AUTHOR]

UVCANDELS is a Hubble Space Telescope Cycle-26 Treasury Program awarded 164 orbits of primary ultraviolet (UV) F275W imaging and coordinated parallel optical F435W imaging in four CANDELS fields—GOODS-N, GOODS-S, EGS, and COSMOS—covering a total area of ∼426 arcmin2. This is ∼2.7 times larger than the area covered by previous deep-field space UV data combined, reaching a depth of about 27 and 28 ABmag (5 σ in 0."2 apertures) for F275W and F435W, respectively. Along with new photometric catalogs, we present an analysis of the rest-frame UV luminosity function (LF), relying on our UV-optimized aperture photometry method, yielding a factor of 1.5 increase over H-isophot aperture photometry in the signal-to-noise ratios of galaxies in our F275W imaging. Using well-tested photometric redshift measurements, we identify 5810 galaxies at redshifts 0.6 < z < 1, down to an absolute magnitude of M UV = −14.2. In order to minimize the effect of uncertainties in estimating the completeness function, especially at the faint end, we restrict our analysis to sources above 30% completeness, which provides a final sample of 4726 galaxies at −21.5 < M UV < −15.5. We performed a maximum likelihood estimate to derive the best-fit parameters of the UV LF. We report a best-fit faint-end slope of α = − 1.359 − 0.041 + 0.041 at z ∼ 0.8. Creating subsamples at z ∼ 0.7 and z ∼ 0.9, we observe a possible evolution of α with redshift. The unobscured UV luminosity density at M UV < −10 is derived as ρ U V = 1.339 − 0.030 + 0.027 (× 10 26 e r g s − 1 H z − 1 M p c − 3) using our best-fit LF parameters. The new F275W and F435 photometric catalogs from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes. [ABSTRACT FROM AUTHOR]

We present volume-averaged neutral hydrogen fractions x H I and ionized bubble radii R b measured with Ly α damping wing absorption of galaxies at the epoch of reionization. We combine JWST/NIRSpec spectra taken by the CEERS, GO-1433, DDT-2750, and JADES programs and obtain a sample containing 27 bright UV-continuum (M UV < −18.5 mag) galaxies at 7 < z < 12. We construct four composite spectra binned by redshift and find the clear evolution of the softening break toward high redshift at rest-frame 1216 Å, suggesting an increase of Ly α damping wing absorption. We estimate Ly α damping wing absorption in the galaxy spectra with realistic templates including Ly α emission and circumgalactic medium absorptions. Assuming the standard inside-out reionization picture having an ionized bubble with radius R b around a galaxy embedded in the intergalactic medium with x H I, we obtain x H I (R b) values generally increasing (decreasing) from x H I = 0.53 − 0.47 + 0.18 to 0.92 − 0.10 + 0.08 ( log R b = 1.67 − 0.16 + 0.14 to − 0.69 − 0.24 + 0.89 comoving Mpc) at redshift 7.12 − 0.08 + 0.06 to 9.91 − 1.15 + 1.49 . The redshift evolution of x H I indicates a moderately late reionization history consistent with the one previously suggested from the electron scattering of cosmic microwave background and the evolution of UV luminosity function with an escape fraction f esc ∼ 0.2. Our R b measurements suggest that bubble sizes could be up to a few dex larger than the cosmic average values estimated by analytic calculations for a given x H I, while our R b measurements are roughly comparable with the values for merged ionized bubbles around bright galaxies predicted by recent numerical simulations. [ABSTRACT FROM AUTHOR]

We analyzed images of every northern hemisphere Sd galaxy listed in the Third Reference Catalogue of Bright Galaxies with a relatively face-on inclination (θ ≤ 30°). Specifically, we measured the spiral arms' winding angle, ϕ, in 85 galaxies. We applied a novel black hole mass planar scaling relation involving the rotational velocities (from the literature) and pitch angles of each galaxy to predict central black hole masses. This yielded 23 galaxies, each having at least a 50% chance of hosting a central intermediate-mass black hole (IMBH), 102 < M • ≤ 105 M ☉. These 23 nearby (≲50 Mpc) targets may be suitable for an array of follow-up observations to check for active nuclei. Based on our full sample of 85 Sd galaxies, we estimate that the typical Sd galaxy (which tends to be bulgeless) harbors a black hole with log (M • / M ☉) = 6.00 ± 0.14 , but with a 27.7% chance of hosting an IMBH, making this morphological type of galaxy fertile ground for hunting elusive IMBHs. Thus, we find that a ∼106 M ☉ black hole corresponds roughly to the onset of bulge development and serves as a conspicuous waypoint along the galaxy–supermassive black hole coevolution journey. Our survey suggests that >1.22% of bright galaxies (B T ≲ 15.5 mag) in the local Universe host an IMBH (i.e., the "occupation fraction"), which implies a number density >4.96 × 10−6 Mpc−3 for central IMBHs. Finally, we observe that Sd galaxies exhibit an unexpected diversity of properties that resemble the general population of spiral galaxies, albeit with an enhanced signature of the eponymous prototypical traits (i.e., low masses, loosely wound spiral arms, and smaller rotational velocities). [ABSTRACT FROM AUTHOR]

We present statistics of z ∼ 6–9 galaxy outflows indicated by spatially extended gas emission and broad lines. With a total of 61 spectroscopically confirmed galaxies at z ∼ 6–9 in the JWST CEERS, GLASS, and ERO data, we find four galaxies with [O iii ] + H β ionized gas emission significantly extended beyond the kiloparsec-scale stellar components on the basis of the emission line images constructed by the subtraction of NIRCam broadband (line on/off-band) images. By comparison with low- z galaxies, the fraction of galaxies with spatially extended gas, 4/18, at z ∼ 6–9 is an order of magnitude higher than that at z ∼ 0–1, which can be explained by events triggered by frequent major mergers at high redshift. We also investigate medium- and high-resolution NIRSpec spectra of 30 galaxies at z ∼ 6–9, and identify five galaxies with broad (140−800 km s−1) lines in the [O iii ] forbidden line emission, suggestive of galaxy outflows. One galaxy at z = 6.38 shows both spatially extended gas emission and the broad lines, while none of the galaxies with spatially extended gas emission or broad lines present a clear signature of active galactic nuclei (AGN) in either the line diagnostics or type 1 AGN line broadening (>1000 km s−1), which suggests that outflows are mainly driven by stellar feedback. The existence of galaxies with/without spatially extended gas emission or broad lines may be explained by different viewing angles toward outflows or the fact that these are galaxies in the early, late, or post phases of galaxy outflows at high redshift, where relatively large fractions of such galaxies indicate longer-duration and/or more frequent outflows at early cosmic epochs. [ABSTRACT FROM AUTHOR]

We characterize the earliest galaxy population in the JADES Origins Field, the deepest imaging field observed with JWST. We make use of ancillary Hubble Space Telescope optical images (five filters spanning 0.4–0.9 μ m) and novel JWST images with 14 filters spanning 0.8−5 μ m, including seven medium-band filters, and reaching total exposure times of up to 46 hr per filter. We combine all our data at >2.3 μ m to construct an ultradeep image, reaching as deep as ≈31.4 AB mag in the stack and 30.3–31.0 AB mag (5 σ, r = 0.″1 circular aperture) in individual filters. We measure photometric redshifts and use robust selection criteria to identify a sample of eight galaxy candidates at redshifts z = 11.5−15. These objects show compact half-light radii of R 1/2 ∼ 50−200 pc, stellar masses of M ⋆ ∼ 107−108 M ☉, and star formation rates ∼ 0.1−1 M ☉ yr−1. Our search finds no candidates at 15 < z < 20, placing upper limits at these redshifts. We develop a forward-modeling approach to infer the properties of the evolving luminosity function without binning in redshift or luminosity that marginalizes over the photometric redshift uncertainty of our candidate galaxies and incorporates the impact of nondetections. We find a z = 12 luminosity function in good agreement with prior results, and that the luminosity function normalization and UV luminosity density decline by a factor of ∼2.5 from z = 12 to z = 14. We discuss the possible implications of our results in the context of theoretical models for evolution of the dark matter halo mass function. [ABSTRACT FROM AUTHOR]

The measurements of morphological indicators of galaxies are often influenced by a series of observational effects. In this study, we utilize a sample of 832 TNG50 simulated galaxies with log(M */ M ⊙)> 9 at 0.5 < z < 3 to investigate the differences in nonparametric morphological indicators (C, A, Gini, M 20, A O, and D O) derived from noise-free and high-resolution TNG50 images and mock images simulated to have the same observational conditions as JWST/NIRCam. We quantify the relationship between intrinsic and observed values of the morphological indicators and accordingly apply this calibration to 4733 galaxies in the same stellar mass and redshift ranges observed in JWST CEERS and JADES surveys. We find a significant evolution of morphological indicators with rest-frame wavelength (λ rf) at λ rf < 1 μ m, while essentially no obvious variations occur at λ rf > 1 μ m. The morphological indicators of star-forming galaxies (SFGs) and quiescent galaxies (QGs) are significantly different. The morphologies of QGs exhibit a higher sensitivity to rest-frame wavelength than SFGs. After analyzing the evolution of morphological indicators in the rest-frame V band (0.5–0.7 μ m) and rest-frame J band (1.1–1.4 μ m), we find that the morphologies of QGs evolve substantially with both redshift and stellar mass. For SFGs, the C, Gini, and M 20 show a rapid evolution with stellar mass at log(M */ M ⊙) ≥ 10.5, while the A O, D O, and A evolve with both redshift and stellar mass. Our comparison shows that TNG50 simulations effectively reproduce the morphological indicators we measured from JWST observations when the impact of dust attenuation is considered. [ABSTRACT FROM AUTHOR]

By using an ultradeep JWST/MIRI image at 5.6 μ m in the Hubble eXtreme Deep Field, we constrain the role of strong H α emitters (HAEs) during "cosmic reionization" at z ≃ 7–8. Our sample of HAEs is comprised of young (<35 Myr) galaxies, except for one single galaxy (≈300 Myr), with low stellar masses (≲109 M ⊙). These HAEs show a wide range of rest-frame UV continuum slopes (β), with a median value of β = −2.15 ± 0.21, which broadly correlates with stellar mass. We estimate the ionizing photon production efficiency (ξ ion,0) of these sources (assuming f esc,LyC = 0%), which yields a median value log 10 (ξ ion , 0 / (Hz erg − 1)) = 25.50 − 0.12 + 0.10 . We show that ξ ion,0 positively correlates with H α equivalent width and specific star formation rate. Instead ξ ion,0 weakly anticorrelates with stellar mass and β. Based on the β values, we predict f esc , LyC = 4 % − 2 + 3 , which results in log 10 (ξ ion / (Hz erg − 1)) = 25.55 − 0.13 + 0.11 . Considering this and related findings from the literature, we find a mild evolution of ξ ion with redshift. Additionally, our results suggest that these HAEs require only modest escape fractions (f esc,rel) of 6%–15% to reionize their surrounding intergalactic medium. By only considering the contribution of these HAEs, we estimated their total ionizing emissivity ( N ̇ ion ) as N ̇ ion = 10 50.53 ± 0.45 s − 1 Mpc − 3 . When comparing their N ̇ ion with non-HAE galaxies across the same redshift range, we find that that strong, young, and low-mass emitters may have played an important role during cosmic reionization. [ABSTRACT FROM AUTHOR]

Quiescent galaxies generally possess denser cores than star-forming galaxies with similar mass. As a measurement of the core density, the central stellar mass surface density within a radius of 1 kpc (Σ1) was thus suggested to be closely related to galaxy quenching. Massive star-forming galaxies with high Σ1 do not fit into this picture. To understand the origin of such galaxies, we compare the spatially resolved stellar population and star formation properties of massive (>1010.5 M ⊙) blue spiral galaxies with high and low Σ1, divided by Σ1 = 109.4 M ⊙ kpc−2, based on the final release of MaNGA integral field unit data. We find that both high-Σ1 and low-Σ1 blue spirals show large diversities in stellar population and star formation properties. Despite the diversities, high-Σ1 blue spirals are statistically different from the low-Σ1 ones. Specifically, the radial profiles of the luminosity-weighted age and Mgb/〈Fe〉 show that high-Σ1 blue spirals consist of a larger fraction of galaxies with younger and less α -element-enhanced centers than their low-Σ1 counterparts, ∼55% versus ∼30%. The galaxies with younger centers mostly have higher central specific star formation rates, which still follow the spaxel-based star formation main-sequence relation. Examinations of the H α velocity field and the optical structures suggest that galactic bars or galaxy interactions should be responsible for the rejuvenation of these galaxies. The remaining ∼45% of high-Σ1 blue spirals are consistent with the inside-out growth scenario. [ABSTRACT FROM AUTHOR]

The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use smoothed particle hydrodynamics simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L* galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume Romulus25, with stellar masses between log(M */ M ⊙) = 9.5–11.5. We measure the fraction of metals remaining in the interstellar medium (ISM) and CGM of each galaxy and calculate the expected mass of each SMBH based on the M BH– σ relation (Kormendy & Ho 2013). The deviation of each SMBH from its expected mass, Δ M BH, is compared to the potential of its host via σ. We find that SMBHs with accreted mass above M BH– σ are more effective at removing metals from the ISM than undermassive SMBHs in star-forming galaxies. Overall, overmassive SMBHs suppress the total star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little to no evacuation of gas from the CGM out of their halos, in contrast with other simulations. Finally, we predict that C iv column densities in the CGM of L* galaxies are unlikely to depend on host galaxy SMBH mass. Our results show that the scatter in the low-mass end of the M BH– σ relation may indicate how effective an SMBH is in the local redistribution of mass in its host galaxy. [ABSTRACT FROM AUTHOR]

The large-scale environment of the cosmic web is believed to impact galaxy evolution, but there is still no consensus regarding the mechanisms. We use a semi-analytic model (SAM) galaxy catalog to study the star formation and dust content of local galaxies in different cosmic environments of the cosmic web—namely voids, filaments, walls, and nodes. We find a strong impact of the environment only for galaxies with M stars ≲ 1010.8 M ⊙: the less dense the environment, the larger the star formation rate and dust content at fixed stellar mass. This is attributed to the fact that galaxies in less dense environments typically feature younger stellar populations, a slower evolution of their stellar mass, and delayed star formation compared to galaxies in denser environments. As for galaxies with M stars ≳ 1010.8 M ⊙, the differences among environments are milder, due to the disk-instability-driven supermassive black hole (SMBH) growth implemented in the SAM, which makes SMBH growth, and thus galaxy quenching, environment-insensitive. We qualitatively test our predictions against observations by identifying environments in Sloan Digital Sky Survey Data Release 16 using dust masses derived from the GAMA survey. The agreement is encouraging, particularly at log M stars / M ⊙ ≳ 10.5 – 11 , where the specific star formation rates and dust masses appear quite environment-insensitive. This result confirms the importance of in situ growth channels of SMBHs. [ABSTRACT FROM AUTHOR]

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

Galaxy interactions and mergers can enhance or reduce star formation, but a complete understanding of the involved processes is still lacking. The effect of dwarf galaxy mergers is even less clear than their massive counterpart. We present a study on a dwarf merger remnant in the Virgo cluster, VCC322, which might form a triple system with VCC334 and VCC319. We identify a prominent long and straight tail-like substructure that has a size comparable to its host galaxy VCC322. By comparing the color–color (g − r versus r − H) distribution with simple stellar population models, we infer that the metallicity and stellar age of this tail are Z ⋆ ∼ 0.02 Z ⊙ and t ⋆ ∼ 10 Gyr, respectively. In VCC319, we find a sign of isophotal twisting. This suggests that VCC319 may be subject to tidal interaction. An analysis of the Sloan Digital Sky Survey optical spectra of VCC322 indicates mass- and light-weighted ages of about 109.8 yr and 107.5 yr, respectively, indicating an ongoing star formation activity. However, the star formation in VCC322 seems suppressed when compared to other star-forming dwarfs of comparable stellar masses. Our finding of shock excitation of optical emission lines indicates that interaction-induced shock may contribute to the heating of cold gas and suppression of star formation. [ABSTRACT FROM AUTHOR]

We present a multiple emission line study of ∼1300 H α emitters (HAEs) at z ∼ 2.3 in the ZFOURGE survey. In contrast to the traditional spectroscopic method, our sample is selected based on the flux excess in the ZFOURGE K s broadband data relative to the best-fit stellar continuum. Using the same method, we also extract the strong diagnostic emission lines for these individual HAEs: [O iii ] λ λ 4959, 5007 and [O ii ] λ λ 3726, 3729. Our measurements demonstrate good consistency with those obtained from spectroscopic surveys. We investigate the relationship between the equivalent widths (EWs) of these emission lines and various galaxy properties, including stellar mass, stellar age, star formation rate, specific star formation rate, and ionization state (O32). We have identified a discrepancy between HAEs at z ∼ 2.3 and typical local star-forming galaxies observed in the Sloan Digital Sky Survey, suggesting the evolution of lower gas-phase metallicity (Z) and higher ionization parameters (U) with redshift. Notably, we have observed a significant number of low-mass HAEs exhibiting exceptionally high EW[O iii ]. Their galaxy properties are comparable to those of extreme objects, such as extreme O3 emitters and Ly α emitters at z ≃ 2–3. Considering that these characteristics may indicate potential strong Lyman continuum leakage, higher-redshift analogs of the low-mass HAEs could be significant contributors to the cosmic reionization. Further investigations of this particular population are required to gain a clearer understanding of galaxy evolution and cosmic reionization. [ABSTRACT FROM AUTHOR]

We have investigated the properties (e.g., age, metallicity) of the stellar populations of a ringlike tidal stellar stream (or streams) around the edge-on galaxy SPRC047 (z = 0.031) using spectral energy distribution (SED) fits to integrated broadband aperture flux densities. We used visual images in six different bands and Spitzer/IRAC 3.6 μ m data. We have attempted to derive best-fit stellar population parameters (metallicity, age) in three noncontiguous segments of the stream. Due to the very low surface brightness of the stream, we have performed a deconvolution with a Richardson–Lucy–type algorithm of the low spatial resolution 3.6 μ m IRAC image, thereby reducing the effect of the point-spread function aliased emission from the bright edge-on central galaxy at the locations of our three stream segments. Our SED fits that used several different star formation (SF) history priors, from an exponentially decaying SF burst to continuous SF, indicate that the age–metallicity–dust degeneracy is not resolved, most likely because of inadequate wavelength coverage and low signal-to-noise ratios of the low surface brightness features. We also discuss how future deep visual–near-infrared observations, combined with absolute flux calibration uncertainties at or below the 1% level, complemented by equally well absolute flux-calibrated observations in ultraviolet and mid-infrared bands, would improve the accuracy of broadband SED fitting results for low surface brightness targets, such as stellar streams around nearby galaxies that are not resolved into stars. [ABSTRACT FROM AUTHOR]

We present an analysis of the quenching of star formation in massive galaxies (M * > 109.5 M ⊙) within the first 0.5–3 Gyr of the Universe's history utilizing JWST-CEERS data. We utilize a combination of advanced statistical methods to accurately constrain the intrinsic dependence of quenching in a multidimensional and intercorrelated parameter space. Specifically, we apply random forest classification, area statistics, and a partial correlation analysis to the JWST-CEERS data. First, we identify the key testable predictions from two state-of-the-art cosmological simulations (IllustrisTNG and EAGLE). Both simulations predict that quenching should be regulated by supermassive black hole mass in the early Universe. Furthermore, both simulations identify the stellar potential (ϕ *) as the optimal proxy for black hole mass in photometric data. In photometric observations, where we have no direct constraints on black hole masses, we find that the stellar potential is the most predictive parameter of massive galaxy quenching at all epochs from z = 0–8, exactly as predicted by simulations for this sample. The stellar potential outperforms stellar mass, galaxy size, galaxy density, and Sérsic index as a predictor of quiescence at all epochs probed in JWST-CEERS. Collectively, these results strongly imply a stable quenching mechanism operating throughout cosmic history, which is closely connected to the central gravitational potential in galaxies. This connection is explained in cosmological models via massive black holes forming and growing in deep potential wells, and subsequently quenching galaxies through a mix of ejective and preventative active galactic nucleus feedback. [ABSTRACT FROM AUTHOR]

We present a star formation rate function (SFRF) at z ∼ 6 based on star formation rates (SFRs) derived by spectral energy distribution fitting on data from rest-frame UV to optical wavelengths of galaxies in the CANDELS GOODS-South and North fields. The resulting SFRF shows an excess compared to the previous estimations by using rest-frame UV luminosity functions (LFs) corrected for the dust attenuation and is comparable to that estimated from a far-infrared LF. This suggests that the number density of dust-obscured intensively star-forming galaxies at z ∼ 6 has been underestimated in the previous approach based only on rest-frame UV observations. We parameterize the SFRF using the Schechter function and obtain the best-fit parameter of the characteristic SFR (SFR*) when the faint-end slope and characteristic number density are fixed. The best-fit SFR* at z ∼ 6 is comparable to that at z ∼ 2, when the cosmic star formation activity reaches its peak. Together with SFRF estimations with a similar approach using rest-frame UV to optical data, the SFR* is roughly constant from z ∼ 2 to ∼6 and may decrease above z ∼ 6. Since the SFR* is sensitive to the high-SFR end of the SFRF, this evolution of SFR* suggests that the high-SFR end of the SFRF grows rapidly during the epoch of reionization and reaches a similar level observed at z ∼ 2. [ABSTRACT FROM AUTHOR]

We present the first results from a high-resolution simulation with a focus on galactic wind driving for an isolated galaxy with a halo mass of ∼1011 M ⊙ (similar to the Large Magellanic Cloud) and a total gas mass of ∼6 × 108 M ⊙, resulting in ∼108 gas cells at ∼4 M ⊙ mass resolution. We adopt a resolved stellar feedback model with nonequilibrium cooling and heating, including photoelectric heating and photoionizing radiation, as well as supernovae, coupled to the second-order meshless finite-mass method for hydrodynamics. These features make this the largest resolved interstellar medium (ISM) galaxy model run to date. We find mean star formation rates around 0.05 M ⊙ yr−1 and evaluate typical time-averaged loading factors for mass (η M ∼ 1.0, in good agreement with recent observations) and energy (η E ∼ 0.01). The bulk of the mass of the wind is transported by the warm (T < 5 × 105 K) phase, while there is a similar amount of energy transported in the warm and the hot phases (T > 5 × 105 K). We find an average opening angle of 30° for the wind, decreasing with higher altitude above the midplane. The wind mass loading is decreasing (flat) for the warm (hot) phase as a function of the star formation surface rate density ΣSFR, while the energy loading shows inverted trends with ΣSFR, decreasing for the warm wind and increasing for the hot wind, although with very shallow slopes. These scalings are in good agreement with previous simulations of resolved wind driving in the multiphase ISM. [ABSTRACT FROM AUTHOR]

Quenching of star formation in the central galaxies of cosmological halos is thought to result from energy released as gas accretes onto a supermassive black hole. The same energy source also appears to lower the central density and raise the cooling time of baryonic atmospheres in massive halos, thereby limiting both star formation and black hole growth, by lifting the baryons in those halos to greater altitudes. One predicted signature of that feedback mechanism is a nearly linear relationship between the central black hole's mass (M BH) and the original binding energy of the halo's baryons. We present the increasingly strong observational evidence supporting a such a relationship, showing that it extends up to halos of mass M halo ∼ 1014 M ⊙. We then compare current observational constraints on the M BH– M halo relation with numerical simulations, finding that black hole masses in IllustrisTNG appear to exceed those constraints at M halo < 1013 M ⊙ and that black hole masses in EAGLE fall short of observations at M halo ∼ 1014 M ⊙. A closer look at IllustrisTNG shows that quenching of star formation and suppression of black hole growth do indeed coincide with black hole energy input that lifts the halo's baryons. However, IllustrisTNG does not reproduce the observed M BH– M halo relation because its black holes gain mass primarily through accretion that does not contribute to baryon lifting. We suggest adjustments to some of the parameters in the IllustrisTNG feedback algorithm that may allow the resulting black hole masses to reflect the inherent links between black hole growth, baryon lifting, and star formation among the massive galaxies in those simulations. [ABSTRACT FROM AUTHOR]

A key question in galaxy evolution has been the importance of the apparent "clumpiness" of high-redshift galaxies. Until now, this property has been primarily investigated in rest-frame UV, limiting our understanding of their relevance. Are they short-lived or are they associated with more long-lived massive structures that are part of the underlying stellar disks? We use JWST/NIRCam imaging from the Cosmic Evolution and Epoch of Reionization Survey to explore the connection between the presence of these "clumps" in a galaxy and its overall stellar morphology, in a mass-complete ( log M * / M ⊙ > 10.0) sample of galaxies at 1.0 < z < 2.0. Exploiting the uninterrupted access to rest-frame optical and near-IR light, we simultaneously map the clumps in galactic disks across our wavelength coverage, along with measuring the distribution of stars among their bulges and disks. First, we find that the clumps are not limited to the rest-frame UV and optical, but are also apparent in near-IR with ∼60% spatial overlap. This rest-frame near-IR detection indicates that clumps would also feature in the stellar-mass distribution of the galaxy. A secondary consequence is that these will hence be expected to increase the dynamical friction within galactic disks leading to gas inflow. We find a strong negative correlation between how clumpy a galaxy is and strength of the bulge. This firmly suggests an evolutionary connection, either through clumps driving bulge growth or the bulge stabilizing the galaxy against clump formation, or a combination of the two. Finally, we find evidence of this correlation differing from rest-frame optical to near-IR, which could suggest a combination of varying formation modes for the clumps. [ABSTRACT FROM AUTHOR]

The nearby dwarf galaxy POX 52 at z = 0.021 hosts an active galactic nucleus (AGN) with a black hole (BH) mass of M BH ∼ 105–6 M ⊙ and an Eddington ratio of ∼0.1–1. This object provides the rare opportunity to study both AGN and host-galaxy properties in a low-mass highly accreting system. To do so, we collected its multiwavelength data from X-ray to radio. First, we construct a spectral energy distribution, and by fitting it with AGN and host-galaxy components, we constrain AGN-disk and dust-torus components. Then, while considering the AGN-disk emission, we decompose optical Hubble Space Telescope images. As a result, it is found that a classical bulge component is probably present, and its mass (M bulge) is consistent with an expected value from a local relation. Lastly, we analyze new quasi-simultaneous X-ray (0.2–30 keV) data obtained by the Nuclear Spectroscopic Telescope Array and XMM-Newton. The X-ray spectrum can be reproduced by multicolor blackbody, warm and hot coronae, and disk and torus reflection components. Based on this, the spin is estimated to be a spin = 0.998−0.814, which could suggest that most of the current BH mass was achieved by prolonged mass accretion. Given the presence of the bulge, POX 52 would have undergone a galaxy merger, while the M BH– M bulge relation and the inferred prolonged accretion could suggest that AGN feedback occurred. Regarding the AGN structure, the spectral slope of the hot corona, its relative strength to the bolometric emission, and the torus structure are found to be consistent with Eddington-ratio dependencies found for nearby AGNs. [ABSTRACT FROM AUTHOR]

Submillimeter observations reveal the star formation activity obscured by dust in the young Universe. It still remains unclear how galaxies detected at submillimeter wavelengths are related to ultraviolet/optical-selected galaxies in terms of their observed quantities, physical properties, and evolutionary stages. Deep near- and mid-infrared observational data are crucial to characterize the stellar properties of galaxies detected with submillimeter emission. In this study, we make use of a galaxy catalog from the Spitzer Matching survey of the UltraVISTA ultra-deep Stripes. By crossmatching with a submillimeter source catalog constructed with archival data of the Atacama Large Millimeter/submillimeter Array (ALMA), we search for galaxies at z > 2 with a submillimeter detection in our galaxy catalog. We find that the ALMA-detected galaxies at z > 2 are systematically massive and have redder K s –[4.5] colors than the nondetected galaxies. The redder colors are consistent with the larger dust reddening values of the ALMA-detected galaxies obtained from spectral energy distribution (SED) fitting. We also find that the ALMA-detected galaxies tend to have brighter 4.5 μ m magnitudes. This may suggest that they tend to have smaller mass-to-light ratios and thus to be younger than star-forming galaxies fainter at submillimeter wavelengths with similar stellar masses. We identify starburst galaxies with high specific star formation rates among both ALMA-detected and nondetected SMUVS sources. Irrespective of their brightness at submillimeter wavelengths, these populations have similar dust reddening values, which may suggest a variety of dust SED shapes among the starburst galaxies at z > 2. [ABSTRACT FROM AUTHOR]

The assembly of galaxies over cosmic time is tightly connected to the assembly of their host dark matter halos. We investigate the stellar mass growth history and the chemical enrichment history of central galaxies in the Sloan Digital Sky Survey Mapping Nearby Galaxies at Apache Point Observatory. We find that the derived stellar metallicity of passive central galaxies is always higher than that of the star-forming ones. This stellar metallicity enhancement becomes progressively larger toward low-mass galaxies (at a given epoch) and earlier epochs (at a given stellar mass), which suggests strangulation as the primary mechanism for star formation quenching in central galaxies not only in the local Universe but also very likely at higher redshifts up to z ∼ 3. We show that at the same present-day stellar mass, passive central galaxies assembled half of their final stellar mass ∼2 Gyr earlier than star-forming central galaxies, which agrees well with the semi-analytic model. Exploring the semi-analytic model, we find that this is because passive central galaxies reside in, on average, more massive halos with a higher halo mass increase rate across cosmic time. As a consequence, passive central galaxies are assembled faster and also quenched earlier than their star-forming counterparts. While at the same present-day halo mass, different halo assembly history also produces a very different final stellar mass of the central galaxy within, and halos assembled earlier host more massive centrals with a higher quenched fraction, in particular around the "golden halo mass" at 1012 M ⊙. Our results call attention back to the dark matter halo as a key driver of galaxy evolution. [ABSTRACT FROM AUTHOR]

Dual active galaxy nuclei (dAGNs) trace the population of post-merger galaxies and are the precursors to massive black hole (MBH) mergers, an important source of gravitational waves that may be observed by the Laser Interferometer Space Antenna (LISA). In Paper I of this series, we used the population of ≈2000 galaxy mergers predicted by the TNG50-3 simulation to seed semi-analytic models of the orbital evolution and coalescence of MBH pairs with initial separations of ≈1 kpc. Here, we calculate the dAGN luminosities and separations of these pairs as they evolve in post-merger galaxies, and show how the coalescence fraction of dAGNs changes with redshift. We find that because of the several gigayear-long dynamical friction timescale for orbital evolution, the fraction of dAGNs that eventually end in an MBH merger grows with redshift and exceeds 50% beyond z dAGN ≈ 1. Dual AGNs in galaxies with bulge masses ≲1010 M ⊙, or consisting of near-equal-mass MBHs, evolve more quickly and have higher than average coalescence fractions. At any redshift, dAGNs observed with small separations (≲0.7 kpc) have a higher probability of merging before z = 0 than more widely separated systems. Radiation feedback effects can significantly reduce the number of MBH mergers, and this could be manifested as a larger than expected number of widely separated dAGNs. We present a method to estimate the MBH coalescence rate as well as the potential LISA detection rate given a survey of dAGNs. Comparing these rates to the eventual LISA measurements will help determine the efficiency of dynamical friction in post-merger galaxies. [ABSTRACT FROM AUTHOR]

The "tuning-fork" (TF) analysis of CO and H α emission has been used to estimate the lifetimes of molecular clouds in nearby galaxies. With simple model calculations, we show that this analysis does not necessarily estimate cloud lifetimes, but instead captures a duration of the cloud evolutionary cycle, from dormant to star-forming, and then back to a dormant phase. We adopt a hypothetical setup in which molecular clouds (e.g., traced in CO) live forever and form stars (e.g., H ii regions) at some frequency, which then drift away from the clouds. The TF analysis still returns a timescale for the immortal clouds. This model requires drifting motion to separate the newborn stars from the clouds, and we discuss its origin. We also discuss the physical origin of the characteristic spatial separation term in the TF analysis and a bias due to systematic error in the determination of the reference timescale. [ABSTRACT FROM AUTHOR]

Using N -body simulations, we explore the effects of growing a supermassive black hole (SMBH) prior to or during the formation of a stellar bar. Keeping the final mass and growth rate of the SMBH fixed, we show that if it is introduced before or while the bar is still growing, the SMBH does not cause a decrease in bar amplitude. Rather, in most cases, it is strengthened. In addition, an early-growing SMBH always either decreases the buckling amplitude, delays buckling, or both. This weakening of buckling is caused by an increase in the disk vertical velocity dispersion at radii well beyond the nominal black hole sphere of influence. While we find considerable stochasticity and sensitivity to initial conditions, the only case where the SMBH causes a decrease in bar amplitude is when it is introduced after the bar has attained a steady state. In this case, we confirm previous findings that the decrease in bar strength is a result of scattering of bar-supporting orbits with small pericenter radii. By heating the inner disk both radially and vertically, an early-growing SMBH increases the fraction of stars that can be captured by the inner Lindblad resonance (ILR) and the vertical ILR, thereby strengthening both the bar and the boxy-peanut-shaped bulge. Using orbital frequency analysis of star particles, we show that when an SMBH is introduced early and the bar forms around it, the bar is populated by different families of regular bar-supporting orbits than when the bar forms without an SMBH. [ABSTRACT FROM AUTHOR]

Recent cosmological hydrodynamical simulations have produced populations of numerical galaxies whose global star-forming properties are in good agreement with those of observed galaxies. Proper modeling of energetic feedback from supernovae and active galactic nuclei is critical to the ability of simulations to reproduce observed galaxy properties, and historically, such modeling has proven to be a challenge. Here, we analyze the local properties of central and satellite galaxies in the z = 0 snapshot of the TNG100 simulation as a test of feedback models. We generate a face-on projection of stellar particles in TNG100 galaxies, from which we demonstrate the existence of a resolved star-forming main sequence (ΣSFR–Σ* relation) with a slope and normalization that is in reasonable agreement with previous studies. We also present radial profiles of various galaxy populations for two parameters: the distance from the resolved main-sequence line (ΔΣSFR) and the luminosity-weighted stellar age (Age L ). We find that, on average, high-mass central and satellite galaxies quench from the inside out, while low-mass central and satellite galaxies have similar, flatter profiles. Overall, we find that, with the exception of the starburst population, the TNG100 feedback models yield simulated galaxies whose radial distributions of Age L and ΔΣSFR agree with those of observed galaxies. [ABSTRACT FROM AUTHOR]

As the progenitors of present-day galaxy clusters, protoclusters are excellent laboratories to study galaxy evolution. Since existing observations of protoclusters are limited to the detected constituent galaxies at UV and/or infrared wavelengths, the details of how typical galaxies grow in these young, pre-virialized structures remain uncertain. We measure the total stellar mass and star formation within protoclusters, including the contribution from faint undetected members by performing a stacking analysis of 211 z = 2–4 protoclusters selected as Planck cold sources. We stack Wide-field Infrared Survey Explorer and Herschel/SPIRE images to measure the angular size and the spectral energy distribution of the integrated light from the protoclusters. The fluxes of protoclusters selected as Planck cold sources can be contaminated by line-of-sight interlopers. Using the WebSky simulation, we estimate that a single protocluster contributes 33% ± 15% of the flux of a Planck cold source on average. After this correction, we obtain a total star formation rate of 7.3 ± 3.2 × 103 M ⊙yr−1 and a total stellar mass of 4.9 ± 2.2 × 1012 M ⊙. Our results indicate that protoclusters have, on average, 2× more star formation and 4× more stellar mass than the total contribution from individually detected galaxies in spectroscopically confirmed protoclusters. This suggests that much of the total flux within z = 2–4 protoclusters comes from galaxies with luminosities lower than the detection limit of SPIRE (L IR < 3 × 1012 L ⊙). Lastly, we find that protoclusters subtend a half-light radius of 2.′8 (4.2–5.8 cMpc), which is consistent with simulations. [ABSTRACT FROM AUTHOR]

Mapping of multiple lines such as the fine-structure emission from [C ii ] (157.7 μ m), [O iii ] (52 and 88.4 μ m), and rotational emission lines from CO are of particular interest for upcoming line intensity mapping (LIM) experiments at millimeter wavelengths, due to their brightness features. Several upcoming experiments aim to cover a broad range of scientific goals, from detecting signatures of the epoch of reionization to the physics of star formation and its role in galaxy evolution. In this paper, we develop a semianalytic approach to modeling line strengths as functions of the star formation rate (SFR) or infrared luminosity based on observations of local and high- z galaxies. This package, LIMpy (Line Intensity Mapping in Python), estimates the intensity and power spectra of [C ii ], [O iii ], and CO rotational transition lines up to the J levels (1–0) to (13–12) based both on analytic formalism and on simulations. We develop a relation among halo mass, SFR, and multiline intensities that permits us to construct a generic formula for the evolution of several line strengths up to z ∼ 10. We implement a variety of star formation models and multiline luminosity relations to estimate the astrophysical uncertainties on the intensity power spectrum of these lines. As a demonstration, we predict the signal-to-noise ratio of [C ii ] detection for an EoR-Spec-like instrument on the Fred Young Submillimeter Telescope. Furthermore, the ability to use any halo catalog allows the LIMpy code to be easily integrated into existing simulation pipelines, providing a flexible tool to study intensity mapping in the context of complex galaxy formation physics. [ABSTRACT FROM AUTHOR]

Using deep near-infrared Keck/MOSFIRE observations, we analyze the rest-optical spectra of eight star-forming galaxies in the COSMOS and GOODS-N fields. We reach integration times of ∼10 hr in the deepest bands, pushing the limits on current ground-based observational capabilities. The targets fall into two redshift bins, of five galaxies at z ∼ 1.7 and three galaxies at z ∼ 2.5, and were selected as likely to yield significant auroral-line detections. Even with long integration times, detection of the auroral lines remains challenging. We stack the spectra together into subsets based on redshift, improving the signal-to-noise ratio on the [O iii ] λ 4364 auroral emission line and, in turn, enabling a direct measurement of the oxygen abundance for each stack. We compare these measurements to commonly employed strong-line ratios alongside measurements from the literature. We find that the stacks fall within the distribution of z > 1 literature measurements, but a larger sample size is needed to robustly constrain the relationships between strong-line ratios and oxygen abundance at high redshift. We additionally report detections of [O i ] λ 6302 for nine individual galaxies and composite spectra of 21 targets in the MOSFIRE pointings. We plot their line ratios on the [O iii ] λ 5008/H β versus [O i ] λ 6302/H α diagnostic diagram, comparing our targets to local galaxies and H ii regions. We find that the [O i ]/H α ratios in our sample of galaxies are consistent with being produced in gas ionized by α -enhanced massive stars, as has been previously inferred for rapidly forming galaxies at early cosmic times. [ABSTRACT FROM AUTHOR]

This study examines the impact of cluster environments on galaxy properties using data from the Hyper Suprime-Cam Subaru Strategic Program and an optically selected CAMIRA cluster sample. Specifically, the study analyzes the fractions of quiescent and green valley galaxies with stellar masses above 108.6 M ⊙ at z ∼ 0.2 and 109.8 M ⊙ at z ∼ 1.1, investigating their trends in radius and density. The results indicate that a slow quenching mechanism is at work, as evidenced by a radially independent specific star formation rate reduction of 0.1 dex for star-forming galaxies in a cluster environment. The study also finds that slow quenching dominates fast quenching only for low-mass galaxies (<109.2 M ⊙) near the cluster edge, based on their contributions to the quiescent fraction. After controlling for M *, z, and local overdensity, the study still finds a significant radial gradient in the quiescent fraction, indicating active ram pressure stripping in dense environments. That said, analyzing the density trend of the quiescent fraction with other fixed parameters suggests that radial and density-related quenching processes are equally crucial for low-mass cluster galaxies. The study concludes that ram pressure stripping is the primary environmental quenching mechanism for high stellar mass galaxies in clusters. By contrast, ram pressure stripping and density-related quenching processes act comparably for low-mass cluster galaxies around the center. Near the cluster boundary, starvation and harassment become the leading quenching processes for low stellar mass galaxies. [ABSTRACT FROM AUTHOR]