The study of the spatially resolved Star Formation Rate-Mass (Sigma_SFR-Sigma_M) relation gives important insights on how galaxies assemble at different spatial scales. Here we present the analysis of the Sigma_SFR-Sigma_M of 40 local cluster galaxies undergoing ram pressure stripping drawn from the GAs Stripping Phenomena in galaxies (GASP) sample. Considering their integrated properties, these galaxies show a SFR enhancement with respect to undisturbed galaxies of similar stellar mass; we now exploit spatially resolved data to investigate the origin and location of the excess. Even on ~1kpc scales, stripping galaxies present a systematic enhancement of Sigma_SFR (~0.35 dex at Sigma_M =108^M_sun/kpc^2) at any given Sigma_M compared to their undisturbed counterparts. The excess is independent on the degree of stripping and of the amount of star formation in the tails and it is visible at all galactocentric distances within the disks, suggesting that the star formation is most likely induced by compression waves from ram pressure. Such excess is larger for less massive galaxies and decreases with increasing mass. As stripping galaxies are characterised by ionised gas beyond the stellar disk, we also investigate the properties of 411 star forming clumps found in the galaxy tails. At any given stellar mass density, these clumps are systematically forming stars at a higher rate than in the disk, but differences are reconciled when we just consider the mass formed in the last few 10^8yr ago, suggesting that on these timescales the local mode of star formation is similar in the tails and in the disks.

Using MUSE observations from the GASP survey, we study 54 galaxies undergoing ram-pressure stripping (RPS) spanning a wide range in galaxy mass and host cluster mass. We use this rich sample to study how the star formation rate (SFR) in the tails of stripped gas depends on the properties of the galaxy and its host cluster. We show that the interplay between all the parameters involved is complex and that there is not a single, dominant one in shaping the observed amount of SFR. Hence, we develop a simple analytical approach to describe the mass fraction of stripped gas and the SFR in the tail, as a function of the cluster velocity dispersion, galaxy stellar mass, clustercentric distance and speed in the intracluster medium. Our model provides a good description of the observed gas truncation radius and of the fraction of star-formation rate (SFR) observed in the stripped tails, once we take into account the fact that the star formation efficiency in the tails is a factor $\sim 5$ lower than in the galaxy disc, in agreement with GASP ongoing HI and CO observations. We finally estimate the contribution of RPS to the intracluster light (ICL) and find that the average SFR in the tails of ram-pressure stripped gas is $\sim 0.22 M_{\odot}\,\mathrm{yr}^{-1} $ per cluster. By extrapolating this result to evaluate the contribution to the ICL at different epochs, we compute an integrated average value per cluster of $\sim 4 \times 10^9 M_\odot$ of stars formed in the tails of RPS galaxies since $z\sim 1$.

We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.

Feedback is indispensable in galaxy formation. However, lacking resolutions, cosmological simulations often use ad hoc feedback parameters. Conversely, small-box simulations, while better resolving the feedback, cannot capture gas evolution beyond the simulation domain. We aim to bridge the gap by implementing small-box results of supernovae-driven outflows into dark matter halo-scale simulations and studying their impact on large scales. Galactic outflows are multiphase, but small-box simulations show that the hot phase (T$\approx$ 10$^{6-7}$ K) carries the majority of energy and metals. We implement hot outflows in idealized simulations of the Milky Way halo, and examine how they impact the circumgalactic medium (CGM). In this paper, we discuss the case when the star formation surface density is low and therefore the emerging hot outflows are gravitationally bound by the halo. We find that outflows form a large-scale, metal-enriched atmosphere with fountain motions. As hot gas accumulates, the inner atmosphere becomes "saturated". Cool gas condenses, with a rate balancing the injection of the hot outflows. This balance leads to a universal density profile of the hot atmosphere, independent of mass outflow rate. The atmosphere has a radially-decreasing temperature, naturally producing the observed X-ray luminosity and column densities of O VI, O VII, O VIII. The self-regulated atmosphere has a baryon and a metal mass of $(0.5-1.2)\times 10^{10}M_\odot$ and $(0.6-1.4)\times 10^8 M_\odot$, respectively, small compared to the ``missing" baryons and metals from the halo. We conjecture that the missing materials reside at even larger radii, ejected by more powerful outflows in the past.

We present HI observations of the jellyfish galaxy, JO201. This massive galaxy (M$_{\ast} = 3.5 \times 10^{10}$ M$_\odot$) is falling along the line-of-sight towards the centre of a rich cluster (M$_{200} \sim 1.6 \times 10^{15}$ M$_\odot$, $\sigma_{cl} \sim 982$ km/s) at a high velocity $\geq$3363 km/s. Its H$\alpha$ emission shows a $\sim$40 kpc tail confined closely to its stellar disc and a $\sim$100 kpc tail extending further out. We find HI emission coinciding only with the shorter clumpy H$\alpha$ tail. In total, we measure an HI mass of M$_{\rm HI} = 1.65 \times 10^{9}$ M$_\odot$, which is about 60% lower than expected based on its stellar mass and stellar surface density. We compared JO201 to another jellyfish in the GASP sample, JO206 (of similar mass but residing in a 10$\times$ less massive cluster), and find that they are similarly HI-deficient. Of the total HI mass in JO201, about 30% lies outside the galaxy disc in projection. This HI fraction is probably a lower limit since most of the HI is redshifted relative to the stellar disc and could be outside the disc. The global star formation rate (SFR) analysis of JO201 suggests that its observed SFR would be expected if it had 10$\times$ its current HI mass. The disc is the main contributor of the high star formation efficiency at a given HI gas density for both galaxies, but their tails also show higher star formation efficiencies compared to the outer regions of field galaxies. Generally, we find that JO201 and JO206 are similar based on their HI content, stellar mass and star formation rate. This finding is unexpected considering their different environments. A toy model comparing the ram pressure of the ICM versus the restoring forces of these galaxies suggests that the ram pressure strength exerted on them could be comparable if we consider their 3D orbital velocities and radial distances relative to the clusters.

