In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from observations of stellar velocity dispersion, gas rotation curves, weak and strong gravitational lensing, and/or X-ray maps. We make our comparison in terms of the maximal surface density of haloes, circumventing the need for semi-analytic or parametric models for dark matter density profiles. We find that the maximal surface density as a function of halo mass is well reproduced by simulations of CDM that include baryons, and that a SIDM cross-section of 1 cm$^2/$g is marginally excluded. This constraint is driven mainly by galaxy clusters, with less massive systems being relatively less affected by DM self-interactions with velocity-independent cross-section.

Using a set of cosmological hydro-dynamical simulations, we constrained the properties of primordial magnetic fields by studying their impact on the formation and evolution of dwarf galaxies. We performed a large set of simulations (8 dark matter only and 72 chemo-hydrodynamical) including primordial magnetic fields through the extra density fluctuations they induce at small length scales ($k \geq 10\,h\,\rm{Mpc^{-1}}$) in the matter power spectrum. We explored a large variety of primordial magnetic fields with strength $B_\lambda$ ranging from $0.05$ to $0.50\,\textrm{nG}$ and magnetic energy spectrum slopes $n_B$ from $-2.9$ to $-2.1$. Strong magnetic fields characterized by a high amplitude ($B_\lambda=0.50,\,0.20\,\textrm{nG}$ with $n_B=-2.9$) or by a steep initial power spectrum slope ($n_B=-2.1,-2.4$, with $B_\lambda=0.05\,\textrm{nG}$) induce perturbations in the mass scales from $10^7$ to $10^9\,\rm{M}_{\odot}$. In this context emerging galaxies see their star formation rate strongly boosted. They become more luminous and metal rich than their counterparts without primordial magnetic fields. Such strong fields are ruled out by their inability to reproduce the observed scaling relations of dwarf galaxies. They predict dwarf galaxies to be at the origin of an unrealistically early reionization of the Universe and also overproduce luminous satellites in the Local Group. Weaker magnetic fields impacting the primordial density field at corresponding masses $\lesssim 10^6\,\rm{M}_{\odot}$, produce a large number of mini dark halos orbiting the dwarfs, however out of reach for current lensing observations. This study allows for the first time to constrain the properties of primordial magnetic fields based on realistic cosmological simulations of dwarf galaxies.

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

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.

In the Local Group, quenched gas-poor dwarfs galaxies are most often found close to the Milky Way and Andromeda, while star forming gas-rich ones are located at greater distances. This so-called morphology-density relation is often interpreted as the consequence of the ram pressure stripping of the satellites during their interaction with the Milky Way hot halo gas. While this process has been often investigated, self-consistent high resolution simulations were still missing. In this study, we have analysed the impact of both the ram pressure and tidal forces induced by a host galaxy on dwarf models as realistic as possible emerging from cosmological simulations. These models were re-simulated using both a wind tunnel and a moving box technique. The secular mass growth of the central host galaxy, as well as the gas density and temperature profiles of its hot halo have been taken into account. We show that while ram pressure is very efficient at stripping the hot and diffuse gas of the dwarf galaxies, it can remove their cold gas ($T < 10^3$~[K]) only in very specific conditions. Depending on the infall time of the satellites relatively to the build-up stage of the massive host, star formation can thus be prolonged instead of being quenched. This is the direct consequence of the clumpy nature of the cold gas and the thermal pressure the hot gas exerts onto it. We discuss the possibility that the variety in satellite populations among Milky Way-like galaxies reflects their accretion histories.

We modify the chemo-dynamical code GEAR to simulate the impact of self-interacting dark matter on the observable quantities of 19 low mass dwarf galaxies with a variety star forming properties. We employ a relatively high, velocity independent cross-section of $\sigma/m = 10$cm$^2$/g and extract, in addition to integrated quantities, the total mass density profile, the luminosity profile, the line-of-sight velocities, the chemical abundance and the star formation history. We find that despite the creation of large cores at the centre of the dark matter haloes, the impact of SIDM on the observable quantities of quenched galaxies is indiscernible, dominated mostly by the stochastic build up of the stellar matter. As such we conclude that it is impossible to make global statements on the density profile of dwarf galaxies from single or small samples. Although based mostly on quenched galaxies, this finding supports other recent work putting into question the reliability of inferred cored density profiles that are derived from observed line-of-sight velocities.

