CyberSec Research

Aims: This paper investigates the coevolution of metals and dust for 173 galaxies at $4.0<z<11.4$ observed with JWST/NIRSpec in the CEERS project. We focus on galaxies with extremely low dust attenuation to understand the physical mechanisms at play. Methods: We developed a new version of the \texttt{CIGALE} code that integrates spectroscopic and photometric data. By statistically comparing observations with modeled spectra, we derive physical parameters to constrain these mechanisms. Results: Our analysis reveals a population of 49 extremely low dust attenuation galaxies (GELDAs), consistent with $A_{FUV} = 0.0$ within $2\sigma$ and $M_{\star} < 10^9 M_\odot$. The stacked spectrum of GELDAs shows a very blue UV slope $\beta_{FUV} = -2.451 \pm 0.066$ and a Balmer decrement H$\alpha$/H$\beta = 2.932 \pm 0.660$, consistent with no dust and Case B recombination with minimal underlying absorption. Notably, GELDAs are more prevalent at $z > 8.8$ (83.3\%) than at lower redshifts (26.3\%), suggesting they could dominate in the early Universe. Using a far-infrared dust spectrum from the ALPINE sample, we study $M_{dust}$ vs. $M_{\star}$ trends. These exhibit upper and lower sequences connected by transitional galaxies. Our comparison with models indicates a critical transition around $M_{\star} \approx 10^{8.5}\,M_\odot$, from dust dominated by stellar sources (SNe and AGB stars) to dust growth via gas accretion. This corresponds to a metallicity of $12 + \log_{10}(O/H) = 7.60$ ($Z/Z_\odot \approx 0.1$), aligning with the point where ISM dust growth matches stellar dust production. The sample has a high gas fraction ($f_{\mathrm{gas}} \gtrsim 0.9$), with no significant gas expulsion, and high surface gas densities. This leads to low star formation efficiencies compared to sub-millimeter galaxies. GELDAs may help explain the observed excess of bright galaxies at $z \gtrsim 9$.

Adding to the RESOLVE and ECO Gas in Galaxy Groups (G3) initiative, we examine possible drivers of group-integrated HI-to-halo mass ratios ($M_{\rm HI,grp}/M_{\rm halo}$) and group X-ray emission, including group halo mass ($M_{\rm halo}$), virialization as probed by crossing time ($t_{\rm cross}$), presence of active galactic nuclei (AGN), and group-integrated fractional stellar mass growth rate (FSMGR$_{\rm grp}$). G3 groups span $M_{\rm halo}=10^{11-14.5}\,M_\odot$ with comprehensive HI and AGN information, which we combine with X-ray stacking of ROSAT All-Sky data. We detect hot gas emission exceeding AGN and X-ray binary backgrounds confidently for $M_{\rm halo}=10^{12.6-14}\,M_\odot$ and unambiguously for $M_{\rm halo}>10^{14}\,M_\odot$, reflecting an inverse dependence of $M_{\rm\,HI,grp}/M_{\rm halo}$ and hot gas emission on halo mass. At fixed halo mass, $M_{\rm\,HI,grp}/M_{\rm halo}$ transitions to greater spread below $t_{\rm cross}\sim2$ Gyr. Dividing groups across this transition, lower-$t_{\rm cross}$ groups show elevated X-ray emission compared to higher-$t_{\rm cross}$ groups for $M_{\rm halo}>10^{13.3}\,M_\odot$, but this trend reverses for $M_{\rm halo}=10^{12.6-13.3}\,M_\odot$. Additionally, AGN-hosting halos below $M_{\rm halo}\sim10^{12.1}\,M_\odot$ exhibit a broad, $\sim$0.25 dex deep valley in $M_{\rm HI,grp}/M_{\rm halo}$ compared to non-AGN-hosting halos with correspondingly reduced FSMGR$_{\rm grp}$. When diluted by non-AGN-hosting halos, this valley becomes shallower and narrower, falling roughly between $M_{\rm halo}=10^{11.5}\,M_\odot$ and $M_{\rm halo}=10^{12.1}\,M_\odot$ in the overall $M_{\rm\,HI,grp}/M_{\rm\,halo}$ vs. $M_{\rm halo}$ relation. We may also detect a second, less easily interpreted valley at $M_{\rm halo}\sim10^{13}\,M_\odot$. Neither valley matches theoretical predictions of a deeper valley at or above $M_{\rm halo}=10^{12.1}\,M_\odot$.