We present JVLA-C observations of the HI gas in JO204, one of the most striking jellyfish galaxies from the GASP survey. JO204 is a massive galaxy in the low-mass cluster Abell 957 at z=0.04243. The HI map reveals an extended 90 kpc long ram-pressure stripped tail of neutral gas, stretching beyond the 30 kpc long ionized gas tail and pointing away from the cluster center. The HI mass seen in emission is (1.32 $ \pm 0.13) \times 10^{9} \rm M_{\odot}$, mostly located in the tail. The northern part of the galaxy disk has retained some HI gas, while the southern part has already been completely stripped and displaced into an extended unilateral tail. Comparing the distribution and kinematics of the neutral and ionized gas in the tail indicates a highly turbulent medium. Moreover, we observe associated HI absorption against the 11 mJy central radio continuum source with an estimated HI absorption column density of 3.2 $\times 10^{20}$ cm$^{-2}$. The absorption profile is significantly asymmetric with a wing towards higher velocities. We modelled the HI absorption by assuming that the HI and ionized gas disks have the same kinematics in front of the central continuum source, and deduced a wider absorption profile than observed. The observed asymmetric absorption profile can therefore be explained by a clumpy, rotating HI gas disk seen partially in front of the central continuum source, or by ram-pressure pushing the neutral gas towards the center of the continuum source, triggering the AGN activity.

We introduce a collection of primarily centrally star-forming galaxies that are selected by disk color to have truncated disk star formation. We show that common explanations for centrally-concentrated star formation -- low stellar mass, bars, and high-density environments, do not universally apply to this sample. To gain insight into our sample, we compare these galaxies to a parent sample of strongly star-forming galaxies and to a parent sample of galaxies with low specific star formation rates. We find that in star formation and color space from ultraviolet to the infrared these galaxies either fall between the two samples or agree more closely with galaxies with high-specific star formation rates. Their morphological characteristics also lie between high- and low-specific star formation rate galaxies, although their Petrosian radii agree well with that of the low-specific star formation rate parent sample. We discuss whether this sample is likely to be quenching or showing an unusual star-formation distribution while continuing to grow through star formation. Future detailed studies of these galaxies will give us insights into how the local conditions within a galaxy balance environmental influence to govern the distribution of star formation. In this first paper in a series, we describe the global properties that identify this sample as separate from more average spiral galaxies, and identify paths forward to explore the underlying causes of their differences.

With MUSE, Chandra, VLA, ALMA and UVIT data from the GASP programme we study the multiphase baryonic components in a jellyfish galaxy (JW100) with a stellar mass 3.2 X 10^{11} M_sun hosting an AGN. We present its spectacular extraplanar tails of ionized and molecular gas, UV stellar light, X-ray and radio continuum emission. This galaxy represents an excellent laboratory to study the interplay between different gas phases and star formation, and the influence of gas stripping, gas heating, and AGN. We analyze the physical origin of the emission at different wavelengths in the tail, in particular in-situ star formation (related to Halpha, CO and UV emission), synchrotron emission from relativistic electrons (producing the radio continuum) and heating of the stripped interstellar medium (ISM) (responsible for the X-ray emission). We show the similarities and differences of the spatial distributions of ionized gas, molecular gas and UV light, and argue that the mismatch on small scales (1kpc) is due to different stages of the star formation process. We present the relation Halpha--X-ray surface brightness, which is steeper for star-forming regions than for diffuse ionised gas regions with high [OI]/Halpha ratio. We propose that ISM heating due to interaction with the intracluster medium (either for mixing, thermal conduction or shocks) is responsible for the X-ray tail, the observed [OI]-excess and the lack of star formation in the northern part of the tail. We also report the tentative discovery in the tail of the most distant (and among the brightest) currently known ULX, a point-like ultraluminous X-ray source commonly originating in a binary stellar system powered either by an intermediate-mass black hole or a magnetized neutron star.

