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

We present a new set of galaxy scaling relations for the relative mass content of atomic gas, molecular gas and stars. Such relations are driven by disc gravitational instability, and originate from the low galaxy-to-galaxy variance of Toomre's $Q$ stability parameter. We test such relations using more than 100 galaxies, from massive spirals to dwarf irregulars, thus spanning several orders of magnitude in stellar mass ($M_{\star}\approx10^{6\mbox{-}11}\,\mbox{M}_{\odot}$) and atomic gas mass ($M_{\mathrm{HI}}\approx10^{7\mbox{-}10.5}\,\mbox{M}_{\odot}$). Such tests demonstrate (i) that our scaling relations are physically motivated and tightly constrained, (ii) that the mass-averaged gravitational instability properties of galaxy discs are remarkably uniform across the sequence Sa-dIrr, and (iii) that specific angular momentum plays an important role in such a scenario. Besides providing new insights into a very important topic in galaxy evolution, this work provides a simple formula (Eq. 5) that one can use for generating other galaxy relations driven by disc instability. We explain how to do that, mention a few possible applications, and stress the importance of testing our approach further.

The stellar velocity dispersion, $\sigma$, is a quantity of crucial importance for spiral galaxies, where it enters fundamental dynamical processes such as gravitational instability and disc heating. Here we analyse a sample of 34 nearby spirals from the Calar Alto Legacy Integral Field Area (CALIFA) spectroscopic survey, deproject the line-of-sight $\sigma$ to $\sigma_{R}$ and present reliable radial profiles of $\sigma_{R}$ as well as accurate measurements of $\langle\sigma_{R}\rangle$, the radial average of $\sigma_{R}$ over one effective (half-light) radius. We show that there is a trend for $\sigma_{R}$ to increase with decreasing $R$, that $\langle\sigma_{R}\rangle$ correlates with stellar mass ($M_{\star}$) and tested correlations with other galaxy properties. The most significant and strongest correlation is the one with $M_{\star}$: $\langle\sigma_{R}\rangle \propto M_{\star}^{0.5}$. This tight scaling relation is applicable to spiral galaxies of type Sa $\mbox{--}$ Sd and stellar mass $M_{\star}\approx10^{9.5}\mbox{--}10^{11.5}\ \mbox{M}_{\odot}$. Simple models that relate $\sigma_{R}$ to the stellar surface density and disc scale length roughly reproduce that scaling, but overestimate $\langle\sigma_{R}\rangle$ significantly.

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

We introduce a new diagnostic for exploring the link between angular momentum and local gravitational instability in galaxy discs. Our diagnostic incorporates the latest developments in disc instability research, is fully consistent with approximations that are widely used for measuring the stellar specific angular momentum, $j_{\star}=J_{\star}/M_{\star}$, and is also very simple. We show that such a disc instability diagnostic hardly correlates with $j_{\star}$ or $M_{\star}$, and is remarkably constant across spiral galaxies of any given type (Sa$\!-\!$Sd), stellar mass ($M_{\star}=10^{9.5}\!-\!10^{11.5}\,\mathrm{M}_{\odot}$) and velocity dispersion anisotropy ($\sigma_{z\star}/\sigma_{R\star}=0\!-\!1$). The fact that $M_{\star}$ is tightly correlated with star formation rate ($\mathrm{SFR}$), molecular gas mass ($M_{\mathrm{mol}}$), metallicity ($12+\log\mathrm{O/H}$) and other fundamental galaxy properties thus implies that nearby star-forming spirals self-regulate to a quasi-universal disc stability level. This proves the existence of the self-regulation process postulated by several star formation models, but also raises important caveats.

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

We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a subhalo abundance matching technique coupled to merger trees extracted from a N-body simulation. Our approach consist of forcing the model to match the observed SMF at redshift $z \sim 2.3$, and let it evolve down to $z \sim 0.3$ according to a $\tau$ model, an exponentially declining functional form which describes the star formation rate decay of both satellite and central galaxies. In this study, we use three different sets of SMFs: ZFOURGE data from Tomczak et al.; UltraVISTA data from Ilbert et al. and COSMOS data from Davidzon et al. We also build a mock survey combining UltraVISTA with ZFOURGE. Our modelling of quenching timescales is consistent with the evolution of the SMF down to $z \sim 0.3$, with different accuracy depending on the particular survey used for calibration. We tested our model against the observed SMFs at low redshift and it predicts residuals (observation versus model) within $1\sigma$ observed scatter along most of the stellar mass range investigated, and with mean residuals below 0.1 dex in the range $\sim [10^{8.7}-10^{11.7}] M_{\odot}$. We then compare the SFR-$M_*$ relation predicted by the model with the observed one at different redshifts. The predicted SFR-$M_*$ relation underpredicts the median SFR at fixed stellar mass relative to observations at all redshifts. Nevertheless, the shapes are consistent with the observed relations up to intermediate-mass galaxies, followed by a rapid decline for massive galaxies.

We introduce simulations aimed at assessing how well weak gravitational lensing of 21cm radiation from the Epoch of Reionization ($z \sim 8$) can be measured by an SKA-like radio telescope. A simulation pipeline has been implemented to study the performance of lensing reconstruction techniques. We show how well the lensing signal can be reconstructed using the three-dimensional quadratic lensing estimator in Fourier space assuming different survey strategies. The numerical code introduced in this work is capable of dealing with issues that can not be treated analytically such as the discreteness of visibility measurements and the inclusion of a realistic model for the antennae distribution. This paves the way for future numerical studies implementing more realistic reionization models, foreground subtraction schemes, and testing the performance of lensing estimators that take into account the non-Gaussian distribution of HI after reionization. If multiple frequency channels covering $z \sim 7-11.6$ are combined, Phase 1 of SKA-Low should be able to obtain good quality images of the lensing potential with a total resolution of $\sim 1.6$ arcmin. The SKA-Low Phase 2 should be capable of providing images with high-fidelity even using data from $z\sim 7.7 - 8.3$. We perform tests aimed at evaluating the numerical implementation of the mapping reconstruction. We also discuss the possibility of measuring an accurate lensing power spectrum. Combining data from $z \sim 7-11.6$ using the SKA2-Low telescope model, we find constraints comparable to sample variance in the range $L<1000$, even for survey areas as small as $25\mbox{ deg}^2$.

