Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

Within the next few years, the Square Kilometer Array (SKA) or one of its pathfinders will hopefully provide a detection of the 21-cm signal fluctuations from the Epoch of Reionization (EoR). Then, the main goal will be to accurately constrain the underlying astrophysical parameters. Currently, this is mainly done with Bayesian inference using Markov Chain Monte Carlo sampling. Recently, studies using neural networks trained to performed inverse modelling have shown interesting results. We build on these by improving the accuracy of the predictions using neural network and exploring other supervised learning methods: the kernel and ridge regressions. Based on a large training set of 21-cm power spectra, we compare the performances of these supervised learning methods. When using an un-noised signal as input, we improve on previous neural network accuracy by one order of magnitude and, using local ridge kernel regression, we gain another factor of a few. We then reach a rms prediction error of a few percents of the 1-sigma confidence level due to SKA thermal noise (as estimated with Bayesian inference). This last performance level requires optimizing the hyper-parameters of the method: how to do that perfectly in the case of an unknown signal remains an open question. For an input signal altered by a SKA-type thermal noise, our neural network recovers the astrophysical parameter values with an error within half of the 1$\sigma$ confidence level due to the SKA thermal noise. This accuracy improves to 10$\%$ of the 1$\sigma$ level when using the local ridge kernel regression (with optimized hyper-parameters). We are thus reaching a performance level where supervised learning methods are a viable alternative to determine the best-fit parameters values.

With a statistical detection of the 21 cm signal fluctuations from the Epoch of Reionization (EoR) expected in the next few years, there is an interest in developing robust and precise techniques to constrain the underlying astrophysical parameters. Bayesian inference with Markov Chain Monte Carlo, or different types of supervised learning for backward modelling (from signal to parameters) are examples of such techniques. They usually require many instances of forward modelling (from parameters to signal) in sampling the parameters space, either when performing the steps of the Markov Chain or when building a training sample for supervised learning. As forward modelling can be costly (if performed with numerical simulations for example), we should attempt to perform an optimal sampling according to some principle. With this goal in mind, we present an approach based on defining a metric on the space of observables, induced by the manner through which the modelling creates a mapping from the parameter space onto the space of observables. This metric bears a close connection to Jeffreys' prior from information theory. It is used to generate a homogeneous and isotropic sampling of the signal space with two different methods. We show that when the resulting optimized samplings, created with 21cmFAST, are used to train a neural network we obtain a modest reduction of the error on parameter reconstruction of ~10% (compared to a na\"ive sampling of the same size). Excluding the borders of the parameter space region, the improvement is more substantial, on the order of 30-40%.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

Quasars contribute to the 21-cm signal from the Epoch of Reionization (EoR) primarily through their ionizing UV and X-ray emission. However, their radio continuum and Lyman-band emission also regulates the 21-cm signal in their direct environment, potentially leaving the imprint of their duty cycle. We develop a model for the radio and UV luminosity functions of quasars from the EoR, and constrain it using recent observations. Our model is consistent with the z~7.5 quasar from Banados et al 2017, and also predicts only a few quasars suitable for 21-cm forest observations (10mJy) in the sky. We exhibit a new effect on the 21-cm signal observed against the CMB: a radio-loud quasar can leave the imprint of its duty cycle on the 21-cm tomography. We apply this effect in a cosmological simulation and conclude that the effect of typical radio-loud quasars is most likely negligible in an SKA field of view. For a 1-10mJy quasar the effect is stronger though hardly observable at SKA resolution. Then we study the contribution of the lyman band Ly-alpha to Ly-beta) emission of quasars to the Wouthuisen-Field coupling. The collective effect of quasars on the 21-cm power spectrum is larger than the thermal noise at low k, though featureless. However, a distinctive pattern around the brightest quasars in an SKA field of view may be observable in the tomography, encoding the duration of their duty cycle. This pattern has a high signal-to-noise ratio for the brightest quasar in a typical SKA shallow survey.