Exploiting the sample of 30 local star-forming, undisturbed late-type galaxies in different environments drawn from the GAs Stripping Phenomena in galaxies with MUSE (GASP), we investigate the spatially resolved Star Formation Rate-Mass ({\Sigma}SFR-{\Sigma}_star) relation. Our analysis includes also the galaxy outskirts (up to >4 effective radii, re), a regime poorly explored by other Integral Field Spectrograph surveys. Our observational strategy allows us to detect H{\alpha} out to more than 2.7re for 75% of the sample. Considering all galaxies together, the correlation between the {\Sigma}SFR and {\Sigma}_star is quite broad, with a scatter of 0.3 dex. It gets steeper and shifts to higher {\Sigma}_star values when external spaxels are excluded and moving from less to more massive galaxies. The broadness of the overall relation suggests galaxy-by-galaxy variations. Indeed, each object is characterized by a distinct {\Sigma}SFR-{\Sigma}_star relation and in some cases the correlation is very loose. The scatter of the relation mainly arises from the existence of bright off-center star-forming knots whose {\Sigma}SFR-{\Sigma}_star relation is systematically broader than that of the diffuse component. The {\Sigma}SFR-{\Sigma}tot gas (total gas surface density) relation is as broad as the {\Sigma}SFR-{\Sigma}_star relation, indicating that the surface gas density is not a primary driver of the relation. Even though a large galaxy-by-galaxy variation exists, mean {\Sigma}SFR and {\Sigma}_star values vary of at most 0.7 dex across galaxies. We investigate the relationship between the local and global SFR-M_star relation, finding that the latter is driven by the existence of the size-mass relation.

We present VLA HI observations of JO206, a prototypical ram-pressure stripped galaxy in the GASP sample. This massive galaxy (M$_{\ast} =$ 8.5 $\times$ 10$^{10}$ M$_{\odot}$) is located at a redshift of $z =$ 0.0513, near the centre of the low-mass galaxy cluster, IIZw108 ($\sigma \sim575$ km/s). JO206 is characterised by a long tail ($\geq$90 kpc) of ionised gas stripped away by ram-pressure. We find a similarly long HI tail in the same direction as the ionised gas tail and measure a total HI mass of $3.2 \times 10^{9}$ M$_{\odot}$. This is about half the expected HI mass given the stellar mass and surface density of JO206. A total of $1.8 \times 10^{9}$ M$_{\odot}$ (60%) of the detected HI is in the gas stripped tail. An analysis of the star formation rate shows that the galaxy is forming more stars compared to galaxies with the same stellar and HI mass. On average we find a HI gas depletion time of $\sim$0.5 Gyr which is about four times shorter than that of "normal" spiral galaxies. We performed a spatially resolved analysis of the relation between star formation rate density and gas density in the disc and tail of the galaxy at the resolution of our HI data. The star formation efficiency of the disc is about 10 times higher than that of the tail at fixed HI surface densities. Both the inner and outer parts of JO206 show an enhanced star formation compared to regions of similar HI surface density in field galaxies. The enhanced star formation is due to ram-pressure stripping during the galaxy's first infall into the cluster.

We investigate the importance of varying the ram pressure to more realistically mimic the infall of a cluster satellite galaxy when comparing ram pressure stripping simulations to observations. We examine the gas disk and tail properties of stripped cluster galaxies in eight "wind-tunnel" hydrodynamical simulations with either varying or constant ram pressure strength. In simulations without radiative cooling, applying a varying wind leads to significantly different density and velocity structure in the tail than found when applying a constant wind, although the stripping rate, disk mass, and disk radius remain consistent in both scenarios. In simulations with radiative cooling, the differences between a constant and varying wind are even more pronounced. Not only is there a difference in morphology and velocity structure in the tails, but a varying wind leads to a much lower stripping rate, even after the varying wind has reached the ram pressure strength of the constant wind. Also, galaxies in constant and varying wind simulations with the same gas disk mass do not have in the same gas disk radius. A constant wind cannot appropriately model the ram pressure stripping of a galaxy entering a cluster. We conclude that simulations attempting detailed comparisons with observations must take the variation of the ram pressure profile due to a galaxy's orbit into consideration.

In spatially resolved galaxy observations, star formation rate radial profiles are found to correlate with total specific star formation rates. A central depletion in star formation is thought to correlate with the globally depressed star formation rates of, for example, galaxies within the Green Valley. We present, for the first time, radial specific star formation rate profiles for a statistical sample of simulated galaxies from the Illustris and EAGLE large cosmological simulations. For galaxies on the star-forming sequence, simulated specific star formation rate profiles are in loose agreement with observations, although galaxies from the EAGLE simulation are too centrally peaked, and galaxies from Illustris have a steeper decline at large radii. However, both galaxy samples show centrally concentrated star formation for galaxies in the Green Valley at all galaxy stellar masses, indicating that quenching occurs from the outside-in, in strong conflict with observations of inside-out quenching. These results appear in spite of the different feedback models in these two simulations. We conclude that the distribution of star formation within galaxies is a strong additional constraint for star formation and feedback models in simulations, in particular related to the quenching of star formation.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

We present a study of the physical properties of JO201, a unique disk galaxy with extended tails undergoing extreme ram-pressure stripping as it moves through the massive cluster Abell 85 at supersonic speeds mostly along the line of sight. JO201 was observed with MUSE as part of the GASP programme. In a previous paper (GASP II) we studied the stellar and gas kinematics. In this paper we present emission-line ratios, gas-phase metallicities and ages of the stellar populations across the galaxy disk and tails. We find that while the emission at the core of the galaxy is dominated by an active galactic nucleus (AGN), the disk is composed of star-forming knots surrounded by excited diffuse gas. The collection of star-forming knots presents a metallicity gradient steadily decreasing from the centre of the galaxy outwards, and the ages of the stars across the galaxy show that the tails formed <10^9 yr ago. This result is consistent with an estimate of the stripping timescale (1 Gyr), obtained from a toy orbital model. Overall, our results independently and consistently support a scenario in which a recent or ongoing event of intense ram-pressure stripping acting from the outer disk inwards, causes removal and compression of gas, thus altering the AGN and star-formation activity within and around the galaxy.