We present the results of a set of high resolution chemo-dynamical simulations of dwarf galaxies in a $\Lambda$CDM cosmology. Out of an original 3.4 Mpc$^3$/h$^3$ cosmological box, a sample of 27 systems are zoomed-in from z=70 to z=0. Gas and stellar properties are confronted to the observations in the greatest details: in addition to the galaxy global properties, we investigate the model galaxy velocity dispersion profiles, half-light radii, star formation histories, metallicity distributions, and [Mg/Fe] abundance ratios. The formation and sustainability of the metallicity gradients and kinetically distinct stellar populations are also tackled. We show how the properties of six Local Group dwarf galaxies, NGC 6622, Andromeda II, Sculptor, Sextans, Ursa Minor and Draco are reproduced, and how they pertain to three main galaxy build-up modes. Our results indicate that the interaction with a massive central galaxy could be needed for a handful of Local Group dwarf spheroidal galaxies only, the vast majority of the systems and their variety of star formation histories arising naturally from a $\Lambda$CDM framework. We find that models fitting well the local Group dwarf galaxies are embedded in dark haloes of mass between $5\times 10^8$ to a few $10^9\,\rm{M_\odot}$, without any missing satellite problem. We confirm the failure of the abundance matching approach at the mass scale of dwarf galaxies. Some of the observed faint however gas-rich galaxies with residual star formation, such as Leo T and Leo P, remain challenging. They point out the need of a better understanding of the UV-background heating.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

Chemo-dynamical N-body simulations are an essential tool for understanding the formation and evolution of galaxies. As the number of observationally determined stellar abundances continues to climb, these simulations are able to provide new constraints on the early star formaton history and chemical evolution inside both the Milky Way and Local Group dwarf galaxies. Here, we aim to reproduce the low $\alpha$-element scatter observed in metal-poor stars. We first demonstrate that as stellar particles inside simulations drop below a mass threshold, increases in the resolution produce an unacceptably large scatter as one particle is no longer a good approximation of an entire stellar population. This threshold occurs at around $10^3\,\rm{M_\odot}$, a mass limit easily reached in current (and future) simulations. By simulating the Sextans and Fornax dwarf spheroidal galaxies we show that this increase in scatter at high resolutions arises from stochastic supernovae explosions. In order to reduce this scatter down to the observed value, we show the necessity of introducing a metal mixing scheme into particle-based simulations. The impact of the method used to inject the metals into the surrounding gas is also discussed. We finally summarise the best approach for accurately reproducing the scatter in simulations of both Local Group dwarf galaxies and in the Milky Way.

We present the high resolution spectroscopic study of five -3.9<=[Fe/H]<=-2.5 stars in the Local Group dwarf spheroidal, Sculptor, thereby doubling the number of stars with comparable observations in this metallicity range. We carry out a detailed analysis of the chemical abundances of alpha, iron peak, light and heavy elements, and draw comparisons with the Milky Way halo and the ultra faint dwarf stellar populations. We show that the bulk of the Sculptor metal-poor stars follows the same trends in abundance ratios versus metallicity as the Milky Way stars. This suggests similar early conditions of star formation and a high degree of homogeneity of the interstellar medium. We find an outlier to this main regime, which seems to miss the products of the most massive of the TypeII supernovae. In addition to its value to help refining galaxy formation models, this star provides clues to the production of cobalt and zinc. Two of our sample stars have low odd-to-even barium isotope abundance ratios, suggestive of a fair proportion of s-process; we discuss the implication for the nucleosynthetic origin of the neutron capture elements.

Tidal dwarf galaxies (TDGs) are recycled objects that form within the collisional debris of interacting/merging galaxies. They are expected to be devoid of non-baryonic dark matter, since they can form only from dissipative material ejected from the discs of the progenitor galaxies. We investigate the gas dynamics in a sample of six bona-fide TDGs around three interacting and post-interacting systems: NGC 4694, NGC 5291, and NGC 7252 ("Atoms for Peace"). For NGC 4694 and NGC 5291 we analyse existing HI data from the Very Large Array (VLA), while for NGC 7252 we present new HI observations from the Jansky VLA together with long-slit and integral-field optical spectroscopy. For all six TDGs, the HI emission can be described by rotating disc models. These HI discs, however, have undergone less than a full rotation since the time of the interaction/merger event, raising the question of whether they are in dynamical equilibrium. Assuming that these discs are in equilibrium, the inferred dynamical masses are consistent with the observed baryonic masses, implying that TDGs are devoid of dark matter. This puts constraints on putative "dark discs" (either baryonic or non-baryonic) in the progenitor galaxies. Moreover, TDGs seem to systematically deviate from the baryonic Tully-Fisher relation. These results provide a challenging test for alternative theories like MOND.

