Search

Context. The UVJ color-color diagram is a widely used diagnostic to separate star-forming and quiescent galaxies. Observational data from photometric surveys reveal a strong stellar mass trend, with higher-mass star-forming galaxies being systematically more dust-reddened. Aims. We analyze the UVJ diagram in the TNG100 cosmological simulation at cosmic noon ($z\approx2$). Specifically, we focus on the trend between UVJ colors and mass which has so far not been reproduced in any cosmological simulation. Methods. We apply the SKIRT dust radiative transfer code to the TNG100 simulation to generate rest-frame UVJ fluxes. These UVJ colors are then compared to observational data from several well-studied extragalactic fields from the CANDELS/3D-HST programs, augmented by recent JWST/NIRCam photometry. Results. Quiescent and low-mass ($M_\star\lesssim10^{10.5}\,\mathrm{M}_\odot$) galaxies at cosmic noon do not require significant levels of dust reddening, as opposed to massive ($M_\star\gtrsim10^{11}\,\mathrm{M}_\odot$) star-forming galaxies. An extensive range of possible dust models fall short of the required dust reddening in V-J color for massive star-forming galaxies, with the simulated galaxies being too blue by $\approx0.9\,\mathrm{mag}$. Conclusions. We find that only variations in the star-to-dust geometries of the simulated galaxies can yield V-J colors that are red enough to match the observations. A toy model with isolated dust screens around younger stellar populations (with ages below $\sim1\,\mathrm{Gyr}$) can reproduce the observational data, while all conventional dust radiative transfer models (where the dust distribution follows the metals in the interstellar medium) fail to achieve the required V-J colors., Comment: Main text 17 pages, 12 figures. Resubmitted to A&A after first referee report. Comments warmly welcome!

The James Webb Space Telescope (JWST) is unveiling the rest-frame near-IR structure of galaxies. We measure the evolution with redshift of the rest-frame optical and near-IR S\'ersic index ($n$), and examine the dependence on stellar mass and star-formation activity across the redshift range $0.5\leq z\leq2.5$. We infer rest-frame near-IR S\'ersic profiles for $\approx 15.000$ galaxies in publicly available NIRCam imaging mosaics from the COSMOS-Web and PRIMER surveys. We augment these with rest-frame optical S\'ersic indices, previously measured from HST imaging mosaics. The median S\'ersic index evolves slowly or not at all with redshift, except for very high-mass galaxies ($M_\star > 10^{11}~{\text{M}}_\odot$), which show an increase from $n\approx 2.5$ to $n\approx 4$ at $z<1$. High-mass galaxies have higher $n$ than lower-mass galaxies ($M_\star=10^{9.5}~{\text{M}}_\odot$) at all redshifts, with a stronger dependence in the rest-frame near-IR than in the rest-frame optical at $z>1$. This wavelength dependence is caused by star-forming galaxies that have lower optical than near-IR $n$ at z>1 (but not at z<1). Both at optical and near-IR wavelengths, star-forming galaxies have lower $n$ than quiescent galaxies, fortifying the connection between star-formation activity and radial stellar mass distribution. At $z>1$ the median near-IR $n$ varies strongly with star formation activity, but not with stellar mass. The scatter in near-IR $n$ is higher in the green valley (0.25 dex) than on the star-forming sequence and among quiescent galaxies (0.18 dex) -- this trend is not seen in the optical because dust and young stars contribute to the variety in optical light profiles. Our newly measured rest-frame near-IR radial light profiles motivate future comparisons with radial stellar mass profiles of simulated galaxies as a stringent constraint on processes that govern galaxy formation., Comment: Accepted for publication on A&A

Current galaxy formation simulations often approximate star-formation, necessitating models of star-forming regions to produce observables. In the first paper of the series, we introduced TODDLERS, a time-resolved model of UV-mm emission from star-forming regions implemented in the radiative transfer code SKIRT. This work uses SKIRT-TODDLERS to produce synthetic observations, demonstrating its potential through observables related to star-formation and comparing results with existing models in SKIRT. We calculate broadband and line emission maps for 30 Milky Way-like galaxies from the Auriga simulation at z=0. Analyzing FUV and IR data, we calculate kpc-resolved IR correction factors (k_IR), quantifying the ratio of FUV luminosity absorbed by dust to reprocessed IR luminosity. We use IR maps to calculate kpc-scale MIR (8 micron / 24 micron) and FIR (70 micron / 500 micron) colors. H alpha and H beta line maps are used to study the Balmer decrement and dust correction. We also verify the fidelity of our model FIR fine structure lines as SFR indicators. We find that the Auriga integrated UV-mm SEDs show higher FUV/NUV attenuation and lower 24 micron emission when using TODDLERS instead of the existing models in SKIRT, alleviating tensions with observations. The light-weighted mean k_IR increases with aperture and inclination, while its correlation with kpc-resolved sSFR is weaker than literature values. Kpc-scale MIR-FIR colors agree excellently with local observations, with anti-correlation degree varying by galaxy morphology. The Balmer decrement effectively corrects for dust, with attenuation law varying with dust amount. H alpha emission attenuation levels are comparable to high-density regions in state-of-the-art simulations. FIR fine-structure line emission-based luminosity-SFR relations align with global observations, [CII] showing best agreement., Comment: 22 pages, 20 figures. Accepted for publication in Astronomy and Astrophysics on October 21, 2024