Euclid will image ~14000 deg^2 of the extragalactic sky at visible and NIR wavelengths, providing a dataset of unprecedented size and richness that will facilitate a multitude of studies into the evolution of galaxies. In the vast majority of cases the main source of information will come from broad-band images and data products thereof. Therefore, there is a pressing need to identify or develop scalable yet reliable methodologies to estimate the redshift and physical properties of galaxies using broad-band photometry from Euclid, optionally including ground-based optical photometry also. To address this need, we present a novel method to estimate the redshift, stellar mass, star-formation rate, specific star-formation rate, E(B-V), and age of galaxies, using mock Euclid and ground-based photometry. The main novelty of our property-estimation pipeline is its use of the CatBoost implementation of gradient-boosted regression-trees, together with chained regression and an intelligent, automatic optimization of the training data. The pipeline also includes a computationally-efficient method to estimate prediction uncertainties, and, in the absence of ground-truth labels, provides accurate predictions for metrics of model performance up to z~2. We apply our pipeline to several datasets consisting of mock Euclid broad-band photometry and mock ground-based ugriz photometry, to evaluate the performance of our methodology for estimating the redshift and physical properties of galaxies detected in the Euclid Wide Survey. The quality of our photometric redshift and physical property estimates are highly competitive overall, validating our modeling approach. We find that the inclusion of ground-based optical photometry significantly improves the quality of the property estimation, highlighting the importance of combining Euclid data with ancillary ground-based optical data. (Abridged)

Self-interacting dark matter (SIDM) theories predict that dark matter halos experience core-collapse in late-stage evolution, a process where the halo's inner region rapidly increases in density and decreases in size. This process can be modeled by treating the dark matter as a gravothermal fluid, and solving the fluid equations to predict the density profile evolution. This model is incomplete without calibration to N-body simulations, through a constant factor $\beta$ included in the thermal conductivity for the long-mean-free-path limit. The value of $\beta$ employed in the gravothermal fluid formalism has varied between studies, with no clear universal value in the literature. In this work, we use the N-body code Arepo to conduct a series of isolated core-collapse simulations across a range of scattering cross-sections, halo concentrations, and halo masses to calibrate the heat transfer parameter $\beta$. We find that $\beta$ is independent of cross-section, halo concentration, and halo mass for velocity independent elastic scattering cross-sections. We present a model for an effective $\beta$ as a function of a dimensionless cross-section, to describe halo evolution in the long mean free path limit, and show that it accurately captures halo evolution as long as the cross section is not too large. This effective model facilitates comparisons between simulations and the gravothermal model, and enables fast predictions of the dark matter density profile at any given time without running N-body simulations.

The study of molecular line emission is crucial to unveil the kinematics and the physical conditions of gas in star-forming regions. Our aim is to quantify the reliability of using individual molecular transitions to derive physical properties of the bulk of the H2 gas, looking at morphological correlations in their overall integrated molecular line emission with the cold dust. For this study we selected transitions of H2CO, CH$_3$OH, DCN, HC$_3$N, CH$_3$CN, CH$_3$OCHO, SO, and SiO and compared them with the 1.38 mm dust continuum emission at different spatial scales in the ALMAGAL sample, that observed a total of 1013 targets covering all evolutionary stages of the high-mass star-formation process and different conditions of clump fragmentation. We used the method of the histogram of oriented gradients (HOG) implemented in the tool astroHOG to compare the morphology of integrated line emission with maps of the 1.38 mm dust continuum emission. Moreover, we calculated the Spearman's correlation coefficient, and compared it with our astroHOG results. Only H$_2$CO, CH$_3$OH, and SO show emission on spatial scales comparable with the diffuse continuum emission. However, from the HOG method, the median correlation of the emission of each of these species with the continuum is only $\sim$24-29%. In comparison with the dense fragments these molecular species still have low values of correlation. On the other hand DCN, HC$_3$N, CH$_3$CN, and CH$_3$OCHO show a good correlation with the dense dust fragments, above 60%. The worst correlation is seen with SiO, both with the extended continuum emission and with compact sources. From the comparison of the results of the HOG method and the Spearman's correlation coefficient, the HOG method gives much more reliable results than the intensity-based coefficient in estimating the level of similarity of the emission morphology.

