We present the first spatially resolved study of molecular gas in the vicinity of a Gamma Ray Burst, using CO(2-1) emission line observations with the Atacama Large Millimetre Array (ALMA) at ~50 pc scales. The host galaxy of GRB 980425 contains a ring of high column density HI gas which is likely to have formed due to a collision between the GRB host and its companion galaxy, within which the GRB is located. We detect eleven molecular gas clumps in the galaxy, seven of which are within the gas ring. The clump closest to the GRB position is at a projected separation of ~280 pc. Although it is plausible that the GRB progenitor was ejected from clusters formed in this clump, we argue that the in situ formation of the GRB progenitor is the most likely scenario. We measure the molecular gas masses of the clumps and find them to be sufficient for forming massive star clusters. The molecular gas depletion times of the clumps show a variation of ~2 dex, comparable with the large variation in depletion times found in starburst galaxies in the nearby Universe. This demonstrates the presence of starburst modes of star formation on local scales in the galaxy, even while the galaxy as a whole cannot be categorised as a starburst based on its global properties. Our findings suggest that the progenitor of GRB 9802425 was originated in a young massive star cluster formed in the starburst mode of star formation.

The stellar-to-halo mass relation (SHMR) is one of the main sources of information we have on the connection between galaxies and their dark matter haloes. Here we analyse in detail two popular forms of the SHMR, M*/Mh vs Mh and M*/Mh vs M*, and compare them with another physically motivated scaling relation, M*/Mh vs GMh/j*sigmahat*. Although this relation cannot predict the halo mass explicitely, it connects the stellar mass fraction to fundamental galaxy properties such as specific angular momentum (j*) and velocity dispersion (sigmahat*) via disc gravitational instability. Our detailed comparative analysis is based on one of the largest sample of galaxies with both high-quality rotation curves and near-infrared surface photometry, and leads to the following results: (i) M*/Mh vs Mh and M*/Mh vs M* are not just two alternative parametrizations of the same relation, but two significantly different relations; (ii) M*/Mh vs GMh/j*sigmahat* outperforms the two popular relations in terms of tightness, correlation strength and significance; (iii) j* and sigmahat* play an equally important role in our scaling relation, and it is their interplay that constrains M*/Mh so tightly; (iv) the evolution of M*/Mh, j* and sigmahat* is regulated by disc gravitational instability: when M*/Mh varies, j* and sigmahat* also vary as predicted by our scaling relation, thus erasing the memory of such evolution. This implies that the process of disc gravitational instability is intriguingly uniform across disc galaxies of all morphological types: from lenticulars to blue compact dwarfs. In particular, the cosmic variance of Toomre's Q is 0.2 dex, a universal value for both stars and atomic gas.

Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[$\alpha$/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fueled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer gas disk, inclined with respect to the main galactic plane, and with a significantly poorer chemical content. The first passage of the last major merger galaxy triggers tidal compression in the outer disk, which increases the gas density and eventually leads to star formation, at a metallicity 0.75 dex lower than the inner galaxy. This forms the first stars of the low-[$\alpha$/Fe] sequence. These in situ stars have halo-like kinematics, similarly to what is observed in the Milky Way, due to the inclination of the outer disk which eventually aligns with the inner one via gravitational torques. We show that this tilting disk scenario is likely to be common in Milky-Way like galaxies. This process implies that the low-[$\alpha$/Fe] sequence is populated in situ, simultaneously from two formation channels, in the inner and the outer galaxy, with distinct metallicities. This contrasts with purely sequential scenarios for the assembly of the Milky Way disk and could be tested observationally.

Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot disk populated by high-[$\alpha$/Fe] stars, formed both in situ and in accreted satellite galaxies. We find that the diversity of evolutionary tracks followed by the simulated galaxy and its progenitors leads to very little overlap of the in situ and accreted populations for any given chemical composition. At a given age, the spread in [$\alpha$/Fe] abundance ratio results from the diversity of physical conditions in VINTERGATAN and its satellites, with an enhancement in [$\alpha$/Fe] found in stars formed during starburst episodes. Later, the cessation of the merger activity promotes the in situ formation of stars in the low-[$\alpha$/Fe] regime, in a radially extended, thin and overall kinematically colder disk, thus establishing chemically bimodal thin and thick disks, in line with observations. We draw links between notable features in the [Fe/H] - [$\alpha$/Fe] plane with their physical causes, and propose a comprehensive formation scenario explaining self-consistently, in the cosmological context, the main observed properties of the Milky Way.

Spectroscopic surveys of the Milky Way's stars have revealed spatial, chemical and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at $z\sim 1.5$, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[$\alpha$/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[$\alpha$/Fe] stars. This event leads to a bimodality in [$\alpha$/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since $z\sim 1$ is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scaleheights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H]-profiles above the midplane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, $\alpha$-enhancements, and ages of stars, e.g. the most metal-poor stars in the low-[$\alpha$/Fe] sequence are found to have a scaleheight comparable to old high-[$\alpha$/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [$\alpha$/Fe]-[Fe/H] cannot be captured.

Roughly ten per cent of OB stars are kicked out of their natal clusters before ending their life as supernovae. These so-called runaway stars can travel hundreds of parsecs into the low-density interstellar medium, where momentum and energy from stellar feedback is efficiently deposited. In this work we explore how this mechanism affects large scale properties of the galaxy, such as outflows. To do so we use a new model which treats OB stars and their associated feedback processes on a star-by-star basis. With this model we compare two hydrodynamical simulations of Milky Way-like galaxies, one where we include runaways, and one where we ignore them. Including runaway stars leads to twice as many supernovae explosions in regions with gas densities ranging from 1e-5 cm^-3 to 1e-3 cm^-3. This results in more efficient heating of the inter-arm regions, and drives strong galactic winds with mass loading factors boosted by up to one order of magnitude. These outflows produce a more massive and extended multi-phase circumgalactic medium, as well as a population of dense clouds in the halo. Conversely, since less energy and momentum is released in the dense star forming regions, the cold phase of the interstellar medium is less disturbed by feedback effects.

We present a detailed study of the composition of 20 M giants in the Galactic center with 15 of them confirmed to be in the Nuclear Star Cluster. As a control sample we have also observed 7 M giants in the Milky Way Disk with similar stellar parameters. All 27 stars are observed using the NIRSPEC spectograph on the KECK II telescope in the K-band at a resolving power of R=23,000. We report the first silicon abundance trends versus [Fe/H] for stars in the Galactic center. While finding a disk/bulge like trend at subsolar metallicities, we find that [Si/Fe] is enhanced at supersolar metallicities. We speculate on possible enrichment scenarios to explain such a trend. However, the sample size is modest and the result needs to be confirmed by additional measurements of silicon and other \textalpha-elements. We also derive a new distribution of [Fe/H] and find the most metal rich stars at [Fe/H]=+0.5 dex, confirming our earlier conclusions that the Galactic center hosts no stars with extreme chemical composition.

We report that the density profiles of globular cluster (GC) systems in a sample of 17 early-type galaxies (ETGs) show breaks at the radii where the gravitational acceleration exerted by the stars equals the galactic acceleration scale $a_0$ known from the radial acceleration relation or MOND. The match with the other characteristic radii in the galaxy is not that close. We propose possible explanations in the frameworks of the $\Lambda$CDM model and MOND. We find tentative evidence that in the $\Lambda$CDM context, GCs reveal not only the masses of the dark halos through the richness of the GC systems but also the concentrations through the break radii of the GC systems.

Dense stellar systems in general and star clusters in particular have recently regained the interest of the extragalactic and even cosmology communities, due to the role they could play as actors and probes of re-ionization, galactic archeology and the dark matter content of galaxies, among many others. In the era of the exploitation and the preparation of large stellar surveys (Gaia, APOGEE, 4MOST, WEAVE), of the detection of gravitational waves mostly originating from dense regions like the cores of clusters (Ligo, LISA), and in an always more holistic view of galaxy formation (HARMONI, Euclid, LSST, soon to be known as the Vera Rubin Survey Telescope, VRST), a complete theory on the formation and evolution of clusters is needed to interpret the on-going and forthcoming data avalanche. In this context, the community carries an effort to model the aspects of star cluster formation and evolution in galactic and even cosmological context. However, it is not always easy to understand the caveats and the shortcuts taken in theories and simulations, and their implications on the conclusions drawn. I take the opportunity of this document to highlight three of these topics and discuss why some shortcuts taken by the community are or could be misleading.

The formation of globular clusters remains an open debate. Dwarf starburst galaxies are efficient at forming young massive clusters with similar masses as globular clusters and may hold the key to understanding their formation. We study star cluster formation in a tidal debris - including the vicinity of three tidal dwarf galaxies - in a massive gas dominated collisional ring around NGC~5291. These dwarfs have physical parameters which differ significantly from local starbursting dwarfs. They are gas-rich, highly turbulent, have a gas metallicity already enriched up to half-solar, and are expected to be free of dark matter. The aim is to study massive star cluster formation in this as yet unexplored type of environment. We use imaging from the Hubble Space Telescope using broadband filters covering the wavelength range from the near-ultraviolet to the near-infrared. We determine the masses and ages of the cluster candidates by using the spectral energy distribution-fitting code CIGALE, carefully considering age-extinction degeneracy effects on the estimation of the physical parameters. We find that the tidal dwarf galaxies in the ring of NGC 5291 are forming star clusters with an average efficiency of $\sim40\%$, comparable to blue compact dwarf galaxies. We also find massive star clusters for which the photometry suggests that they were formed at the very birth of the tidal dwarf galaxies and have survived for several hundred million years. Therefore our study shows that extended tidal dwarf galaxies and compact clusters may be formed simultaneously. In the specific case observed here, the young star clusters are not massive enough to survive for a Hubble time. However one may speculate that similar objects at higher redshift, with higher star formation rate, might form some of the long lived globular clusters.

We report the Atacama Large Millimeter/submillimeter Array (ALMA) observations of the metal rich host galaxy of superluminous supernova (SLSN) PTF10tpz, a barred spiral galaxy at z=0.03994. We find the CO(1-0) emission to be confined within the bar of the galaxy. The distribution and kinematics of molecular gas in the host galaxy resemble gas flows along two lanes running from the tips of the bar towards the galaxy center. These gas lanes end in a gaseous structure in the inner region of the galaxy, likely associated with an inner Lindblad resonance. The interaction between the large-scale gas flows in the bar and the gas in the inner region plausibly leads to the formation of massive molecular clouds and consequently massive clusters. This in turn can result in formation of massive stars, and thus the likely progenitor of the SLSN in a young, massive cluster. This picture is consistent with SLSN PTF10tpz being located near the inner structure. We find the molecular gas in the vicinity of the SLSN to have high surface densities, comparable with those in interacting galaxies or starburst regions in nearby galaxies. This lends support to high densities being favorable conditions for formation of SLSNe progenitors.

Observations find a median star formation efficiency per free-fall time in Milky Way Giant Molecular Clouds (GMCs) on the order of $\epsilon_{\rm ff}\sim 1\%$ with dispersions of $\sim0.5\,{\rm dex}$. The origin of this scatter in $\epsilon_{\rm ff}$ is still debated and difficult to reproduce with analytical models. We track the formation, evolution and destruction of GMCs in a hydrodynamical simulation of a Milky Way-like galaxy and by deriving cloud properties in an observationally motivated way, measure the distribution of star formation efficiencies which are in excellent agreement with observations. We find no significant link between $\epsilon_{\rm ff}$ and any measured global property of GMCs (e.g. gas mass, velocity dispersion). Instead, a wide range of efficiencies exist in the entire parameter space. From the cloud evolutionary tracks, we find that each cloud follow a \emph{unique} evolutionary path which gives rise to wide diversity in all properties. We argue that it is this diversity in cloud properties, above all else, that results in the dispersion of $\epsilon_{\rm ff}$.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

(abridged) Context. Gravitational fields at the outskirts of early-type galaxies (ETGs) are difficult to constrain observationally. It thus remains poorly explored how well the LCDM and MOND hypotheses agree with ETGs. Aims. This led us to gather a large sample of ETGs and examine homogeneously which dark matter halos they occupy, whether the halos follow the theoretically predicted stellar-to-halo mass relation (SHMR) and the halo mass-concentration relation (HMCR), whether ETGs obey MOND and the radial acceleration relation (RAR) observed for late-type galaxies (LTGs), and finally whether LCDM or MOND perform better in ETGs. Methods. We employed Jeans analysis of radial velocities of globular clusters (GCs). We analysed nearly all ETGs having more than about 100 archival GC radial velocity measurements. The GC systems of our 17 ETGs extend mostly over ten effective radii. A LCDM simulation of GC formation helped us to interpret the results. Results. Successful LCDM fits are found for all galaxies, but compared to the theoretical HMCR and SHMR, the best-fit halos usually have concentrations that are too low and stellar masses that are too high for their masses. This might be because of tidal stripping of the halos or because ETGs and LTGs occupy different halos. Most galaxies can be fitted by the MOND models successfully as well, but for some of the galaxies, especially those in centers of galaxy clusters, the observed GC velocity dispersions are too high. This might be a manifestation of the additional dark matter that MOND requires in galaxy clusters. Additionally, we find many signs that the GC systems were perturbed by galaxy interactions. Formal statistical criteria prefer the best-fit LCDM models over the MOND models, but this might be due to the higher flexibility of the LCDM models. The MOND approach can predict the GC velocity dispersion profiles better.

Modelling the propagation of supernova (SN) bubbles, in terms of energy, momentum and spatial extent, is critical for simulations of galaxy evolution which do not capture these scales. To date, small scale models of SN feedback predict that the evolution of above-mentioned quantities can be solely parameterised by average quantities of the surrounding gas, such as density. However, most of these studies neglect the turbulent motions of this medium. In this paper, we study the propagation and evolution of SNe in turbulent environments. We confirm that the time evolution of injected energy and momentum can be characterised by the average density. However, the details of the density structure of the interstellar medium play a crucial role in the spatial extent of the bubble, even at a given average density. We demonstrate that spherically symmetric models of SN bubbles do not model well their spatial extent, and therefore cannot not be used to design sub-grid models of SNe feedback at galactic and cosmological scales.

The physical origin of enhanced star formation activity in interacting galaxies remains an open question. Knowing whether starbursts are triggered by an increase of the quantity of dense gas or an increase of the star formation efficiency would improve our understanding of galaxy evolution and allow to transpose the results obtained in the local Universe to high redshift galaxies. In this paper, we analyze a parsec-resolution simulation of an Antennae-like model of interacting galaxies. We find that the interplay of physical processes has complex and important variations in time and space, through different combinations of mechanisms like tides, shear and turbulence. These can have similar imprints on observables like depletion time and CO emission. The densest gas only constitutes the tail of the density distribution of some clouds, but exists in large excess in others. The super-linearity of the star formation rate dependence on gas density implies that this excess translates into a reduction of depletion times, and thus leads to a deviation from the classical star formation regime, visible up to galactic scales. These clouds are found in all parts of the galaxies, but their number density varies from one region to the next, due to different cloud assembly mechanisms. Therefore, the dependence of cloud and star formation-related quantities (like CO flux and depletion time) on the scale at which they are measured also varies across the galaxies. We find that the $\alpha_{\rm CO}$ conversion factor between the CO luminosity and molecular gas mass has even stronger spatial than temporal variations in a system like the Antennae. These results raise a number of cautionary notes for the interpretation of observations of unresolved star-forming regions, but also predict that the diversity of environments for star formation will be better captured by the future generations of instruments.

We use state-of-the-art cosmological zoom simulations to explore the distribution of neutral gas in and around galaxies that gives rise to high column density \ion{H}{i} \mbox{Ly-$\alpha$} absorption (formally, sub-DLAs and DLAs) in the spectra of background quasars. Previous cosmological hydrodynamic simulations under-predict the mean projected separations $(b)$ of these absorbers relative to the host, and invoke selection effects to bridge the gap with observations. On the other hand, single lines of sight (LOS) in absorption cannot uniquely constrain the galactic origin. Our simulations match all observational data, with DLA and sub-DLA LOS existing over the entire probed parameter space ($-4\lesssim $[M/H]$\lesssim 0.5$, $b<50$ kpc) at all redshifts ($z\sim 0.4 - 3.0$). We demonstrate how the existence of DLA LOS at $b\gtrsim 20-30$ kpc from a massive host galaxy require high numerical resolution, and that these LOS are associated with dwarf satellites in the main halo, stripped metal-rich gas and outflows. Separating the galaxy into interstellar ("\ion{H}{i} disc") and circumgalactic ("halo") components, we find that both components significantly contribute to damped \ion{H}{i} absorption LOS. Above the sub-DLA (DLA) limits, the disc and halo contribute with $\sim 60 (80)$ and $\sim 40 (20)$ per cent, respectively. Our simulations confirm analytical model-predictions of the DLA-distribution at $z\lesssim 1$. At high redshift ($z\sim 2-3$) sub-DLA and DLAs occupy similar spatial scales, but on average separate by a factor of two by $z\sim 0.5$. On whether sub-DLA and DLA LOS sample different stellar-mass galaxies, such a correlation can be driven by a differential covering-fraction of sub-DLA to DLA LOS with stellar mass. This preferentially selects sub-DLA LOS in more massive galaxies in the low-$z$ universe.

We present the first spatially resolved (~0.5 kpc) measurements of the molecular gas depletion time $\tau_{depl}$ across the disk of the interacting spiral galaxy NGC\,2276, a system with an asymmetric morphology in various SFR tracers. To estimate $\tau_{depl}$, we use new NOEMA observations of the $^{12}$CO(1-0) emission tracing the bulk molecular gas reservoir in NGC 2276, and extinction-corrected H$\alpha$ measurements obtained with the PMAS/PPaK integral field unit for robust estimates of the SFR. We find a systematic decrease in $\tau_{depl}$ of 1-1.5 dex across the disk of NGC 2276, with a further, abrupt drop in $\tau_{depl}$ of ~1 dex along the galaxy's western edge. The global $\tau_{depl}$ in NGC 2776 is $\tau_{depl}=0.55$ Gyr, insistent with literature measurements for the nearby galaxy population. Such a large range in $\tau_{depl}$ on sub-kpc scales has never previously been observed within an individual isolated or pre-merger system. When using a metallicity-dependent molecular gas conversion factor X$\rm_{CO}$ the variation decreases by 0.5 dex. We attribute the variation in $\tau_{depl}$ to the influence of galactic-scale tidal forces and ram pressure on NGC 2276's molecular interstellar medium (ISM). Our observations add to the growing body of numerical and observational evidence that galaxy-galaxy interactions significantly modify the molecular gas properties and star-forming activity within galactic disks throughout the interaction, and not just during the final merger phase.

The conversion factor $\alpha_{\rm CO}$ from the observable CO(1-0) luminosity to the mass of molecular gas is known to vary between isolated galaxies and some mergers, but the underlying reasons are not clearly understood. Thus, the value(s) of $\alpha_{\rm CO}$ to be adopted remain highly uncertain. To provide better constraints, we apply the large velocity gradient method to a series of hydrodynamical simulations of galaxies and derive the evolution of $\alpha_{\rm CO}$. We report significant variations of $\alpha_{\rm CO}$, and identify three distinct regimes: disk galaxies, starbursts and post-burst phases. We show that estimating the star formation rate over 20 Myr smooths out some of these differences, but still maintains a distinction between disks and starbursts. We find a tighter correlation of $\alpha_{\rm CO}$ with the gas depletion time than with star formation rate, yet with deviations induced by the transitions to and from the starburst episodes. We conclude that $\alpha_{\rm CO}$ fluctuates because of both the feedback energy and the velocity dispersion. Identifying the phase of an interaction by classical means (e.g. morphology, luminosity) could then help selecting the relevant conversion factor to be used and get more accurate estimates of the molecular masses of galaxies.

Using hydrodynamical simulations of entire galactic discs similar to the Milky Way, reaching 4.6pc resolution, we study the origins of observed physical properties of giant molecular clouds (GMCs). We find that efficient stellar feedback is a necessary ingredient in order to develop a realistic interstellar medium (ISM), leading to molecular cloud masses, sizes, velocity dispersions and virial parameters in excellent agreement with Milky Way observations. GMC scaling relations observed in the Milky Way, such as the mass-size ($M$--$R$), velocity dispersion-size ($\sigma$--$R$), and the $\sigma$--$R\Sigma$ relations, are reproduced in a feedback driven ISM when observed in projection, with $M\propto R^{2.3}$ and $\sigma\propto R^{0.56}$. When analysed in 3D, GMC scaling relations steepen significantly, indicating potential limitations of our understanding of molecular cloud 3D structure from observations. Furthermore, we demonstrate how a GMC population's underlying distribution of virial parameters can strongly influence the scatter in derived scaling relations. Finally, we show that GMCs with nearly identical global properties exist in different evolutionary stages, where a majority of clouds being either gravitationally bound or expanding, but with a significant fraction being compressed by external ISM pressure, at all times.

Massive star clusters are observed in a broad range of galaxy luminosity and types, and are assumed to form in dense gas-rich environments. Using a parsec-resolution hydrodynamical simulation of an isolated gas-rich low mass galaxy, we discuss here the non-linear effects of stellar feedback on the properties of star clusters with a focus on the progenitors of nuclear clusters. Our simulation shows two categories of star clusters: those for which feedback expels gas leftovers associated with their formation sites, and those, in a denser environment around which feedback fails at totally clearing the gas. We confirm that radiation feedback (photo-ionization and radiative pressure) plays a more important role than type-II supernovae in destroying dense gas structures, and altering or quenching the subsequent cluster formation. It also disturbs the cluster mass growth, by increasing the internal energy of the gas component to the point when radiation pressure overcomes the cluster gravity. We discuss how these effects may depend on the local properties of the interstellar medium, and also on the details of the subgrid recipes, which can affect the available cluster gas reservoirs, the evolution of potential nuclear clusters progenitors, and the overall galaxy morphology.

Their ubiquity and extreme densities make star clusters probes of prime importance of galaxy evolution. Old globular clusters keep imprints of the physical conditions of their assembly in the early Universe, and younger stellar objects, observationally resolved, tell us about the mechanisms at stake in their formation. Yet, we still do not understand the diversity involved: why is star cluster formation limited to 1e5 Msun objects in the Milky Way, while some dwarf galaxies like NGC 1705 are able to produce clusters 10 times more massive? Why do dwarfs generally host a higher specific frequency of clusters than larger galaxies? How to connect the present-day, often resolved, stellar systems to the formation of globular clusters at high redshift? And how do these links depend on the galactic and cosmological environments of these clusters? In this review, I present recent advances on star cluster formation and evolution, in galactic and cosmological context. The emphasis is put on the theory, formation scenarios and the effects of the environment on the evolution of the global properties of clusters. A few open questions are identified.

Stephan's Quintet, a compact group of galaxies, is often used as a laboratory to study a number of phenomena, including physical processes in the interstellar medium, star formation, galaxy evolution, and the formation of fossil groups. As such, it has been subject to intensive multi-wavelength observation campaigns. Yet, models lack constrains to pin down the role of each galaxy in the assembly of the group. We revisit here this system with multi-band deep optical images obtained with MegaCam on the Canada-France-Hawaii Telescope (CFHT), focusing on the detection of low surface brightness (LSB) structures. They reveal a number of extended LSB features, some new, and some already visible in published images but not discussed before. An extended diffuse, reddish, lopsided, halo is detected towards the early-type galaxy NGC 7317, the role of which had so far been ignored in models. The presence of this halo made of old stars may indicate that the group formed earlier than previously thought. Finally, a number of additional diffuse filaments are visible, some close to the foreground galaxy NGC 7331 located in the same field. Their structure and association with mid--IR emission suggest contamination by emission from Galactic cirrus.

We use hydrodynamical simulations of a Cartwheel-like ring galaxy, modelled as a nearly head-on collision of a small companion with a larger disc galaxy, to probe the evolution of the gaseous structures and flows, and to explore the physical conditions setting the star formation activity. Star formation is first quenched by tides as the companion approaches, before being enhanced shortly after the collision. The ring ploughs the disc material as it radially extends, and almost simultaneously depletes its stellar and gaseous reservoir into the central region, through the spokes, and finally dissolve 200 Myr after the collision. Most of star formation first occurs in the ring before this activity is transferred to the spokes and then the nucleus. We thus propose that the location of star formation traces the dynamical stage of ring galaxies, and could help constrain their star formation histories. The ring hosts tidal compression associated with strong turbulence. This compression yields an azimuthal asymmetry, with maxima reached in the side furthest away from the nucleus, which matches the star formation activity distribution in our models and in observed ring systems. The interaction triggers the formation of star clusters significantly more massive than before the collision, but less numerous than in more classical galaxy interactions. The peculiar geometry of Cartwheel-like objects thus yields a star (cluster) formation activity comparable to other interacting objects, but with notable second order differences in the nature of turbulence, the enhancement of the star formation rate, and the number of massive clusters formed.

Kinematically cold tidal streams of globular clusters (GC) are excellent tracers of the Galactic gravitational potential at moderate Galactocentric distances, and can also be used as probes of the law of gravity on Galactic scales. Here, we compare for the first time the generation of such streams in Newtonian and Milgromian gravity (MOND). We first compute analytical results to investigate the expected shape of the GC gravitational potential in both frameworks, and we then run N-body simulations with the Phantom of Ramses code. We find that the GCs tend to become lopsided in MOND. This is a consequence of the external field effect which breaks the strong equivalence principle. When the GC is filling its tidal radius the lopsidedness generates a strongly asymmetric tidal stream. In Newtonian dynamics, such markedly asymmetric streams can in general only be the consequence of interactions with dark matter subhalos, giant molecular clouds, or interaction with the Galactic bar. In these Newtonian cases, the asymmetry is the consequence of a very large gap in the stream, whilst in MOND it is a true asymmetry. This should thus allow us in the future to distinguish these different scenarios by making deep observations of the environment of the asymmetric stellar stream of Palomar 5. Moreover, our simulations indicate that the high internal velocity dispersion of Palomar 5 for its small stellar mass would be natural in MOND.

We study the impact of stellar feedback in shaping the density and velocity structure of neutral hydrogen (HI) in disc galaxies. For our analysis, we carry out $\sim 4.6$pc resolution $N$-body+adaptive mesh refinement (AMR) hydrodynamic simulations of isolated galaxies, set up to mimic a Milky Way (MW), and a Large and Small Magellanic Cloud (LMC, SMC). We quantify the density and velocity structure of the interstellar medium using power spectra and compare the simulated galaxies to observed HI in local spiral galaxies from THINGS (The HI Nearby Galaxy Survey). Our models with stellar feedback give an excellent match to the observed THINGS HI density power spectra. We find that kinetic energy power spectra in feedback regulated galaxies, regardless of galaxy mass and size, show scalings in excellent agreement with super-sonic turbulence ($E(k)\propto k^{-2}$) on scales below the thickness of the HI layer. We show that feedback influences the gas density field, and drives gas turbulence, up to large (kpc) scales. This is in stark contrast to density fields generated by large scale gravity-only driven turbulence. We conclude that the neutral gas content of galaxies carries signatures of stellar feedback on all scales.

We present a cosmological zoom-in simulation of a Milky Way-like galaxy used to explore the formation and evolution of star clusters. We investigate in particular the origin of the bimodality observed in the colour and metallicity of globular clusters, and the environmental evolution through cosmic times in the form of tidal tensors. Our results self-consistently confirm previous findings that the blue, metal-poor clusters form in satellite galaxies which are accreted onto the Milky Way, while the red, metal-rich clusters form mostly in situ or, to a lower extent in massive, self-enriched galaxies merging with the Milky Way. By monitoring the tidal fields these populations experience, we find that clusters formed in situ (generally centrally concentrated) feel significantly stronger tides than the accreted ones, both in the present-day, and when averaged over their entire life. Furthermore, we note that the tidal field experienced by Milky Way clusters is significantly weaker in the past than at present-day, confirming that it is unlikely that a power-law cluster initial mass function like that of young massive clusters, is transformed into the observed peaked distribution in the Milky Way with relaxation-driven evaporation in a tidal field.

Galaxy mergers are believed to trigger strong starbursts. This is well assessed by observations in the local Universe. However the efficiency of this mechanism has poorly been tested so far for high redshift, actively star forming, galaxies. We present a suite of pc-resolution hydrodynamical numerical simulations to compare the star formation process along a merging sequence of high and low z galaxies, by varying the gas mass fraction between the two models. We show that, for the same orbit, high-redshift gas-rich mergers are less efficient than low-redshift ones at producing starbursts: the star formation rate excess induced by the merger and its duration are both around 10 times lower than in the low gas fraction case. The mechanisms that account for the star formation triggering at low redshift - the increased compressive turbulence, gas fragmentation, and central gas inflows - are only mildly, if not at all, enhanced for high gas fraction galaxy encounters. Furthermore, we show that the strong stellar feedback from the initially high star formation rate in high redshift galaxies does not prevent an increase of the star formation during the merger. Our results are consistent with the observed increase of the number of major mergers with increasing redshift being faster than the respective increase in the number of starburst galaxies.

The radii of young (<100 Myr) star clusters correlate only weakly with their masses. This shallow relation has been used to argue that impulsive tidal perturbations, or `shocks', by passing giant molecular clouds (GMCs) preferentially disrupt low-mass clusters. We show that this mass-radius relation is in fact the result of the combined effect of two-body relaxation and repeated tidal shocks. Clusters in a broad range of environments including those like the solar neighbourhood evolve towards a typical radius of a few parsecs, as observed, independent of the initial radius. This equilibrium mass-radius relation is the result of a competition between expansion by relaxation and shrinking due to shocks. Interactions with GMCs are more disruptive for low-mass clusters, which helps to evolve the globular cluster mass function (GCMF). However, the properties of the interstellar medium in high-redshift galaxies required to establish a universal GCMF shape are more extreme than previously derived, challenging the idea that all GCs formed with the same power-law mass function.

Together with interstellar turbulence, gravitation is one key player in star formation. It acts both at galactic scales in the assembly of gas into dense clouds, and inside those structures for their collapse and the formation of pre-stellar cores. To understand to what extent the large scale dynamics govern the star formation activity of galaxies, we present hydrodynamical simulations in which we generalise the behaviour of gravity to make it differ from Newtonian dynamics in the low acceleration regime. We focus on the extreme cases of interacting galaxies, and compare the evolution of galaxy pairs in the dark matter paradigm to that in the Milgromian Dynamics (MOND) framework. Following up on the seminal work by Tiret & Combes, this paper documents the first simulations of galaxy encounters in MOND with a detailed Eulerian hydrodynamical treatment of baryonic physics, including star formation and stellar feedback. We show that similar morphologies of the interacting systems can be produced by both the dark matter and MOND formalisms, but require a much slower orbital velocity in the MOND case. Furthermore, we find that the star formation activity and history are significantly more extended in space and time in MOND interactions, in particular in the tidal debris. Such differences could be used as observational diagnostics and make interacting galaxies prime objects in the study of the nature of gravitation at galactic scales.

Nuclear Clusters (NCs) are common stellar systems in the centres of galaxies. Yet, the physical mechanisms involved in their formation are still debated. Using a parsec-resolution hydrodynamical simulation of a dwarf galaxy, we propose an updated formation scenario for NCs. In this 'wet migration scenario', a massive star cluster forms in the gas-rich disc, keeping a gas reservoir, and growing further while it migrates to the centre via a combination of interactions with other substructures and dynamical friction. A wet merger with another dense cluster and its own gas reservoir can occur, although this is not a pre-requisite for the actual formation of the NC. The merging process does significantly alter the properties of the NC (mass, morphology, star formation history), also quenching the on-going local star formation activity, thus leading to interesting observational diagnostics for the physical origin of NCs. A population of lower mass clusters co-exist during the simulation, but these are either destroyed via tidal forces, or have high angular momentum preventing them to interact with the NC and contribute to its growth. The proposed updated scenario emphasises the role of gas reservoirs associated with the densest star clusters formed in a gas-rich low-mass galaxy.

The strong time-dependence of the dynamics of galactic bars yields a complex and rapidly evolving distribution of dense gas and star forming regions. Although bars mainly host regions void of any star formation activity, their extremities can gather the physical conditions for the formation of molecular complexes and mini-starbursts. Using a sub-parsec resolution hydrodynamical simulation of a Milky Way-like galaxy, we probe these conditions to explore how and where bar (hydro-)dynamics favours the formation or destruction of molecular clouds and stars. The interplay between the kpc-scale dynamics (gas flows, shear) and the parsec-scale (turbulence) is key to this problem. We find a strong dichotomy between the leading and trailing sides of the bar, in term of cloud fragmentation and in the age distribution of the young stars. After orbiting along the bar edge, these young structures slow down at the extremities of the bar, where orbital crowding increases the probability of cloud-cloud collision. We find that such events increase the Mach number of the cloud, leading to an enhanced star formation efficiency and finally the formation of massive stellar associations, in a fashion similar to galaxy-galaxy interactions. We highlight the role of bar dynamics in decoupling young stars from the clouds in which they form, and discuss the implications on the injection of feedback into the interstellar medium, in particular in the context of galaxy formation.

When the cosmic star formation history peaks (z ~ 2), galaxies vigorously fed by cosmic reservoirs are gas dominated and contain massive star-forming clumps, thought to form by violent gravitational instabilities in highly turbulent gas-rich disks. However, a clump formation event has not been witnessed yet, and it is debated whether clumps survive energetic feedback from young stars, thus migrating inwards to form galaxy bulges. Here we report spatially resolved spectroscopy of a bright off-nuclear emission line region in a galaxy at z = 1.987. Although this region dominates the star formation in the galaxy disk, its stellar continuum remains undetected in deep imaging, revealing an extremely young (age < 10 Myr) massive clump, forming through the gravitational collapse of > 10$^9$ M$_{\odot}$ of gas. Gas consumption in this young clump is > 10 times faster than in the host galaxy, displaying high star formation efficiency during this phase, in agreement with our hydrodynamic simulations. The frequency of older clumps with similar masses coupled with our initial estimate of their formation rate (~ 2.5 Gyr$^{-1}$) supports long lifetimes (~ 500 Myr), favouring scenarios where clumps survive feedback and grow the bulges of present-day galaxies.

In the fundamental quest of the rotation curve of the Milky Way, the tangent-point (TP) method has long been the simplest way to infer velocities for the inner, low latitude regions of the Galactic disk from observations of the gas component. We test the validity of the method on realistic gas distribution and kinematics of the Milky Way, using a numerical simulation of the Galaxy. We show that the resulting velocity profile strongly deviates from the true rotation curve of the simulation, as it overstimates it in the central regions, and underestimates it around the bar corotation. Also, its shape strongly depends on the orientation of the stellar bar. The discrepancies are caused by highly non-uniform azimuthal velocities, and the systematic selection by the TP method of high-velocity gas along the bar and spiral arms, or low-velocity gas in less dense regions. The velocity profile is in good agreement with the rotation curve only beyond corotation, far from massive asymmetric structures. Therefore the observed velocity profile of the Milky Way inferred by the TP method is expected to be very close to the true Galactic rotation curve for 4.5<R<8 kpc. Another consequence is that the Galactic velocity profile for R<4-4.5 kpc is very likely flawed by the non-uniform azimuthal velocities, and does not represent the true Galactic rotation curve, but instead local motions. The real shape of the innermost rotation curve is probably shallower than previously thought. Using a wrong rotation curve has a dramatic impact on the modelling of the mass distribution, in particular for the bulge component of which derived enclosed mass within the central kpc and scale radius are, respectively, twice and half of the actual values. We thus strongly argue against using terminal velocities or the velocity curve from the TP method for modelling the mass distribution of the Milky Way. (abridged)

The growth of galaxies through adiabatic accretion of dark matter is one of the main drivers of galaxy evolution. By isolating it from other processes like mergers, we analyse how it affects the evolution of star clusters. Our study comprises a fast and approximate exploration of the orbital and intrinsic cluster parameter space, and more detailed monitoring of their evolution, through N-body simulations for a handful of cases. We find that the properties of present-day star clusters and their tidal tails differ very little, whether the clusters are embedded in a growing galactic halo for 12 Gyr, or in a static one.

We introduce a numerical method to integrate tidal effects on collisional systems, using any definition of the external potential as a function of space and time. Rather than using a linearisation of the tidal field, this new method follows a differential technique to numerically evaluate the tidal acceleration and its time derivative. Theses are then used to integrate the motions of the components of the collisional systems, like stars in star clusters, using a predictor-corrector scheme. The versatility of this approach allows the study of star clusters, including their tidal tails, in complex, multi-components, time-evolving external potentials. The method is implemented in the code nbody6 (Aarseth 2003).

We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with very high resolutions (3pc and 10^3 Msun) and detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes and galactic winds. The simulations are analyzed with a Large Velocity Gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer and opacity effects, and the intensity emerging from galaxies, to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially-concentrated stellar feedback sources, although CO(1-0) to CO(3-2) emission is overall dominated by low-excitation gas around the densest clumps. These results overall highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn impact the CO emission properties.

The interstellar dust content in galaxies can be traced in extinction at optical wavelengths, or in emission in the far-infrared. Several studies have found that radiative transfer models that successfully explain the optical extinction in edge-on spiral galaxies generally underestimate the observed FIR/submm fluxes by a factor of about three. In order to investigate this so-called dust energy balance problem, we use two Milky Way-like galaxies produced by high-resolution hydrodynamical simulations. We create mock optical edge-on views of these simulated galaxies (using the radiative transfer code SKIRT), and we then fit the parameters of a basic spiral galaxy model to these images (using the fitting code FitSKIRT). The basic model includes smooth axisymmetric distributions along a S\'ersic bulge and exponential disc for the stars, and a second exponential disc for the dust. We find that the dust mass recovered by the fitted models is about three times smaller than the known dust mass of the hydrodynamical input models. This factor is in agreement with previous energy balance studies of real edge-on spiral galaxies. On the other hand, fitting the same basic model to less complex input models (e.g. a smooth exponential disc with a spiral perturbation or with random clumps), does recover the dust mass of the input model almost perfectly. Thus it seems that the complex asymmetries and the inhomogeneous structure of real and hydrodynamically simulated galaxies are a lot more efficient at hiding dust than the rather contrived geometries in typical quasi-analytical models. This effect may help explain the discrepancy between the dust emission predicted by radiative transfer models and the observed emission in energy balance studies for edge-on spiral galaxies.

The distinction between globular clusters and dwarf galaxies has been progressively blurred by the recent discoveries of several extended star clusters, with size (20-30 pc) and luminosity (-6 < Mv < -2) comparable to the one of faint dwarf spheroidals. In order to explain their sparse structure, it has been suggested that they formed as star clusters in dwarf galaxy satellites that later accreted onto the Milky Way. If these clusters form in the centre of dwarf galaxies, they evolve in a tidally-compressive environment where the contribution of the tides to the virial balance can become significant, and lead to a super-virial state and subsequent expansion of the cluster, once removed. Using N-body simulations, we show that a cluster formed in such an extreme environment undergoes a sizable expansion, during the drastic variation of the external tidal field due to the accretion process. However, we show that the expansion due to the removal of the compressive tides is not enough to explain the observed extended structure, since the stellar systems resulting from this process are always more compact than the corresponding clusters that expand in isolation due to two-body relaxation. We conclude that an accreted origin of extended globular clusters is unlikely to explain their large spatial extent, and rather favor the hypothesis that such clusters are already extended at the stage of their formation.

We study the connection between the large-scale dynamics and the gas fueling toward a central black hole via the analysis of a Milky Way-like simulation at sub-parsec resolution. This allows us to follow a set of processes at various scales (e.g., the triggering of inward gas motion towards inner resonances via the large-scale bar, the connection to the central black hole via mini spirals) in a self-consistent manner. This simulation provides further insights on the role of shear for the inhibition of star formation within the bar in regions with significant amount of gas. We also witness the decoupling of the central gas and nuclear cluster from the large-scale disc, via interactions with the black hole. This break of symmetry in the mass distribution triggers the formation of gas clumps organised in a time-varying 250 pc ring-like structure, the black hole being offset by about 70 pc from its centre. Some clumps form stars, while most get disrupted or merge. Supernovae feedback further creates bubbles and filaments, some of the gas being expelled to 100 pc or higher above the galaxy plane. This helps remove angular momentum from the gas, which gets closer to the central dark mass. Part of the gas raining down is being accreted, forming a 10~pc polar disc-like structure around the black hole, leading to an episode of star formation. This gives rise to multiple stellar populations with significantly different angular momentum vectors, and may lead to a natural intermittence in the fueling of the black hole.

We present a hydrodynamical simulation of an Antennae-like galaxy merger at parsec resolution, including a multi-component model for stellar feedback and reaching numerical convergence in the global star formation rate for the first time. We analyse the properties of the dense stellar objects formed during the different stages of the interaction. Each galactic encounter triggers a starburst activity, but the varying physical conditions change the triggering mechanism of each starburst. During the first two pericenter passages, the starburst is spatially extended and forms many star clusters. However, the starburst associated to the third, final passage is more centrally concentrated: stars form almost exclusively in the galactic nucleus and no new star cluster is formed. The maximum mass of stars clusters in this merger is more than 30 times higher than those in a simulation of an isolated Milky Way-like galaxy. Antennae-like mergers are therefore a formation channel of young massive clusters possibly leading to globular clusters. Monitoring the evolution of a few clusters reveals the diversity of formation scenarios including the gathering and merger of gas clumps, the monolithic formation and the hierarchical formation in sub-structures inside a single cloud. Two stellar objects formed in the simulation yield the same properties as ultra-compact dwarf galaxies. They share the same formation scenario than the most massive clusters, but have a larger radius either since birth, or get it after a violent interaction with the galactic center. The diversity of environments across space and time in a galaxy merger can account for the diversity of the stellar objects formed, both in terms of mass and size.

Using parsec-resolution simulations of a typical galaxy merger, we study the triggering of starbursts by connecting the (inter-)galactic dynamics to the structure of the interstellar medium. The gravitational encounter between two galaxies enhances tidal compression over large volumes, which increases and modifies the turbulence, in particular its compressive mode with respect to the solenoidal one. This generates an excess of dense gas leading to intense star formation activity. Along the interaction, the compressive turbulence modifies the efficiency of gas-to-star conversion which, in the Schmidt-Kennicutt diagram, drives the galaxies from the sequence of discs to that of starbursts.

The Schmidt-Kennicutt relation links the surface densities of gas to the star formation rate in galaxies. The physical origin of this relation, and in particular its break, i.e. the transition between an inefficient regime at low gas surface densities and a main regime at higher densities, remains debated. Here, we study the physical origin of the star formation relations and breaks in several low-redshift galaxies, from dwarf irregulars to massive spirals. We use numerical simulations representative of the Milky Way, the Large and the Small Magellanic Clouds with parsec up to subparsec resolution, and which reproduce the observed star formation relations and the relative variations of the star formation thresholds. We analyze the role of interstellar turbulence, gas cooling, and geometry in drawing these relations, at 100 pc scale. We suggest in particular that the existence of a break in the Schmidt- Kennicutt relation could be linked to the transition from subsonic to supersonic turbulence and is independent of self-shielding effects. This transition being connected to the gas thermal properties and thus to the metallicity, the break is shifted toward high surface densities in metal-poor galaxies, as observed in dwarf galaxies. Our results suggest that together with the collapse of clouds under self-gravity, turbulence (injected at galactic scale) can induce the compression of gas and regulate star formation.

We present the hierarchical structure of the gas and its rapid evolution in the central region of a simulation of the entire Milky Way, run at subparsec resolution. We emphasize the coupling between the kpc-scale dynamics, the molecular ring and the central 5 pc disk feeding the super massive black hole.

We aim at addressing the questions related to galaxy mass assembly through major and minor wet merging processes in the redshift range 1<z<2. A consequent fraction of Milky Way like galaxies are thought to have undergone an unstable clumpy phase at this early stage. Using the adaptive mesh refinement code RAMSES, with a recent physically-motivated implementation of stellar feedback, we build the Merging and Isolated high-Redshift Adaptive mesh refinement Galaxies (MIRAGE) sample. It is composed of 20 mergers and 3 isolated idealized disks simulations with global physical properties in accordance with the 1<z<2 mass complete sample MASSIV. The numerical hydrodynamical resolution reaches 7 parsecs in the smallest Eulerian cells. Our simulations include: star formation, metal line cooling, metallicity advection, and a recent implementation of stellar feedback which encompasses OB-type stars radiative pressure, photo-ionization heating, and supernovae. The initial conditions are set to match the z~2 observations, thanks to a new public code DICE. The numerical resolution allows us to follow the formation and evolution of giant clumps formed in-situ from Jeans instabilities triggered by high initial gas fraction. The star formation history of isolated disks shows stochastic star formation rate, which proceeds from the complex behavior of the giant clumps. Our minor and major gas-rich merger simulations do not trigger starbursts, suggesting a saturation of the star formation in a turbulent and clumpy interstellar medium fed by substantial accretion from the circum-galactic medium. Our simulations are close to the normal regime of the disk-like star formation on a Schmidt-Kennicutt diagram. The mass-size relation and its rate of evolution matches observations, suggesting that the inside-out growth mechanisms of the stellar disk do not necessarily require to be achieved through a cold accretion.

The organization of the interstellar medium in disk galaxies obeys the large scale dynamics set by kpc-size structures. Improving our knowledge of how the dense, molecular gas is distributed in a disk is an important step in our understanding of star formation at galactic scale. Using a recently published simulation of a Milky Way-like galaxy, we explore the formation and dynamical organization of the star forming gas in a proto-typical disk. Along spiral arms, we report the formation of regularly spaced clouds, called beads on a string and spurs. The former form through gravitational instabilities while the later originate from Kelvin-Helmholtz instabilities. We propose that the co-existence of both structures in the same galaxy can be explained by a different role of the disk dynamics, depending on the location within the disk. In particular, we highlight the impact of the pitch angle of the spiral arm in the development of either type of structure.

We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy, at the resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through photo-ionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at subparsec resolution for a few cloud lifetimes, and at 0.05 pc for about a cloud crossing time. Turbulence cascade and gravitation from the kpc scales are de facto included in smaller structures like molecular clouds. We show that the formation of a bar influences the dynamics of the central ~100 pc by creating resonances. At larger radii, the spiral arms host the formation of regularly spaced clouds: beads on a string and spurs. These instabilities pump turbulent energy into the gas, generally in the supersonic regime. Because of asymmetric drift, the supernovae explode outside of their gaseous nursery, which diminishes the effect of feedback on the structure of clouds. The evolution of clouds is thus mostly due to fragmentation and gas consumption, regulated mainly by supersonic turbulence. The transition from turbulence supported to self-gravitating gas is detected in the gas density probability distribution function at ~2000 cm^-3. The power spectrum density suggests that gravitation governs the hierarchical organisation of structures from the galactic scale down to a few parsecs.

Star-forming disk galaxies at high redshift are often subject to violent disk instability, characterized by giant clumps whose fate is yet to be understood. The main question is whether the clumps disrupt within their dynamical timescale (<50Myr), like molecular clouds in today's galaxies, or whether they survive stellar feedback for more than a disk orbital time (~300Myr) in which case they can migrate inward and help building the central bulge. We present 3.5-7pc resolution AMR simulations of high-redshift disks including photo-ionization, radiation pressure, and supernovae feedback (Renaud et al. 2013, and Perret et al., this astro-ph issue). Our modeling of radiation pressure determines the mass loading and initial velocity of winds from basic physical principles. We find that the giant clumps produce steady outflow rates comparable to and sometimes somewhat larger than their star formation rate, with velocities largely sufficient to escape galaxy. The clumps also lose mass, especially old stars, by tidal stripping, and the stellar populations contained in the clumps hence remain relatively young (<=200Myr), as observed. The clumps survive gaseous outflows and stellar loss, because they are wandering in gas-rich turbulent disks from which they can re-accrete gas at high rates compensating for outflows and tidal stripping, overall keeping realistic and self-regulated gaseous and stellar masses. Our simulations produce gaseous outflows with velocities, densities and mass loading consistent with observations, and at the same time suggest that the giant clumps survive for hundreds of Myr and complete their migration to the center of high-redshift galaxies, without rapid dispersion and reformation of clumps.

Interacting galaxies favor the formation of star clusters but are also suspected to affect their evolution through an intense and rapidly varying tidal field. Treating this complex behaviour remains out-of-reach of (semi-)analytical studies. By computing the tidal field from galactic models and including it into star-by-star N-body simulations of star clusters, we monitor the structure and mass evolution of a population of clusters in a galaxy major merger, taking the Antennae galaxies (NGC 4038/39) as a prototype. On the long timescale (~ 10^9 yr), the merger only indirectly affects the evolution of clusters by modifying their orbits in or around the galaxies: the mass-loss of clusters in the merger remnant is faster, while clusters ejected in the tidal debris survive much longer, compared to in an isolated galaxy. The tidal perturbations of the galactic collisions themselves are too short lived and not strong enough to significantly influence the structure and dissolution of realistically dense/massive star clusters.

We present an analytical model of the relation between the surface density of gas and star formation rate in galaxies and clouds, as a function of the presence of supersonic turbulence and the associated structure of the interstellar medium. The model predicts a power-law relation of index 3/2, flattened under the effects of stellar feedback at high densities or in very turbulent media, and a break at low surface densities when ISM turbulence becomes too weak to induce strong compression. This model explains the diversity of star formation laws and thresholds observed in nearby spirals and their resolved regions, the Small Magellanic Cloud, high-redshift disks and starbursting mergers, as well as Galactic molecular clouds. While other models have proposed interstellar dust content and molecule formation to be key ingredients to the observed variations of the star formation efficiency, we demonstrate instead that these variations can be explained by interstellar medium turbulence and structure in various types of galaxies.

Long tails and streams of stars are the most noticeable upshots of galaxy collisions. Their origin as gravitational, tidal, disturbances has however been recognized only less than fifty years ago and more than ten years after their first observations. This Review describes how the idea of galactic tides emerged, in particular thanks to the advances in numerical simulations, from the first ones that included tens of particles to the most sophisticated ones with tens of millions of them and state-of-the-art hydrodynamical prescriptions. Theoretical aspects pertaining to the formation of tidal tails are then presented. The third part of the review turns to observations and underlines the need for collecting deep multi-wavelength data to tackle the variety of physical processes exhibited by collisional debris. Tidal tails are not just stellar structures, but turn out to contain all the components usually found in galactic disks, in particular atomic / molecular gas and dust. They host star-forming complexes and are able to form star-clusters or even second-generation dwarf galaxies. The final part of the review discusses what tidal tails can tell us (or not) about the structure and content of present-day galaxies, including their dark components, and explains how tidal tails may be used to probe the past evolution of galaxies and their mass assembly history. On-going deep wide-field surveys disclose many new low-surface brightness structures in the nearby Universe, offering great opportunities for attempting galactic archeology with tidal tails.

We present smoothed particle hydrodynamic models of the interactions in the compact galaxy group, Stephan's Quintet. This work is extension of the earlier collisionless N-body simulations of Renaud et al. in which the large-scale stellar morphology of the group was modeled with a series of galaxy-galaxy interactions in the simulations. Including thermohydrodynamic effects in this work, we further investigate the dynamical interaction history and evolution of the intergalactic gas of Stephan's Quintet. The major features of the group, such as the extended tidal features and the group-wide shock, enabled us to constrain the models reasonably well, while trying to reproduce multiple features of the system. We found that reconstructing the two long tails extending from NGC 7319 toward NGC 7320c one after the other in two separate encounters is very difficult and unlikely, because the second encounter usually destroys or distorts the already-generated tidal structure. Our models suggest the two long tails may be formed simultaneously from a single encounter between NGC 7319 and 7320c, resulting in a thinner and denser inner tail than the outer one. The tails then also run parallel to each other as observed. The model results support the ideas that the group-wide shock detected in multi-wavelength observations between NGC 7319 and 7318b and the starburst region north of NGC 7318b are triggered by the high-speed collision between NGC 7318b and the intergalactic gas. Our models show that a gas bridge is formed by the high-speed collision and clouds in the bridge continue to interact for some tens of millions of years after the impact. This produces many small shocks in that region, resulting a much longer cooling time than that of a single impact shock.

We present a novel and flexible tensor approach to computing the effect of a time-dependent tidal field acting on a stellar system. The tidal forces are recovered from the tensor by polynomial interpolation in time. The method has been implemented in a direct-summation stellar dynamics integrator (NBODY6) and test-proved through a set of reference calculations: heating, dissolution time and structural evolution of model star clusters are all recovered accurately. The tensor method is applicable to arbitrary configurations, including the important situation where the background potential is a strong function of time. This opens up new perspectives in stellar population studies reaching to the formation epoch of the host galaxy or galaxy cluster, as well as for star-burst events taking place during the merger of large galaxies. A pilot application to a star cluster in the merging galaxies NGC 4038/39 (the Antennae) is presented.

The evolution of compact groups of galaxies may represent one of the few places in the nearby universe in which massive galaxies are being forged through a complex set of processes involving tidal interaction, ram-pressure stripping, and perhaps finally "dry-mergers" of galaxies stripped of their cool gas. Using collisionless N-body simulations, we propose a possible scenario for the formation of one of the best studied compact groups: Stephan's Quintet. We define a serial approach which allows us to consider the history of the group as sequence of galaxy-galaxy interactions seen as relatively separate events in time, but chained together in such a way as to provide a plausible scenario that ends in the current configuration of the galaxies. By covering a large set of parameters, we claim that it is very unlikely that both major tidal tails of the group have been created by the interaction between the main galaxy and a single intruder. We propose instead a scenario based on two satellites orbiting the main disk, plus the recent involvement of an additional interloper, coming from the background at high speed. This purely N-body study of the quintet will provide a parameter-space exploration of the basic dynamics of the group that can be used as a basis for a more sophisticated N-body/hydrodynamic study of the group that is necessary to explain the giant shock structure and other purely gaseous phenomena observed in both the cold, warm and hot gas in the group.

In interacting galaxies, strong tidal forces disturb the global morphology of the progenitors and give birth to the long stellar, gaseous and dusty tails often observed. In addition to this destructive effect, tidal forces can morph into a transient, protective setting called compressive mode. Such modes then shelter the matter in their midst by increasing its gravitational binding energy. This thesis focuses on the study of this poorly known regime by quantifying its properties thanks to numerical and analytical tools applied to a spectacular merging system of two galaxies, commonly known as the Antennae galaxies. N-body simulations of this pair yield compressive modes in the regions where observations reveal a burst of star formation. Furthermore, characteristic time- and energy scales of these modes match well those of self-gravitating substructures such as star clusters and tidal dwarf galaxies. These results suggest that the compressive modes of tidal fields plays an important role in the formation and evolution of young clusters, at least in a statistical sense, over a lapse of ~10 million years. Preliminary results from simulations of stellar associations highlight the importance of embedding the clusters in the evolving background galaxies to account precisely for their morphology and internal evolution.

We present a new high-resolution N-body/SPH simulation of an encounter of two gas-rich disk galaxies which closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation and feedback from SNII. The large-scale morphology and kinematics are determined by the internal structure and the orbit of the progenitor disks. The properties of the central region, in particular the starburst in the overlap region, only match the observations for a very short time interval of ~20 Myr after the second encounter. This indicates that the Antennae galaxies are in a special phase only about 40 Myr after the second encounter and 50 Myr before their final collision. This is the only phase in the simulation when a gas-rich overlap region between the nuclei is forming accompanied by enhanced star formation. The star formation rate as well as the recent star formation history in the central region agree well with observational estimates. For the first time this new model explains the distributed extra-nuclear star formation in the Antennae galaxies as a consequence of the recent second encounter. The proposed model predicts that the Antennae are in a later merger stage than the Mice (NGC 4676) and would therefore lose their first place in the classical Toomre sequence.

The disruptive effect of galactic tides is a textbook example of gravitational dynamics. However, depending on the shape of the potential, tides can also become fully compressive. When that is the case, they might trigger or strengthen the formation of galactic substructures (star clusters, tidal dwarf galaxies), instead of destroying them. We perform N-body simulations of interacting galaxies to quantify this effect. We demonstrate that tidal compression occurs repeatedly during a galaxy merger, independently of the specific choice of parameterization. With a model tailored to the Antennae galaxies, we show that the distribution of compressive tides matches the locations and timescales of observed substructures. After extending our study to a broad range of parameters, we conclude that neither the importance of the compressive tides (~15% of the stellar mass) nor their duration (~ 10 Myr) are strongly affected by changes in the progenitors' configurations and orbits. Moreover, we show that individual clumps of matter can enter compressive regions several times in the course of a simulation. We speculate that this may spawn multiple star formation episodes in some star clusters, through e.g., enhanced gas retention.

We point out a strong time-evolution of the mass-to-light conversion factor eta commonly used to estimate masses of unresolved star clusters from observed cluster spectro-photometric measures. We present a series of gas-dynamical models coupled with the Cambridge stellar evolution tracks to compute line-of-sight velocity dispersions and half-light radii weighted by the luminosity. We explore a range of initial conditions, varying in turn the cluster mass and/or density, and the stellar population's IMF. We find that eta, and hence the estimated cluster mass, may increase by factors as large as 3 over time-scales of 50 million years. We apply these results to an hypothetic cluster mass distribution function (d.f.) and show that the d.f. shape may be strongly affected at the low-mass end by this effect. Fitting truncated isothermal (Michie-King) models to the projected light profile leads to over-estimates of the concentration parameter c of delta c ~ 0.3 compared to the same functional fit applied to the projected mass density.

We use smoothed particle hydrodynamics (SPH) models to study the large-scale morphology and dynamical evolution of the intergalactic gas in Stephan's Quintet, and compare to multiwavelength observations. Specifically, we model the formation of the hot X-ray gas, the large-scale shock, and emission line gas as the result of NGC 7318b colliding with the group. We also reproduce the N-body model of Renaud and Appleton for the tidal structures in the group.

Physics lectures always refer to the tides as a disruptive effect. However, tides can also be compressive. When the potential of two galaxies overlap, as happens during a merger, fully compressive tides can develop and have a strong impact on the dynamics of substructures such as star clusters or tidal dwarf galaxies. Using N-body simulations of a large set of mergers, we noticed the importance of these tidal modes at cluster scale. With a model of the Antennae galaxies, we conclude that the positions and timescales of these tidal modes match the actual distribution of young clusters. A detailed study of the statistics of the compressive tides shows a stunning correlation between this purely gravitational effect and the observed properties of the star clusters. In this contribution, we introduce the concept of compressive tide and show its relevance in the special case of the Antennae galaxies. We extend our conclusions to a broad range of parameters and discuss their implications on several critical points such as the infant mortality, multiple star formation epochs in clusters or the age distribution.

The tidal field of galaxies is known generally to be disruptive. However, in the case of galaxy mergers, a compressive mode of tidal wave may develop and last long enough to cocoon the formation of star clusters. Using an N-body simulation of the Antennae galaxies, we derive the positions of these compressive regions and the statistics of their duration. Excellent agreement between the spatial distribution of tides and observed young clusters is found, while the characteristic e-folding times of 10 to 30 Myrs derived for the tidal field compare well with cluster formation time-scales.

Gravitational tides are widely understood to strip and destroy galactic substructures. In the course of a galaxy merger, however, transient totally compressive tides may develop and prevent star forming regions from dissolving, after they condensed to form clusters of stars. We study the statistics of such compressive modes in an N-body model of the galaxy merger NGC 4038/39 (the Antennae) and show that ~15% of the disc material undergoes compressive tides at pericentre. The spatial distribution of observed young clusters in the overlap and nuclear regions of the Antennae matches surprisingly well the location of compressive tides obtained from simulation data. Furthermore, the statistics of time intervals spent by individual particles embedded in a compressive tide yields a log-normal distribution of characteristic time ~10 Myr, comparable to star cluster formation timescales. We argue that this generic process is operative in galaxy mergers at all redshifts and possibly enhances the formation of star clusters. We show with a model calculation that this process will prevent the dissolution of a star cluster during the formation phase, even for a star formation efficiency as low as ~10%. The transient nature of compressive tides implies that clusters may dissolve rapidly once the tidal field switches to the usual disruptive mode.