Individual Galaxies Research Articles

Gas metallicity distributions in SDSS-IV MaNGA galaxies: what drives gradients and local trends?

ABSTRACTThe gas metallicity distributions across individual galaxies and across galaxy samples can teach us much about how galaxies evolve. Massive galaxies typically possess negative metallicity gradients, and mass and metallicity are tightly correlated on local scales over a wide range of galaxy masses; however, the precise origins of such trends remain elusive. Here, we employ data from SDSS-IV MaNGA to explore how gas metallicity depends on the local stellar mass density and on galactocentric radius within individual galaxies. We also consider how the strengths of these dependencies vary across the galaxy mass-size plane. We find that radius is more predictive of local metallicity than stellar mass density in extended lower-mass galaxies, while we find density and radius to be almost equally predictive in higher-mass and more compact galaxies. Consistent with previous work, we find a mild connection between metallicity gradients and large-scale environment; however, this is insufficient to explain variations in gas metallicity behaviour across the mass-size plane. We argue our results to be consistent with a scenario in which extended galaxies have experienced smooth gas accretion histories, producing negative metallicity gradients over time. We further argue that more compact and more massive systems have experienced increased merging activity that disrupts this process, leading to flatter metallicity gradients and more dominant density-metallicity correlations within individual galaxies.

The contribution of magnetized galactic outflows to extragalactic Faraday rotation

ABSTRACT Galactic outflows driven by star formation and active galactic nuclei blow bubbles into their local environments, causing galactic magnetic fields to be carried into intergalactic space. We explore the redshift-dependent effect of these magnetized bubbles on the Faraday rotation measure (RM) of extragalactic radio sources. Using the IllustrisTNG cosmological simulations, we separate the contribution from magnetic bubbles from that of the volume-filling magnetic component expected to be due to the seed field originating in the early universe. We use this separation to extract the redshift dependence of each component and to compare IllustrisTNG model predictions with observation measurements of the NRAO VLA sky survey (NVSS). We find that magnetized bubbles provide a sizeable contribution to the extragalactic RM, with redshift-independent 〈|RM|〉 ≃ 13 rad/m2 for sources at redshifts z ≥ 2. This is close to the mean residual RM of 16 rad/m2 found from NVSS data in this redshift range. Using the IllustrisTNG simulations, we also evaluate a simple model for the contribution to residual RM from individual host galaxies and show that this contribution is negligible at high-redshift. While the contribution from magnetic bubbles in the IllustrisTNG model is currently compatible with observational measurements of residual RM, the next-generation RM sky surveys, which will be free from the wrapping uncertainty, have larger statistics and better sensitivity should be able to observe predicted flat contribution from magnetic bubbles at large redshifts. This should allow to experimentally probe magnetic bubbles and check models of galaxy feedback in cosmological simulations.

EDIG-CHANGES I: extended Hα emission from the extraplanar diffuse ionized gas (eDIG) around CHANG-ES galaxies

ABSTRACT The extraplanar diffuse ionized gas (eDIG) represents the cool/warm ionized gas reservoir around galaxies. We present spatial analysis of the Hα images of 22 nearby edge-on spiral galaxies taken with the Apache Point Observatory 3.5-m telescope (eDIG-CHANGES). We conduct an exponential fit to the vertical Hα intensity profiles of the galaxies, of which 16 can be decomposed into thin + thick disk components. The median value of the Hα scale height of the thick disk is $\langle h_{\rm H\alpha }\rangle =1.13\pm 0.14\rm ~kpc$. We further examine the dependence of hHα on the stellar mass, SFR, and SFR surface density (SFRSD) of the galaxies. We find a tight sublinear correlation between hHα and SFR, expressed in hHα ∝ SFRα, where α ≈ 0.29. Moreover, the offset of individual galaxies from the best-fit SFR-hHα relation, expressed in hHα/SFRα, shows significant anti-correlation with SFRSD. We further compare the vertical extension of the eDIG to multi-wavelength measurements of other CGM phases. We find the eDIG slightly more extended than the neutral gas. This indicates the existence of some extended ionizing sources, in addition to the leaking photons from the disk star formation regions. Most galaxies have an X-ray scale height smaller than Hα, suggesting the majority of the X-ray photons are actually from the thick disk instead of the extended CGM. hHα is comparable to the L-band radio continuum scale height. This indicates that the thermal and non-thermal electrons have similar spatial distributions, a natural result if both are transported outwards by a galactic wind. This further indicates the thermal gas, cosmic rays, and magnetic field may be close to energy equipartition.