Based on MUSE data from the GASP survey, we study the Halpha-emitting extraplanar tails of 16 cluster galaxies at z~0.05 undergoing ram pressure stripping. We demonstrate that the dominating ionization mechanism of this gas (between 64% and 94% of the Halpha emission in the tails depending on the diagnostic diagram used) is photoionization by young massive stars due to ongoing star formation (SF) taking place in the stripped tails. This SF occurs in dynamically quite cold HII clumps with a median Halpha velocity dispersion sigma = 27 km s^-1. We study the characteristics of over 500 star-forming clumps in the tails and find median values of Halpha luminosity L_{Halpha} = 4 X 10^38 erg s^-1, dust extinction A_V=0.5 mag, star formation rate SFR=0.003 M_sun yr^-1, ionized gas density n_e =52 cm^-3, ionized gas mass M_gas = 4 X 10^4 Msun, and stellar mass M_{*} = 3 X 10^6 Msun. The tail clumps follow scaling relations (M_gas-M_{*}, L_{Halpha} -sigma, SFR-M_gas) similar to disk clumps, and their stellar masses are comparable to Ultra Compact Dwarfs and Globular Clusters.The diffuse gas component in the tails is ionized by a combination of SF and composite/LINER-like emission likely due to thermal conduction or turbulence. The stellar photoionization component of the diffuse gas can be due either to leakage of ionizing photons from the HII clumps with an average escape fraction of 18%, or lower luminosity HII regions that we cannot individually identify.

Exploiting the data from the GAs Stripping Phenomena in galaxies with MUSE (GASP) program, we compare the integrated Star Formation Rate- Mass relation (SFR-M_ast) relation of 42 cluster galaxies undergoing ram pressure stripping ("stripping galaxies") to that of 32 field and cluster undisturbed galaxies. Theoretical predictions have so far led to contradictory conclusions about whether ram pressure can enhance the star formation in the gas disks and tails or not and until now a statistically significant observed sample of stripping galaxies was lacking. We find that stripping galaxies occupy the upper envelope of the control sample SFR-M_ast relation, showing a systematic enhancement of the SFR at any given mass. The star formation enhancement occurs in the disk (0.2 dex), and additional star formation takes place in the tails. Our results suggest that strong ram pressure stripping events can moderately enhance the star formation also in the disk prior to gas removal.

GAs Stripping Phenomena in galaxies with MUSE (GASP) is a program aimed at studying gas removal processes in nearby galaxies in different environments. We present the study of four galaxies that are part of the same group (z= 0.06359) and highlight the multitude of mechanisms affecting the spatially resolved properties of the group members. One galaxy is passive and shows a regular stellar kinematics. The analysis of its star formation history indicates that the quenching process lasted for a few Gyr and that the star formation declined throughout the disk in a similar way, consistent with strangulation. Another galaxy is characterised by a two-component stellar disk with an extended gas disk that formed a few 10^8 yr ago, most likely as a consequence of gas accretion. The third member is a spiral galaxy at the edges of the group, but embedded in a filament. We hypothesise that the compression exerted by the sparse intergalactic medium on the dense circumgalactic gas switches on star formation in a number of clouds surrounding the galaxy ("cosmic web enhancement"). Alternatively, also ram pressure stripping might be effective. Finally, the fourth galaxy is a spiral with a truncated ionised gas disk and an undisturbed stellar kinematics. An analytical model of the galaxy's restoring pressure, and its location and velocity within the cluster, suggest ram pressure is the most likely physical mechanism in action. This is the first optical evidence for stripping in groups.

Recent hydrodynamical simulations predict that stellar feedback in intermediate mass galaxies (IMGs) can drive strong fluctuations in structure (e.g., half-light radius, $R_e$). This process operates on timescales of only a few hundred Myr and persists even at late cosmic times. One prediction of this quasi-periodic, galactic-scale "breathing" is an $anti$-correlation between star formation rate (SFR) and half-light radius as central gas overdensities lead to starbursts whose feedback drags stars to larger radii while star formation dwindles. We test this prediction with a sample of 322 $isolated$ IMGs with stellar masses of $10^{9.0} \leq M/M_{\odot} \leq 10^{9.5}$ at $0.3<z<0.4$ in the HST $I_{814}$ COSMOS footprint. We find that IMGs with higher specific SFRs (SSFR $>10^{-10}$ yr$^{-1}$) are the most extended with median sizes of $R_e \sim 3-3.4$ kpc and are mostly disk-dominated systems. In contrast, IMGs with lower SSFRs ($<10^{-10}$ yr$^{-1}$) are a factor of $\sim 2-3$ more compact with median sizes of $R_e \sim 0.9-1.6$ kpc and have more significant bulge contributions to their light. These observed trends are opposite the predictions for stellar feedback that operate via the "breathing" process described above. We discuss various paths to reconcile the observations and simulations, all of which likely require a different implementation of stellar feedback in IMGs that drastically changes their predicted formation history.

