Four expressions involving sums of position and velocity coordinates bounding the total angular momentum of particle systems, and by extension of any continuous or discontinuous material systems, are derived which are tighter for any particle configuration than similar inequalities derived by Sundman (1913), Saari (2005) and Scheeres (2012). Eight distinct inequalities can thus be ordered according to their tightness to angular momentum.

The cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme submillimetre galaxies at high redshift. Here we search for GMCs in a typical Milky Way progenitor at z = 1.036. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we mapped the CO(4-3) emission of this gravitationally lensed galaxy at high resolution, reading down to 30 parsecs, which is comparable to the resolution of CO observations of nearby galaxies. We identify 17 molecular clouds, characterized by masses, surface densities and supersonic turbulence all of which are 10-100 times higher than present-day analogues. These properties question the universality of GMCs and suggest that GMCs inherit their properties from ambient interstellar medium. The measured cloud gas masses are similar to the masses of stellar clumps seen in the galaxy in comparable numbers. This corroborates the formation of molecular clouds by fragmentation of distant turbulent galactic gas disks, which then turn into stellar clumps ubiquitously observed in galaxies at cosmic noon.

Based on a high quality $N$-body simulation of a double bar galaxy model, we investigate the evolution of the bar properties, including their size, strength and instantaneous pattern speed derived by using three distinct methods: the Fourier, Jacobi integral, and moment of inertia methods. The interaction of the two bars, which rotate at distinct speeds, primarily affects the size, strength and pattern speed of the inner bar. When the two bars are perpendicular to each other, the size and the pattern speed of the inner bar decrease and its strength increases. The emergence of a strong Fourier $m=1$ mode increases the oscillation amplitude of the size, strength and pattern speed of the inner bar. On the other hand, the characteristics of the outer bar are substantially influenced by its adjacent spiral structure. When the spiral structure disappears, the size of the outer bar increases and its strength and pattern speed decrease. Consequently, the ratio of the pattern speed of the outer bar with respect to the inner bar is not constant and increases with time. Overall, the double bar and disk system displays substantial high frequency semi-chaotic fluctuations of the pattern strengths and speeds both in space and time, superposed on the slow secular evolution, which invalidates the assumption that the actions of individual stars should be well conserved in barred galaxies, such as the Milky Way.

Condensation of H 2 in the interstellar medium (ISM) has long been seen as a possibility, either by deposition on dust grains or thanks to a phase transition combined with self-gravity. H 2 condensation might explain the observed low efficiency of star formation and might help to hide baryons in spiral galaxies. Our aim is to quantify the solid fraction of H 2 in the ISM due to a phase transition including self-gravity for different densities and temperatures in order to use the results in more complex simulations of the ISM as subgrid physics. We used molecular dynamics simulations of fluids at different temperatures and densities to study the formation of solids. Once the simulations reached a steady state, we calculated the solid mass fraction, energy increase, and timescales. By determining the power laws measured over several orders of magnitude, we extrapolated to lower densities the higher density fluids that can be simulated with current computers. The solid fraction and energy increase of fluids in a phase transition are above 0.1 and do not follow a power law. Fluids out of a phase transition are still forming a small amount of solids due to chance encounters of molecules. The solid mass fraction and energy increase of these fluids are linearly dependent on density and can easily be extrapolated. The timescale is below one second, the condensation can be considered instantaneous. The presence of solid H 2 grains has important dynamic implications on the ISM as they may be the building blocks for larger solid bodies when gravity is included. We provide the solid mass fraction, energy increase, and timescales for high density fluids and extrapolation laws for lower densities.

The effective potential neighboring the corotation resonance region in barred galaxies is shown to be strongly time-dependent in any rotating frame because of the competition of nearby perturbations of similar strengths with differing rotation speeds. Contrary to the generally adopted assumption, that in the bar rotating frame the corotation region should possess four stationary equilibrium points (Lagrange points), with high quality N-body simulations we localize the instantaneous equilibrium points and find that they circulate or oscillate broadly in azimuth with respect to the pattern speeds of the inner or outer perturbations. This implies that at the particle level the Jacobi integral is not well conserved around the corotation radius. That is, angular momentum exchanges decouple from energy exchanges, enhancing the chaotic diffusion of stars through the corotation region.

