Studies of galaxy populations classified according to their kinematic behaviours and dynamical state using the Projected Phase Space Diagram (PPSD) are affected by misclassification and contamination, leading to systematic errors in determining the characteristics of the different galaxy classes. We propose a method to statistically correct the determination of galaxy properties' distributions accounting for the contamination caused by misclassified galaxies from other classes. Using a sample of massive clusters and galaxies in their surroundings taken from the MultiDark Planck 2 simulation combined with the semi-analytic model of galaxy formation SAG, we compute the confusion matrix associated to a classification scheme in the PPSD. Based on positions in the PPSD, galaxies are classified as cluster members, backsplash galaxies, recent infallers, infalling galaxies, and interlopers. This classification is determined using probabilities calculated by the code ROGER, along with a threshold criterion. By inverting the confusion matrix, we are able to get better determinations of distributions of galaxy properties such as colour. Compared to a direct estimation based solely on the predicted galaxy classes, our method provides better estimates of the mass-dependent colour distribution for the galaxy classes most affected by misclassification: cluster members, backsplash galaxies, and recent infallers. We apply the method to a sample of observed X-ray clusters and galaxies. Our method can be applied to any classification of galaxies in the PPSD, and to any other galaxy property besides colour, provided an estimation of the confusion matrix. Blue, low-mass galaxies in clusters are almost exclusively recent infaller galaxies that have not yet been quenched by the environmental action of the cluster. Backsplash galaxies are on average redder than expected.

We investigated the population of low-mass central galaxies with M-star=10(9.5)-10(10) h(-1) M-circle dot, inhabiting regions near massive groups and clusters of galaxies using the IllustrisTNG300 and MDPL2-SAG simulations. We set out to study their evolutionary histories, aiming to find hints about the large-scale conformity signal they produce. We also used a control sample of central galaxies with the same stellar mass range located far away from massive structures. For both samples, we find a subpopulation of galaxies accreted by another halo in the past, but now considered central galaxies; we refer to these objects as former satellites. The number of former satellites is higher for quenched central galaxies near massive systems, with fractions of 45% and 17% in IllustrisTNG300 and MDPL2-SAG, respectively. The differences in the numerical resolution of each simulation lead to the different fractions of former satellites. Our results in TNG300 show that former satellites "pollute" the sample of central galaxies because they suffered environmental processes when they were satellites hosted typically by massive dark matter halos (M-200 >= 10(13) h(-1) M-circle dot) since z less than or similar to 0.5. After removing former satellites, the evolutionary trends for quenched central galaxies near massive structures are fairly similar to those of the quenched control galaxies, showing small differences at low redshift. For MDPL2-SAG instead, former satellites were hosted by less massive halos, with a mean halo mass around 10(11.4) h(-1) M-circle dot, and the evolutionary trends remain equal before and after removing former satellite galaxies. We also measured the two-halo conformity, that is, the correlation in the specific star formation rate between low-mass central galaxies and their neighbors at megaparsec scales, and how former satellites contribute to the signal at three different redshifts: z=0, 0.3, and 1. The time evolution of the conformity signal in the simulations presents apparent contradictory results: it decreases from z=0 to z=1 in MDPL2-SAG, while it increases in TNG300. However, after removing former satellites in the latter, the signal is strongly reduced, but practically does not change at z <= 0.3, and it disappears at z = 1. We compare our findings with recent literature data and discuss the conformity measurements, as different approaches can lead to varying results.

Boxy/peanut (b/p) or X-shaped bulges have been extensively explored with theory and numerical simulations of isolated galaxies. However, it is only recently that advances in hydrodynamical cosmological simulations have made it possible to explore b/p bulges in a cosmological setting, with much remaining to be understood about their formation and evolution. By using the Auriga magnetohydrodynamical cosmological zoom-in simulations, we characterize the structural parameters of b/p bulges and how they form and evolve throughout cosmic history. We develop a method for estimating the b/p strength that allows us to identify the formation time and size of these structures. We find that b/p bulges in Auriga form between similar to 1.1 and 1.6 Gyr after bar formation, following a 'buckling' episode; some galaxies undergo multiple bucklings and events of b/p growth, with some b/p structures 'dissolving' between buckling events. We find that at z=0, the b/p bulges have an extent of almost half the bar length. Finally, we analyse the evolution of the b/p fraction over redshift, finding that at z=0, two-thirds of galaxies host a bar, and of these, 45 per cent have a b/p. This b/p fraction is within the observed range at z=0, although on the low end as compared to some observational studies. The b/p fraction decreases to 20 per cent at z=0.5, and falls to zero at z similar to 1; this is in line with the observed trend of declining b/p fraction with redshift. We discuss possible culprits for the apparent mismatch in b/p occurrence between observations and cosmological simulations, what causes them to form (or not) in these simulations, and what this might reveal about models of galaxy formation and evolution.

