Showing 1–50 of 57 results for author: Quilis, V

Vortex-p: a Helmholtz-Hodge and Reynolds decomposition algorithm for particle-based simulations

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis, Frederick Groth, Tirso Marin-Gilabert, Klaus Dolag

Abstract: Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent velocity fields are the Helmholtz-Hodge decomposition (compressive+solenoidal; HHD) and the Reynolds decomposition (bulk+turbulent; RD). These problems are relativ… ▽ More Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent velocity fields are the Helmholtz-Hodge decomposition (compressive+solenoidal; HHD) and the Reynolds decomposition (bulk+turbulent; RD). These problems are relatively simple to perform numerically for uniformly-sampled data, such as the one emerging from Eulerian, fix-grid simulations; but their computation is remarkably more complex in the case of non-uniformly sampled data, such as the one stemming from particle-based or meshless simulations. In this paper, we describe, implement and test vortex-p, a publicly available tool evolved from the vortex code, to perform both these decompositions upon the velocity fields of particle-based simulations, either from smoothed particle hydrodynamics (SPH), moving-mesh or meshless codes. The algorithm relies on the creation of an ad-hoc adaptive mesh refinement (AMR) set of grids, on which the input velocity field is represented. HHD is then addressed by means of elliptic solvers, while for the RD we adapt an iterative, multi-scale filter. We perform a series of idealised tests to assess the accuracy, convergence and scaling of the code. Finally, we present some applications of the code to various SPH and meshless finite-mass (MFM) simulations of galaxy clusters performed with OpenGadget3, with different resolutions and physics, to showcase the capabilities of the code. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Computer Physics Communications. 20 pages, 13 figures + appendices. The code is available in https://github.com/dvallesp/vortex-p and the documentation can be accessed in https://vortex-particles.github.io

Journal ref: Computer Physics Communications, Volume 304, November 2024, 109305

Cosmic accretion shocks as a tool to measure the dark matter mass of galaxy clusters

Authors: David Vallés-Pérez, Vicent Quilis, Susana Planelles

Abstract: Cosmological accretion shocks created during the formation of galaxy clusters are a ubiquitous phenomenon all around the Universe. These shocks, and their features, are intimately related with the gravitational energy put into play during galaxy cluster formation. Studying a sample of simulated galaxy clusters and their associated accretion shocks, we show that objects in our sample sit in a plane… ▽ More Cosmological accretion shocks created during the formation of galaxy clusters are a ubiquitous phenomenon all around the Universe. These shocks, and their features, are intimately related with the gravitational energy put into play during galaxy cluster formation. Studying a sample of simulated galaxy clusters and their associated accretion shocks, we show that objects in our sample sit in a plane within the three dimensional-space of cluster total mass, shock radius, and Mach number (a measure of shock intensity). Using this relation, and considering that forthcoming new observations will be able to measure shock radii and intensities, we put forward the idea that the dark matter content of galaxy clusters could be indirectly measured with an error up to around 30 per cent at the $1σ$ confidence level. This procedure would be a new and independent method to measure the dark matter mass in cosmic structures, and a novel constraint to the accepted $Λ$CDM paradigm. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Nature Astronomy. This version corresponds to the preprint, prior to any formal peer review. The published version, including peer-review improvements, is freely accessible through the SharedIt platform in https://rdcu.be/dMm5v. This version: 13 pages, 3 figures (main paper + methods), plus supplementary material (9 pages, 10 figures)

CAVITY, Calar Alto Void Integral-field Treasury surveY and project extension

Authors: I. Pérez, S. Verley, L. Sánchez-Menguiano, T. Ruiz-Lara, R. García-Benito, S. Duarte Puertas, A. Jiménez, J. Domínguez-Gómez, D. Espada, R. F. Peletier, J. Román, M. I. Rodríguez, P. Sánchez Alarcón, M. Argudo-Fernández, G. Torres-Ríos, B. Bidaran, M. Alcázar-Laynez, R. van de Weygaert, S. F. Sánchez, U. Lisenfeld, A. Zurita, E. Florido, J. M. van der Hulst, G. Blázquez-Calero, P. Villalba-González , et al. (36 additional authors not shown)

Abstract: We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evoluti… ▽ More We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evolution as they are largely unaffected by the complex physical processes that transform galaxies in high-density environments. These void galaxies can hold the key as well to answer current challenges to the $Λ$CDM paradigm. The Calar Alto Void Integral-field Treasury surveY (CAVITY) is a Legacy project approved by the Calar Alto Observatory to obtain spatially resolved spectroscopic information of $\sim300$ void galaxies in the Local Universe (0.005 < z < 0.050) covering from -17.0 to -21.5 in $\rm r$ band absolute magnitude. It officially started in January 2021 and has been awarded 110 useful dark observing nights at the 3.5 m telescope using the PMAS spectrograph. Complementary follow-up projects including deep optical imaging, integrated, as well as resolved CO data, and integrated HI spectra, have joint the PMAS observations and naturally complete the scientific aim of characterising galaxies in cosmic voids. The extension data has been denominated CAVITY+. The data will be available to the whole community in different data releases, the first of which is planned for July 2024, and it will provide the community with PMAS data cubes for around 100 void galaxies through a user friendly, and well documented, database platform. We present here the survey, sample selection, data reduction, quality control schemes, science goals, and some examples of the scientific power of the CAVITY and CAVITY+ data. △ Less