We have proposed a method for measuring weak lensing using the Lyman-alpha forest. Here we estimate the noise expected in weak lensing maps and power spectra for different sets of observational parameters. We find that surveys of the size and quality of the ones being done today and ones planned for the future will be able to measure the lensing power spectrum at a source redshift of z~2.5 with high precision and even be able to image the distribution of foreground matter with high fidelity on degree scales. For example, we predict that Lyman-alpha forest lensing measurement from the Dark Energy Spectroscopic Instrument survey should yield the mass fluctuation amplitude with statistical errors of 1.5%. By dividing the redshift range into multiple bins some tomographic lensing information should be accessible as well. This would allow for cosmological lensing measurements at higher redshift than are accessible with galaxy shear surveys and correspondingly better constraints on the evolution of dark energy at relatively early times.

The angular positions of quasars are deflected by the gravitational lensing effect of foreground matter. The Lyman-alpha forest seen in the spectra of these quasars is therefore also lensed. We propose that the signature of weak gravitational lensing of the forest could be measured using similar techniques that have been applied to the lensed Cosmic Microwave Background, and which have also been proposed for application to spectral data from 21cm radio telescopes. As with 21cm data, the forest has the advantage of spectral information, potentially yielding many lensed "slices" at different redshifts. We perform an illustrative idealized test, generating a high resolution angular grid of quasars (of order arcminute separation), and lensing the Lyman-alphaforest spectra at redshifts z=2-3 using a foreground density field. We find that standard quadratic estimators can be used to reconstruct images of the foreground mass distribution at z~1. There currently exists a wealth of Lya forest data from quasar and galaxy spectral surveys, with smaller sightline separations expected in the future. Lyman-alpha forest lensing is sensitive to the foreground mass distribution at redshifts intermediate between CMB lensing and galaxy shear, and avoids the difficulties of shape measurement associated with the latter. With further refinement and application of mass reconstruction techniques, weak gravitational lensing of the high redshift Lya forest may become a useful new cosmological probe.

The velocity dispersion of cold interstellar gas, sigma, is one of the quantities that most radically affect the onset of gravitational instabilities in galaxy discs, and the quantity that is most drastically approximated in stability analyses. Here we analyse the stability of a large sample of nearby star-forming spirals treating molecular gas, atomic gas and stars as three distinct components, and using radial profiles of sigma_CO and sigma_HI derived from HERACLES and THINGS observations. We show that the radial variations of sigma_CO and sigma_HI have a weak effect on the local stability level of galaxy discs, which remains remarkably flat and well above unity, but is low enough to ensure (marginal) instability against non-axisymmetric perturbations and gas dissipation. More importantly, the radial variation of sigma_CO has a strong impact on the size of the regions over which gravitational instabilities develop, and results in a characteristic instability scale that is one order of magnitude larger than the Toomre length of molecular gas. Disc instabilities are driven, in fact, by the self-gravity of stars at kpc scales. This is true across the entire optical disc of every galaxy in the sample, with few exceptions. In the linear phase of the disc instability process, stars and molecular gas are strongly coupled, and it is such a coupling that ultimately triggers local gravitational collapse/fragmentation in the molecular gas.

Jets from Active Galactic Nuclei (AGN) inflate large cavities in the hot gas environment around galaxies and galaxy clusters. The large-scale gas circulation promoted within such cavities by the jet itself gives rise to backflows that propagate back to the center of the jet-cocoon system, spanning all the physical scales relevant for the AGN. Using an Adaptive Mesh Refinement code, we study these backflows through a series of numerical experiments, aiming at understanding how their global properties depend on jet parameters. We are able to characterize their mass flux down to a scale of a few kiloparsecs to about $0.5\,\mathrm{M_\odot/y}$for as long as $15$ or $20$ Myr, depending on jet power. We find that backflows are both spatially coherent and temporally textbf{intermittent}, independently of jet power in the range $10^{43-45}$ erg/s. Using the mass flux thus measured, we model analytically the effect of backflows on the central accretion region, where a Magnetically Arrested Disk lies at the center of a thin circumnuclear disk. Backflow accretion onto the disk modifies its density profile, producing a flat core and tail. We use this analytic model to predict how accretion beyond the BH magnetopause is modified, and thus how the jet power is temporally modulated. Under the assumption that the magnetic flux stays frozen in the accreting matter, and that the jets are always launched via the Blandford-Znajek (1977) mechanism, we find that backflows are capable of boosting the jet power up to tenfold during relatively short time episodes (a few Myr).

We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a Subhalo Abundance Matching technique coupled to merger trees extracted from a N-body simulation. Our approach, which considers both galaxy mergers and stellar stripping, is to force the model to match the observed SMF at redshift $z>2$, and let it evolve down to the present time according to the observed (SFR-$M_*$) relation. In this study, we use two different sets of SMFs and two SFR-$M_*$ relations: a simple power law and a relation with a mass-dependent slope. Our analysis shows that the evolution of the SMF is more consistent with a SFR-$M_*$ relation with a mass-dependent slope, in agreement with predictions from other models of galaxy evolution and recent observations. In order to fully and realistically describe the evolution of the SMF, both mergers and stellar stripping must be considered, and we find that both have almost equal effects on the evolution of SMF at the massive end. Taking into account the systematic uncertainties in the observed data, the high-mass end of the SMF obtained by considering stellar stripping results in good agreement with recent observational data from the Sloan Digital Sky Survey (SDSS). At $\log M_* < 11.2$, our prediction at z=0.1 is close to \citet{li-white09} data, but the high-mass end ($\log M_* > 11.2$) is in better agreement with \citet{dsouza15} data which account for more massive galaxies.

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