We present a novel Lyman alpha (Ly$\alpha$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$\alpha$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$\alpha$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$\alpha$ photons in the interstellar medium of young star-forming galaxies like Ly$\alpha$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$\alpha$ emission, which can be compared to the observations of LAEs.

We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spatial, i.e. overlapping, disc populations with different scaleheights, kinematics and metallicities. In this model the thick, metal-poor, and centrally concentrated disc populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient, which is mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced -- as shown already in previous works in the literature. However, as we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions of the Bulge seen in the data. On the other hand, the model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p. For this scenario to reproduce the observations, the $\alpha$-enhanced, metal-poor, thick disc populations must have a significant mass contribution in the inner regions -- as has been suggested for the Milky Way.

The 21-cm signal from the Epoch of Reionization (EoR) is expected to be detected in the next few years, either with existing instruments or by the upcoming SKA and HERA projects. In this context there is a pressing need for publicly available high-quality templates covering a wide range of possible signals. These are needed both for end-to-end simulations of the up-coming instruments, as well as to develop signal analysis methods. In this work we present such a set of templates, publicly available, for download at https://21ssd.obspm.fr/. The database contains 21-cm brightness temperature lightcones at high and low resolution, and several derived statistical quantities for 45 models spanning our choice of 3D parameter space. These data are the result of fully coupled radiative hydrodynamic high resolution ($1024^3$) simulations performed with the LICORICE code. Both X-ray and Lyman line transfer is performed to account for heating and Wouthuysen-Field coupling fluctuations. We also present a first exploitation of the data using the power spectrum and the Pixel Distribution Function (PDF) as functions of redshifts, computed from lightcone data. We analyse how these two quantities behave when varying the model parameters while taking into account the thermal noise expected of a typical SKA survey. Finally, we show that the power spectrum and the PDF have different -- and to some extent complementary -- abilities to distinguish between different models. This opens the door to formulating an optimal sampling of the parameter space, dependant on the chosen diagnostics.

We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio almost half that of the thick disc bar. The thin disc b/p bulge has a pronounced X-shape, while the thick disc b/p is weaker with a rather boxy shape. This leads to the signature of the b/p bulge in the thick disc to be weaker and further away from the plane than in the thin disc. Regarding the kinematics, we find that the los velocity of thick disc stars in the outer parts of the b/p bulge can be \emph{larger} than that of thin disc stars, by up to 40\% and 20\% for side-on and Milky Way-like orientations of the bar respectively. This is due to the different orbits followed by thin and thick disc stars in the bar-b/p region, which are affected by the fact that: i) thin disc stars are trapped more efficiently in the bar - b/p instability and thus lose more angular momentum than their thick disc counterparts and ii) thick disc stars have large radial excursions and therefore stars from large radii with high angular momenta can be found in the bar region. We also find that the difference between the los velocities of the thin and thick disc in the b/p bulge ($\Delta v_{los}$) correlates with the initial difference between the radial velocity dispersions of the two discs ($\Delta \sigma$) . We therefore conclude that stars in the bar - b/p bulge will have considerably different morphologies and kinematics depending on the kinematic properties of the disc population they originate from.

The 21 cm signal from the Epoch of Reionization should be observed within the next decade. While a simple statistical detection is expected with SKA pathfinders, the SKA will hopefully produce a full 3D mapping of the signal. To extract from the observed data constraints on the parameters describing the underlying astrophysical processes, inversion methods must be developed. For example, the Markov Chain Monte Carlo method has been successfully applied. Here we test another possible inversion method: artificial neural networks (ANN). We produce a training set which consists of 70 individual sample. Each sample is made of the 21 cm power spectrum at different redshifts produced with the 21cmFast code plus the value of three parameters used in the semi-numerical simulations that describe astrophysical processes. Using this set we train the network to minimize the error between the parameter values it produces as an output and the true values. We explore the impact of the architecture of the network on the quality of the training. Then we test the trained network on the new set of 54 test samples with different values of the parameters. We find that the quality of the parameter reconstruction depends on the sensitivity of the power spectrum to the different parameters at a given redshift, that including thermal noise and sample variance decreases the quality of the reconstruction and that using the power spectrum at several redshifts as an input to the ANN improves the quality of the reconstruction. We conclude that ANNs are a viable inversion method whose main strength is that they require a sparse exploration of the parameter space and thus should be usable with full numerical simulations.

The $\Lambda$CDM cosmological scenario predicts that our Galaxy should contain hundreds of stellar streams at the solar vicinity, fossil relics of the merging history of the Milky Way and more generally of the hierarchical growth of galaxies. Because of the mixing time scales in the inner Galaxy, it has been claimed that these streams should be difficult to detect in configuration space but can still be identifiable in kinematic-related spaces like the energy/angular momenta spaces, E-Lz and Lperp-Lz, or spaces of orbital/velocity parameters. By means of high-resolution, dissipationless N-body simulations, containing between 25$\times10^6$ and 35$\times10^6$ particles, we model the accretion of a series of up to four 1:10 mass ratio satellites then up to eight 1:100 satellites and we search systematically for the signature of these accretions in these spaces. In all spaces considered (1) each satellite gives origin to several independent overdensities; (2) overdensities of multiple satellites overlap; (3) satellites of different masses can produce similar substructures; (4) the overlap between the in-situ and the accreted population is considerable everywhere; (5) in-situ stars also form substructures in response to the satellite(s) accretion. These points are valid even if the search is restricted to kinematically-selected halo stars only. As we are now entering the 'Gaia era', our results warn that an extreme caution must be employed before interpreting overdensities in any of those spaces as evidence of relics of accreted satellites. Reconstructing the accretion history of our Galaxy will require a substantial amount of accurate spectroscopic data, that, complemented by the kinematic information, will possibly allow us to (chemically) identify accreted streams and measure their orbital properties. (abridged)