As galaxy simulations increase in resolution more attention is being paid towards the evolution of dwarf galaxies and how the simulations compare to observations. Despite this increasing resolution we are however, far away from resolving the interactions of satellite dwarf galaxies and the hot coronae which surround host galaxies. We describe a new method which focuses only on the local region surrounding an infalling dwarf in an effort to understand how the hot baryonic halo will alter the chemodynamical evolution of dwarf galaxies. Using this method we examine how a dwarf, similar to Sextans dwarf spheroidal, evolves in the corona of a Milky Way like galaxy. We find that even at high perigalacticons the synergistic interaction between ram pressure and tidal forces transform a dwarf into a stream, suggesting that Sextans was much more massive in the past in order survive its perigalacticon passage. In addition the large confining pressure of the hot corona allows gas that was originally at the outskirts to begin forming stars, initially forming stars of low metallicity compared to the dwarf evolved in isolation. This increase in star formation eventually allows a dwarf galaxy to form more metal rich stars compared to one in isolation, but only if the dwarf retains gas for a sufficiently long period of time. In addition, dwarfs which formed substantial numbers of stars post-infall will have a slightly elevated [Mg/Fe] at high metallicity ([Fe/H] -1.5).

The damping on the fluctuation spectrum and the presence of thermal velocities as properties of warm dark matter particles like sterile neutrinos imprint a distinct signature found from the structure formation mechanisms to the internal structures of halos. Using warm dark matter simulations we explore these effects on the structure formation for different particle energies and we find that the formation of structure is more complex than originally assumed, a combination of top-down collapse and hierarchical (bottom-up) clustering on multiple scales. The degree on which one scenario is more prominent with respect to the other depends globally on the energy of the particle and locally on the morphology and architecture of the analyzed region. The presence of shells and caustics in warm dark matter halos is another important effect seen in simulations. Furthermore, we discuss the impact of thermal velocities on the structure formation from theoretical considerations as well as from the analysis of the simulations. We re-examine the assumptions considered when estimating the velocity dispersion for warm dark matter particles that have been adopted in previous works for more than a decade and we give an independent estimation for the velocities. We identify some inconsistencies in previous published results. The relation between the warm dark matter particle mass and its corresponding velocity dispersion is strongly model dependent, hence the constraints on particle mass from simulation results are weak. Finally, we review the technical difficulties that arise in warm dark matter simulations along with possible improvements of the methods.

The effect of the local environment on the evolution of dwarf spheroidal galaxies is poorly understood. We have undertaken a suite of simulations to investigate the tidal impact of the Milky Way on the chemodynamical evolution of dwarf spheroidals that resemble present day classical dwarfs using the SPH code GEAR. After simulating the models through a large parameter space of potential orbits the resulting properties are compared with observations from both a dynamical point of view, but also from the, often neglected, chemical point of view. In general, we find that tidal effects quench the star formation even inside gas-endowed dwarfs. Such quenching, may produce the radial distribution of dwarf spheroidals from the orbits seen within large cosmological simulations. We also find that the metallicity gradient within a dwarf is gradually erased through tidal interactions as stellar orbits move to higher radii. The model dwarfs also shift to higher $\langle$[Fe/H]$\rangle$/L ratios, but only when losing $>$$20\%$ of stellar mass.

Using a suite of detailed numerical simulations we estimate the level of anisotropy generated by the time evolution along the light cone of the 21cm signal from the epoch of reionization. Our simulations include the physics necessary to model the signal during both the late emission regime and the early absorption regime, namely X-ray and Lyman-band 3D radiative transfer in addition to the usual dynamics and ionizing UV transfer. The signal is analysed using correlation functions perpendicular and parallel to the line of sight (LOS). We reproduce general findings from previous theoretical studies: the overall amplitude of the correlations and the fact that the light cone anisotropy is visible only on large scales (100 cMpc). However, the detailed behaviour is different. At 3 different epochs, the amplitude of the correlations along and perpendicular to the LOS differ from each other, indicating anisotropy. These 3 epochs are associated with 3 events of the global reionization history: the overlap of ionized bubbles, the onset of mild heating by X-rays in regions around the sources, and the onset of efficient Lyman-alpha coupling in regions around the sources. A 20x20 deg^2 survey area may be necessary to mitigate sample variance when we use the directional correlation functions. On a 100 cMpc scale the light cone anisotropy dominates over the anisotropy generated by peculiar velocity gradients computed in the linear regime. By modelling instrumental noise and limited resolution, we find that the anisotropy should be easily detectable by the SKA, assuming perfect foreground removal, the limiting factor being a large enough survey size. In the case of the LOFAR, it is likely that only first anisotropy episode will fall in the observing frequency range and will be detectable only if sample variance is much reduced (i.e. a larger than 20x20 deg^2 survey, which is not presently planned).