We study the evolution of the Red Sequence (RS) scatter in galaxy clusters and groups simultaneously using predictions from our simulations (cosmological hydrodynamic + semi-analytical) as well as observational data from the HAWK-I Cluster Survey (HCS), a sample of galaxy clusters at redshifts $0.8 < z < 1.5$. We analyze the intrinsic scatter of the RS to investigate whether the stellar age can be its main effective driver, at the same time assessing the role of metallicity variations in shaping the RS building at around epoch $z\sim$1 and beyond. To this purpose we rely on various methods to derive the average age and age spread from the RS colour scatter of the HCS sample, with the aid of population synthesis models. The RS scatter predicted by the models at $z< 0.7$ is found to not depend on the star formation history adopted, whilst at $z>the correlation found between age and rest-frame colour scatters is quite robust, although all age scatter estimations ultimately depend on the definition of RS as well as on the completeness limits adopted. We find that the age spread of RS galaxies predicted by both hydrodynamical simulations and SAM increases with cosmic epoch, while the ratio between the age spread and the average age remains approximately constant. Both trends are in agreement with observational results from both the HCS and other literature samples.

The mass-metallicity relation shows that the galaxies with the lowest mass have the lowest metallicities. As most dwarf galaxies are in group environments, interaction effects such as tides could contribute to this trend. We perform a series of smoothed particle hydrodynamics (SPH) simulations of dwarf galaxies in external tidal fields to examine the effects of tides on their metallicities and metallicity gradients. In our simulated galaxies, gravitational instabilities drive gas inwards and produce centralized star formation and a significant metallicity gradient. Strong tides can contribute to these instabilities, but their primary effect is to strip the outer low-metallicity gas, producing a truncated gas disk with a large metallicity. This suggests that the role of tides on the mass-metallicity relation is to move dwarf galaxies to higher metallicities.

We explore the role that gravitational instability plays in NGC 1068, a nearby Seyfert galaxy that exhibits unusually vigorous starburst activity. For this purpose, we use the Romeo-Falstad disc instability diagnostics and data from BIMA SONG, SDSS and SAURON. Our analysis illustrates that NGC 1068 is a gravitationally unstable "monster". Its starburst disc is subject to unusually powerful instabilities. Several processes, including AGN/stellar feedback, try to quench such instabilities from inside out by depressing the surface density of molecular gas across the central kpc, but they do not succeed. Gravitational instability "wins" because it is driven by the stars via their much higher surface density. In this process, stars and molecular gas are strongly coupled, and it is such a coupling that ultimately triggers local gravitational collapse/fragmentation in the molecular gas.

By means of our own cosmological-hydrodynamical simulation and semi-analytical model we studied galaxy population properties in clusters and groups, spanning over 10 different bands from UV to NIR, and their evolution since redshift z=2. We compare our results in terms of galaxy red/blue fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with recent observational data reaching beyond z=1.5. Different selection criteria were tested in order to retrieve galaxies belonging to the RS: either by their quiescence degree measured from their specific SFR ("Dead Sequence"), or by their position in a colour-colour plane which is also a function of sSFR. In both cases, the colour cut and the limiting magnitude threshold were let evolving with redshift, in order to follow the natural shift of the characteristic luminosity in the LF. We find that the Butcher-Oemler effect is wavelength-dependent, with the fraction of blue galaxies increasing steeper in optical colours than in NIR. Besides, only when applying a lower limit in terms of fixed absolute magnitude, a steep BO effect can be reproduced, while the blue fraction results less evolving when selecting samples by stellar mass or an evolving magnitude limit. We then find that also the RS-LFR behaviour, highly debated in the literature, is strongly dependent on the galaxy selection function: in particular its very mild evolution recovered when measured in terms of stellar mass, is in agreement with values reported for some of the highest redshift confirmed (proto)clusters. As to differences through environments, we find that normal groups and (to a lesser extent) cluster outskirts present the highest values of both star forming fraction and LFR at low z, while fossil groups and cluster cores the lowest: this separation among groups begins after z~0.5, while earlier all group star forming properties are undistinguishable.

We present an unprecedented measurement of the disc stability and local instability scales in the luminous infrared Seyfert 1 host, NGC7469, based on ALMA observations of dense gas tracers and with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. Our derived disc stability parameter accounts for a thick disc structure and its value falls below unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic instability scale and find a striking agreement between our measured scale of ~ 180 pc, and the typical sizes of individual complexes of young and massive star clusters seen in high-resolution images.

The late-type spiral galaxy NGC 6946 is a prime example of molecular gas dynamics driven by "bars within bars". Here we use data from the BIMA SONG and HERACLES surveys to analyse the structure and stability of its molecular disc. Our radial profiles exhibit a clear transition at distance R ~ 1 kpc from the galaxy centre. In particular, the surface density profile breaks at R ~ 0.8 kpc and is well fitted by a double exponential distribution with scale lengths R_1 ~ 200 pc and R_2 ~ 3 kpc, while the 1D velocity dispersion sigma decreases steeply in the central kpc and is approximately constant at larger radii. The fact that we derive and use the full radial profile of sigma rather than a constant value is perhaps the most novel feature of our stability analysis. We show that the profile of the Q stability parameter traced by CO emission is remarkably flat and well above unity, while the characteristic instability wavelength exhibits clear signatures of the nuclear starburst and inner bar within bar. We also show that CO-dark molecular gas, stars and other factors can play a significant role in the stability scenario of NGC 6946. Our results provide strong evidence that gravitational instability, radial inflow and disc heating have driven the formation of the inner structures and the dynamics of molecular gas in the central kpc.