The 21-cm forest is a promising probe of the Epoch of Reionization. The local state of the intergalactic medium (IGM) is encoded in the spectrum of a background source (radio-loud quasars or gamma ray burst afterglow) by absorption at the local 21-cm wavelength, resulting in a continuous and fluctuating absorption level. Small-scale structures (filaments and minihaloes) in the IGM are responsible for the strongest absorption features. The absorption can also be modulated on large scales by inhomogeneous heating and Wouthuysen-Field coupling. We present the results from a simulation that attempts to preserve the cosmological environment while resolving some of the small-scale structures (a few kpc resolution in a 50 Mpc/h box). The simulation couples the dynamics and the ionizing radiative transfer and includes X-ray and Lyman lines radiative transfer for a detailed physical modelling. As a result we find that soft X-ray self-shielding, Lyman-alpha self-shielding and shock heating all have an impact on the predicted values of the 21-cm optical depth of moderately overdense structures like filaments. An correct treatment of the peculiar velocities is also critical. Modelling these processes seems necessary for accurate predictions and can be done only at high enough resolution. As a result, based on our fiducial model, we estimate that LOFAR should be able to detect a few (strong) absorptions features in a frequency range of a few tens of MHz for a 20 mJy source located at z=10, while the SKA would extract a large fraction of the absorption information for the same source.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

The local intensity of the 21 cm signal emitted during the Epoch of Reionization that will be mapped by the SKA is modulated by the amount of neutral hydrogen. Consequently, understanding the process of reionization of the intergalactic medium (IGM) is crucial for predicting and interpreting the upcoming observations. After presenting the basic physics and most meaningful quantities pertaining to the process of reionization, we will review recent progress in our understanding of the production and escape of ionizing photons in primordial galaxies and of their absorption in the IGM especially in so-called minihalos and Lyman Limit Systems.

By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as some studies have recently suggested. On the basis of several arguments emphasized in this paper, we propose that the bulge population which, in the Milky Way, is observed not to be part of the peanut structure corresponds to the old galactic thick disk, thus implying that the Milky Way is a pure thin+thick disk galaxy, with only a possible limited contribution of a classical bulge.

By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013).

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 Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

We explore the impacts of ultraviolet (UV) radiation feedback on galaxies during the epoch of reionisation by cosmological simulations in which hydrodynamics and the transfer of the H and He ionising photons are consistently coupled. Moreover we take into account H_2 non-equilibrium chemistry, including photo-dissociation. The most striking feature of the simulations is a high spatial resolution for the radiative transfer (RT) calculation which enables us to start considering not only external UV feedback processes but also internal UV feedback processes in each galaxy. We find that the star formation is significantly suppressed due to the internal UV and supernova (SN) feedback. In low mass galaxies with M<10^9Msun, a large amount of gas is evacuated by photo-evaporation as previous studies have shown, which results in the suppression of star formation. Surprisingly, star formation in massive halos is also strongly suppressed despite the fact that these halos hardly lose any gas by photo-evaporation. The suppression of star formation in massive halos is mainly caused by following two factors; (i) small scale clumpy structures in the galaxies are smoothened by the internal feedback, (ii) although the dense gas in the galaxies is mostly neutral, the H_2 formation and cooling processes are disturbed by mild photo-heating. Photo-dissociating radiations actually suppress star formation, but the magnitude of the effect is not so large in massive galaxies. Even though our simulation volume is too small to be a representative patch of the Universe during reionisation, we find that our simulated star formation rate densities and HI fractions at z=6-7 are consistent with those found in observations.

By means of N-body simulations, we show that radial migration in galaxy disks, induced by bar and spiral arms, leads to significant azimuthal variations in the metallicity distribution of old stars at a given distance from the galaxy center. Metals do not show an axisymmetric distribution during phases of strong migration. Azimuthal variations are visible during the whole phase of strong bar phase, and tend to disappear as the effect of radial migration diminishes, together with a reduction in the bar strength. These results suggest that the presence of inhomogeneities in the metallicity distribution of old stars in a galaxy disk can be a probe of ongoing strong migration. Such signatures may be detected in the Milky Way by Gaia (and complementary spectroscopic data), as well as in external galaxies, by IFU surveys like CALIFA and ATLAS3D. Mixing - defined as the tendency toward a homogeneous, azimuthally symmetric, stellar distribution in the disk - and migration turns out to be two distinct processes, the effects of mixing starting to be visible when strong migration is over.

First, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter haloes of ~10^4M_sun can produce first luminous objects with the aid of dark matter cusps. Therefore, the mass of first objects is smaller by roughly two orders of magnitude than in the previous prediction. This implies that the number of Pop III stars formed in the early universe could be significantly larger than hitherto thought. Secondly, the feedback by photo-ionization and photo-dissociation photons in the first objects is explored with radiation hydrodynamic simulations, and it is demonstrated that multiple stars can form in a 10^5M_sun halo. Thirdly, the fragmentation of an accretion disk around a primordial protostar is explored with photo-dissociation feedback. As a result, it is found that the photo-dissociation can reduce the mass accretion rate onto protostars. Also, protostars as small as 0.8M_sun may be ejected and evolve with keeping their mass, which might be detected as "real first stars" in the Galactic halo. Finally, state-of-the-art radiation hydrodynamic simulations are performed to investigate the internal ionization of first galaxies and the escape of ionizing photons. We find that UV feedback by forming massive stars enhances the escape fraction even in a halo as massive as > 6* 10^9M_sun, while it reduces the star formation rate significantly. This may have a momentous impact on the cosmic reionization.

Galaxies accrete their mass by means of both smooth accretion from the cosmic web, and the mergers of smaller entities. We wish to quantify the respective role of these two modes of accretion, which could determine the morphological types of galaxies observed today. Multi-zoom cosmological simulations are used to estimate as a function of time the evolution of mass in bound systems, for dark matter as well as baryons. The baryonic contents of dark matter haloes are studied. Merger histories are followed as a function of external density, and the different ways in which mass is assembled in galaxies and the stellar component accumulated are quantified. We find that most galaxies assemble their mass through smooth accretion, and only the most massive galaxies also grow significantly through mergers. The mean fraction of mass assembled by accretion is 77 %, and by mergers 23 %. We present typical accretion histories of hundreds of galaxies: masses of the most massive galaxies increase monotonically in time, mainly through accretion, many intermediate-mass objects also experience mass-loss events such as tidal stripping and evaporation. However, our simulations suffer from the overcooling of massive galaxies caused by the neglect of active galaxy nuclei (AGN) feedback. The time by which half of the galay mass has assembled, both in dark matter and baryons, is a decreasing function of mass, which is compatible with the observations of a so-called downsizing. At every epoch in the universe, there are low-mass galaxies actively forming stars, while more massive galaxies form their stars over a shorter period of time within half the age of the universe.

We present a novel method to investigate cosmic reionization, using joint spectral information on high redshift Lyman Alpha Emitters (LAE) and quasars (QSOs). Although LAEs have been proposed as reionization probes, their use is hampered by the fact their Ly{\alpha} line is damped not only by intergalactic HI but also internally by dust. Our method allows to overcome such degeneracy. First, we carefully calibrate a reionization simulation with QSO absorption line experiments. Then we identify LAEs in two simulation boxes at z=5.7 and z=6.6 and we build synthetic images/spectra of a prototypical LAE. At redshift 5.7, we find that the Ly{\alpha} transmissivity (T_LAE) ~ 0.25, almost independent of the halo mass. This constancy arises from the conspiracy of two effects: (i) the intrinsic Ly{\alpha} line width and (ii) the infall peculiar velocity. At higher redshift, z=6.6, where the transmissivity is instead largely set by the local HI abundance and LAE transmissivity consequently increases with halo mass from 0.15 to 0.3. Although outflows are present, they are efficiently pressure-confined by infall in a small region around the LAE; hence they only marginally affect transmissivity. Finally, we cast LOS originating from background QSOs passing through foreground LAEs at different impact parameters, and compute the quasar transmissivity (T_QSO). At smaller impact parameters, d < 1 cMpc, a positive correlation between T_QSO and halo mass is found at z = 5.7, which tends to become less pronounced (i.e. flatter) at larger distances. By cross-correlating T_LAE and T_QSO, we can obtain a HI density estimate unaffected by dust. At z= 5.7, the cross-correlation is relatively weak,whereas at z = 6.6 we find a clear positive correlation. We conclude by briefly discussing the perspectives for the application of the method to existing and forthcoming data.

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.

It has been known since many decades that galaxy interactions can induce star formation (hereafter SF) enhancements and that one of the driving mechanisms of this enhancement is related to gas inflows into the central galaxy regions, induced by asymmetries in the stellar component, like bars. In the last years many evidences have been accumulating, showing that interacting pairs have central gas-phase metallicities lower than those of field galaxies, by {\sim} 0.2-0.3 dex on average. These diluted ISM metallicities have been explained as the result of inflows of metal-poor gas from the outer disk to the galaxy central regions. A number of questions arises: What's the timing and the duration of this dilution? How and when does the SF induced by the gas inflow enrich the circumnuclear gas with re-processed material? Is there any correlation between the timing and strength of the dilution and the timing and intensity of the SF? By means of Tree-SPH simulations of galaxy major interactions, we have studied the effect that gas inflows have on the ISM dilution, and the effect that the induced SF has, subsequently, in re-enriching the nuclear gas. In this contribution, we present the main results of this study.