Next-generation photometric and spectroscopic surveys will detect faint galaxies in massive clusters, advancing our understanding of galaxy formation in dense environments. Comparing these observations with theoretical models requires high-resolution cluster simulations. Hydrodynamical simulations resolve galaxy properties in halos, but face challenges simulating low-mass galaxies due to computational limitations. In contrast, dark matter-only (DMO) simulations can provide higher resolution but need models to populate subhalos with galaxies. In this work, we introduce a fast and efficient emulator of hydrodynamical cluster simulations, based on the semi-analytic models (SAMs) SAGE and SAG. The SAMs are calibrated using cluster galaxies from hydrodynamical simulations at intermediate resolution, ensuring consistency in stellar masses and luminosities across redshifts. These SAMs are then applied to DMO simulations from The Three Hundred Project at three resolutions. We show that the SAG model, unlike SAGE, better emulates galaxy properties even at the highest resolution. This improvement is due to the treatment of orphan galaxies, which contribute significantly to the total population. SAG enables the study of dwarf galaxies down to stellar masses of 10^7 solar masses at the highest resolution, an order of magnitude smaller than those in the hydrodynamical simulations, corresponding to approximately four magnitudes fainter. This shows that SAMs can be effectively calibrated to provide fast and accurate predictions of hydrodynamical simulations, offering an efficient alternative to explore galaxy populations in dense environments.

The intrinsic properties of galaxies are influenced by their environments, underscoring the environment's critical role in galaxy formation and evolution. Traditionally, these environments are categorized into four fixed classifications: knots, filaments, walls, and voids, which collectively describe the complex organization of galaxies within large-scale structures. We propose an alternative description that complements the traditional quadripartite categorization by introducing a continuous framework, allowing for a more nuanced examination of the relationship between the intrinsic properties of galaxies and their environments. This complementary description is applied using one of the most prevalent methodologies: categorization using the eigenvalues of the Hessian matrix extracted from the matter density field. We integrated our findings into a semi-analytical model of galaxy formation, combined with cosmological numerical simulations, to analyze how the intrinsic properties of galaxies are influenced by environmental changes. In our study, we find a continuous distribution of eigenvalue ratios, revealing a clear dependence of galaxy properties on their surrounding environments. This method allowed us to identify critical values at which transitions in the behavior of key astrophysical galaxy properties become evident.

Context. Nebular emission lines are powerful diagnostics for the physical processes at play in galaxy formation and evolution. Moreover, emission-line galaxies (ELGs) are one of the main targets of current and forthcoming spectroscopic cosmological surveys. Aims. We investigate the contributions to the line luminosity functions (LFs) of different galaxy populations in the local Universe, providing a benchmark for future surveys of earlier cosmic epochs. Methods. The large statistics of the observations from the SDSS DR7 main galaxy sample and the MPA-JHU spectral catalog enabled us to precisely measure the H α , H β , [O II ], [O III ], and, for the first time, the [N II ], and [S II ] emission-line LFs over ∼2.4 Gyrs in the low- z Universe, 0.02 < z < 0.22. We present a generalized 1/ V max LF estimator capable of simultaneously correcting for spectroscopic, r -band magnitude, and emission-line incompleteness. We studied the contribution to the LF of different types of ELGs classified using two methods: (i) the value of the specific star formation rate (sSFR), and (ii) the line ratios on the Baldwin–Phillips–Terlevich (BPT) and the WHAN (i.e., H α equivalent width, EW H α , versus the [N II ]/H α line ratio) diagrams. Results. The ELGs in our sample are mostly star forming, with 84 percent having sSFR > 10 −11 yr −1 . When classifying ELGs using the BPT+WHAN diagrams, we find that 63.3 percent are star forming, only 0.03 are passively evolving, and 1.3 have nuclear activity (Seyfert). The rest are low-ionization narrow emission-line regions (LINERs) and composite ELGs. We found that a Saunders function is the most appropriate to describe all of the emission-line LFs, both observed and dust-extinction-corrected (i.e., intrinsic). They are dominated by star-forming regions, except for the bright end of the [O III ] and [N II ] LFs (i.e., L [N II ] > 10 42 erg s −1 , L [O III ] > 10 43 erg s −1 ), where the contribution of Seyfert galaxies is not negligible. In addition to the star-forming population, composite galaxies, and LINERs are the ones that contribute the most to the ELG numbers at L < 10 41 erg s −1 . We do not observe significant evolution with redshift of our ELGs at 0.02 < z < 0.22. All of our results, including data points and analytical fits, are publicly available. Conclusions. Local ELGs are dominated by star-forming galaxies, except for the brightest [N II ] and [O III ] emitters, which have a large contribution of Seyfert galaxies. The local line luminosity functions are best described by Saunders functions. We expect these two conclusions to hold up at higher redshifts for the ELG targeted by current cosmological surveys, such as DESI and Euclid .