Intracluster light is already abundant at redshift beyond unity.

Intracluster light (ICL) is diffuse light from stars that are gravitationally bound not to individual member galaxies, but to the halo of galaxyclusters. Leading theories1,2 predict that the ICL fraction, defined by the ratio of the ICL to the total light, rapidly decreases with increasing redshift, to the level of a few per cent at z > 1. However, observational studies have remained inconclusive about the fraction beyond redshift unity because, to date, only two clusters in this redshift regime have been investigated. One shows a much lower fraction than the mean value at low redshift3, whereas the other possesses a fraction similar to the low-redshift value4. Here we report an ICL study of ten galaxy clusters at 1 ≲ z ≲ 2 based on deep infrared imaging data. Contrary to the leading theories, our study finds that ICL is already abundant at z ≳ 1, with a mean ICL fraction of approximately 17%. Moreover, no significant correlation between cluster mass and ICL fraction or between ICL colour and cluster-centric radius is observed. Our findings suggest that gradual stripping can no longer be the dominant mechanism of ICL formation. Instead, our study supports the scenario wherein the dominant ICL production occurs in tandem with the formation and growth of the brightest cluster galaxies and/or through the accretion of preprocessed stray stars.

CLEAR: Spatially Resolved Emission Lines and Active Galactic Nuclei at 0.6 < z < 1.3

We investigate spatially resolved emission-line ratios in a sample of 219 galaxies (0.6 < z < 1.3) detected using the G102 grism on the Hubble Space Telescope Wide Field Camera 3 taken as part of the CANDELS Lyα Emission at Reionization survey to measure ionization profiles and search for low-luminosity active galactic nuclei (AGN). We analyze [O iii] and Hβ emission-line maps, enabling us to spatially resolve the [O iii]/Hβ emission-line ratio across the galaxies in the sample. We compare the [O iii]/Hβ ratio in galaxy centers and outer annular regions to measure ionization differences and investigate the potential of sources with nuclear ionization to host AGN. We investigate some of the individual galaxies that are candidates to host strong nuclear ionization and find that they often have low stellar mass and are undetected in X-rays, as expected for low-luminosity AGN in low-mass galaxies. We do not find evidence for a significant population of off-nuclear AGN or other clumps of off-nuclear ionization. We model the observed distribution of [O iii]/Hβ spatial profiles and find that most galaxies are consistent with a small or zero difference between their nuclear and off-nuclear line ratios, but 6%–16% of galaxies in the sample are likely to host nuclear [O iii]/Hβ that is ∼0.5 dex higher than in their outer regions. This study is limited by large uncertainties in most of the measured [O iii]/Hβ spatial profiles; therefore, deeper data, e.g., from deeper HST/WFC3 programs or from JWST/NIRISS, are needed to more reliably measure the spatially resolved emission-line conditions of individual high-redshift galaxies.

Hierarchical Bayesian Inference of Photometric Redshifts with Stellar Population Synthesis Models

We present a Bayesian hierarchical framework to analyze photometric galaxy survey data with stellar population synthesis (SPS) models. Our method couples robust modeling of spectral energy distributions with a population model and a noise model to characterize the statistical properties of the galaxy populations and real observations, respectively. By self-consistently inferring all model parameters, from high-level hyperparameters to SPS parameters of individual galaxies, one can separate sources of bias and uncertainty in the data. We demonstrate the strengths and flexibility of this approach by deriving accurate photometric redshifts for a sample of spectroscopically confirmed galaxies in the COSMOS field, all with 26-band photometry and spectroscopic redshifts. We achieve a performance competitive with publicly released photometric redshift catalogs based on the same data. Prior to this work, this approach was computationally intractable in practice due to the heavy computational load of SPS model calls; we overcome this challenge by the addition of neural emulators. We find that the largest photometric residuals are associated with poor calibration for emission-line luminosities and thus build a framework to mitigate these effects. This combination of physics-based modeling accelerated with machine learning paves the path toward meeting the stringent requirements on the accuracy of photometric redshift estimation imposed by upcoming cosmological surveys. The approach also has the potential to create new links between cosmology and galaxy evolution through the analysis of photometric data sets.

The role of impact parameter in typical close galaxy flybys