The larger number of models of asteroid shapes and their rotational states derived by the lightcurve inversion give us better insight into both the nature of individual objects and the whole asteroid population. With a larger statistical sample we can study the physical properties of asteroid populations, such as main-belt asteroids or individual asteroid families, in more detail. Shape models can also be used in combination with other types of observational data (IR, adaptive optics images, stellar occultations), e.g., to determine sizes and thermal properties. We use all available photometric data of asteroids to derive their physical models by the lightcurve inversion method and compare the observed pole latitude distributions of all asteroids with known convex shape models with the simulated pole latitude distributions. We used classical dense photometric lightcurves from several sources and sparse-in-time photometry from the U.S. Naval Observatory in Flagstaff, Catalina Sky Survey, and La Palma surveys (IAU codes 689, 703, 950) in the lightcurve inversion method to determine asteroid convex models and their rotational states. We also extended a simple dynamical model for the spin evolution of asteroids used in our previous paper. We present 119 new asteroid models derived from combined dense and sparse-in-time photometry. We discuss the reliability of asteroid shape models derived only from Catalina Sky Survey data (IAU code 703) and present 20 such models. By using different values for a scaling parameter cYORP (corresponds to the magnitude of the YORP momentum) in the dynamical model for the spin evolution and by comparing synthetics and observed pole-latitude distributions, we were able to constrain the typical values of the cYORP parameter as between 0.05 and 0.6.

We present accurate time delays for the quadruply imaged quasar HE 0435-1223. The delays were measured from 575 independent photometric points obtained in the R-band between January 2004 and March 2010. With seven years of data, we clearly show that quasar image A is affected by strong microlensing variations and that the time delays are best expressed relative to quasar image B. We measured Delta_t(BC) = 7.8+/-0.8 days, Delta_t(BD) = -6.5+/-0.7 days and Delta_t_CD = -14.3+/-0.8 days. We spacially deconvolved HST NICMOS2 F160W images to derive accurate astrometry of the quasar images and to infer the light profile of the lensing galaxy. We combined these images with a stellar population fitting of a deep VLT spectrum of the lensing galaxy to estimate the baryonic fraction, $f_b$, in the Einstein radius. We measured f_b = 0.65+0.13-0.10 if the lensing galaxy has a Salpeter IMF and f_b = 0.45+0.04-0.07 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, sigma_ap = 222+/-34 km/s. We used f_b and sigma_ap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solve the Jeans equations numerically for the model and explored the parameter space under the additional requirement that the model must predict the correct astrometry for the quasar images. Given the current error bars on f_b and sigma_ap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with chi^2 < 1. However, narrowing this range is possible, provided a better velocity dispersion measurement becomes available. In addition, increasing the depth of the current HST imaging data of HE 0435-1223 will allow us to combine our constraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object.

We present a fully parallel chemo-dynamical Tree/SPH code, GEAR, which allows to perform high resolution simulations with detailed chemical diagnostics. Starting from the public version of Gadget-2, we included the complex treatment of the baryon physics: gas cooling, star formation law, chemical evolution and supernovae feedback. We qualified the performances of GEAR with the case of dSph galaxies. GEAR conserves the total energy budget of the systems to better than 5% over 14Gyr and proved excellent convergence of the results with numerical resolution. We showed that models of dSphs in a static Euclidean space, where the expansion of the universe is neglected are valid. In addition, we tackled some of the existing open questions in the field, like the stellar mass fraction of dSphs and its link with the predicted dark matter halo mass function, the effect of the supernova feedback, the spatial distribution of the stellar populations, and the origin of the diversity in star formation histories and chemical abundance patterns. Strong supernovae driven winds seem incompatible with the observed metallicities and luminosities. Despite the fact that newly formed stars are preferentially found in the galaxy central parts, turbulent motions in the gas can quickly erase any metallicity gradient. The variety in dSph properties result from a range of total masses as well as from a dispersion in central densities. The latter is also seen in the haloes emerging from a LCDM cosmogony.