ABSTRACT Effective yields, yeff, are defined by fundamental galaxy properties (i.e. stellar mass M⋆, gas mass Mgas, and gas-phase metallicity). For a closed-box model, yeff is constant and equivalent to the mass in metals returned to the gas per unit mass locked in long-lived stars. Deviations from such behaviour have been often considered observational signatures of past feedback events. By analysing eagle simulations with different feedback models, we evaluate the impact of supernovae (SNe) and active galactic nuclei (AGNs) feedback on yeff at redshift z = 0. When removing supermassive black holes (BHs) and, hence, AGN effects, in simulations, galaxies are located around a plane in the M⋆–Mgas–O/H parameter space (being O/H a proxy for gas metallicity, as usual), with such a plane roughly describing a surface of constant yeff. As the ratio between BH mass and M⋆ increases, galaxies deviate from that plane towards lower yeff as a consequence of AGN feedback. For galaxies not strongly affected by AGN feedback, a stronger SN feedback efficiency generates deviations towards lower yeff, while galaxies move towards the opposite side of the plane (i.e. towards higher values of yeff) as SN feedback becomes weaker. Star-forming galaxies observed in the Local Universe are located around a similar 3D plane. Our results suggest that the features of the scatter around the observed plane are related to the different feedback histories of galaxies, which might be traced by yeff.

A significant fraction of local galaxies exhibit stellar bars, non-axisymmetric structures composed of stars, gas, and dust. Identifying key differences between the properties of barred and unbarred galaxies can uncover clues about the conditions for triggering bar formation. We explore the early stages of bar formation in a small sample of disc barred galaxies extracted from the TNG50 cosmological simulation, and compare their properties to those of unbarred galaxies. According to our results, the most important difference between barred and unbarred galaxies is that the former have systematically higher fractions of stellar to dark matter mass in their inner regions, from very early stages and prior to the formation of the bars. They harbour high initial gas content, fostering increased star formation rates and leading to a central mass concentration that grows faster over time compared to unbarred galaxies. Examining the evolution of the halo spin within 10 ckpc reveals that barred galaxies have higher angular momentum transfer from the disc to the halo. Curiously, both barred and unbarred galaxies share similar initial low values of the halo spin, consistent with those proposed in the literature for bar formation. Furthermore, we evaluate existing stability criteria to capture the complexity of the process, and investigate the effects of mergers, flybys, and environment as possible drivers of bar formation. We find no clear link between mergers and disc instabilities resulting in the formation of bars, even though some of the simulated barred galaxies might have been influenced by these events.

The feedback mechanisms triggered by supernova (SN) events and active galactic nuclei (AGN) play a central role in regulating the star formation and shaping galaxy properties. However, quantifying the impact and efficiency of these processes remains a challenge. In this study, we use the EAGLE cosmological hydrodynamics simulations to examine different models of SN and AGN feedback. Our goal is to investigate how variations in these processes impact the properties of simulated galaxy populations. Specifically, we focus on the analysis of effective yields, evaluating their capability to trace the effects of feedback processes on scaling relations. Our work contributes to a deeper understanding of the complex relationship between different feedback scenarios and the evolution of galaxies.