Gravitational instabilities play an important role in galaxy evolution and in shaping the interstellar medium (ISM). The ISM is observed to be highly turbulent, meaning that observables like the gas surface density and velocity dispersion depend on the size of the region over which they are measured. In this work we investigate, using simulations of Milky Way-like disc galaxies with a resolution of $\sim 9$ pc, the nature of turbulence in the ISM and how this affects the gravitational stability of galaxies. By accounting for the measured average turbulent scalings of the density and velocity fields in the stability analysis, we can more robustly characterize the average level of stability of the galaxies as a function of scale, and in a straightforward manner identify scales prone to fragmentation. Furthermore, we find that the stability of a disc with feedback-driven turbulence can be well described by a "Toomre-like" $Q$ stability criterion on all scales, whereas the classical $Q$ can formally lose its meaning on small scales if violent disc instabilities occur in models lacking pressure support from stellar feedback.

We investigate various galaxy population properties of the massive X-ray luminous galaxy cluster XDCP J0044.0-2033 at z=1.58, which constitutes the most extreme matter density peak at this redshift currently known. We analyze deep VLT/HAWK-I NIR data in the J- and Ks-bands, complemented by Subaru imaging in i and V, Spitzer observations at 4.5 micron, and new spectroscopic observations with VLT/FORS2. We detect a cluster-associated excess population of about 90 galaxies, which follows a centrally peaked, compact NFW galaxy surface density profile with a concentration of c200~10. Based on the Spitzer 4.5 micron imaging data, we measure a stellar mass fraction of fstar,500=(3.3+-1.4)% consistent with local values. The total J- and Ks-band galaxy luminosity functions of the core region yield characteristic magnitudes J* and Ks* consistent with expectations from simple z_f=3 burst models. However, a detailed look at the morphologies and color distributions of the spectroscopically confirmed members reveals that the most massive galaxies are undergoing a very active mass assembly epoch through merging processes. Consequently, the bright end of the cluster red-sequence is not in place, while at intermediate magnitudes [Ks*,Ks*+1.6] a red-locus population is present, which is then sharply truncated at magnitudes fainter than Ks*+1.6. The dominant cluster core population comprises post-quenched galaxies transitioning towards the red-sequence at intermediate magnitudes, while additionally a significant blue cloud population of faint star-forming galaxies is present even in the densest central regions. Our observations lend support to the scenario in which the dominant effect of the dense z~1.6 cluster environment is an accelerated mass assembly timescale through merging activity that is responsible for driving core galaxies across the mass quenching threshold of log(Mstar/Msun)~10.4.

Gravitational instabilities play a primary role in shaping the clumpy structure and powering the star formation activity of gas-rich high-redshift galaxies. Here we analyse the stability of such systems, focusing on the size and mass ranges of unstable regions in the disc. Our analysis takes into account the mass-size and linewidth-size scaling relations observed in molecular gas, originally discovered by Larson. We show that such relations can have a strong impact on the size and mass of star-forming clumps, as well as on the stability properties of the disc at all observable scales, making the classical Toomre parameter a highly unreliable indicator of gravitational instability. For instance, a disc with Q=1 can be far from marginal instability, while a disc with Q<<1 can be marginally unstable. Our work raises an important caveat: if clumpy discs at high redshift have scale-dependent surface densities and velocity dispersions, as implied by the observed clump scaling relations, then we cannot thoroughly understand their stability and star formation properties unless we perform multi-scale observations. This will soon be possible thanks to dedicated ALMA surveys, which will explore the physical properties of super-giant molecular clouds at the peak of cosmic star formation and beyond.

We investigate the interplay between jets from Active Galactic Nuclei (AGNs) and the surrounding InterStellar Medium (ISM) through full 3D, high resolution, Adaptive Mesh Refinement simulations performed with the FLASH code. We follow the jet- ISM system for several Myr in its transition from an early, compact source to an extended one including a large cocoon. During the jet evolution, we identify three major evolutionary stages and we find that, contrary to the prediction of popular theoretical models, none of the simulations shows a self-similar behavior. We also follow the evolution of the energy budget, and find that the fraction of input power deposited into the ISM (the AGN coupling constant) is of order of a few percent during the first few Myr. This is in broad agreement with galaxy formation models employing AGN feedback. However, we find that in these early stages, this energy is deposited only in a small fraction (< 1%) of the total ISM volume. Finally we demonstrate the relevance of backflows arising within the extended cocoon generated by a relativistic AGN jet within the ISM of its host galaxy, previously proposed as a mechanism for self-regulating the gas accretion onto the central object. These backflows tend later to be destabilized by the 3D dynamics, rather than by hydrodynamic (Kelvin- Helmholtz) instabilities. Yet, in the first few hundred thousand years, backflows may create a central accretion region of significant extent, and convey there as much as a few millions of solar masses.

In this paper, we propose a Q stability parameter that is more realistic than those commonly used, and is easy to evaluate [see Eq. (19)]. Using our Q_N parameter, you can take into account several stellar and/or gaseous components as well as the stabilizing effect of disc thickness, you can predict which component dominates the local stability level, and you can do all that simply and accurately. To illustrate the strength of Q_N, we analyse the stability of a large sample of spirals from The HI Nearby Galaxy Survey (THINGS), treating stars, HI and H_2 as three distinct components. Our analysis shows that H_2 plays a significant role in disc (in)stability even at distances as large as half the optical radius. This is an important aspect of the problem, which was missed by previous (two-component) analyses of THINGS spirals. We also show that HI plays a negligible role up to the edge of the optical disc; and that the stability level of THINGS spirals is, on average, remarkably flat and well above unity.