It is predicted that sources emitting UV radiation in the Lyman band during the epoch of reionization (EoR) showed a series of discontinuities in their Ly-alpha flux radial profile as a consequence of the thickness of the Lyman line series in the primeval intergalactic medium. Through unsaturated Wouthuysen-Field coupling, these spherical discontinuities are also present in the 21 cm emission of the neutral IGM. In this article, we study the effects these discontinuities have on the differential brightness temperature of the 21 cm signal of neutral hydrogen in a realistic setting including all other sources of fluctuations. We focus on the early phases of the EoR, and we address the question of the detectability by the planned Square Kilometre Array. Such a detection would be of great interest, because these structures could provide an unambiguous diagnostic for the cosmological origin of the signal remaining after the foreground cleaning procedure. Also, they could be used as a new type of standard rulers. We determine the differential brightness temperature of the 21 cm signal in the presence of inhomogeneous Wouthuysen-Field effect using simulations which include (hydro)dynamics and both ionizing and Lyman lines 3D radiative transfer with the code LICORICE. We find that the Lyman horizons are clearly visible on the maps and radial profiles around the first sources of our simulations, but for a limited time interval, typically \Delta z \approx 2 at z \sim 13. Stacking the profiles of the different sources of the simulation at a given redshift results in extending this interval to \Delta z \approx 4. When we take into account the implementation and design planned for the SKA (collecting area, sensitivity, resolution), we find that detection will be challenging. It may be possible with a 10 km diameter for the core, but will be difficult with the currently favored design of a 5 km core.

We present the first detection of CO emission lines in the Halpha filaments at distances as far as 50 kpc from the centre of the galaxy NGC 1275. This gas is probably dense (>=10E3 cm-3). However, it is not possible to accurately determine the density and the kinetic temperature of this relatively warm gas (Tkin~20-500K) with the current data only. The amount of molecular gas in the filaments is large 10E9 Msun (assuming a Galactic N(H2)/Ico ratio). This is 10% of the total mass of molecular gas detected in this cD galaxy. This gas has large-scale velocities comparable to those seen in Halpha. The origin of the filaments is still unclear, but their formation is very likely linked to the AGN positive feedback (Revaz et al., 2008) that regulates the cooling of the surrounding X-ray-emitting gas as suggested by numerical simulations. We also present high-resolution spectra of the galaxy core. The spatial characteristics of the double-peaked profile suggest that the molecular web of filaments and streamers penetrates down to radii of less than 2 kpc from the central AGN and eventually feed the galaxy nucleus. The mass of gas inside the very central region is ~10E^9 Msun, and is similar to the mass of molecular gas found in the filaments.

We present simulations of the 21-cm signal during the epoch of reionization. We focus on properly modeling the absorption regime in the presence of inhomogeneous Wouthuysen-Field effect and X-ray heating. We ran radiative transfer simulations for three bands in the source spectrum (Lyman, UV, and X-ray) to fully account for these processes. We find that the brightness temperature fluctuation of the 21 cm signal has an amplitude greater than 100 mK during the early reionization, up to 10 times greater than the typical amplitude of a few 10 mK obtained during the later emission phase. More importantly, we find that even a rather high contribution from QSO-like sources only damps the absorption regime without erasing it. Heating the IGM with X-ray takes time. Our results show that observations of the early reionization will probably benefit from a higher signal-to-noise value than during later stages. After analyzing the statistical properties of the signal (power spectrum and PDF) we find three diagnostics to constrain the level of X-ray, hence the nature of the first sources.

When the surface angular velocity is above about 70% of the critical angular velocity, many interesting features appear which may be tested by interferometric observations, like significant deformation of stars, variation of the effective temperature with the latitude. Also polar winds become important and equatorial disks may appear. Near the critical limit, convection is also favored in the outer layers. In the present paper, we emphasize the need for a proper estimate of the critical velocity since this is the ratio of the actual velocity of the star to that critical velocity which determines the amplitude of the above effects. We recall the existence of two critical velocities. The first one, also called the classical critical velocity is the one to consider when the star has an Eddington factor inferior to 0.639, while the second one is the one to be considered when the Eddington factor is above 0.639. The features of the star at these two critical limits may be very different.