We study the AGB population of the galaxy M31, based on available HST and Spitzer data, to characterize the individual sources in terms of mass, metallicity and formation epoch of the progenitors. Particular attention is dedicated to the derivation of the dust production rate of the stars, in the attempt of determining the global current dust production rate of the galaxy, divided between the silicates and the carbonaceous dust contributions. We use results from stellar evolution modelling complemented by the description of the dust formation process in the wind, to be used in a population synthesis approach, based on the star formation history and age-metallicity relationship obtained in previous investigations. The comparison between the results from synthetic modelling and the data available are used for the characterization of AGB stars in M31. We find that the bulk of the AGB population of M31 is composed by low-mass stars of different metallicity formed between 6 Gyr and 14 Gyr ago, with an additional, significant contribution from the progeny of 1.7-2.5Msun stars formed during the secondary peak in the star formation, which occurred between 1 and 2 Gyr ago. The dust production rate of the galaxy is mostly provided by carbon stars, whose contribution is of the order of 4x10^{-4} Msun/yr, completed by silicates production from massive AGB stars, occurring at a rate of 6x10^{-5} Msun/yr. The implications of the present results on the reliability of AGB modelling are also commented.

Ram pressure stripping (RPS) has a non-negligible impact on the gas content of cluster galaxies. We use the semi-analytic model GAEA and the hydro-simulation TNG to investigate whether cluster galaxies suffer a strong RPS that is sufficient to remove a significant fraction of their gas during the first pericentric passage. We estimate that a ram pressure of $10^{-10.5}$, $10^{-12} $, $10^{-13.5} $ $g cm^{-1} s^{-2}$ can remove at most $90\%$, $50\%$, and $20\%$ of the cold gas reservoir from low-mass galaxies with $9<\log M_{\star}/{\rm M}_{\odot} <9.5$, assuming the gas can be stripped instantaneously. We then use this information to divide the phase space diagram into `strong', `moderate', `weak', and `no' RPS zones. By tracing the orbit of galaxies since $2.5R_{vir}$, we find in both GAEA and TNG that about half of the galaxies in Virgo-like halos ($\log M_h / M_{\odot} \sim 14 $) did not suffer strong RPS during the first pericentric passage. In Coma-like halos ($\log M_h / M_{\odot} \sim 15$), almost all galaxies have suffered strong RPS during the first pericentric passage, which can remove all gas from low-mass galaxies but is insufficient to significantly reduce the gas content of more massive galaxies. In general, results from TNG and GAEA are consistent, with the RPS being only slightly stronger in TNG than in GAEA. Our findings suggest that most cluster galaxies will maintain a notable fraction of their gas and continue forming stars after the first pericentric passage, except for those with low stellar mass ($\log M_{\star}/{\rm M}_{\odot} <9.5$) in very massive halos ($\log M_{h}/{\rm M}_{\odot} > 15$).

Quasi-periodic eruptions (QPEs), the repeated outbursts observed in soft X-ray bands, have attracted broad interest, but their physical origin is under debate. One of the popular models, the star-disk collision model, suggests that QPEs can be produced through periodic collisions of an orbiting star with the accretion disk of a central black hole (BH). However, previous tests of the star-disk collision model mainly focus on the timing analysis. Other observed properties, such as peak luminosities $L_{\rm{p}}$, durations $t_{\rm{e}}$, and radiation temperatures $T_{\rm{p}}$ of the eruptions, are not systematically investigated. For a sample of six QPE sources and two QPE-like sources, we test the star-disk collision model by using these observables to derive the constraints on the stellar radius $R_*$. We find that, except for two sources (eRo-QPE3 and eRo-QPE4), the rest of the sample either has no allowed $R_*$ to simultaneously reproduce the observed $L_{\rm{p}}$ and $t_{\rm{e}}$, or the required $R_*$ is too large to avoid being disrupted by the central BH. For the two exceptions, a stellar radius of the order of $1\ R_{\rm{\odot}}$ is necessary to satisfy all the constraints. Another issue with the simplest version of this model is that it predicts $k T_{\rm{p}} \sim 10\ \rm{eV}$, one order of magnitude lower than the observed value.

A recent study using weak gravitational lensing reveals that there are some isolated galaxies having almost flat rotation curves at very large distance from the galactic centres. According to the authors of the study this provides a strong challenge the standard cold dark matter model, since the dark haloes are too small to explain their observations, especially for small stellar masses. In this article, we show that improving their model, the virial radius is larger than their estimates. The NFW rotational curve, and especially the pseudo Isothermal one, are in agreement with their flat rotational curves, especially for the larger baryonic mass bins used by the authors.

Classical Cepheids can be used as age indicators due to well-established period-age and period-age-color relations. \citet{Desomma2021} refined these relations by including a metallicity term and different Mass-Luminosity assumptions. In this study, we apply the period-age-metallicity relation for the first time to samples of Classical Cepheids in M31 and M33. For both galaxies, we consider Cepheid coordinates and spatial distributions, along with the metallicity gradients by \citet{Zaritsky1994} and \citet{Magrini2007}, to provide a metallicity estimate for each pulsator. By applying the period-age-metallicity relation, we derive individual ages for each Cepheid. Combining these ages and spatial distributions, we construct detailed age maps for both galaxies. Our analysis confirms a radial age gradient, with younger Cepheids preferentially found toward the galactic centers. In M31, we confirm an outer ring at $\sim 11$ kpc, consistent with previous studies, and identify for the first time an inner ring at $\sim 7$ kpc, possibly associated with star formation episodes. Comparing age gradients at different angles, we find a consistent general trend of ages increasing monotonically with radius. At the same time, we observe smaller-scale differences, particularly in the $90^\circ$-$180^\circ$ quadrant, suggesting asymmetric star formation and possible dynamical influences. In contrast, M33 displays a steeper global age gradient, indicating a higher concentration of young stars toward its center. This study highlights the utility of Cepheids as stellar population tracers, providing insights into the star formation and dynamical evolution of spiral galaxies. Future works will extend this methodology to additional galaxies.

We investigate the unusual emission line luminosity ratios observed in the JADES NIRSpec spectroscopy of GN-z11, which reveal exceptionally strong emission lines and a significant detection of the rarely observed N III] $\lambda1748-1753$\r{A} multiplet. These features suggest an elevated N/O abundance, challenging existing models of stellar populations and nebular emission. To assess whether Wolf-Rayet (WR) stars can account for the observed line ratios, we construct a suite of stellar and nebular models incorporating high-resolution stellar spectral libraries, enabling a more accurate treatment of WR evolution and its influence on the ionising radiation field. We find that the inclusion of WR stars is essential for reproducing the observed position of GN-z11 in the C III]/He II versus C III]/C iv diagnostic plane, resolving discrepancies from previous studies. The model-derived metallicity (0.07$\lesssim$Z/Z$_{\odot}\lesssim$0.15), ionisation parameter ($\log\,U$$\approx$-2) and stellar ages are consistent with the literature estimates. However, our models under-predict the N III/O III] ratio, suggesting that WR stars alone cannot fully explain the nitrogen enrichment. This suggests that additional mechanisms, such as rapid chemical enrichment in a young, metal-poor environment, may be necessary to explain the nitrogen excess. While our models successfully reproduce most observed line ratios, further refinements to the models are needed to fully characterise the stellar populations and the enrichment processes of high-redshift galaxies like GN-z11.

The intergalactic medium (IGM) comprises all the matter that lies between galaxies. Hosting the vast majority ($\gtrsim 90\%$) of the baryons in the Universe, the IGM is a critical reservoir and probe for cosmology and astrophysics, providing insights into large-scale structure formation and galaxy evolution. In this Chapter, we present an overview of the general properties of the IGM, focusing on their dependence on cosmic environment and cosmic time. Emphasis is given to the basic physical principles that allow us to model the density, temperature, and ionization state of the IGM, supported by results from cosmological hydrodynamical simulations. We also cover the foundational principles of quasar spectroscopy used to probe the IGM in absorption, with a particular focus on HI absorption lines. Finally, we briefly discuss future prospects and complementary observational techniques to enhance our understanding of the IGM.

We present the first results from the CAPERS survey, utilizing PRISM observations with the JWST/NIRSpec MSA in the PRIMER-UDS field. With just 14 % of the total planned data volume, we spectroscopically confirm two new bright galaxies ($M_{\rm UV}\sim -20.4$) at redshifts $z = 10.562\pm0.034$ and $z = 11.013\pm0.028$. We examine their physical properties, morphologies, and star formation histories, finding evidence for recent bursting star formation in at least one galaxy thanks to the detection of strong (EW$_0\sim70$ A) H$\gamma$ emission. Combining our findings with previous studies of similarly bright objects at high-$z$, we further assess the role of stochastic star formation processes in shaping early galaxy populations. Our analysis finds that the majority of bright ($M_{\rm UV}\lesssim -20$) spectroscopically-confirmed galaxies at $z>10$ were likely observed during a starburst episode, characterized by a median SFR$_{10}$/SFR$_{100}\sim2$, although with substantial scatter. Our work also finds tentative evidence that $z>10$ galaxies are more preferentially in a bursting phase than similarly bright $z\sim6$ galaxies. We finally discuss the prospects of deeper spectroscopic observations of a statistically significant number of bright galaxies to quantify the true impact of bursting star formation on the evolution of the bright end of the ultraviolet luminosity function at these early epochs.

The next generation of weak gravitational lensing surveys has the potential to place stringent constraints on cosmological parameters. However, their analysis is limited by systematics such as the intrinsic alignments of galaxies, which alter weak lensing convergence and can lead to biases in cosmological parameter estimations. In this work, we investigate the impact of intrinsic alignments on non-Gaussian statistics of the weak lensing field using galaxy shapes derived from the IllustrisTNG hydrodynamical simulation. We create two catalogs of ray-traced convergence maps: one that includes the measured intrinsic shape of each galaxy and another where all galaxies are randomly rotated to eliminate intrinsic alignments. We compare an exhaustive list of weak lensing statistics between the two catalogs, including the shear-shear correlation function, the map-level angular power spectrum, one-point, peak count, minimum distribution functions, and Minkowski functionals. For each statistic, we assess the level of statistical distinguishability between catalogs for a set of future survey angular areas. Our results reveal strong small-scale correlation in the alignment of galaxies and statistically significant boosts in weak lensing convergence in both positive and negative directions for high-significance peaks and minimums, respectively. Weak lensing analyses utilizing non-Gaussian statistics must account for intrinsic alignments to avoid significantly compromised cosmological inferences.

Supernova 2023ixf is a type IIP supernova that was observed in May 2023 in the spiral galaxy Messier 101. This was the closest supernova observed of the decade, making this an exciting discovery. Combining the observed brightness and duration with theoretical scaling relations, we model the lightcurve of this supernova in order to reveal the properties of the progenitor star at the time of explosion, including its mass, radius, and explosion energy. We simulate these explosions using the stellar evolution and radiation-hydrodynamics codes MESA+STELLA. We find that SN2023ixf is not easily explained with "normal" stellar evolution, and only models with a small mass of H-rich ejecta can fit the lightcurve. We also find that the late time properties of the lightcurve are better fit by a higher initial-mass star with substantial mass loss during its lifetime, as compared to models with lower initial mass and less mass loss.

Approximately half of all disc galaxies exhibit appreciable warps in both their stellar and HI discs. The typical warp amplitude is small (a few degrees), becoming noticeable only at the periphery of the galaxy disc. As a result, warps remain a relatively poorly studied phenomenon. In this study, we investigate a large sample of distant edge-on galaxies (approximately 1,000 objects) to examine the frequency and characteristics of stellar disc warps up to a redshift of $z \sim 2$. For the selected galaxies, we utilize HST data from the Cosmic Evolution Survey field and JWST observations from the Cosmic Dawn Center Archive. We measure the properties of disc warps and investigate their evolution as a function of redshift. Our results indicate a potential evolution in the observed frequency of strong S-shaped warps (with an amplitude greater than 4$^\circ$) in stellar discs as a function of redshift. At $z \approx 2$, the frequency of strong warps reaches approximately 50%, while at $z \approx 0$, this fraction decreases to around 10-15%. We attribute the observed evolution in the occurrence of strong warps to the changing frequency of galaxy interactions and mergers. If galaxy interactions represent one of the primary mechanisms responsible for the formation of warps, then the prevalence of vertical disc deformations should increase in tandem with the rising interaction and merger rate.

Understanding the processes that transform star-forming galaxies into quiescent ones is key to unraveling the role of environment in galaxy evolution. We present measurements of the luminosity functions (LFs) and stellar mass functions (SMFs) of passive red-sequence galaxies in four galaxy clusters at $0.8 < z < 1.3$, selected using deep VLT observations complemented with data from the GCLASS and GOGREEN surveys. We find a significant enhancement in the abundance of faint/low-mass passive galaxies in both the LFs and SMFs of all four clusters compared to the field. This is further evidenced by a shallower low-mass slope in the composite passive cluster SMF, which yields a Schechter parameter $\alpha = -0.54^{+\,0.03}_{-0.03}$, compared to $\alpha = 0.12^{+\,0.01}_{-0.01}$ for the field. Our findings indicate that quenching processes that act in clusters are enhanced compared to the field, suggesting that environmental quenching mechanisms may already be active by $z\sim1$. To reproduce the observed passive cluster SMF, we estimate that $25\pm5\%$ of the star-forming field population that falls into the cluster must have been quenched. Our results largely support traditional quenching models but highlight the need for deeper studies of larger cluster samples to better understand the role of environmental quenching in the distant Universe.

Dark matter (DM) particles can interact with particles of the Standard Model. Although there exist constraints from direct and indirect detection experiments, the dynamical evolution of astrophysical objects could provide a promising probe for these interactions. Obtaining astrophysical predictions is challenging and limited by our ability to simulate scatterings between DM and baryonic particles within N-body and hydrodynamics simulations. We develop a novel scheme that allows simulating these interacting dark matter (IDM) models and accurately accounts for their angular and velocity dependence, as well as the mass ratio between the DM and baryonic scattering partners. To describe DM-baryon interactions, we use an N-body code together with its implementation of smoothed-particle hydrodynamics and meshless finite mass. The interaction itself is realised in a pairwise fashion by creating a virtual scattering partner from the baryonic particle and letting it interact with a DM particle using a scattering routine initially developed for self-interacting dark matter. After the interaction, the virtual particle is rejoined with the baryonic particle, fulfilling energy and momentum conservation. Through several test problems, we demonstrate that we can reproduce their analytic solutions with our IDM scheme. We comment on various numerical aspects and challenges as well as describe the limitations of our numerical scheme. Furthermore, we study the impact of IDM on halo formation with a collapsing overdensity. Overall, it is possible to accurately model IDM within N-body and hydrodynamics simulations, commonly used in astrophysics. In consequence, our scheme allows for making novel predictions and obtaining new constraints of DM-baryon scattering.