It is well known that galaxies falling into clusters can experience gas stripping due to ram-pressure by the intra-cluster medium (ICM). The most spectacular examples are galaxies with extended tails of optically-bright stripped material known as "jellyfish". We use the first large homogeneous compilation of jellyfish galaxies in clusters from the WINGS and OmegaWINGS surveys, and follow-up MUSE observations from the GASP MUSE programme to investigate the orbital histories of jellyfish galaxies in clusters and reconstruct their stripping history through position vs. velocity phase- space diagrams. We construct analytic models to define the regions in phase-space where ram-pressure stripping is at play. We then study the distribution of cluster galaxies in phase-space and find that jellyfish galaxies have on average higher peculiar velocities (and higher cluster velocity dispersion) than the overall population of cluster galaxies at all clustercentric radii, which is indicative of recent infall into the cluster and radial orbits. In particular, the jellyfish galaxies with the longest gas tails reside very near the cluster cores (in projection) and are moving at very high speeds, which coincides with the conditions of the most intense ram-pressure. We conclude that many of the jellyfish galaxies seen in clusters likely formed via fast (~1- 2 Gyr), incremental, outside-in ram-pressure stripping during first infall into the cluster in highly radial orbits.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift.

When supermassive black holes at the center of galaxies accrete matter (usually gas), they give rise to highly energetic phenomena named Active Galactic Nuclei (AGN). A number of physical processes have been proposed to account for the funneling of gas towards the galaxy centers to feed the AGN. There are also several physical processes that can strip gas from a galaxy, and one of them is ram pressure stripping in galaxy clusters due to the hot and dense gas filling the space between galaxies. We report the discovery of a strong connection between severe ram pressure stripping and the presence of AGN activity. Searching in galaxy clusters at low redshift, we have selected the most extreme examples of jellyfish galaxies, which are galaxies with long tentacles of material extending for dozens of kpc beyond the galaxy disk. Using the MUSE spectrograph on the ESO Very Large Telescope, we find that 6 out of the 7 galaxies of this sample host a central AGN, and two of them also have galactic-scale AGN ionization cones. The high incidence of AGN among the most striking jellyfishes may be due to ram pressure causing gas to flow towards the center and triggering the AGN activity, or to an enhancement of the stripping caused by AGN energy injection, or both. Our analysis of the galaxy position and velocity relative to the cluster strongly supports the first hypothesis, and puts forward ram pressure as another, yet unforeseen, possible mechanism for feeding the central supermassive black hole with gas.

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the Milky Way. Using Markov-Chain Monte-Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass--metallicity relation of the Milky Way satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass--metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g. "pre-heating") is needed to explain the low stellar mass fraction for a given subhalo mass.

We examine the statistics of the low-redshift Ly-alpha forest in an adaptive mesh refinement hydrodynamic cosmological simulation of sufficient volume to include distinct large-scale environments. We compare our HI column density distribution of absorbers both with recent work and between two highly-refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on scales of approximately 20 Mpc). We recover the average results presented in Kollmeier et al. (2014) using different simulation methods. We further break down these results as a function of environment to examine the detailed dependence of absorber statistics on large-scale density. We find that the slope of the HI column density distribution in the 10$^{12.5}$ $\le$ N$_{HI}$/cm$^{-2}$ $\le$ 10$^{14.5}$ range depends on environment such that the slope becomes steeper for higher environmental density, and this difference reflects distinct physical conditions of the intergalactic medium on these scales. We track this difference to the different temperature structures of filaments in varying environments. Specifically, filaments in the overdensity are hotter and, correspondingly, are composed of gas with lower HI fractions than those in underdense environments. Our results highlight that in order to understand the physics driving the HI CDD, we need not only improved accounting of the sources of ionizing UV photons, but also of the physical conditions of the IGM and how this may vary as a function of large-scale environment.

We present 12CO (2-1) data of three Virgo spirals - NGC 4330, NGC 4402 and NGC 4522 obtained using the Submillimeter Array. These three galaxies show clear evidence of ram pressure stripping due to the cluster medium as found in previous HI imaging studies. Using high-resolution CO data, we investigate how the properties of the inner molecular gas disc change while a galaxy is undergoing HI stripping in the cluster. At given sensitivity limits, we do not find any clear signs of molecular gas stripping. However, both its morphology and kinematics appear to be quite disturbed as those of HI. Morphological peculiarities present in the molecular and atomic gas are closely related with each other, suggesting that molecular gas can be also affected by strong ICM pressure even if it is not stripped. CO is found to be modestly enhanced along the upstream sides in these galaxies, which may change the local star formation activity in the disc. Indeed, the distribution of H$\alpha$ emission, a tracer of recent star formation, well coincides with that of the molecular gas, revealing enhancements near the local CO peak or along the CO compression. FUV and H$\alpha$ share some properties in common, but FUV is always more extended than CO/H$\alpha$ in the three galaxies, implying that the star-forming disc is rapidly shrinking as the molecular gas properties have changed. We discuss how ICM pressure affects dense molecular gas and hence star formation properties while diffuse atomic gas is being removed from a galaxy.

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