In search for a synthetic understanding, a scenario for the evolution of the star formation rate and the chemical abundances in galaxies is proposed, combining gas infall from galactic halos, outflow of gas by supernova explosions, and an oscillatory star formation process. The oscillatory star formation model is a consequence of the modelling of the fractional masses changes of the hot, warm and cold components of the interstellar medium. The observed periods of oscillation vary in the range $(0.1-3.0)\times10^{7}$\,yr depending on various parameters existing from giant to dwarf galaxies. The evolution of metallicity varies in giant and dwarf galaxies and depends on the outflow process. Observed abundances in dwarf galaxies can be reproduced under fast outflow together with slow evaporation of cold gases into hot gas whereas slow outflow and fast evaporation is preferred for giant galaxies. The variation of metallicities in dwarf galaxies supports the fact that low rate of SNII production in dwarf galaxies is responsible for variation in metallicity in dwarf galaxies of similar masses as suggested by various authors.

Molecular clouds typically consist of 3/4 H2, 1/4 He and traces of heavier elements. In an earlier work we showed that at very low temperatures and high densities, H2 can be in a phase transition leading to the formation of ice clumps as large as comets or even planets. However, He has very different chemical properties and no phase transition is expected before H2 in dense interstellar medium (ISM) conditions. The gravitational stability of fluid mixtures has been studied before, but these studies did not include a phase transition. First, we study the gravitational stability of van der Waals fluid mixtures using linearized analysis and examine virial equilibrium conditions using the Lennard-Jones intermolecular potential. Then, combining the Lennard-Jones and gravitational potentials, the non-linear dynamics of fluid mixtures are studied via computer simulations using the molecular dynamics code LAMMPS. Along with the classical, ideal-gas Jeans instability criterion, a fluid mixture is always gravitationally unstable if it is in a phase transition because compression does not increase pressure. However, the condensed phase fraction increases. In unstable situations the species can separate: in some conditions He precipitates faster than H2, while in other conditions the converse occurs. Also, for an initial gas phase collapse the geometry is essential. Contrary to spherical or filamentary collapses, sheet-like collapses starting below 15K easily reach H2 condensation conditions because then they are fastest and both the increase of heating and opacity are limited. Depending on density, temperature and mass, either rocky H2 planetoids, or gaseous He planetoids form. H2 planetoids are favoured by high density, low temperature and low mass, while He planetoids need more mass and can form at temperature well above the critical value.

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 existence of substellar cold H2 globules in planetary nebulae and the mere existence of comets suggest that the physics of cold interstellar gas might be much richer than usually envisioned. We study the case of a cold gaseous medium in ISM conditions which is subject to a gas-liquid/solid phase transition. First the equilibrium of general non-ideal fluids is studied using the virial theorem and linear stability analysis. Then the non-linear dynamics is studied by using simulations to characterize the expected formation of solid bodies analogous to comets. The simulations are run with a state of the art molecular dynamics code (LAMMPS). The long-range gravitational forces can be taken into account with short-range molecular forces with finite limited computational resources by using super-molecules, provided the right scaling is followed. The concept of super-molecule is tested with simulations, allowing us to correctly satisfy the Jeans instability criterion for one-phase fluids. The simulations show that fluids presenting a phase transition are gravitationally unstable as well, independent of the strength of the gravitational potential, producing two distinct kinds of sub-stellar bodies, those dominated by gravity ("planetoids") and those dominated by molecular attractive force ("comets"). Observations, formal analysis and computer simulations suggest the possibility of the formation of substellar H2 clumps in cold molecular clouds due to the combination of phase transition and gravity. Fluids presenting a phase transition are gravitationally unstable, independent of the strength of the gravitational potential. Arbitrarily small H2 clumps may form even at relatively high temperatures up to 400 - 600K, according to virial analysis. The combination of phase transition and gravity may be relevant for a wider range of astrophysical situations, such as proto-planetary disks.