In the submm regime, spectral line scans and line intensity mapping (LIM) are new promising probes for the cold gas content and star formation rate of galaxies across cosmic time. However, both of these two measurements suffer from field-to-field variance. We study the effect of field-to-field variance on the predicted CO and [CII] power spectra from future LIM experiments such as CONCERTO, as well as on the line luminosity functions (LFs) and the cosmic molecular gas mass density that are currently derived from spectral line scans. We combined a 117 $\rm deg^2$ dark matter lightcone from the Uchuu cosmological simulation with the simulated infrared dusty extragalactic sky (SIDES) approach. We find that in order to constrain the CO LF with an uncertainty below 20%, we need survey sizes of at least 0.1 $\rm deg^2$. Furthermore, accounting for the field-to-field variance using only the Poisson variance can underestimate the total variance by up to 80%. The lower the luminosity is and the larger the survey size is, the higher the level of underestimate. At $z$<3, the impact of field-to-field variance on the cosmic molecular gas density can be as high as 40% for the 4.6 arcmin$^2$ field, but drops below 10% for areas larger than 0.2 deg$^2$. However, at $z>3$ the variance decreases more slowly with survey size and for example drops below 10% for 1 deg$^2$ fields. Finally, we find that the CO and [CII] LIM power spectra can vary by up to 50% in $\rm 1 deg^2$ fields. This limits the accuracy of the constraints provided by the first 1 deg$^2$ surveys. The level of the shot noise power is always dominated by the sources that are just below the detection thresholds. We provide an analytical formula to estimate the field-to-field variance of current or future LIM experiments. The code and the full SIDES-Uchuu products (catalogs, cubes, and maps) are publicly available.

Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$\nu^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$\nu^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST).

In this work, we investigate the structural properties, distribution and abundance of ΛCDM dark matter subhaloes using the Phi-4096 and Uchuu suite of N-body cosmological simulations. Thanks to the combination of their large volume, high mass resolution and superb statistics, we are able to quantify – for the first time consistently over more than seven decades in ratio of subhalo-to-host-halo mass – dependencies of subhalo properties on mass, maximum circular velocity, V max , host halo mass, and distance to host halo centre. We also dissect the evolution of these dependencies over cosmic time. We provide accurate fits for the subhalo mass and velocity functions, both exhibiting decreasing power-law slopes and with no significant dependence on redshift. We also find subhalo abundance to depend weakly on host halo mass. Subhalo structural properties are codified via a concentration parameter, c V , that does not depend on any pre-defined density profile and relies only on V max . We derive the c V −V max relation and find an important dependence on distance of the subhalo to the host halo centre. Interestingly, we also find subhaloes of the same mass to be significantly more concentrated when they reside inside more massive hosts. Finally, we investigate the redshift evolution of c V , and provide accurate fits. Our results offer an unprecedented detailed characterization of the subhalo population, consistent over a wide range of subhalo and host halo masses, as well as cosmic times. Thus, we expect our work to be particularly useful for any future research involving dark matter halo substructure.

We study the properties of associations of dwarf galaxies and their dependence on the environment. Associations of dwarf galaxies are extended systems composed exclusively of dwarf galaxies, considering as dwarf galaxies those galaxies less massive than $M_{\star, \rm max} = 10^{9.0}$ ${\rm M}_{\odot}\,h^{-1}$. We identify these particular systems using a semi-analytical model of galaxy formation coupled to a dark matter only simulation in the $\Lambda$ Cold Dark Matter cosmological model. To classify the environment, we estimate eigenvalues from the tidal field of the dark matter particle distribution of the simulation. We find that the majority, two thirds, of associations are located in filaments ($ \sim 67$ per cent), followed by walls ($ \sim 26 $ per cent), while only a small fraction of them are in knots ($ \sim 6 $ per cent) and voids ($ \sim 1 $ per cent). Associations located in more dense environments present significantly higher velocity dispersion than those located in less dense environments, evidencing that the environment plays a fundamental role in their dynamical properties. However, this connection between velocity dispersion and the environment depends exclusively on whether the systems are gravitational bound or unbound, given that it disappears when we consider associations of dwarf galaxies that are gravitationally bound. Although less than a dozen observationally detected associations of dwarf galaxies are currently known, our results are predictions on the eve of forthcoming large surveys of galaxies, which will enable these very particular systems to be identified and studied.