Compact elliptical galaxies are characterized by small sizes and high stellar densities. They are thought to form through tidal stripping of massive progenitors. However, only a handful of them were known, preventing us from understanding the role played by this mechanism in galaxy evolution. We present a population of 21 compact elliptical galaxies gathered with the Virtual Observatory. Follow-up spectroscopy and data mining, using high-resolution images and large databases, show that all the galaxies exhibit old metal-rich stellar populations different from those of dwarf elliptical galaxies of similar masses but similar to those of more massive early-type galaxies, supporting the tidal stripping scenario. Their internal properties are reproduced by numerical simulations, which result in compact dynamically hot remnants resembling the galaxies in our sample.

The near proportionality between HI and dark matter in outer galactic disks prompted us to run N-body simulations of galactic disks in which the observed gas content is supplemented by a dark gas component representing between zero and five times the visible gas content. While adding baryons in the disk of galaxies may solve some issues, it poses the problem of disk stability. We show that the global stability is ensured if the ISM is multiphased, composed of two partially coupled phases, a visible warm gas phase and a weakly collisionless cold dark phase corresponding to a fraction of the unseen baryons. The phases are subject to stellar and UV background heating and gas cooling, and their transformation into each other is studied as a function of the coupling strength. This new model, which still possesses a dark matter halo, fits the rotation curves as well as the classical CDM halos, but is the only one to explain the existence of an open and contrasting spiral structure, as observed in the outer HI disks

We present a large sample of fully self-consistent hydrodynamical Nbody/Tree-SPH simulations of isolated dwarf spheroidal galaxies (dSphs). It has enabled us to identify the key physical parameters and mechanisms at the origin of the observed variety in the Local Group dSph properties. The initial total mass (gas + dark matter) of these galaxies is the main driver of their evolution. Star formation (SF) occurs in series of short bursts. In massive systems, the very short intervals between the SF peaks mimic a continuous star formation rate, while less massive systems exhibit well separated SF bursts, as identified observationally. The delay between the SF events is controlled by the gas cooling time dependence on galaxy mass. The observed global scaling relations, luminosity-mass and luminosity-metallicity, are reproduced with low scatter. We take advantage of the unprecedentedly large sample size and data homogeneity of the ESO Large Programme DART, and add to it a few independent studies, to constrain the star formation history of five Milky Way dSphs, Sextans, LeoII, Carina, Sculptor and Fornax. For the first time, [Mg/Fe] vs [Fe/H] diagrams derived from high-resolution spectroscopy of hundreds of individual stars are confronted with model predictions. We find that the diversity in dSph properties may well result from intrinsic evolution. We note, however, that the presence of gas in the final state of our simulations, of the order of what is observed in dwarf irregulars, calls for removal by external processes.

The 21 cm emission of neutral hydrogen is the most promising probe of the epoch of reionization(EoR). In the next few years, the SKA pathfinders will provide statistical measurements of this signal. Numerical simulations predicting these observations are necessary to optimize the design of the instruments. The main difficulty is the computation of the spin temperature of neutral hydrogen which depends on the gas kinetic temperature and on the level of the local Lyman-alpha flux. A T_s >> T_cmb assumption is usual. However, this assumption does not apply early in the reionization history, or even later in the history as long as the sources of X-rays are too weak to heat the intergalactic medium significantly. This work presents the first EoR numerical simulations including, beside dynamics and ionizing continuum radiative transfer, a self-consistent treatment of the Ly-alpha radiative transfer. This allows us to compute the spin temperature more accurately. We use two different box sizes, 20 Mpc/h and 100 Mpc/h, and a star source model. Using the redshift dependence of average quantities, maps, and power spectra, we quantify the effect of using different assumptions to compute the spin temperature and the influence of the box size. The first effect comes from allowing for a signal in absorption. The magnitude of this effect depends on the amount of heating by hydrodynamic shocks and X-rays in the intergalactic medium(IGM). The second effects comes from using the real, local, Lyman-alpha flux. This effect is important for an average ionization fraction of less than 10%: it changes the overall amplitude of the 21 cm signal, and adds its own fluctuations to the power spectrum.