We investigate the gravitational instability of galactic discs, treating stars and cold interstellar gas as two distinct components, and taking into account the phenomenology of turbulence in the interstellar medium (ISM), i.e. the Larson-type scaling relations observed in the molecular and atomic gas. Besides deriving general properties of such systems, we analyse a large sample of galaxies from The HI Nearby Galaxy Survey (THINGS), and show in detail how interstellar turbulence affects disc instability in star-forming spirals. We find that turbulence has a significant effect on both the inner and the outer regions of the disc. In particular, it drives the inner gas disc to a regime of transition between two instability phases and makes the outer disc more prone to star-dominated instabilities.

The Wang-Silk approximation, 1/Q ~ 1/Q_stars + 1/Q_gas, is frequently used for estimating the effective Q parameter in two-component discs of stars and gas. Here we analyse this approximation in detail, and show how its accuracy depends on the radial velocity dispersions and Toomre parameters of the two components. We then propose a much more accurate but still simple approximation for the effective Q parameter, which further takes into account the stabilizing effect of disc thickness. Our effective Q parameter is a natural generalization of Toomre's Q, and as such can be used in a wide variety of contexts, e.g. for predicting star formation thresholds in galaxies or for measuring the stability level of galactic discs at low and high redshifts.

The age and metallicity gradients for a sample of group and cluster galaxies from N-body+hydrodynamical simulation are analyzed in terms of galaxy stellar mass. Dwarf galaxies show null age gradient with a tail of high and positive values for systems in groups and cluster outskirts. Massive systems have generally zero age gradients which turn to positive for the most massive ones. Metallicity gradients are distributed around zero in dwarf galaxies and become more negative with mass; massive galaxies have steeper negative metallicity gradients, but the trend flatten with mass. In particular, fossil groups are characterized by a tighter distribution of both age and metallicity gradients. We find a good agreement with both local observations and independent simulations. The results are also discussed in terms of the central age and metallicity, as well as the total colour, specific star formation and velocity dispersion.

(Abridged) We apply a very general statistical theorem introduced by Cramer (1936) to study the origin of the deviations of the halo spin PDF from the reference lognormal shape. We find that these deviations originate from correlations between two quantities entering the definition of spin, namely the ratio $J/M^{5/2}$ (which depends only on mass) and the total gravitational binding energy $E$. To reach this conclusion, we have made usage of the results deduced from two high spatial- and mass resolution simulations. Our simulations cover a relatively small volume and produce a sample of more than 16.000 gravitationally bound halos, each traced by at least 300 particles. We verify that our results are stable to different systematics, by comparing our results with those derived by the GIF2 and by a more recent simulation performed by Maccio' et al. We find that the spin probability distribution function shows systematic deviations from a lognormal, at all redshifts z <= 1. These deviations depend on mass and redshift: at small masses they change little with redshift, and also the best lognormal fits are more stable. The J-M relationship is well described by a power law of exponent $\alpha$ very near to the linear theory prediction (alpha=5/3), but systematically lower than this at z<= 0.3. We argue that the fact that deviations from a lognormal PDF are present only for high-spin halos could point to a role of large-scale tidal fields in the evolution of the spin PDF.

We explore the gravitational instability of clumpy and turbulent gas discs, taking into account the Larson-type scaling laws observed in giant molecular clouds (GMCs) and HI, as well as more general scaling relations. This degree of freedom is of special interest in view of the coming high-z ISM surveys, and is thus potentially important for understanding the dynamical effects of turbulence at all epochs of galaxy evolution. Our analysis shows that turbulence has a deep impact on the gravitational instability of the disc. It excites a rich variety of stability regimes, several of which have no classical counterpart. Among other diagnostics, we provide two useful tools for observers and simulators: (1) the stability map of turbulence, which illustrates our stability scenario and relates it to the phenomenology of interstellar turbulence: GMC/HI observations, simulations and models; (2) a Toomre-like stability criterion, $Q\geq\bar{Q}$, which applies to a large class of clumpy/turbulent discs. We make specific predictions about GMC and cold-HI turbulence, and point out the implications of our analysis for high-z galaxy surveys.

We study the feedback from an AGN on stellar formation within its host galaxy, mainly using one high resolution numerical simulation of the jet propagation within the interstellar medium of an early-type galaxy. In particular, we show that in a realistic simulation where the jet propagates into a two-phase ISM, star formation can initially be slightly enhanced and then, on timescales of few million years, rapidly quenched, as a consequence both of the high temperatures attained and of the reduction of cloud mass (mainly due to Kelvin-Helmholtz instabilities). We then introduce a model of (prevalently) {\em negative} AGN feedback, where an exponentially declining star formation is quenched, on a very short time scale, at a time t_AGN, due to AGN feedback. Using the Bruzual & Charlot (2003) population synthesis model and our star formation history, we predict galaxy colours from this model and match them to a sample of nearby early-type galaxies showing signs of recent episodes of star formation (Kaviraj et al. 2007). We find that the quantity t_gal - t_AGN, where t_gal is the galaxy age, is an excellent indicator of the presence of feedback processes, and peaks significantly around t_gal - t_AGN \approx 0.85 Gyr for our sample, consistent with feedback from recent energy injection by AGNs in relatively bright (M_{B} \lsim -19) and massive nearby early-type galaxies. Galaxies that have experienced this recent feedback show an enhancement of 3 magnitudes in NUV(GALEX)-g, with respect to the unperturbed, no-feedback evolution. Hence they can be easily identified in large combined near UV-optical surveys.

Observations of turbulent velocity dispersions in the HI component of galactic disks show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high resolution adaptive mesh simulations of gaseous disks embedded within dark matter haloes to explore the roles of cooling, star-formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disk slowly cools until gravitational and thermal instabilities give rise to a multi-phase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large scale driving of global non-axisymmetric modes as well as cloud-cloud tidal interactions and merging. At low star-formation rates these processes alone can explain the observed HI velocity dispersion profiles and the characteristic value of ~10 km/s observed within a wide range of disk galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star-formation rate per unit area >10^-3 M_sun/yr/kpc^2.

