The CGM around unobscured AGN has received much attention in recent years. Comparatively, nebulae associated with obscured AGN are less studied. Here, we simulate the Ly$\alpha$, H$\alpha$, and HeII nebulae around the two types of AGN at $z=2-3$ with ten massive systems from the FIRE simulations based on the unified model to show their differences and to test if they can be used to constrain the AGN model. We post-process the data with the CLOUDY and the Ly$\alpha$ radiative transfer code, RASCAS. Overall, we find that the Ly$\alpha$ nebulae around the unobscured AGN (type-I nebulae) and obscured AGN (type-II nebulae) do not exhibit significant differences in the luminosity, area, and HeII/Ly$\alpha$ when the simulated cutout is set to the halo virial radius. Whereas, the type-II nebulae exhibit less symmetric morphologies, flatter surface brightness profiles, and larger emission line widths (at $R\geq 10$ kpc) than those of the type-I nebulae. These nebulae properties exhibit complicated correlations with the AGN, indicating that nebulae observations can be applied to constrain the AGN engine. However, independent observations on nebulae in the mentioned emissions are insufficient to test the unified model as a priori in observations is not possible to know the direction and opening angle of the ionization cone. We prompt that the joint observations of Ly$\alpha$ nebulae and radio jets can help to reveal the ionization cone to probe the unified model. Our calculations suggest that this method requires $\geq 75$ type-II Ly$\alpha$ nebulae with current instruments to reach a confidence level of $\geq 95\%$.

JWST has revealed an abundance of supermassive black holes (BHs) in the early Universe, and yet the lowest mass seed black holes that gave rise to these populations remain elusive. Here we present a systematic search for broad-line Active Galactic Nuclei (AGNs) in some of the faintest high-$z$ galaxies surveyed yet by combining ultra-deep JWST/NIRSpec G395M spectroscopy with the strong lensing aid in Abell S1063. By employing the profile of the [OIII]$\lambda 5007$ emission lines as a template for narrow-line components and carefully cross-validating with mock observations, we identify a sample of ten broad-line AGNs at $4.5

We present three cosmological radiation-hydrodynamic zoom simulations of the progenitor of a Milky Way-mass galaxy from the MEGATRON suite. The simulations combine on-the-fly radiative transfer with a detailed non-equilibrium thermochemical network (81 ions and molecules), resolving the cold and warm gas in the circumgalactic medium (CGM) on spatial scales down to 20 pc and on average 200 pc at cosmic noon. Comparing our full non-equilibrium calculation with local radiation to traditional post-processed photoionization equilibrium (PIE) models assuming a uniform UV background (UVB), we find that non-equilibrium physics and local radiation fields fundamentally impact the thermochemistry of the CGM. Recombination lags and local radiation anisotropy shift ions away from their PIE+UVB values and modify covering fractions (for example, HI damped Ly$\alpha$ absorbers differ by up to 40%). In addition, a resolution study with cooling-length refinement allows us to double the resolution in the cold and warm CGM gas, reaching 120 pc on average. When refining on cooling length, the mass of the lightest cold clumps decreases tenfold to $\approx 10^4\,M_\odot$, their boundary layers develop sharper ion stratification, and the warm gas is better resolved, boosting the abundance of warm gas tracers such as CIV and OIII. Together, these results demonstrate that non-equilibrium thermochemistry coupled to radiative transfer, combined with physically motivated resolution criteria, is essential to predict circumgalactic absorption and emission signatures and to guide the design of targeted observations with existing and upcoming facilities.

Dwarf galaxies provide powerful laboratories for studying galaxy formation physics. Their early assembly, shallow gravitational potentials, and bursty, clustered star formation histories make them especially sensitive to the processes that regulate baryons through multi-phase outflows. Using high-resolution, cosmological zoom-in simulations of a dwarf galaxy from \textit{the Pandora suite}, we explore the impact of stellar radiation, magnetic fields, and cosmic ray feedback on star formation, outflows, and metal retention. We find that our purely hydrodynamical model without non-thermal physics - in which supernova feedback is boosted to reproduce realistic stellar mass assembly - drives violent, overly enriched outflows that suppress the metal content of the host galaxy. Including radiation reduces the clustering of star formation and weakens feedback. However, the additional incorporation of cosmic rays produces fast, mass-loaded, multi-phase outflows consisting of both ionized and neutral gas components, in better agreement with observations. These outflows, which entrain a denser, more temperate ISM, exhibit broad metallicity distributions while preserving metals within the galaxy. Furthermore, the star formation history becomes more bursty, in agreement with recent JWST findings. These results highlight the essential role of non-thermal physics in galaxy evolution and the need to incorporate it in future galaxy formation models.