According to the hierarchical model, small galaxies form first and merge together to form bigger objects. In parallel, galaxies assemble their mass through accretion from cosmic filaments. Recently, the increased spatial resolution of the cosmological simulations have emphasised that a large fraction of cold gas can be accreted by galaxies. In order to compare the role of both phenomena and the corresponding star formation history, one has to detect the structures in the numerical simulations and to follow them in time, by building a merger tree.

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.

Strong inflows of gas from the outer disk to the inner kiloparsecs are induced during the interaction of disk galaxies. This inflow of relatively low-metallicity gas dilutes the metallicity of the circumnuclear gas. We have investigated several aspects of the process as the timing and duration of the dilution and its correlation with the induced star formation. We analysed major (1:1) gas-rich interactions and mergers, spanning a range of initial orbital characteristics. Star formation and metal enrichment from SNe are included in our model. Our results show that the strongest trend is between the star formation rate and the dilution of the metals in the nuclear region; i.e., the more intense the central burst of star formation, the more the gas is diluted. This trend comes from strong inflows of relatively metal-poor gas from the outer regions of both disks, which fuels the intense star formation and lowers the overall metallicity for a time. The strong inflows happen on timescales of about 10^8 years or less, and the most intense star formation and lowest gas phase metallicities are seen generally after the first pericentre passage. As the star formation proceeds and the merger advances, the dilution reduces and enrichment becomes dominant - ultimately increasing the metallicity of the circumnuclear gas to a level higher than the initial metallicities of the merging galaxies. The "fly-bys" - pairs that interact but do not merge - also cause some dilution. We even see some dilution early in the merger or in the "fly-bys" and thus do not observe a strong trend between the nuclear metallicities and separation in our simulations until the merger is well advanced. We also analyse the O and Fe enrichment of the ISM, and show that the evolution of the alpha/Fe ratios, as well as the dilution of the central gas metallicity, can be used as a clock for "dating" the interaction.

We present the GalMer database, a library of galaxy merger simulations, made available to users through tools compatible with the Virtual Observatory (VO) standards adapted specially for this theoretical database. To investigate the physics of galaxy formation through hierarchical merging, it is necessary to simulate galaxy interactions varying a large number of parameters: morphological types, mass ratios, orbital configurations, etc. On one side, these simulations have to be run in a cosmological context, able to provide a large number of galaxy pairs, with boundary conditions given by the large-scale simulations, on the other side the resolution has to be high enough at galaxy scales, to provide realistic physics. The GalMer database is a library of thousands simulations of galaxy mergers at moderate spatial resolution and it is a compromise between the diversity of initial conditions and the details of underlying physics. We provide all coordinates and data of simulated particles in FITS binary tables. The main advantages of the database are VO access interfaces and value-added services which allow users to compare the results of the simulations directly to observations: stellar population modelling, dust extinction, spectra, images, visualisation using dedicated VO tools. The GalMer value-added services can be used as virtual telescope producing broadband images, 1D spectra, 3D spectral datacubes, thus making our database oriented towards the usage by observers. We present several examples of the GalMer database scientific usage obtained from the analysis of simulations and modelling their stellar population properties, including: (1) studies of the star formation efficiency in interactions; (2) creation of old counter-rotating components; (3) reshaping metallicity profiles in elliptical galaxies; (4) orbital to internal angular momentum transfer; (5) reproducing observed colour bimodality of galaxies.

Many evidences exist for a connection between galaxy interactions and induced star formation. However, a large range of responses of galaxies to tidal interactions is found, both in observations and in numerical simulations. We will discuss some recent results obtained analysing a large sample (~ 1000) of simulations of interacting pairs and their agreement with the most recent observational works.

Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure supported progenitor galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with v/sigma > 1 outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the outer regions acquire part of the angular momentum first; (3) both the baryonic and the dark matter components of the remnant galaxy acquire part of the angular momentum, the relative fractions depend on the initial concentration of the merging galaxies. For this conversion to occur the initial baryonic component must be sufficiently dense and/or the encounter should take place on a orbit with high angular momentum. Systems with these hybrid properties have been recently observed through a combination of stellar absorption lines and planetary nebulae for kinematic studies of early-type galaxies. Our results are in qualitative agreement with such observations and demonstrate that even mergers composed of non-rotating, pressure-supported progenitor galaxies can produce early-type galaxies with significant rotation at large radii.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

Using a large suite of galaxies with a variety of concentrations and metallicity gradients, we study the evolution of non-dissipative ("dry") equal mass mergers. Our purpose in generating these simulations is to explore how the metallicity gradient in dry mergers depends on the structure and metallicity gradients of the galaxies involved in the merger. Specifically, we would like to answer: Could dry mergers lead to metallicity gradients as observed in elliptical galaxies in the local Universe? Do dry mergers always lead to a flattening of the initial (i.e., pre-merger) gradient? From this modeling, we conclude that: The ratio of the remnant and the initial galaxy slopes span a wide range of values, up to values greater than 1 (with values greater than one resulting only when companions have gradients twice the progenitor). For a merger between two ellipticals having identical initial metallicity slopes (i.e., equal companion and galaxy slopes), the metallicity profile of the remnant flattens, with a final gradient about 0.6 times the initial one. Ellipticals can maintain their original pre-merger metallicity gradient if the companion slope is sufficiently steep. The amount of flattening neither depends on the characteristics of the orbit of the progenitors or on their initial concentration. Given the diversity in outcomes of the mergers, we conclude that dry mergers do not violate any observational constraints on the systematic characteristics of metallicity gradients in local ellipticals. In fact, dry mergers may be important within the context of the results of our simulations and may explain the large scatter in gradients for massive ellipticals and the relative paucity of massive ellipticals with no or shallow metallicity gradients.

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 investigate the intensity enhancement and the duration of starburst episodes, triggered by major galaxy interactions and mergers. To this aim, we analyze two large statistical datasets of numerical simulations. These have been obtained using two independent and different numerical techniques to model baryonic and dark matter evolution, that are extensively compared for the first time. One is a Tree-SPH code, the other one is a grid-based N-body sticky-particles code. We show that, at low redshift, galaxy interactions and mergers in general trigger only moderate star formation enhancements. Strong starbursts where the star formation rate is increased by a factor larger than 5 are rare and found only in about 15% of major galaxy interactions and mergers. Merger-driven starbursts are also rather short-lived, with a typical duration of the activity of a few 10^8 yr. These conclusions are found to be robust, independent from the numerical techniques and star formation models. At higher redshifts where galaxies contain more gas, gas inflow-induced starbursts are neither stronger neither longer than their local counterparts. In turn, the formation of massive gas clumps, results of local Jeans instability that can occur spontaneously in gas-rich disks or be indirectly favored by galaxy interactions, could play a more important role in determining the duration and intensity of star formation episodes.

We investigate, by means of numerical simulations, the kinematics of elliptical-spiral merger remnants. Counterrotation can appear both in coplanar and in non-coplanar retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the external regions of the two interacting galaxies acquire part of the orbital angular momentum, due to the action of tidal forces.

We investigate, by means of numerical simulations, the possibility of forming counter-rotating old stellar components by major mergers between an elliptical and a spiral galaxy. We show that counter-rotation can appear both in dissipative and dissipationless retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the external regions of the two interacting galaxies acquire part of the orbital angular momentum, due to the action of tidal forces exerted on each galaxy by the companion.

We have undertaken a large set of simulations of galaxy interactions and mergers (GalMer Project) in order to study the physical processes related to galaxy encounters. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct and retrograde, according to the initial relative energy and impact parameter. The self gravity of stars, gas and dark matter is taken into account through a tree-code algorithm, the gas hydrodynamics through SPH. Star formation is included adopting a density-dependent Schmidt law. This wide library of galaxy interactions and mergers, containing, at present, about 900 simulations of major encounters, represents an unique tool to investigate statistically the chemodynamical evolution of interacting systems. In the following, we present and discuss some results obtained exploring the dataset, together with some future perspectives.

During the epoch of reionization, Ly-alpha photons emitted by the first stars can couple the neutral hydrogen spin temperature to the kinetic gas temperature, providing the opportunity to observe the gas in emission or absorption in the 21-cm line. Given the bright foregrounds, it is of prime importance to determine precisely the fluctuations signature of the signal, to be able to extract it by its correlation power. LICORICE is a Monte-Carlo radiative transfer code, coupled to the dynamics via an adaptative Tree-SPH code. We present here the Ly-alpha part of the implementation, and validate it through three classical tests. Contrary to previous works, we do not assume that P_alpha, the number of scatterings of Ly-alpha photons per atom per second, is proportional to the Ly-alpha background flux, but take into account the scatterings in the Ly-alpha line wings. The latter have the effect to steepen the radial profile of P_alpha around each source, and re-inforce the contrast of the fluctuations. In the particular geometry of cosmic filaments of baryonic matter, Ly-alpha photons are scattered out of the filament, and the large scale structure of P_alpha is significantly anisotropic. This could have strong implications for the possible detection of the 21-cm signal.