N-body + hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a LambdaCDM cosmology are used in order to follow the building-up of the colour-magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the Red Sequence (RS) as they get aged over times and eventually set upon a ``dead sequence'' (DS) once they have stopped their bulk star formation activity. Fainter galaxies keep having significant star formation out to very recent epochs and lie broader around the RS. Environment plays a role as galaxies in groups and cluster outskirts hold star formation activity longer than the central cluster regions. However galaxies experiencing infall from the outskirts to the central parts keep star formation on until they settle on to the DS of the core galaxies. Merging contributes to mass assembly until z~1, after which major events only involve the brightest cluster galaxies. The emerging scenario is that the evolution of the colour-magnitude properties of galaxies within the hierarchical framework is mainly driven by star formation activity during dark matter halos assembly. Galaxies progressively quenching their star formation settle to a very sharp ``red and dead'' sequence, which turns out to be universal, its slope and scatter being almost independent of the redshift (since at least z~1.5) and environment. Differently from the DS, the operatively defined RS evolves more evidently with z, the epoch when it changes its slope being closely corresponding to that at which the passive galaxies population takes over the star forming one: this goes from z~1 in clusters down to 0.4 in normal groups.

The effects of particle discreteness in N-body simulations of Lambda Cold Dark Matter (LambdaCDM) are still an intensively debated issue. In this paper we explore such effects, taking into account the scatter caused by the randomness of the initial conditions, and focusing on the statistical properties of the cosmological density field. For this purpose, we run large sets of LambdaCDM simulations and analyse them using a large variety of diagnostics, including new and powerful wavelet statistics. Among other facts, we point out (1) that dynamical evolution does not propagate discreteness noise up from the small scales at which it is introduced, and (2) that one should aim to satisfy the condition epsilon ~ 2d, where "epsilon" is the force resolution and "d" is the interparticle distance. We clarify what such a condition means, and how to implement it in modern cosmological codes.

Melott has made pioneering studies of the effects of particle discreteness in N-body simulations, a fundamental point that needs careful thought and analysis since all such simulations suffer from numerical noise arising from the use of finite-mass particles. Melott (arXiv:0804.0589) claims that the conclusions of our paper (arXiv:0804.0294) are essentially equivalent to those of his earlier work. Melott is wrong: he has jumped onto one of our conclusions and interpreted that in his own way. Here we point out the whys and the wherefores.

N--body + hydrodynamical simulations of formation and evolution of galaxy groups in a LambdaCDM cosmology have been performed. The properties of the galaxy populations in 12 groups are here discussed, with focus on the colour-magnitude relation in both normal and fossil groups.

We have performed a series of N-body/hydrodynamical (TreeSPH) simulations of clusters and groups of galaxies, selected from cosmological N-body simulations within a $\Lambda$CDM framework: these objects have been re-simulated at higher resolution to $z$=0, in order to follow also the dynamical, thermal and chemical input on to the ICM from stellar populations within galaxies. The simulations include metal dependent radiative cooling, star formation according to different IMFs, energy feedback as strong starburst-driven galactic super-winds, chemical evolution with non-instantaneous recycling of gas and heavy elements, effects of a meta-galactic UV field and thermal conduction in the ICM. In this Paper I of a series of three, we derive results, mainly at $z=0$, on the temperature and entropy profiles of the ICM, its X-ray luminosity, the cluster cold components (cold fraction as well as mass--to--light ratio) and the metal distribution between ICM and stars. In general, models with efficient super-winds, along with a top-heavy stellar IMF, are able to reproduce fairly well the observed $L_X-T$ relation, the entropy profiles and the cold fraction. Observed radial ICM temperature profiles can be matched, except for the gradual decline in temperature inside of $r\sim$~0.1$R_{\rm{vir}}$. Metal enrichment of the ICM gives rise to somewhat steep inner iron gradients; yet, the global level of enrichment compares well to observational estimates after correcting for the stars formed at late times at the base of the cooling flows; also the metal partition between stars and ICM gets into good agreement with observations.

Cosmological (LCDM) TreeSPH simulations of galaxy clusters have been performed, invoking star-formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star-burst and/or AGN driven super-winds, effects of a meta-galactic UV field and thermal conduction. Results for two clusters, one Virgo-like (T~3 keV) and one (sub) Coma-like (T~6 keV), are presented. At z=0, the stellar contents of both clusters consist of a central dominant (cD) galaxy surrounded by galaxies and intracluster (IC) stars. The IC stars are found to contribute 20-40% of the total cluster B-band light and to form at a mean redshift z_f~3, on average about 0.5 Gyr prior to the stars in galaxies. UBVRIJHK surface brightness profiles of the IC star populations are presented; the profile of the larger cluster matches the observed V-band profile of the cD in Abell 1413 (T~8 keV). The IC stellar populations have B-R=1.4-1.5, comparable to larger E galaxies. The mean [Fe/H]~0.0 in the central part of the cluster (r~100 kpc) decreasing to about -0.3 at r_vir. The IC stars are alpha-element enhanced, with a weak trend of [O/Fe] increasing with r, and an overall [O/Fe]~0.4 dex. The IC stars are kinematically colder than the galaxies: The velocity dispersions of the IC stars are at r~100-500 kpc only about half of those of the galaxies increasing to about 70% at r=1-2 Mpc. The projected velocity dispersion in the Virgo-like cluster at R>~50 kpc is 300-600 km/s depending on orientation and R. Rotation is dynamically insignificant for the IC stars. The velocity distributions of IC stars and galaxies are in one cluster highly radially anisotropic, in the other close to isotropic.