The talk consisted mainly in commenting in a linear way the seminal paper in 1964 by Michel Henon and graduate student Carl Heiles at Princeton University: "The applicability of the third integral of motion: Some numerical experiments" in the field of galactic dynamics. Instead of repeating here the lecture of the paper, we advise the reader interested in dynamical systems to study this "must" reading.

Results from the UV satellite GALEX revealed large extensions of disks in some nearby spiral galaxies, extending out to 3 to 4 times the isophotal radius, r25. M63 is a remarkable example of a spiral galaxy with one of the most extended UV disks, so it offers the opportunity to search for the molecular gas and characterize the star formation in outer disk regions as revealed by the UV emission. We obtained deep CO(1-0) and CO(2-1) observations on the IRAM 30 m telescope along the major axis of the M63 disk from the center out to the galactocentric radius rgal = 1.6 r25 and over a bright UV region at rgal = 1.36 r25. CO(1-0) is detected all along the M63 major axis out to r25, and CO(2-1) is confined to rgal = 0.68 r25, which may betray lower excitation temperatures in the outer disk. CO(1-0) is also detected in the external bright UV region of M63. The radial profiles of the CO emission and of the Halpha, 24 micron, NUV and FUV star formation tracers and HI taken from the literature show a severe drop with the galactocentric radius, such that beyond r25 they are all absent with the exception of a faint UV emission and HI. The CO emission detection in the external UV region, where the UV flux is higher than the UV flux observed beyond r25, highlights a tight correlation between the CO and UV fluxes, namely the amount of molecular gas and the intensity of star formation. This external UV region is dominated by the atomic gas, suggesting that HI is more likely the precursor of H2 rather than the product of UV photodissociation. A broken power law needs to be invoked to describe the Kennicutt-Schmidt (K-S) relation of M63 from the center of the galaxy out to rgal = 1.36 r25. While all along the major axis out to r25 the K-S relation is almost linear, in the external UV region the SFR regime is highly nonlinear and characterized by a steep K-S relation and very low star formation efficiency.

We use the Wilkinson Microwave Anisotropy Probe 7-year data (WMAP7) to further probe point source detection technique in the sky maps of the cosmic microwave background (CMB) radiation. The method by Tegmark et al. for foreground reduced maps and the Kolmogorov parameter as the descriptor are adopted for the analysis of WMAP satellite CMB temperature data. Part of the detected points coincide with point sources already revealed by other methods. However, we have also found 2 source candidates for which still no counterparts are known, and identified 7 point sources listed in Planck Early Release Compact Source Catalogue as high reliability sources.

The structure and evolution of the stellar velocity ellipsoid plays an important role in shaping galaxies undergoing bar driven secular evolution and the eventual formation of a boxy/peanut bulge such as present in the Milky Way. Using collisionless N-body simulations, we show that during the formation of such a boxy/peanut bulge, the meridional shear stress of stars, which can be measured by the meridional tilt of the velocity ellipsoid, reaches a characteristic peak in its time evolution. It is shown that the onset of a bar buckling instability is closely connected to the maximum meridional tilt of the stellar velocity ellipsoid. Our findings bring new insight to this complex gravitational instability of the bar which complements the buckling instability studies based on orbital models. We briefly discuss the observed diagnostics of the stellar velocity ellipsoid during such a phenomenon.

We characterize empirically the radial diffusion of stars in the plane of a typical barred disk galaxy by calculating the local spatial diffusion coefficient and diffusion time-scale for bulge-disk-halo N-body self-consistent systems which initially differ in the Safronov-Toomre-Q_T parameter. We find different diffusion scenarios that depend on the bar strength and on the degree of instability of the disk. Marginally stable disks, with Q_T \sim 1, have two families of bar orbits with different values of angular momentum and energy, which determine a large diffusion in the corotation region. In hot disks, Q_T> 1, stellar diffusion is reduced with respect to the case of marginally stable disks. In cold models, we find that spatial diffusion is not constant in time and strongly depends on the activity of the bar, which can move stars all over the disk recurrently. We conclude that to realistically study the impact of radial migration on the chemical evolution modeling of the Milky Way the role of the bar has to be taken into account.