Submitted 24 May, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 19 pages, 11 figures. Accepted for publication in A&A - Replacement after A&A minor language edition

The imprints of galaxy cluster internal dynamics on the Sunyaev-Zeldovich effect

Authors: Óscar Monllor-Berbegal, David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: Forthcoming measurements of the Sunyaev-Zeldovich (SZ) effect in galaxy clusters will dramatically improve our understanding of the main intra-cluster medium (ICM) properties and how they depend on the particular thermal and dynamical state of the associated clusters. Using a sample of simulated galaxy clusters we assess the impact of the ICM internal dynamics on both the thermal and kinetic SZ ef… ▽ More Forthcoming measurements of the Sunyaev-Zeldovich (SZ) effect in galaxy clusters will dramatically improve our understanding of the main intra-cluster medium (ICM) properties and how they depend on the particular thermal and dynamical state of the associated clusters. Using a sample of simulated galaxy clusters we assess the impact of the ICM internal dynamics on both the thermal and kinetic SZ effects (tSZ and kSZ, respectively). We produce synthetic maps of the SZ effect, for the simulated clusters. For each galaxy cluster in the sample, its dynamical state is estimated by using a combination of well-established indicators. We use the correlations between SZ maps and cluster dynamical state, to look for the imprints of the evolutionary events, mainly mergers, on the SZ signals. The kinetic effect shows a remarkable correlation with the dynamical state: unrelaxed clusters present a higher radial profile and an overall stronger signal at all masses and radii. Furthermore, the kSZ signal is correlated with rotation for relaxed clusters, while for the disturbed systems the effect is dominated by other motions such as bulk flows, turbulence, etc. The kSZ effect shows a dipolar pattern when averaging over cluster dynamical classes, especially for the relaxed population. This feature can be exploited to stack multiple kSZ maps in order to recover a stronger dipole signal that would be correlated with the global rotation properties of the sample. The SZ effect can be used as a tool to estimate the dynamical state of galaxy clusters, especially to segregate those clusters with a quiescent evolution from those with a rich record of recent merger events. Our results suggest that the forthcoming observational data measuring the SZ signal in clusters could be used as a complementary strategy to classify the evolutionary history of galaxy clusters. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 13 pages. 11 figures. Accepted for publication in A&A

On the choice of the most suitable indicator for the assembly state of dark matter haloes through cosmic time

Authors: David Vallés-Pérez, Susana Planelles, Óscar Monllor-Berbegal, Vicent Quilis

Abstract: The dynamical state and morphological features of galaxies and galaxy clusters, and their high-redshift precursors, are tightly connected with their assembly history, encoding crucial information about the formation and evolution of such cosmic structures. As a first step towards finding an optimal indicator of the assembly state of observed structures, we use a cosmological simulation of a modera… ▽ More The dynamical state and morphological features of galaxies and galaxy clusters, and their high-redshift precursors, are tightly connected with their assembly history, encoding crucial information about the formation and evolution of such cosmic structures. As a first step towards finding an optimal indicator of the assembly state of observed structures, we use a cosmological simulation of a moderate volume to critically examine the best definition of an indicator that is able to discriminate dark matter haloes undergoing mergers and/or strong accretion from haloes experimenting a relaxed evolution. Using a combination of centre offset, virial ratio, mean radial velocity, sparsity and ellipticity of the dark matter halo, we study how the thresholds on these parameters, as well as their relative weights, should evolve with redshift to provide the best classification possible. This allows us to split a sample of haloes in a totally relaxed, a marginally relaxed and an unrelaxed subsamples. The resulting classification strongly correlates with the merging activity obtained from the analysis of complete merger trees extracted from whole simulation data. The results on how the different indicators depend on redshift and halo mass, and their optimal combination to better match the true assembly history of haloes, could constitute relevant hints to find a suitable set of indicators applicable to observational data. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 15 pages, 10 figures; accepted for publication in MNRAS

Sociology and hierarchy of voids: A study of seven nearby CAVITY galaxy voids and their dynamical CosmicFlows-3 environment

Authors: H. M. Courtois, R. van de Weygaert, M. Aubert, D. Pomarède, D. Guinet, J. Domínguez-Gómez, E. Florido, L. Galbany, R. García-Benito, J. M. van der Hulst, K. Kreckel, R. E. Miura, I. Pérez, S. Planelles, V. Quilis, J. Román, M. Sánchez-Portal