We study the population of backsplash galaxies at $z=0$ in the outskirts of massive, isolated clusters of galaxies taken from the MDPL2-SAG semi-analytic catalogue. We consider four types of backsplash galaxies according to whether they are forming stars or passive at three stagesin their lifetimes: before entering the cluster, during their first incursion through the cluster, and after they exit the cluster. We analyse several geometric, dynamic, and astrophysical aspects of the four types at the three stages. Galaxies that form stars at all stages account for the majority of the backsplash population ($58\%$) and have stellar masses typically below $M_\star\sim 3\times 10^{10} h^{-1}{\rm M}_\odot$ that avoid the innermost cluster's regions and are only mildly affected by it. In a similar mass range, galaxies that become passive after exiting the cluster ($26\%$) follow orbits characterised by small pericentric distance and a strong deflection by the cluster potential well while suffering a strong loss of both dark matter and gas content. Only a small fraction of our sample ($4\%$) become passive while orbiting inside the cluster. These galaxies have experienced heavy pre-processing and the cluster's tidal stripping and ram pressure provide the final blow to their star formation. Finally, galaxies that are passive before entering the cluster for the first time ($12\%$) are typically massive and are not affected significantly by the cluster. Using the bulge/total mass ratio as a proxy for morphology, we find that a single incursion through a cluster do not result in significant morphological changes in all four types.

ABSTRACT In this work, we combine the semi-analytic model of galaxy formation and evolution sag with the 102 relaxed simulated galaxy clusters from the three hundred project, and we study the link between the quenching of star formation (SF) and the physical processes that galaxies experience through their dynamical history in and around clusters. We classify galaxies in four populations based on their orbital history: recent and ancient infallers, and backsplash and neighbouring galaxies. We find that ∼85 per cent of the current population of quenched galaxies located inside the clusters are ancient infallers with low or null content of hot and cold gas. The fraction of quenched ancient infallers increases strongly between the first and second pericentric passage, due to the removal of hot gas by the action of ram-pressure stripping (RPS). The majority of them quenches after the first pericentric passage, but a non-negligible fraction needs a second passage, specially galaxies with $M_\star \le 10^{10.5} \, {\rm M_\odot }$. Recent infallers represent ∼15 per cent of the quenched galaxies located inside the cluster and, on average, they contain a high proportion of hot and cold gas; moreover, pre-processing effects are responsible for quenching the recent infallers prior to infall onto the main cluster progenitor. The ∼65 per cent of quenched galaxies located around clusters are backsplash galaxies, for which the combination of RPS acting during a pre-processing stage and inside the cluster is necessary for the suppression of SF in this population.

ABSTRACT We analyse the connection between the star formation quenching of galaxies and their location in the outskirts of clusters in the redshift range z = [0, 2] by estimating the fraction of red galaxies. More specifically, we focus on galaxies that infall isotropically from those that are infalling alongside filaments. We use a sample of galaxies obtained from the semi-analytical model of galaxy formation sag applied to the MultiDark simulation, mdpl2. In agreement with observational results, we find that the infall regions show levels of star formation that are intermediate between those of galaxies in clusters and in the field. Moreover, we show that, in the redshift range [0–0.85], the quenching of the star formation is stronger in the filamentary region than in the isotropic infall region. We also study the fraction of red galaxies as a function of the normalized distance to the cluster centre and find that, for radii R/R200 > 3, the fraction of red galaxies in the filamentary region is considerably larger than in the isotropic infall region. From the analysis of properties of the main progenitors of galaxies identified at z = 0, we find that they have different evolutionary behaviours depending on the stellar mass and environment. Our results confirm the observational findings that suggest that the infall regions of clusters play an important role in the pre-processing of galaxies along most of the evolutionary history of galaxies.

We connect galaxy properties with their orbital classification by analysing a sample of galaxies with stellar mass $M_{\star} \geq 10^{8.5}h^{-1}M_\odot$ residing in and around massive and isolated galaxy clusters with mass $M_{200} > 10^{15}h^{-1}M_\odot$ at redshift $z=0$. The galaxy population is generated by applying the semi-analytic model of galaxy formation SAG on the cosmological simulation MultiDark Planck 2. We classify galaxies considering their real orbits (3D) and their projected phase-space position using the ROGER code (2D). We define five categories: cluster galaxies, galaxies that have recently fallen into a cluster, backsplash galaxies, infalling galaxies, and interloper galaxies. For each class, we analyse the $g-r$ colour, the specific star formation rate (sSFR), and the stellar age, as a function of the stellar mass. For the 3D classes, we find that cluster galaxies have the lowest sSFR, and are the reddest and the oldest, as expected from environmental effects. Backsplash galaxies have properties intermediate between the cluster and recent infaller galaxies. For each 2D class, we find an important contamination by other classes. We find it necessary to separate the galaxy populations in red and blue to perform a more realistic analysis of the 2D data. For the red population, the 2D results are in good agreement with the 3D predictions. Nevertheless, when the blue population is considered, the 2D analysis only provides reliable results for recent infallers, infalling galaxies and interloper galaxies.