The connection between dark matter halos and galactic baryons is often not well-constrained nor well-resolved in cosmological hydrodynamical simulations. Thus, Halo Occupation Distribution (HOD) models that assign galaxies to halos based on halo mass are frequently used to interpret clustering observations, even though it is well-known that the assembly history of dark matter halos is related to their clustering. In this paper we use high-resolution hydrodynamical cosmological simulations to compare the halo and stellar mass growth of galaxies in a large-scale overdensity to those in a large-scale underdensity (on scales of about 20 Mpc). The simulation reproduces assembly bias, that halos have earlier formation times in overdense environments than in underdense regions. We find that the stellar mass to halo mass ratio is larger in overdense regions in central galaxies residing in halos with masses between 10$^{11}$-10$^{12.9}$ M$_{\odot}$. When we force the local density (within 2 Mpc) at z=0 to be the same for galaxies in the large-scale over- and underdensities, we find the same results. We posit that this difference can be explained by a combination of earlier formation times, more interactions at early times with neighbors, and more filaments feeding galaxies in overdense regions. This result puts the standard practice of assigning stellar mass to halos based only on their mass, rather than considering their larger environment, into question.

Recent observations have constrained the orbit and structure of the Large Magellanic Cloud (LMC), implying a well-constrained pericentric passage about the Milky Way (MW) ~ 50 Myr ago. In this scenario, the LMC's gaseous disk has recently experienced maximal ram pressure stripping, suggesting the current extent of its HI disk directly probes the medium in which it is moving. From the observed stellar and HI distributions of the system we find evidence of a truncated gas profile along the windward ``leading edge' of the LMC disk, despite a far more extended stellar component. We explore the implications of this ram pressure stripping signature, using both analytic prescriptions and full 3-dimensional hydrodynamic simulations of the LMC. Our simulations subject the system to a headwind whose velocity components correspond directly to the recent orbital history of the LMC. We vary the density of this headwind, using a variety of sampled parameters for a Beta-profile for a theoretical MW circumgalactic medium (CGM), comparing the resulting HI morphology directly to observations of the LMC HI and stellar components. This model can match the radial extent of the LMC's leading (windward) edge only in scenarios where the MW CGM density at pericentric passage is n(R = 48.2 +/- 5 kpc) = 1.1 (+.44/-.45) x 1e-4 cm^-3. The implied pericentric density proves insensitive to both the broader CGM structure and temperature profile, thus providing a model-independent constraint on the local gas density. This result imposes an important constraint on the density profile of the MW's CGM, and thus the total baryon content of the MW. From our work, assuming a Beta-profile valid to ~ Rvir, we infer a total diffuse CGM mass M(300 kpc) = 2.6 +/- 1.4 x 1e10 Msun or approximately 15% of a 1e12 Msun MW's baryonic mass budget.

We present a catalog of 1964 isolated, compact neutral hydrogen clouds from the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI DR1). The clouds were identified by a custom machine-vision algorithm utilizing Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds have velocities typically between |VLSR| = 20-400 km/s, linewidths of 2.5-35 km/s, and column densities ranging from 1 - 35 x 10^18 cm^-2. The distances to the clouds in this catalog may cover several orders of magnitude, so the masses may range from less than a Solar mass for clouds within the Galactic disc, to greater than 10^4 Solar Masses for HVCs at the tip of the Magellanic Stream. To search for trends, we separate the catalog into five populations based on position, velocity, and linewidth: high velocity clouds (HVCs); galaxy candidates; cold low velocity clouds (LVCs); warm, low positive-velocity clouds in the third Galactic Quadrant; and the remaining warm LVCs. The observed HVCs are found to be associated with previously-identified HVC complexes. We do not observe a large population of isolated clouds at high velocities as some models predict. We see evidence for distinct histories at low velocities in detecting populations of clouds corotating with the Galactic disc and a set of clouds that is not corotating.

One process affecting gas-rich cluster galaxies is ram pressure stripping, i.e. the removal of galactic gas through direct interaction with the intracluster medium. Galactic magnetic fields may have an important impact on the stripping rate and tail structure. We run the first magnetohydrodynamic simulations of ram pressure stripping that include a galactic magnetic field, using 159 pc resolution throughout our entire domain in order to resolve mixing throughout the tail. We find very little difference in the total amount of gas removed from the unmagnetized and magnetized galaxies, although a magnetic field with a radial component will initially accelerate stripped gas more quickly. In general, we find that magnetic fields in the disk lead to slower velocities in the stripped gas near the disk and faster velocities farther from the disk. We also find that magnetic fields in the galactic gas lead to larger unmixed structures in the tail. Finally, we discuss whether ram pressure stripped tails can magnetize the ICM. We find that the total magnetic energy density grows as the tail lengthens, likely through turbulence. There are microGauss-strength fields in the tail in all of our MHD runs, which survive to at least 100 kpc from the disk (the edge of our simulated region), indicating that the area-filling factor of magnetized tails in a cluster could be large.