Neutral sodium (Na I) is an alkali metal with a favorable absorption cross section such that tenuous gases are easily illuminated at select transiting exoplanet systems. We examine both the time-averaged and time-series alkali spectral flux individually, over 4 nights at a hot Saturn system on a $\sim$ 2.8 day orbit about a Sun-like star WASP-49 A. Very Large Telescope/ESPRESSO observations are analyzed, providing new constraints. We recover the previously confirmed residual sodium flux uniquely when averaged, whereas night-to-night Na I varies by more than an order of magnitude. On HARPS/3.6-m Epoch II, we report a Doppler redshift at $v_{ \Gamma, \mathrm{NaD}} =$ +9.7 $\pm$ 1.6 km/s with respect to the planet's rest frame. Upon examining the lightcurves, we confirm night-to-night variability, on the order of $\sim$ 1-4 % in NaD rarely coinciding with exoplanet transit, not readily explained by stellar activity, starspots, tellurics, or the interstellar medium. Coincident with the $\sim$+10 km/s Doppler redshift, we detect a transient sodium absorption event dF$_{\mathrm{NaD}}$/F$_{\star}$ = 3.6 $\pm$ 1 % at a relative difference of $\Delta F_{\mathrm{NaD}} (t) \sim$ 4.4 $\pm$ 1 %, enduring $\Delta t_{\mathrm{NaD}} \gtrsim$ 40 minutes. Since exoplanetary alkali signatures are blueshifted due to the natural vector of radiation pressure, estimated here at roughly $\sim$ -5.7 km/s, the radial velocity is rather at +15.4 km/s, far larger than any known exoplanet system. Given that the redshift magnitude v$_{\Gamma}$ is in between the Roche limit and dynamically stable satellite orbits, the transient sodium may be a putative indication of a natural satellite orbiting WASP-49 A b., Comment: Accepted to ApJ Letters (2024 August 2)

We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Tr\v{c}ka et al. (2022) of the higher-resolution TNG50 simulation. We verify that TNG100 reproduces observational luminosity functions at low redshifts to excellent precision, unlike TNG50. Additionally, we test the realism of our TNG100 plus SKIRT fluxes by comparing various flux and color relations to data from the GAMA survey. TNG100 broadly reproduces the observed distributions, but we predict ultraviolet colors that are too blue by $\approx\,$0.4 mag, possibly related to the extinction in the star-forming regions subgrid model not being selective enough. Furthermore, we find that the simulated galaxies exhibit mid-infrared fluxes elevated by up to $\approx\,$0.5 mag that we attribute to overly effective stochastic heating of the diffuse dust. All synthetic broadband fluxes and SEDs are made publicly available in three orientations and four apertures, and can readily be used to study TNG100 galaxies in a mock observational fashion., Comment: Main text 14 pages, 5 figures, accepted to MNRAS. Generated data products are publicly available at www.tng-project.org/gebek24

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

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

We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feedback processes including stellar winds, supernovae, and radiation pressure, as well as the gravitational forces on the gas. The semi-analytical evolution model is coupled with the photoionization code Cloudy to calculate time-dependent UV-mm spectral energy distributions (SEDs) from star-forming regions of varying metallicity, star-formation efficiency, birth-cloud density, and mass. The calculated SEDs include the stellar, nebular, and dust continuum emission along with a wide range of emission lines originating from H ii, photo-dissociation, and molecular gas regimes tabulated at high resolution. The SEDs incorporated in SKIRT are generated by calculating a stellar-mass normalized luminosity, which assumes that each emission source is composed of a power-law population of star-forming clouds. When compared to the previous treatment of star-forming regions in SKIRT, TODDLERS shows a better agreement with low-redshift observational data in the IR wavelength range while offering a more comprehensive line-emission support. This paves the way for a variety of applications using simulated galaxies at low and high redshift., Comment: Accepted for publication in MNRAS on 27/09/23

Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to model the atomic-to-molecular transition in TNG50, the highest-resolution run of the IllustrisTNG suite. Comparing integrated quantities such as the HI mass function, we study to what detail the UV field needs to be modelled in order to calculate realistic cold gas properties. We then evaluate new, spatially resolved comparisons for cold gas in galaxies by exploring synthetic maps of atomic hydrogen at redshift zero and compare them to 21-cm observations of local galaxies from the WHISP survey. In terms of non-parametric morphologies, we find that TNG50 HI maps are less concentrated than their WHISP counterparts (median $\Delta C\approx0.3$), due in part to central HI deficits related to the ejective character of supermassive black hole feedback in TNG. In terms of the HI column density distribution function, we find discrepancies between WHISP and IllustrisTNG that depend on the total HI abundance in these datasets as well as the postprocessing method. To fully exploit the synergy between cosmological simulations and upcoming deep HI/H2 data, we advocate the use of accurate methods to estimate the UV radiation field and to generate mock maps., Comment: 21 pages, 10 figures (main text). Accepted to MNRAS