Numerous high-$z$ galaxies have recently been observed with the James Webb Space Telescope (JWST), providing new insights into early galaxy evolution. Their physical properties are typically derived through spectral energy distribution (SED) fitting, but the reliability of this approach for such early systems remains uncertain. Applying {\sc Bagpipes} on simulated SEDs at $z=6$ from the {\sc Sphinx$^{20}$} cosmological simulation, we examine uncertainties in the recovery of stellar masses, star formation rates (SFR$_{10}$), and stellar metallicities from mock JWST/Near-Infrared Camera photometry. Even without dust or emission lines, fitting the intrinsic stellar continuum overestimates the stellar mass by about 60\% on average (and by up to a factor of five for low-mass galaxies with recent starbursts) and underestimates SFR$_{10}$ by a factor of two, owing to inaccurate star formation histories and age-metallicity degeneracies. The addition of dust and nebular emission further amplifies these biases, yielding offsets of approximately +0.3 and -0.4 dex in stellar mass and SFR$_{10}$, respectively, while leaving stellar metallicities largely unconstrained. Incorporating bands free of strong emission lines, such as F410M, helps mitigate stellar mass overestimation by disentangling line emission from older stellar populations. We also find that best-fit or likelihood-weighted estimates are generally more accurate than median posterior values. Although stellar mass functions are reproduced reasonably well, the slope of the star formation main sequence depends sensitively on the adopted fitting model. Overall, these results underscore the importance of careful modelling when interpreting high-$z$ photometry, particularly for galaxies with recent star formation burst and/or strong emission lines, to minimise systematic biases in derived physical properties.

We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ multi-frequency, on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization as well as the growth of PMFs. We run a series of $(5\,\text{cMpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field ($B_0$) and its spectral index ($n_B$). We find that PMFs that have $n_B > -0.562\log_{10}\left(\frac{B_0}{1{\rm n}G}\right) - 3.35$ produce electron optical depths ($\tau_e$) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For $n_B\geq-2.9$, our constraints are considerably tighter than the $\sim{\rm n}G$ constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Lya forest due to PMFs will be challenging because photoionisation and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21cm signal is a sensitive indicator of the properties of the PMFs, even for those that satisfy our $\tau_e$ constraint. Furthermore, strong PMFs can marginally increase the escape of LyC photons by up to 25\% and shrink the effective radii of galaxies by $\sim44\%$ which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500}=-13}$ can probe the sources that provide the majority of photons for reionization out to $z=12$.

The Euclid mission aims to measure the positions, shapes, and redshifts of over a billion galaxies to provide unprecedented constraints on the nature of dark matter and dark energy. Achieving this goal requires a continuous reassessment of the mission's scientific performance, particularly in terms of its ability to constrain cosmological parameters, as our understanding of how to model large-scale structure observables improves. In this study, we present the first scientific forecasts using CLOE (Cosmology Likelihood for Observables in Euclid), a dedicated Euclid cosmological pipeline developed to support this endeavour. Using advanced Bayesian inference techniques applied to synthetic Euclid-like data, we sample the posterior distribution of cosmological and nuisance parameters across a variety of cosmological models and Euclid primary probes: cosmic shear, angular photometric galaxy clustering, galaxy-galaxy lensing, and spectroscopic galaxy clustering. We validate the capability of CLOE to produce reliable cosmological forecasts, showcasing Euclid's potential to achieve a figure of merit for the dark energy parameters $w_0$ and $w_a$ exceeding 400 when combining all primary probes. Furthermore, we illustrate the behaviour of the posterior probability distribution of the parameters of interest given different priors and scale cuts. Finally, we emphasise the importance of addressing computational challenges, proposing further exploration of innovative data science techniques to efficiently navigate the Euclid high-dimensional parameter space in upcoming cosmological data releases.