ABSTRACT We present a model to track the orbital evolution of ‘unresolved subhaloes’ (USHs) in cosmological simulations. USHs are subhaloes that are no longer distinguished by halo finders as self-bound overdensities within their larger host system due to limited mass resolution. These subhaloes would host ‘orphan galaxies’ in semi-analytic models of galaxy formation and evolution (SAMs). Predicting the evolution of the phase-space components of USHs is crucial for the adequate modelling of environmental processes, interactions, and mergers implemented in SAMs that affect the baryonic properties of orphan satellites. Our model takes into account dynamical friction drag, mass-loss by tidal stripping and merger with the host halo, involving three free parameters. To calibrate this model, we consider two DM-only simulations of different mass resolution (MultiDark simulations). The simulation with higher mass resolution (smdpl; $m_{\rm DM} = 9.6 \times 10^7\, h^{-1}\, \mathrm{M_{\odot }}$) provides information about subhaloes that are not resolved in the lower mass resolution one (mdpl2; $m_{\rm DM} = 1.5 \times 10^9\, h^{-1}\, \mathrm{M_{\odot }}$); the orbit of those USHs is tracked by our model. We use as constraining functions the subhalo mass function (SHMF) and the two-point correlation function (2PCF) obtained from smdpl, being the latter a novel aspect of our approach. While the SHMF fails to put tight constraints on the efficiency of dynamical friction and the merger condition, the addition of clustering information helps us to specify the parameters of the model related to the spatial distribution of subhaloes. Our model allows to achieve good convergence between the results of simulations of different mass resolution, with a precision better than 10 per cent for both SHMF and 2PCF.

We analyse the connection between the star formation quenching of galaxies and their location in theoutskirts of clusters in the redshift range $z=[0,2]$ by estimating the fraction of red galaxies. More specifically, we focus on galaxies that infall isotropically from those that are infalling alongside filaments. We use a sample of galaxies obtained from the semi-analytic model of galaxy formation SAG applied to the MultiDark simulation. {\textsc{mdpl2}}. In agreement with observational results, we find that the infall regions show levels of star formation that are intermediate between those of galaxies in clusters and in the field. Moreover, we show that, in the redshift range [0-0.85], the quenching of the star formation is stronger in the filamentary region than in the isotropic infall region. We also study the fraction of red galaxies as a function of the normalised distance to the cluster centre and find that, for radii $R/R_{200}> 3 $, the fraction of red galaxies in the filamentary region is considerably larger than in the isotropic infall region. From the analysis of properties of the main progenitors of galaxies identified at $z = 0$, we find that they have different evolutionary behaviours depending on the stellar mass and environment. Our results confirm the observational findings that suggest that the infall regions of clusters play an important role in the pre-processing of galaxies along most of the evolutionary history of galaxies.

We present a new analytic fitting profile to model the ram pressure (RP) exerted over satellite galaxies on different environments and epochs. The profile is built using the information of the gas particle distribution in hydrodynamical simulations of groups and clusters of galaxies to measure the RP directly. We show that predictions obtained by a previously introduced beta-profile model cannot consistently reproduce the dependence of the RP on halocentric distance and redshift for a given halo mass. It features a systematic underestimation of the predicted RP at high redshifts (z > 1.5), which increases towards the central regions of the haloes and it is independent of halo mass, reaching differences larger than two decades for satellites at r < 0.4R(vir). This behaviour reverses as redshift decreases, featuring an increasing overestimation with halocentric distance at z = 0. As an alternative, we introduce a new universal analytic model for the profiles that can recover the RP dependence on halo mass, halocentric distance, and redshift. We analyse the impact of our new profile on galaxy properties by applying a semi-analytic model of galaxy formation and evolution on top of the simulations. We show that galaxies experiencing large amounts of cumulative RP stripping typically have low stellar masses (M-* <= 10(9.5 )M(circle dot)). Besides, their specific star formation histories depend on the RP modelling applied, particularly at high redshifts (z > 1.5).