We characterize the radial migration of stars in the disk plane by calculating the diffusion coefficient and the diffusion time-scale for a bulge-disk N-body self-consistent system with a marginally-stable Toomre-Q parameter. We find that diffusion is not constant in time, but follows the evolution of the bar, and becomes maximum near the corotation region and in the external disk region, where asymmetric patterns develop.

We analyse the radial and tangential velocity fields in the Galaxy as seen from the Sun by using as a first approximation a simple axisymmetric model, which we then compare with the corresponding fields in a barred N-body model of the Milky Way. This provides a global description of these quantities missing in the literature, showing where they take large values susceptible to be used in future observations even for sources well beyond the Galactic center. Absolute largest proper motions occur at a distance slightly behind the Galactic Center, which are there 1.5 times larger than the highest local proper motions due to the Galactic differential rotation. Large proper motions well beyond the Galactic center are well within the current astrometric accuracy.

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

After general considerations about limits of theories and models, where small changes may imply large effects, we discuss three cases in galactic astrophysics illustrating how galactic dynamics models may become insufficient when previously neglected effects are taken into account: 1) Like in 3D hydrodynamics, the non-linearity of the Poisson-Boltzmann system may imply dissipation through the growth of discontinuous solutions. 2) The relationship between the microscopic exponential sensitivity of N-body systems and the stability of mean field galaxy models. 3) The role of quantum physics in the dynamics of structure formation, considering that cosmological neutrinos are massive and semi-degenerate fermions.

$\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.

We report a new search for 12CO(1-0) emission in high-velocity clouds (HVCs) performed with the IRAM 30 m telescope. This search was motivated by the recent detection of cold dust emission in the HVCs of Complex C. Despite a spatial resolution which is three times better and sensitivity twice as good compared to previous studies, no CO emission is detected in the HVCs of Complex C down to a best 5 sigma limit of 0.16 K km/s at a 22'' resolution. The CO emission non-detection does not provide any evidence in favor of large amounts of molecular gas in these HVCs and hence in favor of the infrared findings. We discuss different configurations which, however, allow us to reconcile the negative CO result with the presence of molecular gas and cold dust emission. H2 column densities higher than our detection limit, N(H2) = 3x10^{19} cm^{-2}, are expected to be confined in very small and dense clumps with 20 times smaller sizes than the 0.5 pc clumps resolved in our observations according to the results obtained in cirrus clouds, and might thus still be highly diluted. As a consequence, the inter-clump gas at the 1 pc scale has a volume density lower than 20 cm^{-3} and already appears as too diffuse to excite the CO molecules. The observed physical conditions in the HVCs of Complex C also play an important role against CO emission detection. It has been shown that the CO-to-H2 conversion factor in low metallicity media is 60 times higher than at the solar metallicity, leading for a given H2 column density to a 60 times weaker integrated CO intensity. And the very low dust temperature estimated in these HVCs implies the possible presence of gas cold enough (< 20 K) to cause CO condensation onto dust grains under interstellar medium pressure conditions and thus CO depletion in gas-phase observations.

Context: The widespread view that cosmological neutrinos, even if massive, are well described since the decoupling redshift z~10^10 down to the present epoch by an almost perfectly collisionless fluid of classical point particles is re-examined. Aims: In view of the likely sub-eV rest mass of neutrinos, the main effects due to their fermionic nature are studied. Methods: By numerical means we calculate the accurate entropy, fugacity and pressure of cosmological neutrinos in the Universe expansion. By solving the Schroedinger equation we derive how and how fast semi-degenerate identical free fermions become entangled. Results: We find that for sub-eV neutrinos the exchange degeneracy has significantly increased during the relativistic to non-relativistic transition epoch at z~10^4-10^5. At all times neutrinos become entangled in less than 10^-6 s, much faster than any plausible decoherence time. The total pressure is increased by quantum effect from 5% at high redshifts to 68% at low redshifts with respect to a collisionless classical fluid. Conclusions: The quantum overpressure has no dynamical consequences in the homogeneous regime at high redshifts, but must be significant for neutrino clustering during the non-linear structure formation epoch at low redshifts.