We report statistically significant detection of H I 21-cm emission from intermediate-redshift ($z\approx0.2$-0.6) galaxies. By leveraging multi-sightline galaxy survey data from the Cosmic Ultraviolet Baryon Survey (CUBS) and deep radio observations from the MeerKAT Absorption Line Survey (MALS), we have established a sample of $\approx6000$ spectroscopically identified galaxies in 11 distinct fields to constrain the neutral gas content at intermediate redshifts. The galaxies sample a broad range in stellar mass -- 8\lesssim\log{M_\rm{star}/\rm{M}_\odot}\lesssim11 with a median of \langle\log{M_\rm{star}/\rm{M}_\odot}\rangle_\rm{med}\approx10 -- and a wide range in redshift -- $0.24\lesssim z\lesssim0.63$ with a median of \langle z\rangle_\rm{med}=0.44. Our detected emission-line signal exceeds $4\,\sigma$ significance in the stacked spectra of all subsamples, and the observed total H I 21-cm line flux translates to a H I mass M_\rm{H\;I}\approx10^{10}\rm{M}_\odot. We find a high H I-to-stellar mass ratio of M_\mathrm{H\;I}/M_\rm{star}\approx6 for low-mass galaxies with \langle\log{M_\rm{star}/\rm{M}_\odot}\rangle \approx9.3 ($>3.7\,\sigma$). For galaxies with \langle\log{M_\rm{star}/\rm{M}_\odot}\rangle\approx10.6, we find M_\mathrm{H\;I}/M_\rm{star}\approx0.3 ($>4.7\,\sigma$). Additionally, the redshift evolution of H I mass in both low- and high-mass field galaxies, inferred from the stacked emission-line signal, aligns well with the expectation from the cosmic star formation history. This suggests that the overall decline in the cosmic star formation activity across the general galaxy population may be connected to a decreasing supply of neutral hydrogen. Finally, our analysis has revealed significant 21-cm signals at distances greater than 75 kpc from these intermediate-redshift galaxies, indicating a substantial reservoir of H I gas in their extended surroundings.

Jellyfish galaxies provide direct evidence of ram pressure stripping in cluster environments. We investigate the role of magnetic fields in the formation of jellyfish galaxies with a multiphase interstellar medium (ISM) using radiation magneto-hydrodynamic simulations. We impose magnetized (MHD) and non-magnetized (HD) winds on the gas-rich dwarf galaxies containing the magnetized or non-magnetized ISM. The MHD winds strip the disk gas more effectively than the HD winds because of the magnetic force acting against the local density gradient, which results in remarkably different ram pressure stripped features. The magnetic fields induced by the MHD winds generate a strong magnetic pressure, which forms smoothed disks and tail gas features. Since the stripped ISM in MHD wind cases travels while being nearly isolated from the intracluster medium (ICM), the stripped ISM mostly forms stars within 20~kpc of the galactic disks. In contrast, non-magnetized winds facilitate the efficient mixing of the stripped ISM with the ICM, resulting in the formation of abundant warm clouds that cool and collapse in the distant ($\sim50-100\,$kpc) tails at times of a few hundred Myr. Consequently, distant tail star formation occurs only in the HD wind runs. Finally, despite the different tail features, the star formation rates in the disk remain similar owing to the interplay between the increased gas stripping and the gas density increase in the disks of the MHD wind runs. These results suggest that the magnetized ICM may have a significant influence on jellyfish galaxies, whereas the magnetized ISM play a minor role.

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation.

Observationally mapping the relation between galaxies and the intergalactic medium (IGM) is of key interest for studies of cosmic reionization. Diffuse hydrogen gas has typically been observed in HI Lyman-$\alpha$ (Ly$\alpha$) absorption in the spectra of bright background quasars. However, it is important to extend these measurements to background galaxies as quasars become increasingly rare at high redshift and rarely probe closely separated sight-lines. Here we use deep integral field spectroscopy in the MUSE eXtremely Deep Field to demonstrate the measurement of the Ly$\alpha$ transmission at $z\approx4$ in absorption to a background galaxy at $z=4.77$. The HI transmission is consistent with independent quasar sight-lines at similar redshifts. Exploiting the high number of spectroscopic redshifts of faint galaxies (500 between $z=4.0-4.7$ within a radius of 8 arcmin) that are tracers of the density field, we show that Ly$\alpha$ transmission is inversely correlated with galaxy density, i.e. transparent regions in the Ly$\alpha$ forest mark under-dense regions at $z\approx4$. Due to large-scale clustering, galaxies are surrounded by excess HI absorption over the cosmic mean out to 4 cMpc/h. We also find that redshifts from the peak of the Ly$\alpha$ line are typically offset from the systemic redshift by +170 km/s. This work extends results from $z\approx 2 - 3$ to higher redshifts and demonstrates the power of deep integral field spectroscopy to simultaneously measure the ionization structure of the IGM and the large-scale density field in the early Universe.