We investigate the enhancement of star formation efficiency in galaxy interactions and mergers, by numerical simulations of several hundred galaxy collisions. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct and retrograde, according to the initial relative energy and impact parameter, and the resulting star formation history is compared to that occuring in the two galaxies when they are isolated. Our principal results are: (1) retrograde encounters have a larger star formation efficiency (SFE) than direct encounters; (2) the amount of gas available in the galaxy is not the main parameter governing the SFE in the burst phase; (3) there is an anticorrelation between the amplitude of the star forming burst and the tidal forces exerted per unit of time, which is due to the large amount of gas dragged outside the galaxy by tidal tails in strong interactions; (4) globally, the Kennicutt-Schmidt law is retrieved statistically for isolated galaxies, interacting pairs and mergers; (5) the enhanced star formation is essentially occurring in nuclear starbursts, triggered by inward gas flows driven by non-axisymmetries in the galaxy disks. Direct encounters develop more pronounced asymmetries than retrograde ones. Based on these statistical results, we derive general laws for the enhancement of star formation in galaxy interactions and mergers, as a function of the main parameters of the encounter.

In this work we focus on the properties of accretion onto galaxies. Through numerical simulations we investigate the geometrical properties of accretion. To span the scale range required in these simulations we have developed a new numerical technique: the multi-zoom method. We run a series of Tree-SPH simulations in smaller and smaller boxes at higher and higher mass resolution, using data recorded at the previous level to account for the matter inflow and the tidal field from outside matter. The code is parallelized using OpenMP. We present a validation test to evaluate the robustness of the method: the pancake collapse. We apply this new multizoom method to study the accretion properties. Zooming in onto galaxies from a cosmological simulation, we select a sample of 10 well resolved galaxies (5000 baryonic particles or more). We sum up their basic properties and plot a Tully-Fisher relation. We find that smooth accretion of intergalactic cold gas dominates mergers for the mass growth of galaxies at z < 2. Next we study the baryonic accretion rate which shows different behaviours depending on the galaxy mass. The bias is also computed at different radii and epochs. Then we present galactocentric angular maps for the accretion integrated between z=2 and z=0, which reveal that accretion is highly anisotropic. Average accretion rates plotted against galactocentric latitude show a variety of behaviours. In general, accretion in the galactic plane is favored, even more for baryonic matter than for dark matter.

Bars in gas-rich spiral galaxies are short-lived. They drive gas inflows through their gravity torques, and at the same time self-regulate their strength. Their robustness has been subject of debate, since it was thought that only the resulting central mass concentrations (CMCs) were weakening bars, and only relatively rare massive CMCs were able to completely destroy them. Through numerical simulations including gas dynamics, we find that with the gas parameters of normal spiral galaxies, the CMC is not sufficient to fully dissolve the bar. But another overlooked mechanism, the transfer of angular momentum from the infalling gas to the stellar bar, can also strongly weaken the bar. In addition, we show that gravity torques are correctly reproduced in simulations, and conclude that bars are transient features, with life-time of 1-2 Gyr in typical Sb-Sc galaxies, because of the combined effects of CMCs and gravity torques, while most existing works had focussed on the CMC effects alone.

We study through numerical simulation the spherical collapse of isothermal gas in Newtonian gravity. We observe a critical behavior which occurs at the threshold of gravitational instability leading to core formation. For a given initial density profile, we find a critical temperature, which is of the same order as the virial temperature of the initial configuration. For the exact critical temperature, the collapse converges to a self-similar form, the first member in Hunter's family of self-similar solutions. For a temperature close to the critical value, the collapse first approaches this critical solution. Later on, in the supercritical case, the collapse converges to another self-similar solution, which is called the Larson-Penston solution. In the subcritical case, the gas bounces and disperses to infinity. We find two scaling laws: one for the collapsed mass in the supercritical case and the other for the maximum density reached before dispersal in the subcritical case. The value of the critical exponent is measured to be $\simeq 0.11$ in the supercritical case, which agrees well with the predicted value $\simeq 0.10567$. These critical properties are quite similar to those observed in the collapse of a radiation fluid in general relativity. We study the response of the system to temperature fluctuation and discuss astrophysical implications for the insterstellar medium structure and for the star formation process. Newtonian critical behavior is important not only because it provides a simple model for general relativity but also because it is relevant for astrophysical systems such as molecular clouds.

The formation and evolution of galactic disks are complex phenomena, where gas and star dynamics are coupled through star formation and the related feedback. The physical processes are so numerous and intricate that numerical models focus, in general, on one or a few of them only. We propose here a numerical model with particular attention to the multiphase nature of the interstellar medium; we consider a warm gas phase (> 10^4 K), treated as a continuous fluid by an SPH algorithm, and a cold gas phase (down to 10K), fragmented in clouds, treated by a low-dissipation sticky particles component. The two gas phases do not have the same dynamics, nor the same spatial distribution. In addition to gravity, they are coupled through mass exchanges due to heating/cooling processes, and supernovae feedback. Stars form out of the cold phase, and re-inject mass to the warm phase through SN explosions and stellar winds. The baryons are embedded in a live cold dark matter component. Baryonic disks, initially composed of pure gas, encounter violent instabilities, and a rapid phase of star formation, that slows down exponentially. Stars form in big clumps, that accumulate in the center to build a bulge. Exponential metallicity gradients are obtained. External infall of gas should be included to maintain a star formation rate in the disk comparable to what is observed in present disk galaxies.