Abstract: Context. The present study addresses a key question related to our understanding of the relation between void galaxies and their environment: the relationship between luminous and dark matter in and around voids. Aims. To explore the extent to which local Universe voids are empty of matter, we study the full (dark+luminous) matter content of seven nearby cosmic voids that are fully contained withi… ▽ More Context. The present study addresses a key question related to our understanding of the relation between void galaxies and their environment: the relationship between luminous and dark matter in and around voids. Aims. To explore the extent to which local Universe voids are empty of matter, we study the full (dark+luminous) matter content of seven nearby cosmic voids that are fully contained within the CosmicFlows-3 volume. Methods. We obtained the matter-density profiles of seven cosmic voids using two independent methods. These were built from the galaxy redshift space two-point correlation function in conjunction with peculiar velocity gradients from the CosmicFlows-3 dataset. Results. The results are striking, because when the redshift survey is used, all voids show a radial positive gradient of galaxies, while based on the dynamical analysis, only three of these voids display a clear underdensity of matter in their center. Conclusions. This work constitutes the most detailed observational analysis of voids conducted so far, and shows that void emptiness should be derived from dynamical information. From this limited study, the Hercules void appears to be the best candidate for a local Universe pure "pristine volume", expanding in three directions with no dark matter located in that void. △ Less

Submitted 8 March, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: Submitted A\&A Nov 29, 2022 - AA/2022/45578 / Accepted March 3rd, 2023

Journal ref: A&A 673, A38 (2023)

The halo finding problem revisited: a deep revision of the ASOHF code

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: Context. New-generation cosmological simulations are providing huge amounts of data, whose analysis becomes itself a cutting-edge computational problem. In particular, the identification of gravitationally bound structures, known as halo finding, is one of the main analyses. A handful of codes developed to tackle this task have been presented during the last years. Aims. We present a deep revisi… ▽ More Context. New-generation cosmological simulations are providing huge amounts of data, whose analysis becomes itself a cutting-edge computational problem. In particular, the identification of gravitationally bound structures, known as halo finding, is one of the main analyses. A handful of codes developed to tackle this task have been presented during the last years. Aims. We present a deep revision of the already existing code ASOHF. The algorithm has been throughfully redesigned in order to improve its capabilities to find bound structures and substructures, both using dark matter particles and stars, its parallel performance, and its abilities to handle simulation outputs with vast amounts of particles. This upgraded version of ASOHF is conceived to be a publicly available tool. Methods. A battery of idealised and realistic tests are presented in order to assess the performance of the new version of the halo finder. Results. In the idealised tests, ASOHF produces excellent results, being able to find virtually all the structures and substructures placed within the computational domain. When applied to realistic data from simulations, the performance of our finder is fully consistent with the results from other commonly used halo finders, with remarkable performance in substructure detection. Besides, ASOHF turns out to be extremely efficient in terms of computational cost. Conclusions. We present a public, deeply revised version of the ASOHF halo finder. The new version of the code produces remarkable results finding haloes and subhaloes in cosmological simulations, with an excellent parallel performance and with a contained computational cost. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 21 pages, 15 figures. Accepted for publication in A&A. Find the code repository in https://github.com/dvallesp/ASOHF and the code documentation in https://asohf.github.io/

Journal ref: A&A 664, A42 (2022)

Long-term FRII jet evolution in dense environments

Authors: Manel Perucho, José María Martí, Vicent Quilis

Abstract: We present long-term numerical three-dimensional simulations of a relativistic outflow propagating through a galactic ambient medium and environment, up to distances $\sim 100$~kpc. Our aim is to study the role of dense media in the global dynamics of the radio source. We use a relativistic gas equation of state, and a basic description of thermal cooling terms. In previous work, we showed that a… ▽ More We present long-term numerical three-dimensional simulations of a relativistic outflow propagating through a galactic ambient medium and environment, up to distances $\sim 100$~kpc. Our aim is to study the role of dense media in the global dynamics of the radio source. We use a relativistic gas equation of state, and a basic description of thermal cooling terms. In previous work, we showed that a linear perturbation could enhance the jet propagation during the early phases of evolution, by introducing obliquity to the jet reverse shock. Here, we show that this effect is reduced in denser media. We find that the dentist-drill effect acts earlier, due to slower jet propagation and an increased growth of the helical instability. The global morphology of the jet is less elongated, with more prominent lobes. The fundamental physical parameters of the jet generated structure derived from our simulations fall within the estimated values derived for FRII jets in the 3C sample. In agreement with previous axisymmetric and three dimensional simulations in lower density media, we conclude that shock heating of the interstellar and intergalactic media is very efficient in the case of powerful, relativistic jets. △ Less

Submitted 7 December, 2021; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Void replenishment: how voids accrete matter over cosmic history

Authors: David Vallés-Pérez, Vicent Quilis, Susana Planelles