We have recently detected CO lines in the well-known filaments around NGC 1275, the galaxy at the centre of the Perseus cluster of galaxies. These previous observations, with the HERA multi-beam array at the IRAM 30m telescope enabled us to make a large map of the CO(2-1) line and to see hints of molecular gas far away from the cluster centre. To confirm the presence of CO emission lines in the outer filaments and to study the CO(2-1)/CO(1-0) line ratio, we observed seven regions of interest again with the 30m telescope in both CO(1-0) and CO(2-1). The regions we observed were: the eastern filament, the horseshoe, the northern filament and a southern extension, all selected from Halpha emission line mapping. Molecular gas is detected in all the observed regions. This result confirms the large extent of the cold molecular gas filaments. We discuss the CO(2-1)/CO(1-0) ratios in the filaments. The eastern filament has optically thick gas, whereas further away, the line ratio increases close to values expected for a warmer optically thin medium. We also show CO(1-0) and CO(2-1) lines in 9 regions closer to the centre. The kinematics of the CO is studied here in more detail and confirms that it follows the motions of the warm H_2 gas found in the near-infrared. Finally, we searched for dense gas tracers around 3C84 and claim here the first detection of HCN(3-2).

We recently found extended CO(2-1) emission from cold molecular gas embedded in the network of Halpha filaments surrounding the galaxy NGC 1275 (Salome et al. 2006). We now present CO(2-1) interferometer maps of the eastern filaments, at high spatial and spectral resolutions. The cold molecular gas is detected by the Plateau de Bure Interferometer along the eastern filaments over an extent of 15'', or with a projected length of 5kpc. In our 2.5'' beam, the main CO filament is mostly unresolved along its minor axis. The multiple peaks along the CO filaments and the low values of the observed CO brightness temperatures imply further unresolved structures that may be giant molecular clouds. These clouds have very narrow line-width emission lines (~30 km/s). The CO emission is optically thick. It very likely traces cold clouds bound under their own self-gravity that may be falling back in the gravitational potential well of the galaxy. Such a picture would agree with current models of ``positive feedback'' in which some of the hot gas around NGC 1275 (a) is trapped by buoyantly rising bubbles inflated by the energy input of the 3C84 AGN, (b) subsequently cools efficiently at a larger radius around the edges of the hot bubbles, and (c) then falls back in self-gravitating clouds of molecular gas toward the center of the galaxy.

$\Lambda$CDM numerical simulations predict that the "missing baryons" reside in a Warm-Hot gas phase in the over-dense cosmic filaments. However, there are now several theoretical and observational arguments that support the fact that galactic disks may be more massive than usually thought, containing a substantial fraction of the "missing baryons". Hereafter, we present new N-body simulations of galactic disks, where the gas content has been multiplied by a factor 5. The stability of the disk is ensured by assuming that the ISM is composed out of two partially coupled phases, a warm phase, corresponding the observed CO and HI gas and a cold collisionless phase corresponding to the unseen baryons.

Emission-lines in the form of filamentary structures is common in bright clusters characterized by short cooling times. In the Perseus cluster, cold molecular gas, tightly linked to the H$\alpha$ filaments, has been recently revealed by CO observations. In order to understand the origin of these filamentary structures, we have investigated the evolution of the hot ICM gas perturbed by the AGN central activity in a Perseus like cluster. Using very-high resolution TreeSPH simulations combined with a multiphase model and a model of plasma bubbles, we have been able to follow the density and temperature evolution of the disturbed ICM gas around the bubbles. Our simulations show that a fraction of the $1-2 \rm{keV}$ gas present at the center of clusters is trapped and entrained by the rising buoyant bubble to higher radius where the AGN heating is less efficient. The radiative cooling makes it cool in a few tens of Myr below $10^4 \rm{K}$, forming cold filamentary structures in the wake and in the rim of the bubbles.

New CO observations at the center of cooling flow in Perseus cluster \cite{salome06} show a clear correlation of the molecular gas with the previously detected H-$\alpha$ filaments \cite{conselice01}. In this poster, we present high resolution multi-phase simulations of the Perseus Cluster, taking into account the AGN feedback in form of hot buoyant bubbles. These simulations show that significant amount of gas can cool far from the center. The AGN feedback provides some heating, but also trigger the hot gas compression, that favours cooling (positive feedback), even at high radius ($R>30 \rm{kpc}$). The cooled gas flows into the cluster core forming the observed filaments.