Recent large surveys have found a reversal of the star formation rate (SFR)-density relation at z=1 from that at z=0 (e.g. Elbaz et al.; Cooper et al.), while the sign of the slope of the color-density relation remains unchanged (e.g. Cucciati et al.; Quadri et al.). We use state-of-the-art adaptive mesh refinement cosmological hydrodynamic simulations of a 21x24x20 (Mpc/h)$^3$ region centered on a cluster to examine the SFR-density and color-density relations of galaxies at z=0 and z=1. The local environmental density is defined by the dark matter mass in spheres of radius 1 Mpc/h, and we probe two decades of environmental densities. Our simulations produce a large increase of SFR with density at z=1, as in the observations of Elbaz et al. We also find a significant evolution to z=0, where the SFR-density relation is much flatter. The color-density relation in our simulations is consistent from z=1 to z=0, in agreement with observations. We find that the increase in the median SFR with local density at z=1 is due to a growing population of star-forming galaxies in higher-density environments. At z=0 and z=1 both the SFR and cold gas mass are tightly correlated with the galaxy halo mass, and therefore the correlation between median halo mass and local density is an important cause of the SFR-density relation at both redshifts. We also show that the local density on 1 Mpc/h scales affects galaxy SFRs as much as halo mass at z=0. Finally, we find indications that the role of the 1 Mpc/h scale environment reverses from z=0 to z=1: at z=0 high-density environments depress galaxy SFRs, while at z=1 high-density environments tend to increase SFRs.

In a hierarchical universe, mergers may be an important mechanism not only in increasing the mass of galaxies but also in driving the color and morphological evolution of galaxies. We use a large sample of ~1000 simulated galaxies of stellar mass greater than 10^9.6 solar masses (for ~4800 observations at multiple redshifts) from a high-res (0.46 h^{-1} kpc) cosmological simulation to determine under what circumstances being a member of a pair influences galaxy properties at z <= 0.2. We identify gravitationally bound pairs, and find a relative fraction of blue-blue, red-red, and blue-red pairs that agrees with observations (Lin et al. 2010). Pairs tend to avoid the extreme environments of clusters and void centres. While pairs in groups can include galaxies that are both blue, both red, or one of each color, in the field it is rare for pair galaxies to both be red. We find that physically bound pairs closer than 250 h^{-1} kpc tend to have higher sSFRs than the galaxy population as a whole. However, the sSFR of a bound galaxy relative to galaxies in a comparable local density environment (determined by the distance to the fifth nearest neighbor, rho_5), differs depending on the local density. In regions of high rho_5 the bound population has a higher fraction of star-forming (bluer) galaxies, whereas there is very little difference between bound and unbound galaxies in low rho_5 regions. This effect on the star-forming fraction may be driven by the higher fraction of bound HI-rich galaxies compared to unbound galaxies, particularly at high local densities. It appears that being in a pair has an incremental, but not overwhelming, effect on the star formation rate of the paired galaxies, compared to the more pronounced trend where galaxies overall have low sSFR (are red) in clusters and higher sSFR (blue) at the centre of voids. This trend depends most strongly on rho_5.(abridged)

We investigate the impact of star formation and feedback on ram pressure stripping using high-resolution adaptive mesh simulations, building on a previous series of papers that systematically investigated stripping using a realistic model for the interstellar medium, but without star formation. We find that star formation does not significantly affect the rate at which stripping occurs, and only has a slight impact on the density and temperature distribution of the stripped gas, indicating that our previous (gas-only) results are unaffected. For our chosen (moderate) ram pressure strength, stripping acts to truncate star formation in the disk over a few hundred million years, and does not lead to a burst of star formation. Star formation in the bulge is slightly enhanced, but the resulting change in the bulge-to-disk ratio is insignificant. We find that stars do form in the tail, primarily from gas that is ablated from the disk and the cools and condenses in the turbulent wake. The star formation rate in the tail is low, and any contribution to the intracluster light is likely to be very small. We argue that star formation in the tail depends primarily on the pressure in the intracluster medium, rather than the ram pressure strength. Finally, we compare to observations of star formation in stripped tails, finding that many of the discrepancies between our simulation and observed wakes can be accounted for by different intracluster medium pressures.

Some tails of ram-pressure stripped galaxies are detected in HI, some in Halpha, and some in X-ray (but never all three so far). We use numerical simulations to probe the conditions for the production of X-ray bright tails, demonstrating that the primary requirement is a high pressure intracluster medium (ICM). This is because the stripped tail is mostly in pressure equilibrium with the ICM, but mixing leaves it with densities and temperatures intermediate between the cold gas in the disk and the hot ICM. Given a high enough ICM pressure, this mixed gas lies in the X-ray bright region of the phase diagram. We compare the simulations to observations of the ram pressure stripped tail of ESO 137-001, showing excellent agreement in the total measured X-ray and Halpha emission and non-flaring morphology of the tail, and consistent HI measurements. Using these comparisons we constrain the level of mixing and efficiency of heat conduction in the intracluster medium (ICM)

Galaxies moving through the intracluster medium (ICM) of a cluster of galaxies can lose gas via ram pressure stripping. This stripped gas forms a tail behind the galaxy which is potentially observable. In this paper, we carry out hydrodynamical simulations of a galaxy undergoing stripping with a focus on the gas properties in the wake and their observational signatures. We include radiative cooling in an adaptive hydrocode in order to investigate the impact of a clumpy, multi-phase interstellar medium. We find that including cooling results in a very different morphology for the gas in the tail, with a much wider range of temperatures and densities. The tail is significantly narrower in runs with radiative cooling, in agreement with observed wakes. In addition, we make detailed predictions of H I, Halpha and X-ray emission for the wake, showing that we generally expect detectable H I and Halpha signatures, but no observable X-ray emission (at least for our chosen ram-pressure strength and ICM conditions). We find that the relative strength of the Halpha diagnostic depends somewhat on our adopted minimum temperature floor (below which we set cooling to zero to mimic physics processes not included in the simulation).