Abstract: Cosmic voids are underdense regions filling up most of the volume in the Universe. They are expected to emerge in regions comprising negative initial density fluctuations, and subsequently expand as the matter around them collapses and forms walls, filaments and clusters. We report results from the analysis of a cosmological simulation specially designed to accurately describe low-density regions,… ▽ More Cosmic voids are underdense regions filling up most of the volume in the Universe. They are expected to emerge in regions comprising negative initial density fluctuations, and subsequently expand as the matter around them collapses and forms walls, filaments and clusters. We report results from the analysis of a cosmological simulation specially designed to accurately describe low-density regions, such as cosmic voids. Contrary to the common expectation, we find that voids also experience significant mass inflows over cosmic history. On average, $10\%$ of the mass of voids in the sample at $z \sim 0$ is accreted from overdense regions, reaching values beyond $35\%$ for a significant fraction of voids. More than half of the mass entering the voids lingers on periods of time $\sim 10 \, \mathrm{Gyr}$ well inside them, reaching inner radii. This would imply that part of the gas lying inside voids at a given time proceeds from overdense regions (e.g., clusters or filaments), where it could have been pre-processed, thus challenging the scenario of galaxy formation in voids, and dissenting from the idea of being pristine environments. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: Accepted for publication in The Astrophysical Journal Letters (ApJL). 11 pages (main text: 6 pages), 3 figures and 3 supplementary figures included in the PDF; the movie on Supplementary Figure 3 can be accessed on https://www.uv.es/davape3/voids/figS3_animation.mp4

Journal ref: The Astrophysical Journal Letters, vol 920, L2 (2021)

Exploring the role of cosmological shock waves in the Dianoga simulations of galaxy clusters

Authors: Susana Planelles, Stefano Borgani, Vicent Quilis, Giuseppe Murante, Veronica Biffi, Elena Rasia, Klaus Dolag, Gian Luigi Granato, Cinthia Ragone-Figueroa

Abstract: Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star format… ▽ More Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalising the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume. △ Less

Submitted 22 August, 2021; originally announced August 2021.

Comments: 19 pages, 15 figures, accepted for publication in MNRAS

Troubled cosmic flows: turbulence, enstrophy and helicity from the assembly history of the intracluster medium

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: Both simulations and observations have shown that turbulence is a pervasive phenomenon in cosmic scenarios, yet it is particularly difficult to model numerically due to its intrinsically multiscale character which demands high resolutions. Additionally, turbulence is tightly connected to the dynamical state and the formation history of galaxies and galaxy clusters, producing a diverse phenomenlogy… ▽ More Both simulations and observations have shown that turbulence is a pervasive phenomenon in cosmic scenarios, yet it is particularly difficult to model numerically due to its intrinsically multiscale character which demands high resolutions. Additionally, turbulence is tightly connected to the dynamical state and the formation history of galaxies and galaxy clusters, producing a diverse phenomenlogy which requires large samples of such structures to attain robust conclusions. In this work, we use an adaptive mesh refinement (AMR) cosmological simulation to explore the generation and dissipation of turbulence in galaxy clusters, in connection to its assembly history. We find that major mergers, and more generally accretion of gas, is the main process driving turbulence in the intracluster medium. We have especially focused on solenoidal turbulence, which can be quantified through enstrophy. Our results seem to confirm a scenario for its generation which involves baroclinicity and compression at the external (accretion) and internal (merger) shocks, followed by vortex stretching downstream of them. We have also looked at the infall of mass to the cluster beyond its virial boundary, finding that gas follows trajectories with some degree of helicity, as it has already developed some vorticity in the external shocks. △ Less

Submitted 7 April, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 19 pages, 12 figures; accepted for publication in MNRAS

Journal ref: MNRAS, 504-1 (2021), 510-527

Unravelling cosmic velocity flows: a Helmholtz-Hodge decomposition algorithm for cosmological simulations

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: In the context of intra-cluster medium turbulence, it is essential to be able to split the turbulent velocity field in a compressive and a solenoidal component. We describe and implement a new method for this aim, i.e., performing a Helmholtz-Hodge decomposition, in multi-grid, multi-resolution descriptions, focusing on (but not being restricted to) the outputs of AMR cosmological simulations. The… ▽ More In the context of intra-cluster medium turbulence, it is essential to be able to split the turbulent velocity field in a compressive and a solenoidal component. We describe and implement a new method for this aim, i.e., performing a Helmholtz-Hodge decomposition, in multi-grid, multi-resolution descriptions, focusing on (but not being restricted to) the outputs of AMR cosmological simulations. The method is based on solving elliptic equations for a scalar and a vector potential, from which the compressive and the solenoidal velocity fields, respectively, are derived through differentiation. These equations are addressed using a combination of Fourier (for the base grid) and iterative (for the refinement grids) methods. We present several idealised tests for our implementation, reporting typical median errors in the order of $1\unicode{x2030}$-$1\%$, and with 95-percentile errors below a few percents. Additionally, we also apply the code to the outcomes of a cosmological simulation, achieving similar accuracy at all resolutions, even in the case of highly non-linear velocity fields. We finally take a closer look to the decomposition of the velocity field around a massive galaxy cluster. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 12 pages, 8 figures; accepted for publication in Computer Physics Communications