In a first part the scope of classical thermodynamics and statistical mechanics is discussed in the broader context of formal dynamical systems, including computer programmes. In this context classical thermodynamics appears as a particular theory suited to a subset of all dynamical systems. A statistical mechanics similar to the one derived with the microcanonical ensemble emerges from dynamical systems provided it contains, 1) a finite non-integrable part of its phase space which is, 2) ergodic at a satisfactory degree after a finite time. The integrable part of phase space provides the constraints that shape the particular system macroscopical properties, and the chaotic part provides well behaved statistical properties over a relevant finite time. More generic semi-ergodic systems lead to intermittent behaviour, thus may be unsuited for a statistical description of steady states. Following these lines of thought, in a second part non-extensive hierarchical systems with statistical scale-invariance and power law interactions are explored. Only the virial constraint, consistent with their microdynamics, is included. No assumptions of classical thermodynamics are used, in particular extensivity and local homogeneity. In the limit of a large hierarchical range new constraints emerge in some conditions that depend on the interaction law range. In particular for the gravitational case, a velocity-site scaling relation is derived which is consistant with the ones empirically observed in the fractal interstellar medium.

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.

To understand galaxies and their evolution, it is necessary to describe how the different scales interact: how the microscopic physics, such as star formation, or the large scale physics, such as galaxy interactions may modify the galaxy global shapes. The purpose of this review is to point out some general or recent topics related to such scale interactions, both observational and theoretical, which are relevant in the present understanding of galaxies.

The thermochemistry of H2 and HD in non-collapsed, non-reionized primordial gas up to the end of the dark age is investigated with recent radiation-matter and chemical reaction rates taking into account the efficient coolant HD, and the possibility of a gas-solid phase transition of H2. In the standard big-bang model we find that these molecules can freeze out and lead to the growth of flakes of solid molecular hydrogen at redshifts z ~ 6-12 in the unperturbed medium and under-dense regions. While this freezing caused by the mere adiabatic cooling of the expanding matter is less likely to occur in collapsed regions due to their higher than radiation background temperature, on the other hand the super-adiabatic expansion in voids strongly favors it. Later reionization (at z ~ 5-6) eventually destroys all these H2 flakes. The possible occurrence of H2 flakes is important for the degree of coupling between matter and radiation, as well as for the existence of a gas-grain chemistry at the end of the dark age.

Observed self-gravitating systems reveal often fragmented non-equilibrium structures that feature characteristic long-range correlations. However, models accounting for non-linear structure growth are not always consistent with observations and a better understanding of self-gravitating $N$-body systems appears necessary. Because unstable gravitating systems are sensitive to non-gravitational perturbations we study the effect of different dissipative factors as well as different small and large scale boundary conditions on idealized $N$-body systems. We find, in the interval of negative specific heat, equilibrium properties differing from theoretical predictions made for gravo-thermal systems, substantiating the importance of microscopic physics and the lack of consistent theoretical tools to describe self-gravitating gas. Also, in the interval of negative specific heat, yet outside of equilibrium, unforced systems fragment and establish transient long-range correlations. The strength of these correlations depends on the degree of granularity, suggesting to make the resolution of mass and force coherent. Finally, persistent correlations appear in model systems subject to an energy flow.