The Multi Unit Spectroscopic Explorer (MUSE) is a second-generation VLT panoramic integral-field spectrograph currently in manufacturing, assembly and integration phase. MUSE has a field of 1x1 arcmin2 sampled at 0.2x0.2 arcsec2 and is assisted by the VLT ground layer adaptive optics ESO facility using four laser guide stars. The instrument is a large assembly of 24 identical high performance integral field units, each one composed of an advanced image slicer, a spectrograph and a 4kx4k detector. In this paper we review the progress of the manufacturing and report the performance achieved with the first integral field unit.

Understanding the escape of Lyman continuum (LyC) and Lyman $\alpha$ (Ly$\alpha$) photons from giant molecular clouds (GMCs) is crucial if we are to study the reionization of the Universe and to interpret spectra of observed galaxies at high redshift. To this end, we perform high-resolution, radiation-magneto-hydrodynamic simulations of GMCs with self-consistent star formation and stellar feedback. We find that a significant fraction (15-70%) of ionizing radiation escapes from the simulated GMCs with different masses ($10^5$ and $10^6\,M_\odot$), as the clouds are dispersed within about $2$-$5\,{\rm Myr}$ from the onset of star formation. The fraction of LyC photons leaked is larger when the GMCs are less massive, metal-poor, less turbulent, and less dense. The most efficient leakage of LyC radiation occurs when the total star formation efficiency of a GMC is about 20%. The escape of Ly$\alpha$ shows a trend similar to that of LyC photons, except that the fraction of Ly$\alpha$ photons escaping from the GMCs is larger ($f_{\rm esc}^{\rm Ly\alpha}\approx f_{900}^{0.27}$) and that a GMC with strong turbulence shows larger $f_{\rm esc}^{\rm Ly\alpha}$. The simulated GMCs show a characteristic velocity separation of $\Delta v\approx 120 \,{\rm km\,s^{-1}}$ in the time-averaged emergent Ly$\alpha$ spectra, suggesting that Ly$\alpha$ could be useful to infer the kinematics of the interstellar and circumgalactic medium. We show that Ly$\alpha$ luminosities are a useful indicator of the LyC escape, provided the number of LyC photons can be deduced through stellar population modeling. Finally, we find that the correlations between the escape fractions of Ly$\alpha$, ultraviolet photons at 1500A, and the Balmer $\alpha$ line are weak.

The physical origin of low escape fractions of ionizing radiation derived from massive star-forming galaxies at $z\sim3$-$4$ is not well understood. We perform idealised disc galaxy simulations to understand how galactic properties such as metallicity and gas mass affect the escape of Lyman Continuum (LyC) photons using radiation-hydrodynamic simulations with strong stellar feedback. We find that the luminosity-weighted escape fraction from a metal-poor ($Z=0.002$) galaxy embedded in a halo of mass $M_h\simeq10^{11}\,M_\odot$ is \left\simeq 10\,\%. Roughly half of the LyC photons are absorbed within scales of 100 pc, and the other half is absorbed in the ISM ($\lesssim 2\, {\rm kpc}$). When the metallicity of the gas is increased to $Z=0.02$, the escape fraction is significantly reduced to \left\simeq1\% because young stars are enshrouded by their birth clouds for a longer time. In contrast, increasing the gas mass by a factor of 5 leads to \left\simeq 5\, \% because LyC photons are only moderately absorbed by the thicker disc. Our experiments suggest that high metallicity is likely more responsible for the low escape fractions observed in massive star-forming galaxies, supporting the scenario in which the escape fraction is decreasing with increasing halo mass. Finally, negligible correlation is observed between the escape fraction and surface density of star formation or galactic outflow rates.