Abstract Close galaxy flybys, interactions during which two galaxies inter-penetrate, are frequent and can significantly affect the evolution of individual galaxies. Equal-mass flybys are extremely rare and almost exclusively distant, while frequent flybys have mass ratios $q=0.1$ or lower, with a secondary galaxy penetrating deep into the primary. This can result in comparable strengths of interaction between the two classes of flybys and lead to essentially the same effects. To demonstrate this, emphasise and explore the role of the impact parameter further, we performed a series of N-body simulations of typical flybys with varying relative impact parameters $b/R_{\mathrm{vir},1}$ ranging from $0.114$ to $0.272$ of the virial radius of the primary galaxy. Two-armed spirals form during flybys, with radii of origin correlated with the impact parameter and strengths well approximated with an inverted S-curve. The impact parameter does not affect the shape of induced spirals, and the lifetimes of a distinguished spiral structure appear to be constant, $T_\mathrm{LF} \sim 2$ Gyr. Bars, with strengths anti-correlated with the impact parameter, form after the encounter is over in simulations with $b/R_{\mathrm{vir},1} \leq 0.178$ and interaction strengths $S\geq0.076$ , but they are short-lived except for the stronger interactions with $S\geq0.129$ . We showcase an occurrence of multiple structures (ring-like, double bar) that survives for an exceptionally long time in one of the simulations. Effects on the pre-existing bar instability, that develops much later, are diverse: from an acceleration of bar formation, little to no effect, to even bar suppression. There is no uniform correlation between these effects and the impact parameter, as they are secondary effects, happening later in a post-flyby stage. Classical bulges are resilient to flyby interactions, while dark matter halos can significantly spin up in the amount anti-correlated with the impact parameter. There is an offset angle between the angular momentum vector of the dark matter halo and that of a disc, and it correlates linearly with the impact parameter. Thus, flybys remain an important pathway for structural evolution within galaxies in the local Universe.

Exploring the genesis of spiral galaxies

Context. A tantalizing enigma in extragalactic astronomy concerns the chronology and driving mechanisms of the build-up of late-type galaxies (LTGs). The standard scenario envisages two formation routes, with classical bulges (CBs) assembling first in a quick and violent quasi-monolithic episode followed by gradual disk assembly, and pseudo-bulges (PBs) progressively forming over gigayear-long timescales through gentle gas inflow from the disk and in situ star formation. The expectation from this antagonistic rationale is the segregation of present-day LTG bulges into two evolutionary distinct groups, which is in sharp contrast with recent observations. Aims. The present study aims for a thorough investigation of the star formation history (SFH) of LTGs with its ultimate goal being to outline a coherent framework for the formation and evolution of spiral galaxies and their main stellar components. Methods. Using population spectral synthesis models, we analyse the spatially resolved SFH of bulges and disks of 135 LTGs from the CALIFA survey, covering the relevant range in LTG mass. Complementarily, characteristic physical properties of bulges and disks, such as mean colours, mass- and light-weighted stellar age and metallicity, and EW(Hα), were contrasted with predictions from evolutionary synthesis models, by adopting exponentially declining SFHs with e-folding times τ between 0.1 and 20 Gyr. Results. Analysis of the SFH of roughly half a million spaxels consistently reveals that the main physical and evolutionary properties of both bulges and disks are continuously distributed across present-day total stellar mass ℳ⋆, T. The τ in spiral galaxies with log(ℳ⋆, T) &gt; 10 increases from the centre to the periphery, suggesting that these systems grow in an inside-out fashion. Quite importantly, the radial gradient of τ in an individual galaxy increases with increasing ℳ⋆, T, which is consistent with a high bulge-to-disk age contrast in high-mass spirals, while lower-mass LTGs display roughly the same τ throughout their entire radial extent, with intermediate mass galaxies in between. Predictions obtained through evolutionary synthesis are overall consistent with observed properties. Finally, bulges and disks of higher mass galaxies exhibit shorter formation timescales as compared to their lower mass counterparts. Conclusions. Collectively, the obtained results evince a coherent and unified picture for the formation and evolution of LTGs, in which PBs and CBs denote extremities of a continuous mass sequence. Our findings are consistent with the framework where bulges are assembled jointly with their parent disks by gradual inside-out growth, at a pace that is regulated by the depth of the galactic potential. This postulate is further supported by the fact that the revealed correlations are entirely devoid of a bimodality, as it would be expected if CBs and PBs were to emerge from two distinct formation routes.

Reconstructing the Assembly of Massive Galaxies. II. Galaxies Develop Massive and Dense Stellar Cores as They Evolve and Head toward Quiescence at Cosmic Noon

We use the spectral energy distribution fitting code Prospector to reconstruct the nonparametric star formation history (SFH) of massive () star-forming galaxies (SFGs) and quiescent galaxies (QGs) at redshift z obs ∼ 2 to investigate the joint evolution of star formation activity and structural properties. Compared to extended SFGs, compact SFGs are more likely to have experienced multiple star formation episodes, with the fractional mass formed during the older (≥1 Gyr) episode being larger, suggesting that high-redshift SFGs assembled their central regions earlier and then kept growing in central mass as they become more compact. The SFH of compact QGs does not significantly differ from the average for this category, and shows an early burst followed by a gradual decline of the star formation rate. The SFH of extended QGs, however, is similar to that of post-starburst galaxies, and their morphology is also frequently disturbed. Knowledge of the SFH also enables us to empirically reconstruct the structural evolution of individual galaxies. While the progenitor effect is clearly observed and accounted for in our analysis, it alone is insufficient to explain the observed structural evolution. We show that, as they evolve from star-forming phase to quiescence, galaxies grow massive dense stellar cores. Quenching begins at the center and then propagates outward to the rest of the structure. We discuss possible physical scenarios for the observed evolution and find that our empirical constraints are in good quantitative agreement with the model predictions from dissipative accretion of gas to the center followed by massive starbursts before final quiescence (wet compaction).

Detection of H i 21 cm emission from a strongly lensed galaxy at z ∼ 1.3

ABSTRACT We report the first 5σ detection of H i 21 cm emission from a star-forming galaxy at redshift z ∼ 1.3 (nearly 9 billion years ago) using upgraded Giant Metrewave Radio Telescope (uGMRT). This is the highest redshift H i detection in emission from an individual galaxy to date. The emission is strongly boosted by the gravitational lens, an early-type elliptical galaxy, at redshift z ∼ 0.13. The measured H i mass of the galaxy is $M_{\rm H\, \small {\rm I}} = (0.90 \pm 0.14 \pm 0.05) \times 10^{10}\, \mathrm{M}_{\odot }$, which is almost twice the inferred stellar mass of the galaxy, indicating an extended structure of the H i gas inside the galaxy. By fitting 2D Gaussian to the H i signal at the peak of the spectral line, we find the source to be marginally resolved with the position angle consistent with the emission being tangential to the critical curve of the lens mass distribution. This indicates that the solid angle of the approaching H i line flux comes very close to the inner lens caustic and results in very high magnification. These results, for the first time, demonstrate the feasibility of observing high-redshift H i in a lensed system with the modest amount of telescope time and open up exciting new possibilities for probing the cosmic evolution of neutral gas with existing and upcoming low-frequency radio telescopes in the near future.

Balancing the efficiency and stochasticity of star formation with dust extinction in z ≳ 10 galaxies observed by JWST

ABSTRACT Early observations with JWST indicate an overabundance of bright galaxies at redshifts z ≳ 10 relative to Hubble-calibrated model predictions. More puzzling still is the apparent lack of evolution in the abundance of such objects between z ∼ 9 and the highest redshifts yet probed, z ∼ 13–17. In this study, we first show that, despite a poor match with JWST luminosity functions (LFs), semi-empirical models calibrated to rest-ultraviolet LFs and colours at 4 ≲ z ≲ 8 are largely consistent with constraints on the properties of individual JWST galaxies, including their stellar masses, ages, and spectral slopes. We then show that order-of-magnitude scatter in the star formation rate of galaxies (at fixed halo mass) can indeed boost the abundance of bright galaxies, provided that star formation is more efficient than expected in low-mass haloes. However, this solution to the abundance problem introduces tension elsewhere: because it relies on the upscattering of low-mass haloes into bright magnitude bins, one expects typical ages, masses, and spectral slopes to be much lower than constraints from galaxies observed thus far. This tension can be alleviated by non-negligible reddening, suggesting that – if the first batch of photometrically selected candidates are confirmed – star formation and dust production could be more efficient than expected in galaxies at z ≳ 10.

The circumgalactic medium of Milky Way-like galaxies in the TNG50 simulation – I: halo gas properties and the role of SMBH feedback

ABSTRACT We analyse the physical properties of gas in the circumgalactic medium (CGM) of 132 Milky Way (MW)-like galaxies at z = 0 from the cosmological magneto-hydrodynamical simulation TNG50, part of the IllustrisTNG project. The properties and abundance of CGM gas across the sample are diverse, and the fractional budgets of different phases (cold, warm, and hot), as well as neutral H i mass and metal mass, vary considerably. Over our stellar mass range of $10^{10.5} &amp;lt; M_\star /{\rm M}_\odot &amp;lt; 10^{10.9}$, radial profiles of gas physical properties from $0.15 &amp;lt; R/R_{\rm 200c} &amp;lt; 1.0$ reveal great CGM structural complexity, with significant variations both at fixed distance around individual galaxies, and across different galaxies. CGM gas is multiphase: the distributions of density, temperature, and entropy are all multimodal, while metallicity and thermal pressure distributions are unimodal; all are broad. We present predictions for magnetic fields in MW-like haloes: a median field strength of $|B|\sim 1\,\mu{\rm G}$ in the inner halo decreases rapidly at larger distance, while magnetic pressure dominates over thermal pressure only within ${\sim}0.2 \times R_{\rm 200c}$. Virial temperature gas at ${\sim}10^6\,{\rm K}$ coexists with a subdominant cool, $\lt 10^5\,{\rm K}$, component in approximate pressure equilibrium. Finally, the physical properties of the CGM are tightly connected to the galactic star formation rate, in turn dependent on feedback from supermassive black holes (SMBHs). In TNG50, we find that energy from SMBH-driven kinetic winds generates high-velocity outflows (≳500–2000 km s−1), heats gas to supervirial temperatures (&amp;gt;106.5–7 K), and regulates the net balance of inflows versus outflows in otherwise quasi-static gaseous haloes.

NGC 3314a/b and NGC 3312: Ram pressure stripping in Hydra I cluster substructure

Cluster substructure and ram pressure stripping in individual galaxies are among the primary pieces of evidence for the ongoing growth of galaxy clusters as they accrete galaxies and groups from their surroundings. We present a multiwavelength study of the center of the Hydra I galaxy cluster, including exquisite new MeerKAT H I and DECam Hα imaging which reveal conclusive evidence for ram pressure stripping in NGC 3312, NGC 3314a, and NGC 3314b through compressed H I contours, well-defined H I tails, and ongoing star formation in the stripped gas. In particular, we quantify the stripped material in NGC 3312, and NGC 3314a, which makes up between 8% and 35% of the gas still in the disk, is forming stars at ∼0.5 M⊙ yr−1, and extends ∼30−60 kpc from the main disk. The estimated stellar mass in the tails is an order of magnitude less than the H I mass. A fourth “ring” galaxy at the same velocity does not show signs of ram pressure in H I. In addition, we used the H I and stellar morphologies, combined with a Beta model of the hot intracluster medium, to constrain the real distances of the galaxies to the cluster center, and we used the chance alignment of NGC 3314b behind NGC 3314a to break the degeneracy between whether the galaxies are in front or in back of the cluster. The drag seen in the H I tails supports our preferred scenario that NGC 3312 and NGC 3314a are moving toward us as part of a foreground substructure which has already passed its pericenter and is on “out fall” from the cluster. The high surviving H I content of the galaxies may suggest that the substructure or intragroup medium can protect them from the harshest effects of ram pressure, or that the galaxies are in fact on more tangential orbits.

Properties of Globular Clusters in Galaxy Clusters: Sensitivity from the Formation and Evolution of Globular Clusters

We investigate the properties of globular clusters (GCs) in a galaxy cluster, using the particle tagging method with a semianalytical approach in a cosmological context. We assume GCs form from dark matter halo mergers and their metallicity is assigned based on the stellar mass of the host dark matter halos and the formation redshift of GCs. Dynamical evolution and disruption of GCs are considered using semianalytical approaches, controlled by several free parameters. In this paper, we investigate how our results are changed by the choice of free parameters. We compare our fiducial results with representative observations, including the mass ratio between the GC system and its host galaxy, the GC occupancy, the number fraction of blue GCs, and the metallicity gradient with the GC mass. Because we can know the positions of GCs with time, comparison with additional observations is possible, e.g., the median radii of the GC system in individual galaxies, the mean projected density profiles of intracluster GCs, and the metallicity and age gradients of GCs with a clustercentric radius. We also find that the specific mass of the GC system in each galaxy is different with a clustercentric radius.

Conditional H i Mass Functions and the H i-to-halo Mass Relation in the Local Universe

We present a new H i mass estimator that relates (M HI/M *) to a linear combination of four galaxy properties: stellar surface mass density, color index u − r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H i mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H i detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H i masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H i mass function (HIMF) of local galaxies, as well as the conditional H i mass function in galaxy groups and the H i–halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 109 M ⊙, but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H i mass in a group increases monotonically with halo mass, but for central galaxies, the H i mass shows weak dependence on halo mass when M h ≳ 1012 M ☉. The observed H i–halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

The resolved scaling relations in DustPedia: Zooming in on the local Universe

Aims.We perform a homogeneous analysis of an unprecedented set of spatially resolved scaling relations (SRs) between interstellar medium (ISM) components, that is to say dust, gas, and gas-phase metallicity, and other galaxy properties, such as stellar mass (Mstar), total baryonic content, and star-formation rate (SFR), in a range of physical scales between 0.3 and 3.4 kpc. We also study some ratios between galaxy components: dust-to-stellar, dust-to-gas, and dust-to-metal ratios.Methods.We use a sample of 18 large, spiral, face-on DustPedia galaxies. The sample consists of galaxies with spatially resolved dust maps corresponding to 15Herschel-SPIRE 500 μm resolution elements across the optical radius, with the morphological stage spanning fromT = 2 to 8,Mstarfrom 2 × 109to 1 × 1011M⊙, SFR from 0.2 to 13M⊙yr−1, and oxygen abundance from 12 + log(O/H) = 8.3 to 8.8.Results.All the SRs are moderate or strong correlations except the dust-H ISR that does not exist or is weak for most galaxies. The SRs do not have a universal form but each galaxy is characterized by distinct correlations, affected by local processes and galaxy peculiarities. The SRs hold, on average, starting from the scale of 0.3 kpc, and if a breaking down scale exists it is below 0.3 kpc. By evaluating all galaxies together at the common scale of 3.4 kpc, differences due to peculiarities of individual galaxies are cancelled out and the corresponding SRs are consistent with those of whole galaxies. By comparing subgalactic and global scales, the most striking result emerges from the SRs involving ISM components: the dust-total gas SR is a good correlation at all scales, while the dust-H2and dust-H ISRs are good correlations at subkiloparsec/kiloparsec and total scales, respectively. For the other explored SRs, there is a good agreement between small and global scales and this may support the picture where the main physical processes regulating the properties and evolution of galaxies occur locally. In this scenario, our results are consistent with the hypothesis of self-regulation of the star-formation process. The analysis of subgalactic ratios between galaxy components shows that they are consistent with those derived for whole galaxies, from low to high redshift, supporting the idea that also these ratios could be set by local processes.Conclusions.Our results highlight the heterogeneity of galaxy properties and the importance of resolved studies on local galaxies in the context of galaxy evolution. They also provide fundamental observational constraints to theoretical models and updated references for high-redshift studies.

Spitzer/IRS Full Spectral Modeling to Characterize Mineralogical Properties of Silicate Dust in Heavily Obscured AGNs

Mid-infrared silicate dust bands observed in heavily obscured active galactic nuclei (AGNs) include information on the mineralogical properties of silicate dust. We aim to investigate the mineralogical picture of the circumnuclear region of heavily obscured AGNs to reveal obscured AGN activities through the picture. In our previous study, we investigated the properties of silicate dust in heavily obscured AGNs, focusing on the mineralogical composition and the crystallinity with Spitzer/IRS 5.3–12 μm spectra. In this study, we model the full-range Spitzer/IRS 5–30 μm spectra of 98 heavily obscured AGNs using a one-dimensional radiative transfer calculation with four dust species in order to evaluate wider ranges of the properties of silicate dust more reliably. Comparing fitting results between four dust models with different sizes and porosities, 95 out of the 98 galaxies prefer a porous silicate dust model without micron-sized large grains. The pyroxene mass fraction and the crystallinity are overall consistent with—but significantly different from—the previous results for the individual galaxies. The pyroxene-poor composition, small dust size, and high porosity are similar to newly formed dust around mass-loss stars as seen in our Galaxy, which presumably originates from the recent circumnuclear starburst activity. The high crystallinity on average suggests dust processing induced by AGN activities.

Emission Line Galaxies in the SHARDS Hubble Frontier Fields. II. Limits on Lyman-continuum Escape Fractions of Lensed Emission Line Galaxies at Redshifts 2 < z < 3.5