We use a set of twelve high-resolution N-body/hydrodynamical simulations in the $\Lambda$CDM cosmology to investigate the origin and formation rate of fossil groups (FGs), which are X-ray bright galaxy groups dominated by a large elliptical galaxy, with the second brightest galaxy being at least two magnitudes fainter. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst driven galactic super winds, effects of a meta-galactic UV field and full stellar population synthesis. We find an interesting correlation between the magnitude gap between the first and second brightest galaxy and the formation time of the group. It is found that FGs have assembled half of their final dark matter mass already at $z\ga1$, and subsequently typically grow by minor merging only, wheras non-FGs on average form later. The early assembly of FGs leaves sufficient time for galaxies of $L \sim L_*$ to merge into the central one by dynamical friction, resulting in the large magnitude gap at $z=0$. A fraction of 33$\pm$16% of the groups simulated are found to be fossil, whereas the observational estimate is $\sim$10-20%. The FGs are found to be X-ray over-luminous relative to non-FGs of the same optical luminosity, in qualitative agreement with observations. Finally, from a dynamical friction analysis is found that only because infall of $L \sim L_*$ galaxies happens along filaments with small impact parameters do FGs exist at all.

Cosmological (LambdaCDM) TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. The properties of the galaxy populations in two clusters, one Virgo-like (T~3 keV) and one (sub) Coma-like (T~6 keV), are discussed. The global star formation rates of the cluster galaxies are found to decrease very significantly with time from redshift z=2 to 0, in agreement with observations. The total K-band luminosity of the cluster galaxies correlates tightly with total cluster mass, and for models without additional AGN feedback, the zero point of the relation matches the observed one fairly well. The match to observed galaxy luminosity functions is reasonable, except for a deficiency of bright galaxies (M_B < -20), which becomes increasingly significant with super-wind strength. Results of a high resolution test indicate that this deficiency is not due to ``over--merging''. The redshift evolution of the luminosity functions from z=1 to 0 is mainly driven by luminosity evolution, but also by merging of bright galaxies with the cD. The colour--magnitude relation of the cluster galaxies matches the observed "red sequence" very well and, on average, galaxy metallicity increases with luminosity. As the brighter galaxies are essentially coeval, the colour--magnitude relation results from metallicity rather than age effects, as observed.

We use the mode summation method together with zeta-function regularization to compute the Casimir energy of a dilute dielectric cylinder. The method is very transparent, and sheds light on the reason the resulting energy vanishes.

Cosmological LambdaCDM TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations include: star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. We report results on the temperature and entropy profiles of the ICM, the X-ray luminosity, cold fraction, M/L and IMLR ratios in gas and stars, metal abundances and gradients. Besides, the properties of the galaxy populations in the two richest clusters are discussed: global star formation rates of the cluster galaxies, the total K-band luminosity, the galaxy luminosity functions at z=0 and their redshift evolution, the colour-magnitude relation (``red sequence'') as resulting from metallicity effects, and the role of the IMF in reproducing colours and abundances of the stellar populations. Finally the contribute (20-40 %) to cluster light from the intra-cluster stars and the cD galaxy has been investigated: surface brightness profiles, mean colours and chemical abundances, kinematics (velocity dispersion and distributions) of the IC stellar populations.

Wavelets are a new and powerful mathematical tool, whose most celebrated applications are data compression and de-noising. In Paper I (Romeo, Horellou & Bergh 2003, astro-ph/0302343), we have shown that wavelets can be used for removing noise efficiently from cosmological, galaxy and plasma N-body simulations. The expected two-orders-of-magnitude higher performance means, in terms of the well-known Moore's law, an advance of more than one decade in the future. In this paper, we describe a wavelet add-on code designed for such an application. Our code can be included in common grid-based N-body codes, is written in Fortran, is portable and available on request from the first author. The code can also be applied for removing noise from standard data, such as signals and images.

Noise is a problem of major concern for N-body simulations of structure formation in the early Universe, of galaxies and plasmas. Here for the first time we use wavelets to remove noise from N-body simulations of disc galaxies, and show that they become equivalent to simulations with two orders of magnitude more particles. We expect a comparable improvement in performance for cosmological and plasma simulations. Our wavelet code will be described in a following paper, and will then be available on request.

We use the results of a high resolution N-body simulation to investigate the role of the environment on the formation and evolution of galaxy-sized halos. Starting from a set of constrained initial conditions, we have produced a final configuration hosting a double cluster in one octant and a large void extending over two octants of the simulation box. We present results for two statistics: the relationship between 1-D velocity dispersion and mass and the probability distribution of the spin parameter $P(\lambda)$. The \svm relationship is well reproduced by the Truncated Isothermal Sphere (TIS) model introduced by Shapiro et al. (1999), although the slope is different from the original prediction. A series of \svm relationships for different values of the anisotropy parameter $\beta$, obtained using the theoretical predictions by Lokas and Mamon (2001) for NFW density profiles are found to be only marginally consistent with the data. Using some properties of the equilibrium TIS models, we construct subsamples of {\em fiducial} equilibrium TIS halos from each of the three subregions, and we study their properties. For these halos, we do find an environmental dependence of their properties, in particular of the spin parameter distribution $P(\lambda)$. We study in more detail the TIS model, and we find new relationships between the truncation radius and other structural parameters. No gravitationally bound halo is found having a radius larger than the critical value for gravithermal instability for TIS halos (\rt $\ge 34.2 r_{0}$, where $r_{0}$ is the core radius of the TIS solution). We do however find a dependence of this relationship on the environment, like for the $P(\lambda)$ statistics. These facts hint at a possible r\^{o}le of tidal fields at determining the statistical properties of halos.