Journal ref: Computer Physics Communications, Volume 263, June 2021, 107892

On the accretion history of galaxy clusters: temporal and spatial distribution

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: We analyse the results of an Eulerian AMR cosmological simulation in order to quantify the mass growth of galaxy clusters, exploring the differences between dark matter and baryons. We have determined the mass assembly histories (MAHs) of each of the mass components and computed several proxies for the instantaneous mass accretion rate (MAR). The mass growth of both components is clearly dominated… ▽ More We analyse the results of an Eulerian AMR cosmological simulation in order to quantify the mass growth of galaxy clusters, exploring the differences between dark matter and baryons. We have determined the mass assembly histories (MAHs) of each of the mass components and computed several proxies for the instantaneous mass accretion rate (MAR). The mass growth of both components is clearly dominated by the contribution of major mergers, but high MARs can also occur during smooth accretion periods. We explored the correlations between MARs, merger events and clusters' environments, finding the mean densities in $1 \leq r/R_{200m} \leq 1.5$ to correlate strongly with $Γ_{200m}$ in massive clusters which undergo major mergers through their MAH. From the study of the dark matter velocity profiles, we find a strong anticorrelation between the MAR proxies $Γ_{200m}$ and $α_{200m}$. Last, we present a novel approach to study the angularly-resolved distribution of gas accretion flows in simulations, which allows to extract and interpret the main contributions to the accretion picture and to assess systematic differences between the thermodynamical properties of each of these contributions using multipolar analysis. We have preliminarily applied the method to the best numerically-resolved cluster in our simulation. Amongst the most remarkable results, we find that the gas infalling through the cosmic filaments has systematically lower entropy compared to the isotropic component, but we do not find a clear distinction in temperature. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: accepted for publication in MNRAS

Journal ref: MNRAS, 499-2 (2020), 2303-2318

Cosmic magnetic fields with MASCLET: an application to galaxy clusters

Authors: Vicent Quilis, José Ma. Martí, Susana Planelles

Abstract: 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… ▽ More 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 μ$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. △ Less

Submitted 27 March, 2020; originally announced March 2020.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

Long-term FRII jet evolution: Clues from three-dimensional simulations

Authors: Manel Perucho, José-María Martí, Vicent Quilis

Abstract: 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 r… ▽ More 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. △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: (1) Universitat de València. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Multi-wavelength mock observations of the WHIM in a simulated galaxy cluster

Authors: Susana Planelles, Petar Mimica, Vicent Quilis, Carlos Cuesta-Martínez

Abstract: 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… ▽ More 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. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: Accepted for publication in MNRAS

Radio mode feedback: Does relativity matter?

Authors: Manel Perucho, José María Martí, Vicent Quilis, Marina Borja-Lloret

Abstract: 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… ▽ More 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. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: Accepted for publication in MNRAS Letters

Is ram-pressure stripping an efficient mechanism to remove gas in galaxies?

Authors: Vicent Quilis, Susana Planelles, Elena Ricciardelli

Abstract: 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… ▽ More 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. △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

Relic galaxies: where are they?

Authors: Luis Peralta de Arriba, Vicent Quilis, Ignacio Trujillo, María Cebrián, Marc Balcells

Abstract: 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… ▽ More 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. △ Less

Submitted 21 April, 2017; v1 submitted 8 December, 2016; originally announced December 2016.

Comments: 6 pages, 4 figures. Updated to match the published version in Highlights on Spanish Astrophysics IX, Proceedings of the XII Scientific Meeting of the Spanish Astronomical Society held on July 18-22, 2016, in Bilbao, Spain. S. Arribas, A. Alonso-Herrero, F. Figueras, C. Hernández-Monteagudo, A. Sánchez-Lavega, S. Pérez-Hoyos (eds.)

Journal ref: Highlights on Spanish Astrophysics IX (2017), 187-192

Massive relic galaxies prefer dense environments

Authors: Luis Peralta de Arriba, Vicent Quilis, Ignacio Trujillo, María Cebrián, Marc Balcells

Abstract: 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 ma… ▽ More 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. △ Less

Submitted 8 December, 2016; v1 submitted 20 May, 2016; originally announced May 2016.

Comments: 9 pages, 6 figures. Updated to match the published version in MNRAS

Report number: IAC preprint PP 016112

Journal ref: Monthly Notices of the Royal Astronomical Society 461 (2016), 156-163

The Spanish Square Kilometre Array White Book