We present an investigation of escape fractions of UV photons from a unique sample of lensed low-mass emission line–selected galaxies at z < 3.5 found in the SHARDS Hubble Frontier Fields medium-band survey. We have used this deep imaging survey to locate 42 relatively low-mass galaxies down to log (M */M ⊙) = 7 in the redshift range 2.4 < z < 3.5 that are candidate line emitters. Using deep multiband Hubble UVIS imaging, we investigate the flux of escaping ionizing photons from these systems, obtaining 1σ upper limits of f rel esc ∼ 7% for individual galaxies and <2% for stacked data. We measure potential escaping Lyman-continuum flux for two low-mass line emitters with values at and , both detected at the ∼3.2σ level. A detailed analysis of possible contamination reveals a <0.1% probability that these detections result from line-of-sight contamination. The relatively low Lyman-continuum escape fraction limit and the low fraction of systems detected are an indication that low-mass line-emitting galaxies may not be as important a source of reionization as hoped if these are analogs of reionization sources. We also investigate the structures of our galaxy sample, finding no evidence for a correlation of escape fraction with asymmetric structure.

JWST/NIRCam observations of stars and H ii regions in z ≃ 6–8 galaxies: properties of star-forming complexes on 150 pc scales

ABSTRACT The onset of the JWST-era provides a much-improved opportunity to characterize the resolved structure of early star-forming systems. Previous Spitzer observations of z ≳ 6 galaxies revealed the presence of old stars and luminous H ii regions (via [O iii]+H β emission), but the poor resolution stunted our ability to map their locations with respect to the star-forming regions identified in the rest-UV. In this paper, we investigate the internal structure of 12 of the most luminous z ≃ 6–8 galaxies in the EGS field observed with recent JWST/NIRCam imaging. The systems appear clumpy in the rest-UV, with more than half of the light coming from ≃ 107 to 109 M⊙ star-forming complexes that are ≃150–480 pc in size. The clumps tend to be dominated by young stars (median = 36 Myr), but we also find large variations in clump ages within individual galaxies. The [O iii]+H β EW varies significantly across individual galaxies (reflecting differences in stellar and gas properties), but the H ii regions largely track the UV-bright complexes. Perhaps surprisingly, the rest-optical continuum is just as clumpy as the UV, and we do not find older (and redder) nuclear stellar components that were previously undetected or faint in the UV. The majority of the stellar mass in bright 6 &amp;lt; z &amp;lt; 8 galaxies appears to be contained in the ≳150 pc-scale clumpy star-forming complexes, reflecting the very active phase of assembly that is common in reionization-era galaxies.

The nature of 500 micron risers – II. Multiplicities and environments of sub-mm faint dusty star-forming galaxies

ABSTRACT We present Submillimeter Common User Bolometer Array 2 (SCUBA-2) and Submillimeter Array (SMA) follow-up observations of four candidate high-redshift dusty star-forming galaxies, selected as sources with rising spectral energy distributions in the 250, 350, and 500 μm Herschel SPIRE bands. Previous SMA observations showed no counterparts to these sources, but in our deeper sub-mm observations we detect counterparts to all four 500 μm risers, with three resolving into multiple systems. For these three multiple systems, the SMA 345 GHz (≈870 μm) observations recover $123 \pm 73{{\ \rm per\ cent}}$, $60 \pm 15{{\ \rm per\ cent}}$, and $19 \pm 4{{\ \rm per\ cent}}$ of the integrated 850 μm flux density from SCUBA-2, indicating that there may be additional sources below our SMA detection limit making up a dense, protocluster core. The fourth 500 μm riser was observed at a lower frequency and so we cannot make a similar comparison. We estimate photometric redshifts based on FIR/sub-mm colours, finding that three fourths likely lie at z ≥ 2. This fits with the interpretation that the 500 μm riser selection criterion selects both intrinsically red, individual galaxies at z &amp;gt; 4, and multiple systems at more moderate redshifts, artificially reddened by the effects of blending. We use the SCUBA-2 850 μm maps to investigate the environments of these 500 μm risers. By constructing cumulative number counts and estimating photometric redshifts for surrounding SCUBA-2 detections, we find that one of our 500 μm risers could plausibly reside in a z ≥ 2 protocluster. We infer that bright 500 μm risers with faint 850 μm flux densities are typically multiple systems at z ≥ 2 that may reside in overdensities of bright sub-mm galaxies.