We describe and test a new version of the adaptive mesh refinement (AMR) cosmological code MASCLET. The new version of the code includes all the ingredients of its previous version plus a description of the evolution of the magnetic field under the approximation of the ideal magneto-hydrodynamics (MHD). To preserve the divergence-free condition of MHD, the original divergence cleaning algorithm of Dedner et al. (2002) is implemented. We present a set of well-known 1D and 2D tests, such as several shock-tube problems, the fast rotor and the Orszag-Tang vortex. The performance of the code in all the tests is excellent with estimated median relative errors of $\nabla \cdot {\bf B}$ in the 2D tests smaller than $5 \times 10^{-5}$ for the fast rotor test, and $5 \times 10^{-3}$ for the Orszag-Tang vortex. As an astrophysical application of the code, we present a simulation of a cosmological box of 40 comoving Mpc side length in which a primordial uniform comoving magnetic field of strength 0.1 nG is seeded. The simulation shows how the magnetic field is channelled along the filaments of gas and is concentrated and amplified within galaxy clusters. Comparison with the values expected from pure compression reveals an additional amplification of the magnetic field caused by turbulence in the central region of the cluster. Values of the order of $\sim 1 \mu$G are obtained in clusters at $z\sim 0$ with median relative errors of $\nabla \cdot {\bf B}$ below 0.4\%. The implications of a proper description of the dynamics of the magnetic field and their possible observational counterparts in future facilities are discussed.

We present a long-term numerical three-dimensional simulation of a relativistic outflow designed to be compared with previous results from axisymmetric, two-dimensional simulations, with existing analytical models and state-of-art observations. We follow the jet evolution from 1~kpc to 200~kpc, using a relativistic gas equation of state and a galactic profile for the ambient medium. We also show results from smaller scale simulations aimed to test convergence and different three-dimensional effects. We conclude that jet propagation can be faster than expected from axisymmetric simulations, covering tens of kiloparsecs in a few million years, until the dentist drill effect produced by the growth of helical instabilities slows down the propagation speed of the jet head. A comparison of key physical parameters of the jet structure as obtained from the simulations with values derived from observations of FRII sources reveals good agreement. Our simulations show that shock heating can play a significant role in the feedback from active galaxies, confirming previous 2D results. A proper description of galactic jets as a relativistic scenario, both dynamical and thermodynamical, reveals an extremely fast and efficient feedback process reheating the ICM, and therefore, with dramatic consequences on the galactic evolution. Our results point towards FRII jets as the source of the energetic electrons observed in radio relics.

About half of the expected total baryon budget in the local Universe is `missing'. Hydrodynamical simulations suggest that most of the missing baryons are located in a mildly overdense, warm-hot intergalactic medium (WHIM), which is difficult to be detected at most wavelengths. In this paper we explore multi-wavelength synthetic observations of a massive galaxy cluster developed in a full Eulerian-AMR cosmological simulation. A novel numerical procedure is applied on the outputs of the simulation, which are post-processed with a full-radiative transfer code that allows to compute the change of the intensity at any frequency along the null-geodesic of photons. We compare the emission from the whole inter-galactic medium (IGM) and from the WHIM component (defined as the gas with a temperature in the range $10^5-10^7$ K) at three observational bands associated to thermal X-rays, thermal and kinematic Sunyaev-Zel'dovich effect, and radio emission. The synthetic maps produced by this procedure could be directly compared with existing observational maps and could be used as a guide for future observations with forthcoming instruments. The analysis of the different emissions associated to a high-resolution galaxy cluster is in broad agreement with previous simulated and observational estimates of both gas components.

Radio mode feedback, associated with the propagation of powerful outflows in active galaxies, is a crucial ingredient in galaxy evolution. Extragalactic jets are well collimated and relativistic, both in terms of thermodynamics and kinematics. They generate strong shocks in the ambient medium, associated with observed hotspots, and carve cavities that are filled with the shocked jet flow. In this Letter, we compare the pressure evolution in the hotspot and the cavity generated by relativistic and classical jets. Our results show that the classical approach underestimates the cavity pressure by a factor larger or equal to 2 for a given shocked volume during the whole active phase. The tension between both approaches can only be alleviated by unrealistic jet flow densities or gigantic jet areas in the classical case. As a consequence, the efficiency of a relativistic jet heating the ambient is typically around 20% larger compared with a classical jet, and the heated volume is 2 to 10 times larger during the time evolution. This conflict translates into two substantially disparate manners, both spatially and temporal, of heating the ambient medium. These differences are expected to have relevant implications on the star formation rates of the host galaxies and their evolution.

The finding that massive galaxies grow with cosmic time fired the starting gun for the search of objects which could have survived up to the present day without suffering substantial changes (neither in their structures, neither in their stellar populations). Nevertheless, and despite the community efforts, up to now only one firm candidate to be considered one of these relics is known: NGC 1277. Curiously, this galaxy is located at the centre of one of the most rich near galaxy clusters: Perseus. Is its location a matter of chance? Should relic hunters focus their search on galaxy clusters? In order to reply this question, we have performed a simultaneous and analogous analysis using simulations (Millennium I-WMAP7) and observations (New York University Value-Added Galaxy Catalogue). Our results in both frameworks agree: it is more probable to find relics in high density environments.

We study how the gas in a sample of galaxies (M* > 10e9 Msun) in clusters, obtained in a cosmological simulation, is affected by the interaction with the intra-cluster medium (ICM). The dynamical state of each elemental parcel of gas is studied using the total energy. At z ~ 2, the galaxies in the simulation are evenly distributed within clusters, moving later on towards more central locations. In this process, gas from the ICM is accreted and mixed with the gas in the galactic halo. Simultaneously, the interaction with the environment removes part of the gas. A characteristic stellar mass around M* ~ 10e10 Msun appears as a threshold marking two differentiated behaviours. Below this mass, galaxies are located at the external part of clusters and have eccentric orbits. The effect of the interaction with the environment is marginal. Above, galaxies are mainly located at the inner part of clusters with mostly radial orbits with low velocities. In these massive systems, part of the gas, strongly correlated with the stellar mass of the galaxy, is removed. The amount of removed gas is sub-dominant compared with the quantity of retained gas which is continuously influenced by the hot gas coming from the ICM. The analysis of individual galaxies reveals the existence of a complex pattern of flows, turbulence and a constant fuelling of gas to the hot corona from the ICM that could make the global effect of the interaction of galaxies with their environment to be substantially less dramatic than previously expected.

We study the preferred environments of $z \sim 0$ massive relic galaxies ($M_\star \gtrsim 10^{10}~\mathrm{M_\odot}$ galaxies with little or no growth from star formation or mergers since $z \sim 2$). Significantly, we carry out our analysis on both a large cosmological simulation and an observed galaxy catalogue. Working on the Millennium I-WMAP7 simulation we show that the fraction of today massive objects which have grown less than 10 per cent in mass since $z \sim 2$ is ~0.04 per cent for the whole massive galaxy population with $M_\star > 10^{10}~\mathrm{M_\odot}$. This fraction rises to ~0.18 per cent in galaxy clusters, confirming that clusters help massive galaxies remain unaltered. Simulations also show that massive relic galaxies tend to be closer to cluster centres than other massive galaxies. Using the New York University Value-Added Galaxy Catalogue, and defining relics as $M_\star \gtrsim 10^{10}~\mathrm{M_\odot}$ early-type galaxies with colours compatible with single-stellar population ages older than 10 Gyr, and which occupy the bottom 5-percentile in the stellar mass-size distribution, we find $1.11 \pm 0.05$ per cent of relics among massive galaxies. This fraction rises to $2.4 \pm 0.4$ per cent in high-density environments. Our findings point in the same direction as the works by Poggianti et al. and Stringer et al. Our results may reflect the fact that the cores of the clusters are created very early on, hence the centres host the first cluster members. Near the centres, high-velocity dispersions and harassment help cluster core members avoid the growth of an accreted stellar envelope via mergers, while a hot intracluster medium prevents cold gas from reaching the galaxies, inhibiting star formation.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

Using a sample of cosmic voids identified in the Sloan Digital Sky Survey Data Release 7, we study the star formation activity of void galaxies. The properties of galaxies living in voids are compared with those of galaxies living in the void shells and with a control sample, representing the general galaxy population. Void galaxies appear to form stars more efficiently than shell galaxies and the control sample. This result can not be interpreted as a consequence of the bias towards low masses in underdense regions, as void galaxy subsamples with the same mass distribution as the control sample also show statistically different specific star formation rates. This highlights the fact that galaxy evolution in voids is slower with respect to the evolution of the general population. Nevertheless, when only the star forming galaxies are considered, we find that the star formation rate is insensitive to the environment, as the main sequence is remarkably constant in the three samples under consideration. This fact implies that environmental effects manifest themselves as fast quenching mechanisms, while leaving the non-quenched galaxies almost unaffected, as their star formation activity is largely regulated by the mass of their halo. We also analyse galaxy properties as a function of void-centric distance and find that the enhancement in the star formation activity with respect to the control sample is observable up to a radial distance 1.5 Rvoid. This result can be used as a suitable definition of void shells. Finally, we find that larger voids show an enhanced star formation activity in the shells with respect to their smaller counterparts, that could be related to the different dynamical evolution experienced by voids of different sizes.

The evolution of powerful extragalactic jets is not only interesting by itself, but also for its impact on the evolution of the host galaxy and its surroundings. We have performed long-term axisymmetric numerical simulations of relativistic jets with different powers to study their evolution through an environment with a pressure and density gradient. Our results show key differences in the evolution of jets with different powers in terms of the spatial and temporal scales of energy deposition. According to our results, the observed morphology in X-ray cavities requires that an important fraction of the jet's energetic budget is in the form of internal energy. Thus, light, lepton-dominated jets are favoured. In all cases, heating is mainly produced by shocks. Cavity overpressure is sustained by an important population of thermal particles. Our simulations reproduce the cool-core structure in projected, luminosity-weighted temperature. We have performed an additional simulation of a slow, massive jet and discuss the differences with its relativistic counterparts. Important qualitative and quantitative differences are found between the non-relativistic and the relativistic jets. Our conclusions point towards a dual-mode of AGN kinetic feedback, depending on the jet power.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

The massive exploitation of cosmic voids for precision cosmology in the upcoming dark energy experiments, requires a robust understanding of their internal structure, particularly of their density profile. We show that the void density profile is insensitive to the void radius both in a catalogue of observed voids and in voids from a large cosmological simulation. However, the observed and simulated voids display remarkably different profile shapes, with the former having much steeper profiles than the latter. Sparsity can not be the main reason for this discrepancy, as we demonstrate that the profile can be recovered with reasonable accuracy even with very sparse samples of tracers. On the other hand, the observed profile shows a significant dependence on the galaxy sample used to trace the matter distribution. Samples including low-mass galaxies lead to shallower profiles with respect to the samples where only massive galaxies are used, as faint galaxies live closer to the void centre. We argue that galaxies are biased tracers when used to probe the matter distribution within voids.

We present the results of a numerical adaptive mesh refinement hydrodynamical and N-body simulation in a $\Lambda CDM$ cosmology. We focus on the analysis of the main properties of massive galaxies ($M_* > 10^{11}\,M_{\odot}$) at $z=0$. For all the massive virtual galaxies we carry out a careful study of their one dimensional density, luminosity, velocity dispersion, and stellar population profiles. In order to best compare with observational data, the method to estimate the velocity dispersion is calibrated by using an approach similar to that performed in the observations, based on the stellar populations of the simulated galaxies. With these ingredients, we discuss the different properties of massive galaxies in our sample according to their morphological types, accretion histories and dynamical properties. We find that the galaxy merging history is the leading actor in shaping the massive galaxies that we see nowadays. Indeed, galaxies having experienced a turbulent life are the most massive in the sample and show the steepest metallicity gradients. Beside the importance of merging, only a small fraction of the final stellar mass has been formed ex-situ (10-50%), while the majority of the stars formed within the galaxy. These accreted stars are significantly older and less metallic than the stars formed in-situ and tend to occupy the most external regions of the galaxies.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

The number of present-day massive galaxies that has survived untouched since their formation at high-z is an important observational constraint to the hierarchical galaxy formation models. Using three different semianalytical models based on the Millenium simulation, we quantify the expected fraction and number densities of the massive galaxies form at z>2 which have evolved in stellar mass less than 10% and 30%. We find that only a small fraction of the massive galaxies already form at z~2 have remained almost unaltered since their formation (<2% with Delta_M*/M*<0.1 and <8% with Delta_M*/M*<0.3). These fractions correspond to the following number densities of massive relics in the present-day Universe: ~1.2x10^-6 Mpc^-3 for Delta_M*/M*<0.1 and ~5.7x10^-6 Mpc^-3 for Delta_M*/M*<0.3. The observed number of relic candidates found in the nearby Universe is today pretty uncertain (with uncertainties up to a factor of ~100) preventing to establish a firm conclusion about the goodness of current theoretical expectations to predict such important number.

Eulerian cosmological codes are especially suited to properly describe the low density regions. This property makes this class of codes excellent tools to study the formation and evolution of cosmic voids. Following such ideas, we present the results of an Eulerian adaptive mesh refinement (AMR) hydrodynamical and N-body simulation, that contrary to the common practice, has been designed to refine the computational grid in the underdense regions of the simulated volume. Thus, the void regions are better described due to the combined effect of the Eulerian character of the numerical technique and the use of high numerical resolution from the AMR approach. To analyse the outcome of this simulation, we have constructed a new void finder optimally suited to find the hierarchy of voids in AMR simulations. The algorithm identifies voids starting from the cells with least density and highest velocity divergence and then expanding the underdense volume until reaching the void walls, defined from the steepness of the density gradient. At redshift z=0, in a cosmological box of comoving side length 100 Mpc/h, we identify hundreds of voids with sizes up to 17 Mpc/h and typical density contrast of -0.8, which show a complex morphology and an intricate hierarchy of nested structures. The analysis of their mass density profile leads to the conclusion that a universal density profile can be applied to voids of any size, density, morphology and redshift.

Shock waves developed during the formation and evolution of cosmic structures encode crucial information on the hierarchical formation of the Universe. We analyze an Eulerian AMR hydro + N-body simulation in a $\Lambda$CDM cosmology focused on the study of cosmological shock waves. The combination of a shock-capturing algorithm together with the use of a halo finder allows us to study the morphological structures of the shock patterns, the statistical properties of shocked cells, and the correlations between the cosmological shock waves appearing at different scales and the properties of the haloes harbouring them. The shocks in the simulation can be split into two broad classes: internal weak shocks related with evolutionary events within haloes, and external strong shocks associated with large-scale events. The shock distribution function contains information on the abundances and strength of the different shocks, and it can be fitted by a double power law with a break in the slope around a Mach number of 20. We introduce a generalised scaling relation that correlates the average Mach numbers within the virial radius of haloes and their virial masses. In this plane, haloes occupy different areas according to their early evolutionary histories: those with a quiet evolution have an almost constant Mach number independently of their masses, whereas haloes undergoing significant merger events very early in their evolution show a linear dependence with their masses. At high redshift, the halo distribution in this scaling relation forms a L-like pattern that changes due to the evolution of the haloes. The analysis of the propagation speed and size of the shock waves around haloes could give some hints on the formation time and main features of the haloes. (Abridged)

Minor merging has been postulated as the most likely evolutionary path to produce the increase in size and mass observed in the massive galaxies since z$\sim$2. In this Letter, we test directly this hypothesis comparing the population of satellites around massive galaxies in cosmological simulations versus the observations. We use state-of-the-art, publically available, Millennium I and II simulations and the associated semi-analytical galaxy catalogues to explore the time evolution of the fraction of massive galaxies that have satellites, the number of satellites per galaxy, the projected distance at which the satellite locate from the host galaxy, and the mass ratio between the host galaxies and their satellites. The three virtual galaxy catalogues considered here, overproduce the fraction of galaxies with satellites by a factor ranging between 1.5 and 6 depending on the epoch, whereas the mean projected distance and ratio of the satellite mass over host mass are in closer agreement with data. The larger pull of satellites in the semi-analytical samples could suggest that the size evolution found in previous hydrodynamical simulations is an artifact due to the larger number of infalling satellites compared to the real Universe. These results advise to revise the physical ingredients implemented in the semi-analytical models in order to reconcile the observed and computed fraction of galaxies with satellites, and eventually, it would leave some room to other mechanisms explaining the galaxy size growth not related to the minor merging.

Galactic jets are powerful energy sources reheating the intra-cluster medium in galaxy clusters. Their crucial role in the cosmic puzzle, motivated by observations, has been established by a great number of numerical simulations missing the relativistic nature of these jets. We present the first relativistic simulations of the very long term evolution of realistic galactic jets. Unexpectedly, our results show no buoyant bubbles, but large cocoon regions compatible with the observed X-ray cavities. The reheating is more efficient and faster than in previous scenarios, and it is produced by the shock wave driven by the jet, that survives for several hundreds of Myrs. Therefore, the X-ray cavities in clusters produced by powerful relativistic jets would remain confined by weak shocks for extremely long periods, whose detection could be an observational challenge.

We study whether non-linear gravitational effects of relic neutrinos on the development of clustering and large-scale structure may be observable by weak gravitational lensing. We compute the density profile of relic massive neutrinos in a spherical model of a cluster of galaxies, for several neutrino mass schemes and cluster masses. Relic neutrinos add a small perturbation to the mass profile, making it more extended in the outer parts. In principle, this non-linear neutrino perturbation is detectable in an all-sky weak lensing survey such as EUCLID by averaging the shear profile of a large fraction of the visible massive clusters in the universe, or from its signature in the general weak lensing power spectrum or its cross-spectrum with galaxies. However, correctly modeling the distribution of mass in baryons and cold dark matter and suppressing any systematic errors to the accuracy required for detecting this neutrino perturbation is severely challenging.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

We present and test a new halo finder based on the spherical overdensity (SO) method. This new adaptive spherical overdensity halo finder (ASOHF) is able to identify dark matter haloes and their substructures (subhaloes) down to the scales allowed by the analysed simulations. The code has been especially designed for the adaptive mesh refinement cosmological codes, although it can be used as a stand-alone halo finder for N-body codes. It has been optimised for the purpose of building the merger tree of the haloes. In order to verify the viability of this new tool, we have developed a set of bed tests that allows us to estimate the performance of the finder. Finally, we apply the halo finder to a cosmological simulation and compare the results obtained to those given by other well known publicly available halo finders.

We present the results of an Eulerian adaptive mesh refinement (AMR) hydrodynamical and N-body simulation in a $\Lambda$CDM cosmology. The simulation incorporates common cooling and heating processes for primordial gas. A specific halo finder has been designed and applied in order to extract a sample of galaxy clusters directly obtained from the simulation without considering any resimulating scheme. We have studied the evolutionary history of the cluster halos, and classified them into three categories depending on the merger events they have undergone: major mergers, minor mergers, and relaxed clusters. The main properties of each one of these classes and the differences among them are discussed. The collisions among galaxy clusters are produced naturally by the non-linear evolution in the simulated cosmological volume, no controlled collisions have been considered. We pay special attention to discuss the role of merger events as a source of feedback and reheating, and their effects on the existence of cool cores in galaxy clusters, as well as in the scaling relations.

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

We discuss the emission of gravitational radiation produced in encounters of dark matter galactic halos. To this aim we perform a number of numerical simulations of typical galaxy mergers, computing the associated gravitational radiation waveforms as well as the energy released in the processes. Our simulations yield dimensionless gravitational wave amplitudes of the order of $10^{-13}$ and gravitational wave frequencies of the order of $10^{-16}$ Hz, when the galaxies are located at a distance of 10 Mpc. These values are of the same order as those arising in the gravitational radiation originated by strong variations of the gravitational field in the early Universe, and therefore, such gravitational waves cannot be directly observed by ground-based detectors. We discuss the feasibility of an indirect detection by means of the B-mode polarization of the Cosmic Microwave Background (CMB) induced by such waves. Our results show that the gravitational waves from encounters of dark matter galactic halos leave much too small an imprint on the CMB polarization to be actually observed with ongoing and future missions.

It has been recently argued through numerical work that rotating stars with a high degree of differential rotation are dynamically unstable against bar-mode deformation, even for values of the ratio of rotational kinetic energy to gravitational potential energy as low as O(0.01). This may have implications for gravitational wave astronomy in high-frequency sources such as core collapse supernovae. In this paper we present high-resolution simulations, performed with an adaptive mesh refinement hydrodynamics code, of such low T/|W| bar-mode instability. The complex morphological features involved in the nonlinear dynamics of the instability are revealed in our simulations, which show that the excitation of Kelvin-Helmholtz-like fluid modes outside the corotation radius of the star leads to the saturation of the bar-mode deformation. While the overall trends reported in an earlier investigation are confirmed by our work, we also find that numerical resolution plays an important role during the long-term, nonlinear behaviour of the instability, which has implications on the dynamics of rotating stars and on the attainable amplitudes of the associated gravitational wave signals.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

We present a complete analysis, which includes morphology, kinematics, stellar populations and N-body simulations of CGCG 480-022, the most distant (cz=14317 km/s) isolated galaxy studied so far in such detail. The results all support the hypothesis that this galaxy has suffered a major merger event with a companion of ~0.1 times its mass. Morphology reveals the presence of a circumnuclear ring and possibly further ring debris. The radial velocity curve looks symmetrical, whilst the velocity dispersion increases with radius reaching values that do not correspond to a virialized system. Moreover, this galaxy deviates significantly from the Fundamental Plane and the Faber-Jackson relation. The stellar population analysis show that the ring is younger and more metal-rich, which suggest that it has undergone a fairly recent burst of star formation. Both morphological and dynamical results are in broad agreeement with our N-body simulations.

A method for ray-tracing through n-body simulations has been recently proposed. It is based on a periodic universe covered by simulation boxes. Photons move along appropriate directions to avoid periodicity effects. Here, an improved version of this method is applied to simulate lensed CMB maps and maps of lens deformations. Particle mesh n-body simulations with appropriate boxes and resolutions are used to evolve the nonlinear inhomogeneities until present time. The resulting maps are statistically analyzed to look for deviations from Gaussianity. These deviations are measured --for the first time-- using correlations for configurations of n directions (n being smaller than 7). A wide range of angular scales are considered. Some interesting prospects are outlined.

A new cosmological multidimensional hydrodynamic and N-body code based on an Adaptive Mesh Refinement scheme is described and tested. The hydro part is based on modern high-resolution shock-capturing techniques, whereas N-body approach is based on the Particle Mesh method. The code has been specifically designed for cosmological applications. Tests including shocks, strong gradients, and gravity have been considered. A cosmological test based on Santa Barbara cluster is also presented. The usefulness of the code is discussed. In particular, this powerful tool is expected to be appropriate to describe the evolution of the hot gas component located inside asymmetric cosmological structures.

We work in the framework of an inflationary cold dark matter universe with cosmological constant, in which the cosmological inhomogeneities are considered as gravitational lenses for the CMB photons. This lensing deforms the angular distribution of the CMB maps in such a way that the induced deformations are not Gaussian. Our main goal is the estimation of the deviations with respect to Gaussianity appeared in the distribution of deformations. In the new approach used in this paper, matter is evolved with a particle-mesh N-body code and, then, an useful ray-tracing technique designed to calculate the correlations of the lens deformations induced by nonlinear structures is applied. Our approach is described in detail and tested. Various correlations are estimated at an appropriate angular scale. The resulting values point out both deviations with respect to Gaussian statistics and a low level of correlation in the lens deformations.

We investigate the evolution of angular momentum in SPH simulations of galaxy formation, paying particular attention to artificial numerical effects. We find that a cold gas disc forming in an ambient hot gas halo receives a strong hydrodynamic torque from the hot gas. By splitting the hydrodynamic force into artificial viscosity and pressure gradients, we find that the angular momentum transport is caused not by the artificial viscosity but by the pressure gradients. Using simple test simulations of shear flows, we conclude that the pressure gradient-based viscosity can be divided into two components: one due to the noisiness of SPH and the other to ram pressure. The former is problematic even with very high resolution because increasing resolution does not reduce the noisiness. On the other hand, the ram pressure effect appears only when a cold gas disc or sheet does not contain enough particles. In such a case, holes form in the disc or sheet, and then ram pressure from intra-hole hot gas, causes significant deceleration. In simulations of galactic disc formation, star formation usually decreases the number of cold gas particles, and hole formation leads to the fragmentation of the disc. To circumvent these problem, we modify the SPH algorithm, decoupling the cold from the hot gas phases, i.e. inhibiting the hydrodynamic interaction between cold and hot particles. This, a crude modelling of a multi-phase fluid in SPH cosmological simulations, leads to the formation of smooth extended cold gas discs and to better numerical convergence. The decoupling is applicable in so far as the self-gravitating gas disc with negligible external pressure is a good approximation for a cold gas disc. (abridged)

Non-linear evolution of cosmological energy density fluctuations triggers deviations from Gaussianity in the temperature distribution of the cosmic microwave background. A method to estimate these deviations is proposed. N-body simulations -- in a $\Lambda$CDM cosmology -- are used to simulate the strongly non-linear evolution of cosmological structures. It is proved that these simulations can be combined with the potential approximation to calculate the statistical moments of the CMB anisotropies produced by non-linear gravity. Some of these moments are computed and the resulting values are different from those corresponding to Gaussianity.

We use a three dimensional hydrodynamical code to simulate the effect of energy injection on cooling flows in the intracluster medium. Specifically, we compare a simulation of a 10$^{15}$ $M_\odot$ cluster with radiative cooling only, with a second simulation in which thermal energy is injected 31 $\kpc$ off-centre, over 64 kpc$^{3}$ at a rate of $4.9\times 10^{44} \ergs$ for 50 Myr. The heat injection forms a hot, low density bubble which quickly rises, dragging behind it material from the cluster core. The rising bubble pushes with it a shell of gas which expands and cools. We find the appearance of the bubble in X-ray temperature and luminosity to be in good qualitative agreement with recent {\it Chandra} observations of cluster cores. Toward the end of the simulation, at 600 Myr, the displaced gas begins to fall back toward the core, and the subsequent turbulence is very efficient at mixing the low and high entropy gas. The result is that the cooling flow is disrupted for up to $\sim 50$ Myr after the injection of energy ceases. Thus, this mechanism provides a very efficient method for regulating cooling flows, if the injection events occur with a 1:1 duty cycle.

Recent observations of high velocity clouds (HVCs) have revealed compression fronts and tail shaped features of HI suggesting that they are interacting with external medium. We perform 3-D hydro-dynamical simulations of HVCs moving through a diffuse hot gaseous component, investigating the behaviour of both extra-galactic dark matter dominated HVCs and nearby pure gas clouds that may be accreting onto the Galactic disk via a Galactic fountain. Both scenarios can give rise to similar features as observed if the external medium has a density $>10^{-4}$ cm$^{-3}$. Observations suggest that this may be too high for a hot ionised halo or intergalactic gas and supports the Galactic fountain origin, a model that can also account for the high fraction of HVCs with tails and asymmetrical morphologies.

The observed morphological evolution since z ~ 0.5 within galaxy clusters provides evidence for a hierarchical universe. This evolution is driven by dynamical effects that operate within the cluster environment -- suppression of star-formation by ram-pressure and viscous stripping of the intra-galactic medium and tidal heating of disks by gravitational encounters.

We present the first 3-dimensional high resolution hydro-dynamical simulations of the interaction between the hot ionised intra-cluster medium and the cold interstellar medium of spiral galaxies. Ram pressure and turbulent/viscous stripping removes 100% of the atomic hydrogen content of luminous galaxies like the Milky Way within 100 million years. These mechanisms naturally account for the morphology of S0 galaxies, the rapid truncation of star formation implied by spectroscopic observations, as well as a host of observational data on the HI morphology of galaxies in clusters.

We compare the structure and substructure of dark matter halos in model universes dominated by collisional, strongly self interacting dark matter (SIDM) and collisionless, weakly interacting dark matter (CDM). While SIDM virialised halos are more nearly spherical than CDM halos, they can be rotationally flattened by as much as 20% in their inner regions. Substructure halos suffer ram-pressure truncation and drag which are more rapid and severe than their gravitational counterparts tidal stripping and dynamical friction. Lensing constraints on the size of galactic halos in clusters are a factor of two smaller than predicted by gravitational stripping, and the recent detection of tidal streams of stars escaping from the satellite galaxy Carina suggests that its tidal radius is close to its optical radius of a few hundred parsecs --- an order of magnitude smaller than predicted by CDM models but consistent with SIDM. The orbits of SIDM satellites suffer significant velocity bias $\sigma_{SIDM}/\sigma_{CDM}=0.85$ and are more circular than CDM, $\beta_{SIDM}} \approx 0.5$, in agreement with the inferred orbits of the Galaxy's satellites. In the limit of a short mean free path, SIDM halos have singular isothermal density profiles, thus in its simplest incarnation SIDM is inconsistent with galactic rotation curves.

Galaxy clusters are sources of gravitational radiation. The main aim of this paper is to give numerical estimates and theoretical description of the relevant features of the gravitational radiation coming from an appropriate spatial distributions of galaxy clusters. Since no analytical approaches are currently available to describe the strongly nonlinear regime, our numerical approach - combining numerical simulations with statistical arguments - seems to be an useful way of studying the main features of that radiation. Although far to be detectable with present technology, we advance some ideas about future observational strategies and its cosmological implications.

A rich galaxy cluster showing strong resemblance with the observed ones is simulated. Cold dark matter spectrum, Gaussian statistics, flat universe, and two components -- baryonic gas plus dark matter particles -- are considered. We have calculated the gravitational-wave output during the epoch of the fully nonlinear and nonsymmetric cluster evolution. The amplitudes and frequencies of the resulting gravitational waves are estimated. Since frequencies are very small --of the order of $10^{-17} Hz$ -- a complete pulse cannot be observed during an admissible integration time; nevertheless, it is proved that these waves can produce an interesting secular effect which appears to be observable with current technology.

Numerical simulations of galaxy clusters including two species -- baryonic gas and dark matter particles --are presented. Cold Dark Matter spectrum, Gaussian statistics and flat universe are assumed. The dark matter component is evolved numerically by means of a standard particle mesh method. The evolution of the baryonic component has been studied numerically by using a multidimensional (3D) hydrodynamical code based on modern high resolution shock capturing techniques. These techniques are specially designed for treating accurately complex flows in which shocks appear and interact. With this picture, the role of shock waves in the formation and evolution of rich galaxy clusters is analyzed. Our results display two well differenced morphologies of the shocked baryonic matter: filamentary at early epochs and quasi-spherical at low redshifts.

The applicability of the potential approximation in the case of open universes is tested. Great Attractor-like structures are considered in the test. Previous estimates of the Cosmic Microwave background anisotropies produced by these structures are analyzed and interpreted. The anisotropies corresponding to inhomogeneous ellipsoidal models are also computed. It is proved that, whatever the spatial symmetry may be, Great Attractor-like objects with extended cores (radius $\sim 10h^{-1}$),located at redshift $z=5.9$ in an open universe with density parameter $\Omega_{0}=0.2$, produce secondary gravitational anisotropies of the order of $10^{-5}$ on angular scales of a few degrees. This anisotropy appears to be an integrated effect along the photon geodesics. Its angular scale is much greater than that subtended by the Great Attractor itself. This is understood taking into account that the integrated effect is produced by the variations of the gravitational potential, which seem to be important in large regions subtending angular scales of various degrees. As a result of the large size of these regions, the spatial curvature of the universe becomes important and, consequently, significant errors ($\sim 30$ per cent) arise in estimations based on the potential approximation.

A cosmological multidimensional hydrodynamic code is described and tested. This code is based on modern high-resolution shock-capturing techniques. It can make use of a linear or a parabolic cell reconstruction as well as an approximate Riemann solver. The code has been specifically designed for cosmological applications. Two tests including shocks have been considered: the first one is a standard shock tube and the second test involves a spherically symmetric shock. Various additional cosmological tests are also presented. In this way, the performance of the code is proved. The usefulness of the code is discussed; in particular, this powerful tool is expected to be useful in order to study the evolution of the hot gas component located inside nonsymmetric cosmological structures.

Previous estimates of the microwave background anisotropies produced by freely falling spherical clusters are discussed. These estimates are based on the Swiss-Cheese and Tolman-Bondi models. It is proved that these models give only upper limits to the anisotropies produced by the observed galaxy clusters. By using spherically symmetric codes including pressureless matter and a hot baryonic gas, new upper limits are obtained. The contributions of the hot gas and the pressureless component to the total anisotropy are compared. The effects produced by the pressure are proved to be negligible; hence, estimations of the cluster anisotropies based on N-body simulations are hereafter justified. After the phenomenon of violent relaxation, any realistic rich cluster can only produce small anisotropies with amplitudes of order $10^{-7}$. During the rapid process of violent relaxation, the anisotropies produced by nonlinear clusters are expected to range in the interval $(10^{-6},10^{-5})$. The angular scales of these anisotropies are discussed.