Authors: J. A. Acosta-Pulido, I. Agudo, A. Alberdi, J. Alcolea, E. J. Alfaro, A. Alonso-Herrero, G. Anglada, P. Arnalte-Mur, Y. Ascasibar, B. Ascaso, R. Azulay, R. Bachiller, A. Baez-Rubio, E. Battaner, J. Blasco, C. B. Brook, V. Bujarrabal, G. Busquet, M. D. Caballero-Garcia, C. Carrasco-Gonzalez, J. Casares, A. J. Castro-Tirado, L. Colina, F. Colomer, I. de Gregorio-Monsalvo , et al. (94 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 17 June, 2015; v1 submitted 10 June, 2015; originally announced June 2015.

Comments: Editors: Miguel A. Perez-Torres (Editor-in-chief), L. Verdes-Montenegro, J. C. Guirado, A. Alberdi, J. Martin-Pintado, R. Bachiller, D. Herranz, J. M. Girart, J. Gorgas, C. Hernandez-Monteagudo, S. Migliari and J. M. Rodriguez Espinosa; 289 pages, published in June 2015 by Sociedad Española de Astronomía

The star formation activity in cosmic voids

Authors: Elena Ricciardelli, Antonio Cava, Jesus Varela, Vicent Quilis

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2014; originally announced October 2014.

Comments: 11 pages, 10 figures, 1 table. Accepted for publication in MNRAS

Large-scale jets from active galactic nuclei as a source of ICM heating: cavities and shocks

Authors: Manel Perucho, José María Martí, Vicent Quilis, Elena Ricciardelli

Abstract: 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 evolut… ▽ More 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. △ Less

Submitted 11 September, 2014; originally announced September 2014.

Comments: Accepted for publication in MNRAS. The quality of some figures has been reduced to adjust their size to arxiv requirements

J-PAS: The Javalambre-Physics of the Accelerated Universe Astrophysical Survey

Authors: N. Benitez, R. Dupke, M. Moles, L. Sodre, J. Cenarro, A. Marin-Franch, K. Taylor, D. Cristobal, A. Fernandez-Soto, C. Mendes de Oliveira, J. Cepa-Nogue, L. R. Abramo, J. S. Alcaniz, R. Overzier, C. Hernandez-Monteagudo, E. J. Alfaro, A. Kanaan, J. M. Carvano, R. R. R. Reis, E. Martinez Gonzalez, B. Ascaso, F. Ballesteros, H. S. Xavier, J. Varela, A. Ederoclite , et al. (127 additional authors not shown)

Abstract: 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$… ▽ More 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 $σ_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Å$ width, placed $100Å$ 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. △ Less

Submitted 20 March, 2014; originally announced March 2014.

Comments: 215 pages, 70 figures. J-PAS "Red Book", containing a description of the technical aspects of the Survey and its main scientific goals

On the universality of void density profiles

Authors: E. Ricciardelli, V. Quilis, J. Varela

Abstract: 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 simulat… ▽ More 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. △ Less

Submitted 12 February, 2014; originally announced February 2014.

Comments: 10 pages, 8 figures. Accepted for publication in MNRAS

Fully cosmological virtual massive galaxies at z=0: kinematical, morphological, and stellar population characterisation

Authors: Javier Navarro-González, Elena Ricciardelli, Vicent Quilis, Alexandre Vazdekis

Abstract: We present the results of a numerical adaptive mesh refinement hydrodynamical and N-body simulation in a $Λ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 o… ▽ More We present the results of a numerical adaptive mesh refinement hydrodynamical and N-body simulation in a $Λ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. △ Less

Submitted 27 September, 2013; originally announced September 2013.

Comments: 20 pages, 12 colour figures, accepted for publication in MNRAS

Expected number of massive galaxy relics in the present-day Universe

Authors: Vicent Quilis, Ignacio Trujillo

Abstract: 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 t… ▽ More 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. △ Less

Submitted 11 July, 2013; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: Accepted for publication in ApJ Letters

The structure of cosmic voids in a LCDM Universe

Authors: E. Ricciardelli, V. Quilis, S. Planelles

Abstract: 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… ▽ More 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. △ Less

Submitted 12 June, 2013; originally announced June 2013.

Comments: 14 pages, 13 figures. Accepted for publication in MNRAS

Structure Finding in Cosmological Simulations: The State of Affairs

Authors: Alexander Knebe, Frazer R. Pearce, Hanni Lux, Yago Ascasibar, Peter Behroozi, Javier Casado, Christine Corbett Moran, Juerg Diemand, Klaus Dolag, Rosa Dominguez-Tenreiro, Pascal Elahi, Bridget Falck, Stefan Gottloeber, Jiaxin Han, Anatoly Klypin, Zarija Lukic, Michal Maciejewski, Cameron K. McBride, Manuel E. Merchan, Stuart I. Muldrew, Mark Neyrinck, Julian Onions, Susana Planelles, Doug Potter, Vicent Quilis , et al. (10 additional authors not shown)

Abstract: 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 remain… ▽ More 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. △ Less

Submitted 26 July, 2013; v1 submitted 2 April, 2013; originally announced April 2013.

Comments: 28 pages containing 13 figures & 4 tables + 9 pages appendix containing another 4 tables + 4 pages of references, accepted for publication in MNRAS

Cosmological shock waves: clues to the formation history of haloes

Authors: Susana Planelles, Vicent Quilis

Abstract: 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 $Λ$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… ▽ More 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 $Λ$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) △ Less

Submitted 4 October, 2012; originally announced October 2012.

Comments: 14 pages, 9 figures, accepted for publication in MNRAS

Satellites around massive galaxies since z$\sim$2: confronting the Millennium simulation with observations

Authors: Vicent Quilis, Ignacio Trujillo

Abstract: 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 simulat… ▽ More 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. △ Less

Submitted 16 May, 2012; originally announced May 2012.

Comments: Accepted for publication in ApJ Letters

Intracluster Medium reheating by relativistic jets

Authors: Manel Perucho, Vicent Quilis, José María Martí

Abstract: 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… ▽ More 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. △ Less

Submitted 13 October, 2011; v1 submitted 27 July, 2011; originally announced July 2011.

Comments: Accepted for publication in ApJ

Neutrino Halos in Clusters of Galaxies and their Weak Lensing Signature

Authors: Francisco Villaescusa-Navarro, Jordi Miralda-Escudé, Carlos Peña-Garay, Vicent Quilis

Abstract: 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, ma… ▽ More 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. △ Less

Submitted 25 April, 2011; originally announced April 2011.

Comments: 13 pages, 11 figures. Submitted to JCAP

Journal ref: JCAP06(2011)027

Haloes gone MAD: The Halo-Finder Comparison Project

Authors: Alexander Knebe, Steffen R. Knollmann, Stuart I. Muldrew, Frazer R. Pearce, Miguel Angel Aragon-Calvo, Yago Ascasibar, Peter S. Behroozi, Daniel Ceverino, Stephane Colombi, Juerg Diemand, Klaus Dolag, Bridget L. Falck, Patricia Fasel, Jeff Gardner, Stefan Gottloeber, Chung-Hsing Hsu, Francesca Iannuzzi, Anatoly Klypin, Zarija Lukic, Michal Maciejewski, Cameron McBride, Mark C. Neyrinck, Susana Planelles, Doug Potter, Vicent Quilis , et al. (12 additional authors not shown)

Abstract: [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 develope… ▽ More [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). △ Less

Submitted 5 April, 2011; originally announced April 2011.

Comments: 27 interesting pages, 20 beautiful figures, and 4 informative tables accepted for publication in MNRAS. The high-resolution version of the paper as well as all the test cases and analysis can be found at the web site http://popia.ft.uam.es/HaloesGoingMAD

ASOHF: a new adaptive spherical overdensity halo finder

Authors: Susana Planelles, Vicent Quilis

Abstract: 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 sta… ▽ More 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. △ Less

Submitted 16 June, 2010; originally announced June 2010.

Comments: Latex format, 16 pages, 11 postscript figures, accepted for publication in Astronomy and Astrophysics

Galaxy cluster mergers

Authors: Susana Planelles, Vicent Quilis

Abstract: We present the results of an Eulerian adaptive mesh refinement (AMR) hydrodynamical and N-body simulation in a $Λ$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.… ▽ More We present the results of an Eulerian adaptive mesh refinement (AMR) hydrodynamical and N-body simulation in a $Λ$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. △ Less

Submitted 22 June, 2009; originally announced June 2009.

Comments: 17 pages, 11 figures, accepted for publication in MNRAS

Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo

Authors: Joachim Stadel, Doug Potter, Ben Moore, Jürg Diemand, Piero Madau, Marcel Zemp, Michael Kuhlen, Vicent Quilis

Abstract: 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 l… ▽ More 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. △ Less

Submitted 22 August, 2008; v1 submitted 21 August, 2008; originally announced August 2008.

Comments: 6 pages, 4 figures, submitted to MNRAS Letters, for full sized images, see http://www.itp.uzh.ch/news.html

Gravitational waves from galaxy encounters

Authors: V. Quilis, A. C. González-García, D. Sáez, J. A. Font

Abstract: 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 grav… ▽ More 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. △ Less

Submitted 23 April, 2007; originally announced April 2007.

Comments: 9 pages with revtex style, 3 ps figures; to be published in Physical Review D

Journal ref: Phys.Rev.D75:104008,2007; Erratum-ibid.D75:109901,2007

AMR simulations of the low T/|W| bar-mode instability of neutron stars

Authors: Pablo Cerdá-Durán, Vicent Quilis, José A. Font

Abstract: 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 supernova… ▽ More 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. △ Less

Submitted 3 April, 2007; originally announced April 2007.

Comments: 20 pages, 9 figures, accepted for publication in Computer Physics Communications.

Journal ref: Comput.Phys.Commun.177:288-297,2007

Fundamental differences between SPH and grid methods

Authors: Oscar Agertz, Ben Moore, Joachim Stadel, Doug Potter, Francesco Miniati, Justin Read, Lucio Mayer, Artur Gawryszczak, Andrey Kravtsov, Joe Monaghan, Ake Nordlund, Frazer Pearce, Vincent Quilis, Douglas Rudd, Volker Springel, James Stone, Elizabeth Tasker, Romain Teyssier, James Wadsley, Rolf Walder

Abstract: 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 instabi… ▽ More 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. △ Less

Submitted 2 October, 2006; originally announced October 2006.

Comments: 15 pages, 13 figures, to be submitted to MNRAS. For high-resolution figures, please see http://www-theorie.physik.unizh.ch/~agertz/

Journal ref: Mon.Not.Roy.Astron.Soc.380:963-978,2007

CGCG 480-022: A distant lonesome merger?

Authors: C. Carretero, A. Vazdekis, C. Gonzalez-Garcia, J. E. Beckman, V. Quilis, ;

Abstract: 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… ▽ More 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. △ Less

Submitted 1 August, 2006; originally announced August 2006.

Comments: Accepted for publication in ApJ Letters

CMB maps lensed by cosmological structures: Simulations and statistical analysis

Authors: L. Antón, P. Cerdá-Durán, V. Quilis, D. Sáez

Abstract: 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… ▽ More 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. △ Less

Submitted 20 April, 2005; originally announced April 2005.

Comments: 28 pages, 11 figures, accepted for publication in ApJ

Journal ref: Astrophys.J. 628 (2005) 1-13

A new multidimensional AMR Hydro+Gravity Cosmological code

Authors: Vicent Quilis

Abstract: 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… ▽ More 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. △ Less

Submitted 20 May, 2004; originally announced May 2004.

Comments: 15 latex pages with mn.sty, 8 eps figures; accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 352 (2004) 1426

Non-Gaussian Signatures in the Lens Deformations of the CMB Sky. A New Ray-Tracing Procedure

Authors: P. Cerdá-Durán, V. Quilis, D. Sáez

Abstract: 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 Gaussian… ▽ More 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. △ Less

Submitted 18 November, 2003; originally announced November 2003.

Comments: 16 pages, 7 figures, to appear in Phys. rev. D

Journal ref: Phys.Rev. D69 (2004) 043002

Momentum transfer across shear flows in Smoothed Particle Hydrodynamic simulations of galaxy formation

Authors: T. Okamoto, A. Jenkins, V. R. Eke, V. Quilis, C. S. Frenk

Abstract: 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 transpor… ▽ More 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) △ Less

Submitted 26 June, 2003; originally announced June 2003.

Comments: 20 pages, 22 figures; accepted for publication in the MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.345:429,2003

CMB anisotropy: deviations from Gaussianity due to non-linear gravity

Authors: A. M. Aliaga, V. Quilis, J. V. Arnau, D. Saez

Abstract: 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 $Λ$CDM cosmology -- are used to simulate the strongly non-linear evolution of cosmological structures. It is proved that these simulations can be comb… ▽ More 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 $Λ$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. △ Less

Submitted 30 October, 2001; originally announced October 2001.

Comments: 6 latex pages with mn.sty, 3 eps figures. Accepted in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 330 (2002) 625

Bubbles, feedback and the intra-cluster medium: three-dimensional hydrodynamic simulations

Authors: Vicent Quilis, Richard G. Bower, Michael L. Balogh

Abstract: 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 in… ▽ More 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. △ Less

Submitted 3 September, 2001; originally announced September 2001.

Comments: 8 Latex pages with MNRAS style, 4 postscript figures. Accepted in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 328 (2001) 1091

Where are the high velocity clouds?

Authors: Vicent Quilis, Ben Moore

Abstract: 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 accretin… ▽ More 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. △ Less

Submitted 25 June, 2001; v1 submitted 14 June, 2001; originally announced June 2001.

Comments: 4 Latex pages with emulateapj.sty, 3 postscript figures. Accepted in ApJLetters. One extra reference added

Journal ref: ApJLetters, 2001, 555, L95-L98

Dynamical effects on galaxies in clusters

Authors: Ben Moore, Vicent Quilis, Richard Bower

Abstract: 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. 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. △ Less

Submitted 15 September, 2000; originally announced September 2000.

Comments: Invited review, XVth IAP meeting

Journal ref: ASP Conf Ser Vol 197, eds Combes, Mamon and Charmandaris, 363

Gone with the wind: the origin of S0 galaxies in clusters

Authors: Vicent Quilis, Ben Moore, Richard Bower

Abstract: 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 morp… ▽ More 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. △ Less

Submitted 1 June, 2000; originally announced June 2000.

Comments: 11 pages, 4 figures. Further simulation details and movies at www.nbody.net

Journal ref: Science 288 (2000) 1617-1620