Following Toomre & Kalnajs (1991), local models of slightly dissipative self-gravitating disks show how inhomogeneous structures can be maintained over several galaxy rotations. Their basic physical ingredients are self-gravity, dissipation and differential rotation. In order to explore the structures resulting from these processes on the kpc scale, local simulation of self-gravitating disks are performed in this paper in 2D as well as in 3D. The third dimension becomes a priori important as soon as matter clumping causes a tight coupling of the 3D equations of motion. The physically simple and general framework of the model permits to make conclusions beyond the here considered scales. A time dependent affine coordinate system is used, allowing to calculate the gravitational forces via a particle-mesh FFT-method, increasing the performance with respect to previous direct force calculations. Persistent patterns, formed by transient structures, whose intensity and morphological characteristic depend on the dissipation rate are obtained and described. Some of our simulations reveal first signs of mass-size and velocity dispersion-size power-law relations, but a clear scale invariant behavior will require more powerful computer techniques.

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.

Gravity is scale free. Thus gravity may form similar structures in self-gravitating systems on different scales. Indeed, observations of the interstellar medium, spiral disks and cosmic structures, reveal similar characteristics. The structures in these systems are very lumpy and inhomogeneous. Moreover some of these structures do not seem to be of random nature, but obey certain power laws. Models of slightly dissipative self-gravitating disks show how such inhomogeneous structures can be maintained on the kpc-scale. The basic physical processes in these models are self-gravity, dissipation and differential rotation. In order to explore the structures resulting from these processes, local simulations of self-gravitating disks are performed in 2D and 3D. We observe persistent patterns, formed by transient structures, whose intensity and morphological characteristic depend on the dissipation rate.

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).

A series of arguments is presented for heavy galaxy disks not only in the optical regions, but also in the dark matter dominated regions of spirals. We are testing this possibility with extreme maximum disk N-body models without any conventional spheroidal dark halo.

The interstellar medium (ISM) reveals strongly inhomogeneous structures at every scale. These structures do not seem completely random since they obey certain power laws. Larson's law (\citeyear{Larson81}) $\sigma \propto R^{\delta}$ and the plausible assumption of virial equilibrium justify to consider fractals as a possible description of the ISM. In the following we investigate how self-gravitation, differential rotation and dissipation affect the matter distribution in galaxies. To this end we have performed 3D-simulations for self-gravitating local boxes embedded in a larger disk, extending the 2D-method of Toomre & Kalnajs (\citeyear{Toomre91}) and Wisdom & Tremaine (\citeyear{Wisdom88}). Our simulations lead to the conclusion that gravitation, shearing and dissipation can be dominantly responsible for maintaining an inhomogeneous and eventually a fractal distribution of the matter.

In this contribution, we investigate the scenario of major mergers between galaxies of comparable mass to the formation of massive counter-rotating disks. Such a scenario can be ruled out in general because they produce rather ellipticals. Yet for a narrow range of initial conditions, major mergers can nevertheless produce anyway rare galaxies such as NGC 4550, and so they are viable alternatives for the formation of these counter-rotating systems. The scenario is successful in producing a remarkably axisymmetric disk possessing strongly counter-rotating populations.

We remind that the assumptions almost universally adopted among astronomers concerning the physics to use to describe rarefied cosmic gases remain often without justifications, mainly because the long range of gravitation invalidates the use of classical thermal physics. In turn, without sufficiently good local thermal equilibrium, macroscopic quantities, such as temperature and pressure, are not defined and the fundamental assumption that locally the medium is smoothed by ``molecular chaos'' to justify the use of differential equations is not granted. The highly inhomogeneous fractal state of the interstellar gas is probably a plain symptom of the large discrepancy between the available theoretical tools, predicting local homogeneity after a few sound crossing times, and reality. Such fundamental problems begin to occur in optically thin media such as stellar atmospheres, but become exacerbated in the interstellar medium, in cooling flows, and in the post-recombination gas, particularly when gravitation becomes energetically dominant, i.e., when the medium is Jeans unstable.