We identify a chain of galaxies along an almost straight line in the nearby Universe with a projected length of ~5 Mpc. The galaxies are distributed within projected distances of only 7-105 kpc from the axis of the identified filament. They have redshifts in a very small range of z=0.0361-0.0370 so that their radial velocities are consistent with galaxy proper motions. The filament galaxies are mainly star-forming and have stellar masses in a range of $\rm 10^{9.1}-10^{10.7}\,M_{\odot}$. We search for systems with similar geometrical properties in the full-sky mock galaxy catalogue of the MillenniumTNG simulations and find that although such straight filaments are unusual and rare, they are predicted by $\Lambda$CDM simulations (4% incidence). We study the cold HI gas in a 1.3 Mpc section of the filament through HI-21cm emission line observations and detect eleven HI sources, many more than expected from the HI mass function in a similar volume. They have HI masses $\rm 10^{8.5}-10^{9.5}\,M_{\odot}$ and are mostly within ~120 kpc projected distance from the filament axis. None of these HI sources has a confirmed optical counterpart. Their darkness together with their large HI-21cm line-widths indicate that they contain gas that might not yet be virialized. These clouds must be marking the peaks of the dark matter and HI distributions over large scales within the filament. The presence of such gas clouds around the filament spines is predicted by simulations, but this is the first time that the existence of such clouds in a filament is observationally confirmed.

We introduce the MUSE gAlaxy Groups in COSMOS (MAGIC) survey, which was built to study the impact of environment on galaxy evolution over the last 8 Gyr. It consists of 17 MUSE fields targeting 14 massive structures at intermediate redshift ($0.321.5$. The spectroscopic redshift completeness is high: in the redshift range of [OII] emitters (0.25 \le z < 1.5), where most of the groups are found, it globally reaches a maximum of 80% down to $z_{app}^{++}=25.9$, and locally decreases from $\sim 100$% to $\sim50$% in magnitude bins from $z_{app}^{++}=23-24$ to $z_{app}^{++}=25.5$. We find that the fraction of quiescent galaxies increases with local density and with the time spent in groups. A morphological dichotomy is also found between bulge-dominated quiescent and disk-dominated star-forming galaxies. As environment gets denser, the peak of the stellar mass distribution shifts towards M_*>10^{10}~M_\odot, and the fraction of galaxies with M_*<10^9~M_\odot decreases significantly, even for star-forming galaxies. We also highlight peculiar features such as close groups, extended nebulae, and a gravitational arc. Our results suggest that galaxies are preprocessed in groups of increasing mass before entering rich groups and clusters. We publicly release two catalogs containing the properties of galaxies and groups, respectively.

We present the first results of V\textit{z}-GAL, a high-redshift CO(\textit{J} = 1 $-$ 0) large survey with the Karl G. Jansky Very Large Array, targeting 92 \textit{Herschel}-selected, infrared-luminous, dusty star-forming galaxies (DSFGs). This flux density-limited sample includes 106 DSFGs in total, along with 14 galaxies from a recent pilot study -- altogether doubling the to-date available high-redshift \coonezero observations. These sources cover redshifts 1 to 6 with available mid/high-\textit{J} CO transitions from the Northern Extended Millimeter Array \zgal survey. We detect \coonezero emission in 90/92 targets above signal-to-noise ratio of 2, while two DSFGs remain undetected. \vzgal also covers additional \cotwoone emission lines in 10 of these sources. We find gas masses of the entire \vzgal sample to be $\rm (\alpha_{CO}/{4.0}) \mu {M}_{\rm H_2}$ = $(2-20) \times {10}^{11}~\mathrm{M_{\odot}}$, calibrated using \coonezero line luminosities. We confirm that these DSFGs -- with derived gas depletion timescales of $(50-600)$ Myr -- represent a heterogeneous population, potentially comprising both main-sequence galaxies and starbursts in the early Universe. Using \coonezero luminosities as an anchor, we robustly derive the CO brightness temperature ratios with a good statistical significance up to \textit{J} = 6. Our measurements reveal a broad range of gas excitation conditions across the \vzgal DSFGs. We also explore [CI](1$-$0)/CO(1$-$0) ratios in 23 \vzgal galaxies that have existing detections of atomic carbon ground-state emission, [CI](1$-$0), from the \zgal survey. Our results show similar [CI]/CO ratio values across cosmic time for starbursts and local star-forming galaxy populations, supporting the use of \cionezero as an alternative cold gas tracer.

Processing of raw data from modern astronomical instruments is nowadays often carried out using dedicated software, so-called "pipelines" which are largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at ESO's Paranal observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software, a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data gets transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis.