[abridged] Aims: We test the effects of re-orienting jets from an active galactic nucleus (AGN) on the intracluster medium in a galaxy cluster environment with short central cooling time. We investigate appearance and properties of the resulting cavities, and the efficiency of jets in providing near-isotropic heating to the cooling cluster core. Methods: We use numerical simulations to explore four models of jets over several active/inactive cycles. We keep the jet power and duration fixed, varying only the jet angle prescription. We track the total energy of the intracluster medium (ICM) in the cluster core over time, and the fraction of the jet energy transferred to the ICM, paying attention to where the energy is deposited. We also compare synthetic X-ray images of the simulated cluster to actual observations. Results: Jets whose re-orientation is minimal ($\lesssim 20^{\circ}$) typically produce conical structures of interconnected cavities, with the opening angle of the cones being $\sim 15-20^{\circ}$, extending to $\sim 300$ kpc from the cluster centre. Such jets transfer about $60\%$ of their energy to the ICM, yet they are not very efficient at heating the cluster core, as the jet energy is deposited further out. Jets that re-orient by $\gtrsim 20^{\circ}$ generally produce multiple pairs of detached cavities. Although smaller, these cavities are inflated within the central 50~kpc and are more isotropically distributed, resulting in more effective heating of the core. Such jets, over few hundreds Myr, can deposit up to $80\%$ of their energy where it is required. Consequently, these models come the closest to an heating/cooling balance and to mitigating runaway cooling of the core, even though all models have identical power/duration profiles. Additionally, the corresponding synthetic X-ray images exhibit structures closely resembling those seen in real cool-core clusters.

Many nearby AGNs display a significant short-term variability. In this work we re-analyze photometric data of four active galactic nuclei observed by Kepler in order to study the flickering activity, having as main goal that of searching for multiple components in the power density spectra. We find that all four objects have similar characteristics, with two break frequencies at approximately log(f/Hz)=-5.2 and -4.7. We consider some physical phenomena whose characteristic time-scales are consistent with those observed, in particular mass accretion fluctuations in the inner geometrically thick disc (hot X-ray corona) and unstable relativistic Rayleigh-Taylor modes. The former is supported by detection of the same break frequencies in the Swift X-ray data of ZW229-15. We also discuss rms-flux relations, and we detect a possible typical linear trend at lower flux levels. Our findings support the hypothesis of a multiplicative character of variability, in agreement with the propagating accretion fluctuation model.

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

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 results from SPH-cosmological simulations, including self-consistent modelling of SN feedback and chemical evolution, of galaxies belonging to two clusters and twelve groups. We reproduce the mass-metallicity (ZM) relation of galaxies classified in two samples according to their star-forming activity, as parametrized by their sSFR, across a redshift range up to z=2. Its slope shows irrelevant evolution in the passive sample, being steeper in groups than in clusters. However, the sub-sample of high-mass passive galaxies only is characterized by a steep increase of the slope with redshift, from which it can be inferred that the bulk of the slope evolution of the ZM relation is driven by the more massive passive objects. (...ABRIDGED...) The ZM relation for the star-forming sample reveals an increasing scatter with redshift, indicating that it is still being built at early epochs. The star-forming galaxies make up a tight sequence in the SFR-M_* plane at high redshift, whose scatter increases with time alongside with the consolidation of the passive sequence. We also confirm the anti-correlation between sSFR and stellar mass, pointing at a key role of the former in determining the galaxy downsizing, as the most significant means of diagnostics of the star formation efficiency. Likewise, an anti-correlation between sSFR and metallicity can be established for the star-forming galaxies, while on the contrary more active galaxies in terms of simple SFR are also metal-richer. We discuss these results in terms of the mechanisms driving the evolution within the high- and low-mass regimes at different epochs: mergers, feedback-driven outflows and the intrinsic variation of the star formation efficiency.

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.

The hypervelocity OB stars in the Milky Way Galaxy were ejected from the central regions some 10-100 million years ago. We argue that these stars, {as well as many more abundant bound OB stars in the innermost few parsecs,} were generated by the interactions of an AGN jet from the central black hole with a dense molecular cloud. Considerations of the associated energy and momentum injection have broader implications for the possible origin of the Fermi bubbles and for the enrichment of the intergalactic medium.

We perform relativistic hydrodynamic simulations of the formation and evolution of AGN cocoons produced by very light powerful jets. We calculate the intensity maps of the Sunyaev-Zel'dovich (SZ) effect at high frequencies for the simulated AGN cocoons using the relativistically correct Wright formalism. Our fully relativistic calculations demonstrate that the contribution from the high temperature gas (kb Te ~ 100 keV) to the SZ signal from AGN cocoons at high frequencies is stronger than that from the shocked ambient intercluster medium owing to the fact that the relativistic spectral functions peak at these temperature values. We present simulations of the SZ effect from AGN cocoons at various frequencies, and demonstrate that SZ observations at 217 GHz and at higher frequencies, such as 857 GHz, will provide us with knowledge about the dynamically-dominant component of AGN cocoons.

We present recent Hubble Space Telescope observations of the inner filament of Centaurus A, using the new Wide Field Camera 3 (WFC3) $F225W, F657N$ and $F814W$ filters. We find a young stellar population near the south-west tip of the filament. Combining the WFC3 dataset with archival Advanced Camera for Surveys (ACS) $F606W$ observations, we are able to constrain the ages of these stars to <=10 Myrs, with best-fit ages of 1-4 Myrs. No further recent star-formation is found along the filament. Based on the location and age of this stellar population, and the fact that there is no radio lobe or jet activity near the star-formation, we propose an updated explanation for the origin of the inner filament. Sutherland et al. (1993) suggested that radio jet-induced shocks can drive the observed optical line emission. We argue that such shocks can naturally arise due to a weak cocoon-driven bow shock (rather than from the radio jet directly), propagating through the diffuse interstellar medium from a location near the inner northern radio lobe. The shock can overrun a molecular cloud, triggering star-formation in the dense molecular cores. Ablation and shock heating of the diffuse gas then gives rise to the observed optical line and X-ray emission. Deeper X-ray observations should show more diffuse emission along the filament.

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.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