Galactic warps represent an old unresolved problem, since the discovery, at the end of the fifties, of the HI warp of the Milky Way. In this paper, we propose a new scenario explaining a large fraction of the observed optical warps. Based on N-body simulations, we show that realistic galactic disks, where the dark matter is essentially distributed in a disk, are subject to bending instabilities. S, U-shaped, as well as asymmetric warps are spontaneously generated and in some cases are long-lived. While this scenario presents the advantage of explaining the three observed types of warps, it also brings new constraints on the dark matter distribution in spiral galaxies. Finally, it gives us a unified picture of galaxies where galactic asymmetries, like bars, spirals and warps result from gravitational instabilities.

Based on N-body simulations, we show that realistic galactic disks are subject to bending instabilities of fire-hose type when the disks are substantially self-gravitating, that is, if they contain dark matter distributed in the disk. Depending on the degree of instability, S and U-shaped, as well as asymmetric warps are generated. In some cases, the warp may last several galactic rotations, particularly when the instability is marginal. Since the bending instability is very sensitive to the disk flattening, the fractions of dark matter distributed in the disk and in the dark halo are constrained. For a Milky Way like galaxy the extended dark halo can not exceed 30-40% of the total mass within 35 kpc if the Milky Way warp results from a bending instability. This mode of warping provides a unified picture of spiral galaxies, where bars, spiral arms and warps result all from disk gravitational instabilities, radial or transverse, which are constantly regenerated by the dissipative gas component.

The Baryonic Tully-Fisher relation (BTF) can be substantially improved when considering that the galactic baryonic mass is likely to consist not only from the detected baryons, stars and gas, but also from a dark baryonic component proportional to the HI gas. The BTF relation is optimally improved when the HI mass is multiplied by a factor of about 3, but larger factors up to 11-16 still improve the fit over the original one using only the detected baryons. The strength of this improved relation is quantified with up-to-date statistical tests such as the Akaike Information Criterion or the Bayesian Information Criterion. In particular they allow to show that supposing a variable $M_\star/L$ ratio instead is much less significant. This result reinforces the suggestion made in several recent works that mass within galactic disks must be a multiple of the HI mass, and that galactic disks are substantially, but not necessarily fully, self-gravitating.

Bars and spiral arms have played an important role as constraints on the dynamics and on the distribution of dark matter in the optical parts of disk galaxies. Dynamics linked to the dissipative nature of gas, and its transformation into stars provide clues that spiral galaxies are driven by dissipation close to a state of \textit{marginal stability} with respect to the dynamics in the galaxy plane. Here we present numerical evidences that warps play a similar role but in the transverse direction. N-body simulations show that typical galactic disks are also marginally stable with respect to a bending instability leading to typical observed warps. The frequent occurrence of warps and asymmetries in the outer galactic disks give therefore, like bars in the inner disks, new constraints on the dark matter, but this time in the outer disks. If disks are marginally stable with respect to bending instabilities, our models suggest that the mass within the HI disks must be a multiple of the detected HI and stars, i.e., disks must be heavier than seen. But the models do not rule out a traditional thick halo with a mass within the HI disk radius similar to the total disk mass.

A substantial fraction of the warps in spiral galaxies may result from bending instabilities if the disks are essentially self-gravitating. With N-body simulations, we show that galaxies with self-gravitating disks as thick as HI disks are subject to bending instabilities generating S-shaped, U-shaped or asymmetric warps. S-shaped warps persist during several rotations and keep the line of node straight. The warp amplitudes generated by bending instabilities remain however modest. Other factors must be invoked for extreme warped disks. However, bending instabilities can account for most of the cases reported in optical surveys, where the warp angle is generally less than $5^\circ$. This mode of warping is very sensitive to the disk flattening. It also constrains the fraction of dark matter distributed in the disk and in the dark halo.

It is often assumed that a warped galaxy can be modeled by a set of rings. This paper verifies numerically the validity of this assumption by the study of periodic orbits populating a heavy self-gravitating warped disk. The phase space structure of a warped model reveals that the circular periodic orbits of a flat disk are transformed in quasi annular periodic orbits which conserve their stability. This lets us also explore the problem of the persistence of a large outer warp. In particular, the consistency of its orbits with the density distribution is checked as a function of the pattern speed.

Two methods generating warped galaxies with N-body simulations are presented. One uses an external potential as a disturber while the other is based on material accretion. The results of both methods are compared. A particular attention is given on the shape of the line of node (LON).