The workshop `Astrophysical Dynamics 1999/2000' followed a homonymous advanced research course, and both activities were organized by me. In this opening paper of the proceedings book, I describe them and document their strong impact on the academic life of the local institutions. The advanced research course was open to graduate students, senior researchers, and motivated under-graduate students with good background in physics and mathematics. The course covered several multi-disciplinary issues of modern research on astrophysical dynamics, and thus also of interest to physicists, mathematicians and engineers. The major topic was gas dynamics, viewed in context with stellar dynamics and plasma physics. The course was complemented by parallel seminars on hot topics given by experts in such fields, and open to a wide scientific audience. In particular, I gave a friendly introduction to wavelets, which are becoming an increasingly powerful tool not only for processing signals and images but also for analysing fractals and turbulence, and which promise to have important applications to dynamical modelling of disc galaxies. The workshop was open to a wide scientific audience. The workshop with published proceedings book was, as a matter of fact, the innovative form of exam that I proposed for the advanced research course. The contributions were refereed and their quality is high on average, exceptionally high in a few cases. The advanced research course and the workshop all together produced great enthusiasm in the students and welcomed the participation of a hundred different people, which means an order of magnitude more than an average graduate course at Chalmers University of Technology and G\"oteborg University.

We approach the theoretical problem of compressing a signal dominated by Gaussian noise. We present expressions for the compression ratio which can be reached, under the light of Shannon's noiseless coding theorem, for a linearly quantized stochastic Gaussian signal (noise). The compression ratio decreases logarithmically with the amplitude of the frequency spectrum $P(f)$ of the noise. Entropy values and compression rates are shown to depend on the shape of this power spectrum, given different normalizations. The cases of white noise (w.n.), $f^{n_p}$ power-law noise ---including $1/f$ noise---, (w.n.$+1/f$) noise, and piecewise (w.n.+$1/f |$ w.n.$+1/f^2$) noise are discussed, while quantitative behaviours and useful approximations are provided.

We present a simple way of coding and compressing the data on board the Planck instruments (HFI and LFI) to address the problem of the on board data reduction. This is a critical issue in the Planck mission. The total information that can be downloaded to Earth is severely limited by the telemetry allocation. This limitation could reduce the amount of diagnostics sent on the stability of the radiometers and, as a consequence, curb the final sensitivity of the CMB anisotropy maps. Our proposal to address this problem consists in taking differences of consecutive circles at a given sky pointing. To a good approximation, these differences are independent of the external signal, and are dominated by thermal (white) instrumental noise. Using simulations and analytical predictions we show that high compression rates, $c_r \simeq 10$, can be obtained with minor or zero loss of CMB sensitivity. Possible effects of digital distortion are also analized. The proposed scheme allows for flexibility to optimize the relation with other critical aspects of the mission. Thus, this study constitutes an important step towards a more realistic modeling of the final sensitivity of the CMB temperature anisotropy maps.

Data on board the future PLANCK Low Frequency Instrument (LFI), to measure the Cosmic Microwave Background (CMB) anisotropies, consist of $N$ differential temperature measurements, expanding a range of values we shall call $R$. Preliminary studies and telemetry allocation indicate the need of compressing these data by a ratio of $c_r \simgt 10$. Here we present a study of entropy for (correlated multi-Gaussian discrete) noise, showing how the optimal compression $c_{r,opt}$, for a linearly discretized data set with $N_{bits}=\log_2{N_{max}}$ bits is given by: $c_r \simeq {N_{bits}/\log_2(\sqrt{2\pi e} ~\sigma_e/\Delta)}$, where $\sigma_e\equiv (det C)^{1/2N}$ is some effective noise rms given by the covariance matrix $C$ and $\Delta \equiv R / N_{max}$ is the digital resolution. This $\Delta$ only needs to be as small as the instrumental white noise RMS: $\Delta \simeq \sigma_T \simeq 2 mK$ (the nominal $\mu K$ pixel sensitivity will only be achieved after averaging). Within the currently proposed $N_{bits}=16$ representation, a linear analogue to digital converter (ADC) will allow the digital storage of a large dynamic range of differential temperature $R= N_{max} \Delta $ accounting for possible instrument drifts and instabilities (which could be reduced by proper on-board calibration). A well calibrated signal will be dominated by thermal (white) noise in the instrument: $\sigma_e \simeq \sigma_T$, which could yield large compression rates $c_{r,opt} \simeq 8$. This is the maximum lossless compression possible. In practice, point sources and $1/f$ noise will produce $\sigma_e > \sigma_T$ and $c_{r,opt} < 8$. This strategy seems safer than non-linear ADC or data reduction schemes (which could also be used at some stage).

Modelling gravity is a fundamental problem that must be tackled in N-body simulations of stellar systems, and satisfactory solutions require a deep understanding of the dynamical effects of softening. In a previous paper (Romeo 1997), we have devised a method for exploring such effects, and we have focused on two applications that reveal the dynamical differences between the most representative types of softened gravity. In the present paper we show that our method can be applied in another, more fruitful, way: for developing new ideas about softening. Indeed, it opens a direct route to the discovery of optimal types of softened gravity for given dynamical requirements, and thus to the accomplishment of a physically consistent modelling of disc galaxies, even in the presence of a cold interstellar gaseous component and in situations that demand anisotropic resolution.

Two questions that naturally arise in N-body simulations of stellar systems are: (1) How can we compare experiments that employ different types of softened gravity? (2) Given a particular type of softened gravity, which choices of the softening length optimize the faithfulness of the experiments to the Newtonian dynamics? We devise a method for exploring the dynamical effects of softening, which provides detailed answers in the case of 2-D simulations of disc galaxies and also solves important aspects of the 3-D problem. In the present paper we focus on two applications that reveal the dynamical differences between the most representative types of softened gravity, including certain anisotropic alternatives. Our method is potentially important not only for testing but also for developing new ideas about softening. Indeed, it opens a direct route to the discovery of optimal types of softened gravity for given dynamical requirements, and thus to the accomplishment of a physically consistent modelling.