It is likely that cold molecular clouds form at high redshift, before galaxies. Considering these cold clouds instead of hot homogeneous gas as the main baryonic component of proto-galactic halos may affect several issues during galaxy formation. In particular, the baryonic matter loses angular momentum through dynamical friction on the dark matter halo. In numerical simulations using hot gas for baryonic matter, this gives rise to the so-called angular momentum problem. In this work, we study the dynamical friction exerted on cold fractal gas by a collisionless background (dark matter) through high-resolution numerical simulations. First, we find that, for values of the parameters relevant during galaxy formation, the friction does not alter the morphology of the fractal, which is mainly driven by internal dynamics. Then, we show that the presence of substructures and inhomogeneities in a body has little effect on the global value of the friction. Parameters such as branching ratio and fractal dimension also have little effect. In fact, we find out that the main effect comes from the deformations and fluctuations of the fractal structure in cold gas clumps. If the deformation time is of the same order of magnitude or shorter than the typical build-up time of the friction, the friction is weakened. We argue that this effect is relevant for galaxy formation, and that the angular momentum problem should not be worsened by including the small scale inhomogeneities of the cold gas which occur at a resolution out of reach of present simulations.

This work studies the hydrodynamics of self-gravitating compressible isothermal fluids. We show that the hydrodynamic evolution equations in absence of viscosity are scale covariant. We study the evolution of the time dependent fluctuations around singular and regular isothermal spheres. We linearize the fluid equations around such stationary solutions and apply Laplace transform to solve them. We find that the system is stable below a critical size (X ~ 9.0 in dimensionless variables) and unstable above; this size is the same critical size found in the study of the thermodynamical stability in the canonical ensemble and associated to a center-to-border density ratio of 32.1 . We prove that the value of this critical size is independent of the Reynolds number of the system. Furthermore, we give a detailed description of the series of successive dynamical instabilities that appear at higher and higher sizes following the geometric progression X_n ~ 10.7^n. We turn then to study exact solutions of the hydrodynamic equations without viscosity and we provide analytic and numerical axisymmetric soliton-type solutions. The stability of exact solutions corresponding to a collapsing filament is studied by computing linear fluctuations. Radial fluctuations growing faster than the background are found for all sizes of the system. However, a critical size (X ~ 4.5) appears, separating a weakly from a strongly unstable regime.

The interstellar medium is observed in a hierarchical fractal structure over several orders of magnitude in scale. Aiming to understand the origin of this structure, we carry out numerical simulations of molecular cloud fragmentation, taking into account self-gravity, dissipation and energy input. Self-gravity is computed through a tree code, with fully or quasi periodic boundary conditions. Energy dissipation is introduced through cloud-cloud ineslatic collisions. Several schemes are tested for the energy input. It appears that energy input from galactic shear allows to achieve a stationary clumped state for the gas, avoiding final collapse. When a stationary turbulent cascade is established, it is possible to derive meaningful statistical studies on the data such as the fractal dimension of the mass distribution.

The self-gravitating thermal gas (non-relativistic particles of mass m at temperature T) is exactly equivalent to a field theory with a single scalar field phi(x) and exponential self-interaction. We build up perturbation theory around a space dependent stationary point phi_0(r) in a finite size domain delta \leq r \leq R ,(delta << R), which is relevant for astrophysical applica- tions (interstellar medium,galaxy distributions).We compute the correlations of the gravitational potential (phi) and of the density and find that they scale; the latter scales as 1/r^2. A rich structure emerges in the two-point correl- tors from the phi fluctuations around phi_0(r). The n-point correlators are explicitly computed to the one-loop level.The relevant effective coupling turns out to be lambda=4 pi G m^2 / (T R). The renormalization group equations (RGE) for the n-point correlator are derived and the RG flow for the effective coupling lambda(tau) [tau = ln(R/delta), explicitly obtained.A novel dependence on tau emerges here.lambda(tau) vanishes each time tau approaches discrete values tau=tau_n = 2 pi n/sqrt7-0, n=0,1,2, ...Such RG infrared stable behavior [lambda(tau) decreasing with increasing tau] is here connected with low density self-similar fractal structures fitting one into another.For scales smaller than the points tau_n, ultraviolet unstable behaviour appears which we connect to Jeans' unstable behaviour, growing density and fragmentation. Remarkably, we get a hierarchy of scales and Jeans lengths following the geometric progression R_n=R_0 e^{2 pi n /sqrt7} = R_0 [10.749087...]^n . A hierarchy of this type is expected for non-spherical geometries,with a rate different from e^{2 n/sqrt7}.