Cluster galaxies moving through the intracluster medium (ICM) are expected to lose some of their interstellar medium (ISM) through ISM-ICM interactions. We perform high resolution (40 pc) three-dimensional hydrodynamical simulations of a galaxy undergoing ram pressure stripping including radiative cooling in order to investigate stripping of a multiphase medium. The clumpy, multiphase ISM is self-consistently produced by the inclusion of radiative cooling, and spans six orders of magnitude in gas density. We find no large variations in the amount of gas lost whether or not cooling is involved, although the gas in the multiphase galaxy is stripped more quickly and to a smaller radius. We also see significant differences in the morphology of the stripped disks. This occurs because the multiphase medium naturally includes high density clouds set inside regions of lower density. We find that the lower density gas is stripped quickly from any radius of the galaxy, and the higher density gas can then be ablated. If high density clouds survive, through interaction with the ICM they lose enough angular momentum to drift towards the center of the galaxy where they are no longer stripped. Finally, we find that low ram pressure values compress gas into high density clouds that could lead to enhanced star formation, while high ram pressure leads to a smaller amount of high-density gas.

Cluster galaxies moving through the intracluster medium (ICM) are expected to lose some of their interstellar medium (ISM) through ram pressure stripping and related ISM-ICM interactions. Using high-resolution cosmological simulations of a large galaxy cluster including star formation, we show that the ram pressure a galaxy experiences at a fixed distance from the cluster center can vary by well over an order of magnitude We find that this variation in ram pressure is due in almost equal parts to variation in the ICM density and in the relative velocity between the galaxy and the ICM. We also find that the ICM and galaxy velocities are weakly correlated for in-falling galaxies.

Galaxies in clusters are gas-deficient and a number of possible explanations for this observation have been advanced, including galaxy-cluster tidal interactions, galaxy harassment, and ISM-ICM gas stripping. In this paper, we use a cosmological simulation of cluster formation and evolution in order to examine this issue from a theoretical standpoint. We follow a large number of galaxies over time and track each galaxy's gas and stellar mass changes to discover what mechanism(s) dominate the evolution of the cluster galaxies. We find that while gas is lost due to a wide variety of mechanisms, the most common way is via a gas-only stripping event, and the amount of gas lost correlates with the ram-pressure the galaxy is experiencing. Although this gas-stripping occurs primarily in the central region (r < 1 Mpc), it is an important mechanism out to the virial radius of the cluster. This is due to the wide scatter in ram-pressure strength that a galaxy experiences at fixed radius. We find that the timescale for complete gas removal is > 1 Gyr. In addition, we find that galaxies in the field and in the cluster periphery (r > 2.4 Mpc) often accrete cool gas; the accretion stops between 1-2.4 Mpc, possibly indicating the onset of galaxy starvation.

The robotic ROTSE-III telescope network detected prompt optical emission contemporaneous with the gamma-ray emission of Swift events GRB051109A and GRB051111. Both datasets have continuous coverage at high signal-to-noise levels from the prompt phase onwards, thus the early observations are readily compared to the Swift XRT and BAT high energy detections. In both cases, the optical afterglow is established, declining steadily during the prompt emission. For GRB051111, there is evidence of an excess optical component during the prompt emission. The component is consistent with the flux spectrally extrapolated from the gamma-rays, using the gamma-ray spectral index. A compilation of spectral information from previous prompt detections shows that such a component is unusual. The existence of two prompt optical components - one connected to the high-energy emission, the other to separate afterglow flux, as indicated in GRB051111 - is not compatible with a simple ``external-external'' shock model for the GRB and its afterglow.

We report on four Chandra grating observations of the oblique magnetic rotator theta^1 Ori C (O5.5 V) covering a wide range of viewing angles with respect to the star's 1060 G dipole magnetic field. We employ line-width and centroid analyses to study the dynamics of the X-ray emitting plasma in the circumstellar environment, as well as line-ratio diagnostics to constrain the spatial location, and global spectral modeling to constrain the temperature distribution and abundances of the very hot plasma. We investigate these diagnostics as a function of viewing angle and analyze them in conjunction with new MHD simulations of the magnetically channeled wind shock mechanism on theta^1 Ori C. This model fits all the data surprisingly well, predicting the temperature, luminosity, and occultation of the X-ray emitting plasma with rotation phase.

The launch of high-spectral-resolution x-ray telescopes (Chandra, XMM) has provided a host of new spectral line diagnostics for the astrophysics community. In this paper we discuss Doppler-broadened emission line profiles from highly supersonic outflows of massive stars. These outflows, or winds, are driven by radiation pressure and carry a tremendous amount of kinetic energy, which can be converted to x rays by shock-heating even a small fraction of the wind plasma. The unshocked, cold wind is a source of continuum opacity to the x rays generated in the shock-heated portion of the wind. Thus the emergent line profiles are affected by transport through a two-component, moving, optically thick medium. While complicated, the interactions among these physical effects can provide quantitative information about the spatial distribution and velocity of the x-ray-emitting and absorbing plasma in stellar winds. We present quantitative models of both a spherically-symmetric wind and a wind with hot plasma confined in an equatorial disk by a dipole magnetic field.