Over the last two decades, realistic studies have often concluded that the first bound objects to form can be very small, much smaller than a solar mass. After recombination, the Jeans mass drops rapidly to the order of a Giant Molecular Cloud (GMC) mass ($\sim 10^5$ M$_\odot$), and the H$_2$ cooling can make the collapse quasi-isothermal; this leads to recursive fragmentation, and formation of clumps so dense that 3-body reactions transform the gas almost entirely to the molecular phase. This could lead to star formation in some places, but since star formation is very inefficient, most of the molecular gas could consist of a fractal built on clumpuscules thermalized with the background radiation, and filling a tiny fraction of the volume. The bulk of the gas mass can therefore be trapped in this phase, well before the first stars re-heat and re-ionize the diffuse gaseous medium}. This results in a very contrasted multi-phase baryonic medium, that has partly remained until the present time.

A single merger scenario for making galaxies such as NGC 4550 possessing equal coplanar counter-rotating stellar disks is investigated by collisionless N-body technique. The scenario is successful in producing an axisymmetric disk made of two almost equal counter-rotating populations. The final disk shows a clear bimodal line profile in the outer part, which demonstrates that disk-disk mergers do not always produce ellipticals.

We present ISOCAM images (1.5"/px) in the LW2 (6.75 microns) and LW3 (15 microns) broad band filters of the circumnuclear starburst of NGC 4321 (M 100). A detailed comparison with the emission at other wavelengths (U, B, Halpha, and CO) is also presented. The nature and intensity of the mid-infrared emission is shown to be remarkably different from one star forming region to the next. For instance, star forming regions at the ends of the nuclear bar may contain enshrouded Wolf-Rayet stars.

For a sample of ten nearby barred galaxies, images at 6.75\um and 15\um of circumnuclear starbursts are being taken by ISOCAM on board of the ISO satellite. We first recall the main goals of our project, give its current status, and finally, as an example, preliminary results for NGC 4321 (M100) are presented.

We review here the main direct or indirect ways to detect the possible presence of large amounts of cold molecular hydrogen in the outer parts of disk galaxies, an hypothesis that we have recently developed. Direct ways range from H2 absorption in the UV domain to detection of the radio hyperfine structure: the ortho-H2 molecule has an hyperfine, or ultrafine, structure in its fundamental state, due to the coupling between the rotation-induced magnetic moment, and the nuclear spin. This gives rise to 2 magnetic dipole transitions, at the wavelengths of 0.5 and 5.5 km. Indirect ways are essentially the detection of the HD and LiH transitions, and in some environments like clusters of galaxies, more heavy trace molecules such as CO. We discuss from this point of view the recent discovery by COBE/FIRAS of a very cold Galactic dust component (4-7 K) which could correspond to a dominating gas mass component of the ISM, if interpreted as standard dust emission. Some of the proposed means could be applied to the well-known molecular clouds, to bring some new light to the problem of the H2/CO conversion ratio.

Today we have numerous evidences that spirals evolve dynamically through various secular or episodic processes, such as bar formation and destruction, bulge growth and mergers, sometimes over much shorter periods than the standard galaxy age of 10-15 Gyr. This, coupled to the known properties of the Hubble sequence, leads to a unique sense of evolution: from Sm to Sa. Linking this to the known mass components provides new indications on the nature of dark matter in galaxies. The existence of large amounts of yet undetected dark gas appears as the most natural option. Bounds on the amount of dark stars can be given since their formation is mostly irreversible and requires obviously a same amount of gas.

The ubiquitous clumpy state of the ISM raises a fundamental and open problem of physics, which is the correct statistical treatment of systems dominated by long range interactions. A simple solvable hierarchical model is presented which explains why systems dominated by gravity prefer to adopt a fractal dimension around 2 or less, like the cold ISM and large scale structures. This has direct relation with the general transparency, or blackness, of the Universe.

We have recently proposed a new candidate for baryonic dark matter: very cold molecular gas, in near-isothermal equilibrium with the cosmic background radiation at 2.73 K. The cold gas, of quasi-primordial abundances, is condensed in a fractal structure, resembling the hierarchical structure of the detected interstellar medium. We present some perspectives of detecting this very cold gas, either directly or indirectly. The H$_2$ molecule has an "ultrafine" structure, due to the interaction between the rotation-induced magnetic moment and the nuclear spins. But the lines fall in the km domain, and are very weak. The best opportunity might be the UV absorption of H$_2$ in front of quasars. The unexpected cold dust component, revealed by the COBE/FIRAS submillimetric results, could also be due to this very cold H$_2$ gas, through collision-induced radiation, or solid H$_2$ grains or snowflakes. The $\gamma$-ray distribution, much more radially extended than the supernovae at the origin of cosmic rays acceleration, also points towards and extended gas distribution.

Recent observational constraints restrict the strict applicability of stellar dynamics in spirals to a few rotation periods. However, stellar dynamics concepts such as periodic orbits are invaluable for understanding the various dynamical processes occurring during much more periods. A distinction of two instability types in stellar systems is pointed out, the first one being well illustrated by the bar instability, and the second one by the bar bending instability. In bars the third dimension brings essential dynamical effects which modify the views about the history of bulges and the spiral secular evolution. Bars may grow, bend, thicken, and dissolve into spheroidal bulges, and spirals may evolve along the Hubble sequence in the sense Sd$\to$Sa. This leads to a much more dynamical picture of isolated galaxies than imagined before.

The conjecture that dark matter in galaxies is mostly cold fractal gas is developed in the more general context of the thermodynamics of the ideal isothermal gas subject to gravitational instability. This simple gas model already contains the contrary ingredients able to prevent an asymptotic equilibrium: any growing gravothermal singularity evaporates in a finite time, and any tendency to uniform gas is gravitationally unstable. The paradox is simply resolved by allowing fractal states, which are then scale-free and steady {\it in average}, but non-differentiable and time-dependent. If we apply to clusters the lessons learned with galaxies, we are led to the conclusion that gas in clusters at a temperature much below the virial temperature should also adopt a fractal structure and become inhomogeneous. The same instrumental biases acting at galactic scale and preventing the detection of the smallest and coldest sub-resolution clumps in the fractal are then even more relevant for cluster gas measurements. The large baryonic mass observed in the cluster hot gas and the morphology-density relation suggest also gaseous dark matter in spirals. If this dark gas component remains undetected in spirals, the same instrumental biases should hold in clusters, where cooling hot gas disappears from detection near the centre.

Major progresses have been made this last year towards a better knowledge of the invisible mass. Michel Spiro will talk in details about the micro-lensing experiments and their promising results; the ROSAT satellite has provided extended X-ray maps of the hot gas, which traces dark matter in galaxy clusters: they reveal lower amounts of dark matter in clusters than was previously derived; the dark to visible mass in clusters is not larger than its value in spiral galaxies. It was shown, by X-ray data and gravitational lenses analysis that the dark matter density is highly peaked towards the cluster centers. A new dark matter candidate has also been proposed, in the form of cold and fractal molecular gas that could be present around most late-type spiral galaxies and account for the observed flat rotation curves.

Gas cloud models taking into account the recently disclosed fractal structure of cold gas are set up, showing that large errors in the classical gas mass determination based on smooth cloud models can easily follow if the gas is in reality fractal. Fractal clouds must present both optically thin and optically thick clumps in any single wavelength observations. The observed fractal dimension of the cold ISM suggests that mass underestimates by a factor 10 or more are typical. Due to its low temperature (around 3 K), and its condensed fractal structure, together with its low metallicity, the outer gas would be almost invisible for usual detectors. (A&A paper (in press) and figures available by anonymous ftp at obssd8.unige.ch in /pub/fractal as postscript file: dm_paper_II.ps (251k) and FIG*.ps, or papers I & II + figures as a compressed tar file dm_papers.Z.tar (2.1 Mb)).

Based on dynamical constraints about the Hubble sequence evolution, observational data and a number of "conspiracies", we propose that the dark matter around spiral galaxies is in the form of cold gas, essentially in molecular form and rotationally supported. (full A&A paper (in press) available by anonymous ftp at obssd8.unige.ch in /pub/fractal as postscript file: dm_paper_I.ps (170k), or papers I & II + figures as a compressed tar file dm_papers.Z.tar (2.1 Mb)).