X-ray cavities and extended radio sources (`cocoons') surrounding active galactic nuclei (AGN) have been detected by the Chandra X-ray mission and radio interferometers. A joint analysis of X-ray and radio maps suggests that pressure values of non-thermal radio-emitting particles derived from the radio maps are not sufficient to inflate the X-ray cavities. We propose using the Sunyaev-Zel'dovich (SZ) effect, whose intensity strongly depends on the pressure, to find the hitherto undetected, dynamically-dominant component in the radio cocoons. We demonstrate that the spectral function at a frequency of 217 GHz has an absolute maximum at a temperature higher than $10^9$ K, therefore the measurement of the SZ effect at this frequency is a powerful tool for potentially revealing the dynamically-dominant component inside AGN jet-driven radio cocoons. A new method is proposed for excluding the contribution from the low energy, non-relativistic electrons to the SZ effect by means of observations at two frequencies. We show how one may correct for a possible contribution from the kinematic SZ effect. The intensity maps of the SZ effect are calculated for the self-similar Sedov solution, and application of a predicted ring-like structure on the SZ map at a frequency of 217 GHz is proposed to determine the energy released during the active jet stage. The SZ intensity map for an AGN cocoon in a distant elliptical is calculated using a 2-D numerical simulation and including relativistic corrections to the SZ effect. We show the intensity spectrum of the SZ effect is flat at high frequencies if gas temperature is as high as $k_\mathrm{b} T_{\mathrm{e}}=500$ keV.

VisIVO is an integrated suite of tools and services specifically designed for the Virtual Observatory. This suite constitutes a software framework for effective visual discovery in currently available (and next-generation) very large-scale astrophysical datasets. VisIVO consists of VisiVO Desktop - a stand alone application for interactive visualization on standard PCs, VisIVO Server - a grid-enabled platform for high performance visualization and VisIVO Web - a custom designed web portal supporting services based on the VisIVO Server functionality. The main characteristic of VisIVO is support for high-performance, multidimensional visualization of very large-scale astrophysical datasets. Users can obtain meaningful visualizations rapidly while preserving full and intuitive control of the relevant visualization parameters. This paper focuses on newly developed integrated tools in VisIVO Server allowing intuitive visual discovery with 3D views being created from data tables. VisIVO Server can be installed easily on any web server with a database repository. We discuss briefly aspects of our implementation of VisiVO Server on a computational grid and also outline the functionality of the services offered by VisIVO Web. Finally we conclude with a summary of our work and pointers to future developments.

We study the internal circulation within the cocoon carved out by the relativistic jet emanating from an AGN within the ISM of its host galaxy. Firstly, we develop a model for the origin of the internal flow, noticing that a significant increase of large scale velocity circulation within the cocoon arises as significant gradients in the density and entropy are created near the hot spot (a consequence of Crocco's vorticity generation theorem). We find simple and accurate approximate solutions for the large scale flow,showing that a backflow towards the few inner parsec region develops. We solve the appropriate fluid dynamic equations, and we use these solutions to predict the mass inflow rates towards the central regions. We then perform a series of 2D simulations of the propagation of jets using FLASH 2.5, in order to validate the predictions of our model. The backflows which arise supply the central AGN region with very low angular momentum gas, at average rates of the order of $0.1-0.8 \rm{M}_{\odot} \rm{yr.}^{-1}$, the exact value seen to be strongly dependent on the central ISM density (for fixed input jet power). The time scales of these inflows are apparently weakly dependent on the jet/ISM parameters, and are of the order of $3-5\times 10^{7} \rm{yrs}$. We then argue that these backflows could (at least partially) feed the AGN, and provide a self-regulatory mechanism of AGN activity, that is not directly controlled by, but instead controls, the star formation rate within the central circumnuclear disk.

We study the internal circulation within the cocoon carved out by a relativistic jet emanating from an AGN, first developing model and then validating it using a series of numerical simulations. We notice that a significant increase of in this flow arises because gradients in the density and entropy develop near the hot spot, as a consequence of Crocco vorticity theorem. We find simple solutions for the streamlines, and we use them to predict the mass inflow rates towards the central regions. The 2D simulations we perform span a rather wide range of mechanical jet's input power and Black Hole masses, and we show that the predicted nuclear mass inflows are in good agreement with the theoretical model.

Aims. The Large Binocular Cameras (LBC) are two twin wide field cameras (FOV ~ 23'x 25') mounted at the prime foci of the 8.4m Large Binocular Telescope (LBT). We performed a weak lensing analysis of the z=0.288 cluster Abell 611 on g-band data obtained by the blue-optimized Large Binocular Camera in order to estimate the cluster mass. Methods. Due to the complexity of the PSF of LBC, we decided to use two different approaches, KSB and Shapelets, to measure the shape of background galaxies and to derive the shear signal produced by the cluster. Then we estimated the cluster mass with both aperture densitometry and parametric model fits. Results. The combination of the large aperture of the telescope and the wide field of view allowed us to map a region well beyond the expected virial radius of the cluster and to get a high surface density of background galaxies (23 galaxies/arcmin^2). This made possible to estimate an accurate mass for Abell 611. We find that the mass within 1.5 Mpc is: $(8\pm3)\times 10^{14} M_\odot$ from the aperture mass technique and $(5\pm1)\times 10^{14} M_\odot$ using the model fitting by a NFW mass density profile, for both Shapelets and KSB methods. This analysis demonstrates that LBC is a powerful instrument for weak gravitational lensing studies.

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.

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.

(Abridged) We study the interaction of relativistic jets from AGNs with the ISM in their host galaxy, using a series of Adaptive Mesh Refinement simulations of the evolution of the interaction between the cocoon produced by the jet with a dense cloud, placed very near the cocoon's path. We vary only the jet's input power between P_{jet} = 10^{41}-10^{47} {\rm erg/sec}. The density Probability Distribution Function (PDF) within the cocoon can be described in terms of two distinct components, which are also spatially distinct: a low- and a high-density component. The PDF of the post-shocked region is well approximated by a modified lognormal distribution, for all values of $P_{jet}$. During the active phase, when the jet is fed by the AGN, the cloud is subject both to compression and stripping, which tend to increase its density and diminish its total mass. When the jet is switched off (i.e. during the passive phase) the shocked cloud cools further and tends to become more filamentary, under the action of a back-flow which develops within the cocoon. We study the evolution of the star formation rate within the cloud, assuming that is determined by a Schmidt-Kennicutt law, and we analyze the different physical factors which have an impact on the star formation rate. We show that, although the star formation rate can occasionally increase, on time scales of the order of $10^{5}-10^{6}$ yrs, star formation will be inhibited and the cloud fragments. The cooling time of the environment within which the cloud is embedded is however very long: thus, star formation from the fragmented cloud remains strongly inhibited.

We study the propagation of relativistic jets originating from AGNs within the Interstellar/Intergalactic Medium of their host galaxies, and use it to build a model for the suppression of stellar formation within the expanding cocoon.

Weak lensing applied to deep optical images of clusters of galaxies provides a powerful tool to reconstruct the distribution of the gravitating mass associated to these structures. We use the shear signal extracted by an analysis of deep exposures of a region centered around the galaxy cluster Abell 209, at redshift z=0.2, to derive both a map of the projected mass distribution and an estimate of the total mass within a characteristic radius. We use a series of deep archival R-band images from CFHT-12k, covering an area of 0.3 deg^2. We determine the shear of background galaxy images using a new implementation of the modified Kaiser-Squires-Broadhurst pipeline for shear determination, which we has been tested against the ``Shear TEsting Program 1 and 2'' simulations. We use mass aperture statistics to produce maps of the 2 dimensional density distribution, and parametric fits using both Navarro-Frenk-White (NFW) and singular-isothermal-sphere profiles to constrain the total mass. The projected mass distribution shows a pronounced asymmetry, with an elongated structure extending from the SE to the NW. This is in general agreement with the optical distribution previously found by other authors. A similar elongation was previously detected in the X-ray emission map, and in the distribution of galaxy colours. The circular NFW mass profile fit gives a total mass of M_{200} = 7.7^{+4.3}_{-2.7} 10^{14} solar masses inside the virial radius r_{200} = 1.8\pm 0.3 Mpc. The weak lensing profile reinforces the evidence for an elongated structure of Abell 209, as previously suggested by studies of the galaxy distribution and velocities.

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.

We investigate the environmental dependence of galaxy population properties in a complete volume-limited sample of 91566 galaxies in the redshift range 0.05 <= z <= 0.095 and with M_r <= -20.0, selected from the Sloan Digital Sky Survey (SDSS) Data Release 4. Our aim is to search for systematic variations in the properties of galaxies with the local galaxy density. In particular, we analize how the (u - r) color index and the morphological type of galaxies (the latter evaluated through the SDSS Eclass and FracDev parameters) are related to the environment and to the luminosity of galaxies, in order to find hints that can be related to the presence of a ''void'' galaxy population. Void galaxies are identified through a highly selective criterion, which takes also into account redshift and allows us to exclude from the sample all the galaxies that are not really close to the center of underdense regions. We study the (u - r) color distribution for galaxies in different luminosity bins, and we look for correlations of color with the environment, the luminosity, and the morphological type of the galaxies. We find that galaxies in underdense regions (voids) have lower luminosity (M_r > -21) and are bluer than cluster galaxies. The transition from overdense to underdense environments is smooth, the fraction of late-type galaxies decreases while the fraction of early-type galaxies increases smoothly from underdense to dense environments. We do not find any sudden transition in the galaxy properties with density, which, according to a suggestion by Peebles (2001), should mark the transition to a population of "void" galaxies in LCDM models. On the contrary, our results suggest a continuity of galaxy properties, from voids to clusters.

FLY is a parallel treecode which makes heavy use of the one-sided communication paradigm to handle the management of the tree structure. In its public version the code implements the equations for cosmological evolution, and can be run for different cosmological models. This reference guide describes the actual implementation of the algorithms of the public version of FLY, and suggests how to modify them to implement other types of equations (for instance the Newtonian ones)

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.

We present a code allowing to evolve three-dimensional self-gravitating collisionless systems with a large number of particles N >= 10^7. FLY (Fast Level-based N-bodY code) is a fully parallel code based on a tree algorithm. It adopts periodic boundary conditions implemented by means of the Ewald summation technique. FLY is based on the one-side communication paradigm and was originally developed on a CRAY T3E system using the SHMEM library and it was ported on SGI ORIGIN 2000 and on IBM SP. FLY (http://www.ct.astro.it/fly/) version 1.1 is an open source freely available code. FLY data output can be analysed with AstroMD, an analysis and visualization tool specifically designed to deal with the visualization and analysis of astrophysical data. AstroMD can manage different physical quantities. It can find out structures without well defined shape or symmetries, and perform quantitative calculations on selected regions. AstroMD (http://www.cineca.it/astromd) is a freely available code.

The code we describe (FLY) is a newly written code (using the tree N-body method), for three-dimensional self-gravitating collisionless systems evolution. FLY is a fully parallel code based on the tree Barnes-Hut algorithm and periodical boundary conditions are implemented by means of the Ewald summation technique. FLY is based on the one-side communication paradigm to share data among the processors, that access to remote private data avoiding any kind of synchronism. The code was originally developed on CRAY T3E system using the logically SHared MEMory access routines (SHMEM) but it runs also on SGI ORIGIN systems and on IBM SP by using the Low-Level Application Programming Interface routines (LAPI). This new code reaches good performance in all systems where it has been well-tested. This performance allows us today to consider the code FLY among the most powerful parallel codes for tree N-body simulations. The FLY version 1.1 is freely available on http://www.ct.astro.it/fly/ and it will be maintained and upgraded with new releases.

N-body codes to perform simulations of the origin and evolution of the Large Scale Structure of the Universe have improved significantly over the past decade both in terms of the resolution achieved and of reduction of the CPU time. However, state-of-the-art N-body codes hardly allow one to deal with particle numbers larger than a few 10^7, even on the largest parallel systems. In order to allow simulations with larger resolution, we have first re-considered the grouping strategy as described in Barnes (1990) (hereafter B90) and applied it with some modifications to our WDSH-PT (Work and Data SHaring - Parallel Tree) code. In the first part of this paper we will give a short description of the code adopting the Barnes and Hut algorithm \cite{barh86} (hereafter BH), and in particular of the memory and work distribution strategy applied to describe the {\it data distribution} on a CC-NUMA machine like the CRAY-T3E system. In the second part of the paper we describe the modification to the Barnes grouping strategy we have devised to improve the performance of the WDSH-PT code. We will use the property that nearby particles have similar interaction list. This idea has been checked in B90, where an interaction list is builded which applies everywhere within a cell C_{group} containing a little number of particles N_{crit}. B90 reuses this interaction list for each particle $ p \in C_{group}$ in the cell in turn. We will assume each particle p to have the same interaction list. Thus it has been possible to reduce the CPU time increasing the performances. This leads us to run simulations with a large number of particles (N ~ 10^7/10^9) in non-prohibitive times.

Over the past few years, the role of visualization for scientific purpose has grown up enormously. Astronomy makes an extended use of visualization techniques to analyze data, and scientific visualization has became a fundamental part of modern researches in Astronomy. With the evolution of high performance computers, numerical simulations have assumed a great role in the scientific investigation, allowing the user to run simulation with higher and higher resolution. Data produced in these simulations are often multi-dimensional arrays with several physical quantities. These data are very hard to manage and to analyze efficiently. Consequently the data analysis and visualization tools must follow the new requirements of the research. AstroMD is a tool for data analysis and visualization of astrophysical data and can manage different physical quantities and multi-dimensional data sets. The tool uses virtual reality techniques by which the user has the impression of travelling through a computer-based multi-dimensional model. AstroMD is a freely available tool for the whole astronomical community (http://www.cineca.it/astromd/).

We report the first detailed X-ray and optical observations of the medium-distant cluster A33 obtained with the Beppo-SAX satellite and with the UH 2.2m and Keck II telescopes at Mauna Kea. The information deduced from X-ray and optical imaging and spectroscopic data allowed us to identify the X-ray source 1SAXJ0027.2-1930 as the X-ray counterpart of the A33 cluster. The faint, $F_{2-10 keV} \approx 2.4 \times 10^{-13} \ergscm2$, X-ray source 1SAXJ0027.2-1930, $\sim 2$ arcmin away from the optical position of the cluster as given in the Abell catalogue, is identified with the central region of A33. Based on six cluster galaxy redshifts, we determine the redshift of A33, $z=0.2409$; this is lower than the value derived by Leir and Van Den Bergh (1977). The source X-ray luminosity, $L_{2-10 keV} = 7.7 \times 10^{43} \ergs$, and intracluster gas temperature, $T = 2.9$ keV, make this cluster interesting for cosmological studies of the cluster $L_X-T$ relation at intermediate redshifts. Two other X-ray sources in the A33 field are identified. An AGN at z$=$0.2274, and an M-type star, whose emission are blended to form an extended X-ray emission $\sim 4$ arcmin north of the A33 cluster. A third possibly point-like X-ray source detected $\sim 3$ arcmin north-west of A33 lies close to a spiral galaxy at z$=$0.2863 and to an elliptical galaxy at the same redshift as the cluster.

Using the results of a high mass resolution ($256^{3}$ particles) N-body simulation of a cluster-forming region, we study the statistic $\sigma_{v}-M$ for low mass halos (${\rm M} < 10^{12} {\rm M}_{\sun}$), and we compare it with three theoretical models. At the final redshift halos are well described by the Truncated Isothermal Sphere model recently introduced by Shapiro, Iliev & Raga (1999). We speculate that this is a consequence of the action of tidal fields on galaxy formation within clusters.

We present and discuss a method to identify substructures in combined angular-redshift samples of galaxies within Clusters. The method relies on the use of Discrete Wavelet Transform (hereafter DWT) and has already been applied to the analysis of the Coma cluster (Gambera et al. 1997). The main new ingredient of our method with respect to previous studies lies in the fact that we make use of a 3D data set rather than a 2D. We test the method on mock cluster catalogs with spatially localized substructures and on a N-body simulation. Our main conclusion is that our method is able to identify the existing substructures provided that: a) the subclumps are detached in part or all of the phase space, b) one has a statistically significant number of redshifts, increasing as the distance decreases due to redshift distortions; c) one knows {\it a priori} the scale on which substructures are to be expected. We have found that to allow an accurate recovery we must have both a significant number of galaxies ($\approx 200$ for clusters at z$\geq 0.4$ or about 800 at z$\leq$ 0.4) and a limiting magnitude for completeness $m_B=16$. The only true limitation to our method seems to be the necessity of knowing {\it a priori} the scale on which the substructure is to be found. This is an intrinsic drawback of the method and no improvement in numerical codes based on this technique could make up for it.

N-body algorithms for long-range unscreened interactions like gravity belong to a class of highly irregular problems whose optimal solution is a challenging task for present-day massively parallel computers. In this paper we describe a strategy for optimal memory and work distribution which we have applied to our parallel implementation of the Barnes & Hut (1986) recursive tree scheme on a Cray T3D using the CRAFT programming environment. We have performed a series of tests to find an " optimal data distribution " in the T3D memory, and to identify a strategy for the " Dynamic Load Balance " in order to obtain good performances when running large simulations (more than 10 million particles). The results of tests show that the step duration depends on two main factors: the data locality and the T3D network contention. Increasing data locality we are able to minimize the step duration if the closest bodies (direct interaction) tend to be located in the same PE local memory (contiguous block subdivison, high granularity), whereas the tree properties have a fine grain distribution. In a very large simulation, due to network contention, an unbalanced load arises. To remedy this we have devised an automatic work redistribution mechanism which provided a good Dynamic Load Balance at the price of an insignificant overhead.

Evidence for clustering within the Coma cluster is found by means of a multiscale analysis of the combined angular-redshift distribution. We have compiled a catalogue of 798 galaxy redshifts from published surveys from the region of the Coma cluster. We examine the presence of substructure and of voids at different scales ranging from $\sim 1 to \sim 16 h^{-1}$ Mpc, using subsamples of the catalogue, ranging from $cz=3000$ km/s to $cz=28000$ km/s. Our substructure detection method is based on the wavelet transform and on the segmentation analysis. The wavelet transform allows us to find out structures at different scales and the segmentation method allows us a quantitative statistical and morphological analysis of the sample. From the whole catalogue we select a subset of 320 galaxies, with redshifts between cz=5858 km/s and cz=8168 km/s that we identify as belonging to the central region of Coma and on which we have performed a deeper analysis, on scales ranging from $180 h^{-1}$ kpc to $1.44 h^{-1}$ Mpc. Our results are expressed in terms of the number of structures or voids and their sphericity for different values of the threshold detection and at all the scales investigated. According to our analysis, there is strong evidence for multiple hierarchical substructure, on scales ranging from a few hundreds of kpc to about $4 h^{-1}$ Mpc. The morphology of these substructures is rather spherical. On the scale of $720 h^{-1}$ kpc we find two main subclusters which where also found before, but our wavelet analysis shows even more substructures, whose redshift position is approximatively marked by these bright galaxies: NGC 4934 & 4840, 4889, 4898 & 4864, 4874 & 4839, 4927, 4875.

We solve numerically the equations of motion for the collapse of a shell of baryonic matter falling into the central regions of a cluster of galaxies, taking into account of the presence of the substructure inducing dynamical friction. The evolution of the expansion parameter a(t) of the perturbation is calculated in spherical systems. The effect of dynamical friction is to reduce the binding radius and the total mass accreted by the central regions. Using a peak density profile given by Bardeen et al. (1986) we show how the binding radius of the perturbation is modified by dinamical friction. We show how dynamical friction modifies the collapse parameter of the perturbation slowing down the collapse.

We reconsider the problem of the age of the stellar discs of late-type Low Surface Brightness (LSB) galaxies by making use of a new IMF recently derived from numerical fluid dynamical simulations (Padoan, Nordlund and Jones, 1997). While a Miller-Scalo IMF cannot adequately describe the photometric properties of LSBs, when we apply the new Padoan et al. (1997) IMF to a simple exponential disc model with parameters appropriate to LSBs, we get excellent fits of the colors and color gradients. We then conclude that: a) The star formation history of LSB disc galaxies can be described by an initial burst of a few times $10^7$ yr followed by a quiescent period with only sporadic star formation; b) LSBs' discs are not young. The age of the LSB disc galaxies for which colors have been measured are all larger than about 10 Gyr$.

We describe a new parallel N-body code for cosmological simulations. The code is based on a work- and data sharing scheme, and is implemented within the Cray Research Corporation's CRAFT programming environment. Different data distribution schemes have been adopted for bodies' and tree's structures. Tests performed for two different types of initial distributions show that the performance scales almost ideally as a function of the size of the system and of the number of processors. We discuss the factors affecting the absolute speedup and how it can be increased with a better tree's data distribution scheme.

We derive the conditional probability distribution of the peculiar velocity around a peak of given overdensity in a Gaussian density field. This distribution characterizes the shear field around local, isolated density peaks. We use the moments of the velocity distribution to study the spatial distribution of the escape velocity from protostructures forming in a standard, flat, scale--invariant CDM model. In the case of high ($\nu \sim 2$) peaks with filtering scales $R_f \sim 2 h^{-1}~\mbox{Mpc}$, we find that at distances between $2 \div 4~{\rm \mbox{Mpc}~h^{-1}}$ (or equivalently $\approx 1 \div 2 ~R_f$) a sensitive fraction of the matter enclosed within an extrapolated (up to the present epoch) volume containing a total gravitating mass $M \approx 2.725\times10^{15} h^{-1} M_{\odot}$ (of the same order as that derived from gravitational lensing observations, see Tyson and Fisher~\cite{tf95}) is {\em gravitationally unbound} to the peak. The fraction of unbound matter varies between approximately between $0.05-0.2$ depending on the central height $\nu$ of the density peak, on the filtering radius and on the adopted normalization for the fluctuation spectrum. We show how these results support our previous conclusions regarding the delay of the collapse induced by the presence of small-scale substructure in clusters of galaxies formed in standard CDM cosmological models.

We reconsider the problem of the origin of Low-Surface-Brightness (LSB) galaxies within the standard Cold Dark Matter (CDM) cosmological model ($\Omega =1, h=0.5$). Taking into account the effect of substructure on the collapse times of low overdensity peaks, we show that the abundance of these peaks is larger than previously expected because of the dragging caused by dynamical friction induced by the presence of small-scale substructure. The number density of these shallow, low-density peaks turns out to be in good agreement with the observed overdensity of the blue LSB galaxies found during recent surveys.

We study the effect of the dynamical friction induced by the presence of substructure on the statistics of the collapse of density peaks. Applying the results of a former paper we show that within high density environments, like rich clusters of galaxies, the collapse of smaller peaks is strongly delayed until very late epochs. A bias of dynamical nature thus naturally arises because high density peaks preferentially collapse For a standard CDM model we find that this dynamical bias can account for a substantial part of the total bias required by observations on cluster scales.

We have implemented a parallel version of the Barnes-Hut 3-D N-body tree algorithm under PVM 3.2.5, adopting an SPMD paradigm. We parallelize the problem by decomposing the physical domain by means of the {\bf Orthogonal Recursive Bisection} oct-tree scheme suggested by Salmon (1991), but we modify the original hypercube communication pattern into an incomplete hypercube, which is more suitable for a generic inhomogenous cluster architecture.\\ We address dynamical load balancing by assigning different "weights" to the spawned tasks according to the dynamically changing workloads of each task. The weights are determined by monitoring the local platforms where the tasks are running and estimating the performance of each task. The monitoring scheme is flexible and allows us to address at the same time cluster and intrinsic sources of load imbalance. We then show measurements of the performance of our code on a test case of astrophysical interest in order to test the performance of our implementation.