The ratio of $\alpha$-elements to iron in galaxies holds valuable information about the star formation history since their enrichment occurs on different timescales. The fossil record of stars in galaxies has mostly been excavated for passive galaxies, since the light of star-forming galaxies is dominated by young stars which have much weaker atmospheric absorption features. Here we use the largest reference cosmological simulation of the EAGLE project to investigate the origin of variations in stellar $\alpha$-enhancement among star-forming galaxies at $z=0$, and their impact on integrated spectra. The definition of $\alpha$-enhancement in a composite stellar population is ambiguous. We elucidate two definitions - termed `mean' and `galactic' $\alpha$-enhancement - in more detail. While a star-forming galaxy has a high `mean' $\alpha$-enhancement when its stars formed rapidly, a galaxy with a large `galactic' $\alpha$-enhancement generally had a delayed star formation history. We find that absorption-line strengths of Mg and Fe correlate with variations in $\alpha$-enhancement. These correlations are strongest for the `galactic' $\alpha$-enhancement. However, we show that these are mostly caused by other effects which are cross-correlated with $\alpha$-enhancement, such as variations in the light-weighted age. This severely complicates the retrieval of $\alpha$-enhancements in star-forming galaxies. The ambiguity is not severe for passive galaxies and we confirm that spectral variations in these galaxies are caused by measurable variations in $\alpha$-enhancements. We suggest that this more complex coupling between $\alpha$-enhancement and star formation histories can guide the interpretation of new observations of star-forming galaxies., Comment: 21 pages, 12 figures, 4 tables. Accepted to ApJ

Hydrostatic equilibrium is an excellent approximation for the dense layers of planetary atmospheres where it has been canonically used to interpret transmission spectra of exoplanets. Here we exploit the ability of high-resolution spectrographs to probe tenuous layers of sodium and potassium gas due to their formidable absorption cross-sections. We present an atmosphere-exosphere degeneracy between optically thick and optically thin mediums, raising the question of whether hydrostatic equilibrium is appropriate for Na I lines observed at exoplanets. To this end we simulate three non-hydrostatic, evaporative, density profiles: (i) escaping, (ii) exomoon, and (iii) torus to examine their imprint on an alkaline exosphere in transmission. By analyzing an evaporative curve of growth we find that equivalent widths of $W_{\mathrm{Na D2}} \sim 1- 10$ mA are naturally driven by evaporation rates $\sim 10^3 - 10^5$ kg/s of pure atomic Na. To break the degeneracy between atmospheric and exospheric absorption, we suggest that if the line ratio is $\mathrm{D2/D1} \gtrsim 1.2$ the gas is optically thin on average and roughly indicating a non-hydrostatic structure of the atmosphere/exosphere. We show this is the case for Na I observations at hot Jupiters WASP-49b and HD189733b and also simulate their K I spectra. Lastly, motivated by the slew of metal detections at ultra-hot Jupiters, we suggest a toroidal atmosphere at WASP-76b and WASP-121b is consistent with the Na I data at present., Comment: 23 pages, 21 figures, accepted by MNRAS

Extrasolar satellites are generally too small to be detected by nominal searches. By analogy to the most active body in the Solar System, Io, we describe how sodium (Na I) and potassium (K I) $\textit{gas}$ could be a signature of the geological activity venting from an otherwise hidden exo-Io. Analyzing $\sim$ a dozen close-in gas giants hosting robust alkaline detections, we show that an Io-sized satellite can be stable against orbital decay below a planetary tidal $\mathcal{Q}_p \lesssim 10^{11}$. This tidal energy is focused into the satellite driving a $\sim 10^{5 \pm 2}$ higher mass loss rate than Io's supply to Jupiter's Na exosphere, based on simple atmospheric loss estimates. The remarkable consequence is that several exo-Io column densities are on average $\textit{more than sufficient}$ to provide the $\sim$ 10$^{10 \pm 1}$ Na cm$^{-2}$ required by the equivalent width of exoplanet transmission spectra. Furthermore, the benchmark observations of both Jupiter's extended ($\sim 1000$ R$_J$) Na exosphere and Jupiter's atmosphere in transmission spectroscopy yield similar Na column densities that are purely exogenic in nature. As a proof of concept, we fit the "high-altitude" Na at WASP 49-b with an ionization-limited cloud similar to the observed Na profile about Io. Moving forward, we strongly encourage time-dependent ingress and egress monitoring along with spectroscopic searches for other volcanic volatiles., Comment: 23 pages, 5 figures, accepted to ApJ