-
Neural network reconstruction of density and velocity fields from the 2MASS Redshift Survey
Authors:
Robert Lilow,
Punyakoti Ganeshaiah Veena,
Adi Nusser
Abstract:
We reconstruct the 3D matter density and peculiar velocity fields in the local Universe up to a distance of $200\,h^{-1}\,\mathrm{Mpc}$ from the Two-Micron All-Sky Redshift Survey (2MRS), using a neural network (NN). We employ a NN with U-net autoencoder architecture and a weighted mean squared error loss function, trained separately to output either the density or velocity field for a given input…
▽ More
We reconstruct the 3D matter density and peculiar velocity fields in the local Universe up to a distance of $200\,h^{-1}\,\mathrm{Mpc}$ from the Two-Micron All-Sky Redshift Survey (2MRS), using a neural network (NN). We employ a NN with U-net autoencoder architecture and a weighted mean squared error loss function, trained separately to output either the density or velocity field for a given input grid of galaxy number counts. The NN is trained on mocks derived from the Quijote N-body simulations, incorporating redshift-space distortions (RSD), galaxy bias and selection effects, closely mimicking the characteristics of 2MRS. The trained NN is benchmarked against a standard Wiener filter (WF) on a validation set of mocks, before applying it to 2MRS. The NN reconstructions effectively approximate the mean posterior estimate of the true density and velocity fields conditioned on the observations. They consistently outperform the WF in terms of reconstruction accuracy, and effectively capture the nonlinear relation between velocity and density. The NN-reconstructed bulk flow of the total survey volume exhibits a significant correlation with the true mock bulk flow, demonstrating that the NN is sensitive to information on super-survey scales encoded in the RSD. When applied to 2MRS, the NN successfully recovers the main known clusters, some of which are partially in the Zone of Avoidance. The reconstructed bulk flows in spheres of different radii less than $100\,h^{-1}\,\mathrm{Mpc}$ are in good agreement with a previous 2MRS analysis that required an additional external bulk flow component inferred from directly observed peculiar velocities. The NN-reconstructed peculiar velocity of the Local Group closely matches the observed CMB dipole in amplitude and Galactic latitude, and only deviates by $18^\circ$ in longitude. The NN-reconstructed fields are publicly available.
△ Less
Submitted 2 April, 2024;
originally announced April 2024.
-
Back to the present: A general treatment for the tidal field from the wake of dynamical friction
Authors:
Rain Kipper,
Peeter Tenjes,
María Benito,
Punyakoti Ganeshaiah Veena,
Aikaterini Niovi Triantafyllaki,
Indrek Vurm,
Moorits Mihkel Muru,
Maret Einasto,
Elmo Tempel
Abstract:
Dynamical friction can be a valuable tool for inferring dark matter properties that are difficult to constrain by other methods. Most applications of dynamical friction calculations are concerned with the long-term angular momentum loss and orbital decay of the perturber within its host. This, however, assumes knowledge of the unknown initial conditions of the system. We advance an alternative met…
▽ More
Dynamical friction can be a valuable tool for inferring dark matter properties that are difficult to constrain by other methods. Most applications of dynamical friction calculations are concerned with the long-term angular momentum loss and orbital decay of the perturber within its host. This, however, assumes knowledge of the unknown initial conditions of the system. We advance an alternative methodology to infer the host properties from the perturber's shape distortions induced by the tides of the wake of dynamical friction, which we refer to as the tidal dynamical friction. As the shape distortions rely on the tidal field that has a predominantly local origin, we present a strategy to find the local wake by integrating the stellar orbits back in time along with the perturber, then removing the perturber's potential and re-integrating them back to the present. This provides perturbed and unperturbed coordinates and hence a change in coordinates, density, and acceleration fields, which yields the back-reaction experienced by the perturber. The method successfully recovers the tidal field of the wake based on a comparison with N-body simulations. We show that similar to the tidal field itself, the noise and randomness of the dynamical friction force due to the finite number of stars is also dominated by regions close to the perturber. Stars near the perturber influence it more but are smaller in number, causing a high variance in the acceleration field. These fluctuations are intrinsic to dynamical friction. We show that a stellar density of $0.0014 {\rm M_\odot\, kpc^{-3}}$ yields an inherent variance of 10% to the dynamical friction. The current method extends the family of dynamical friction methods that allow for the inference of host properties from tidal forces of the wake. It can be applied to specific galaxies, such as Magellanic Clouds, with Gaia data.
△ Less
Submitted 15 November, 2023; v1 submitted 7 November, 2023;
originally announced November 2023.
-
Large-scale density and velocity field reconstructions with neural networks
Authors:
Punyakoti Ganeshaiah Veena,
Robert Lilow,
Adi Nusser
Abstract:
We assess a neural network (NN) method for reconstructing 3D cosmological density and velocity fields (target) from discrete and incomplete galaxy distributions (input). We employ second-order Lagrangian Perturbation Theory to generate a large ensemble of mock data to train an autoencoder (AE) architecture with a Mean Squared Error (MSE) loss function. The AE successfully captures nonlinear featur…
▽ More
We assess a neural network (NN) method for reconstructing 3D cosmological density and velocity fields (target) from discrete and incomplete galaxy distributions (input). We employ second-order Lagrangian Perturbation Theory to generate a large ensemble of mock data to train an autoencoder (AE) architecture with a Mean Squared Error (MSE) loss function. The AE successfully captures nonlinear features arising from gravitational dynamics and the discreteness of the galaxy distribution. It preserves the positivity of the reconstructed density field and exhibits a weaker suppression of the power on small scales than the traditional linear Wiener filter (WF), which we use as a benchmark. In the density reconstruction, the reduction of the AE MSE relative to the WF is $\sim 15 \%$, whereas, for the velocity reconstruction, a relative reduction of up to a factor of two can be achieved. The AE is advantageous to the WF at recovering the distribution of the target fields, especially at the tails. In fact, trained with an MSE loss, any NN estimate approaches the unbiased mean of the underlying target given the input. This implies a slope of unity in the linear regression of the true on the NN-reconstructed field. Only for the special case of Gaussian fields, the NN and WF estimates are equivalent. Nonetheless, we also recover a linear regression slope of unity for the WF with non-Gaussian fields.
△ Less
Submitted 1 June, 2023; v1 submitted 13 December, 2022;
originally announced December 2022.
-
The role of stochastic and smooth processes in regulating galaxy quenching
Authors:
Rain Kipper,
Antti Tamm,
Elmo Tempel,
Roberto de Propris,
Punyakoti Ganeshaiah Veena
Abstract:
Galaxies can be classified as passive ellipticals or star-forming discs. Ellipticals dominate at the high end of the mass range, and therefore there must be a mechanism responsible for the quenching of star-forming galaxies. This could either be due to the secular processes linked to the mass and star formation of galaxies or to external processes linked to the surrounding environment. In this pap…
▽ More
Galaxies can be classified as passive ellipticals or star-forming discs. Ellipticals dominate at the high end of the mass range, and therefore there must be a mechanism responsible for the quenching of star-forming galaxies. This could either be due to the secular processes linked to the mass and star formation of galaxies or to external processes linked to the surrounding environment. In this paper, we analytically model the processes that govern galaxy evolution and quantify their contribution. We have specifically studied the effects of mass quenching, gas stripping, and mergers on galaxy quenching. To achieve this, we first assumed a set of differential equations that describe the processes that shape galaxy evolution. We then modelled the parameters of these equations by maximising likelihood. These equations describe the evolution of galaxies individually, but the parameters of the equations are constrained by matching the extrapolated intermediate-redshift galaxies with the low-redshift galaxy population. In this study, we modelled the processes that change star formation and stellar mass in massive galaxies from the GAMA survey between z~0.4 and the present. We identified and quantified the contributions from mass quenching, gas stripping, and mergers to galaxy quenching. The quenching timescale is on average 1.2 Gyr and a closer look reveals support for the slow-then-rapid quenching scenario. The major merging rate of galaxies is about once per 10~Gyr, while the rate of ram pressure stripping is significantly higher. In galaxies with decreasing star formation, we show that star formation is lost to fast quenching mechanisms such as ram pressure stripping and is countered by mergers, at a rate of about 41% Gyr$^{-1}$ and to mass quenching 49% Gyr$^{-1}$. (abridged)
△ Less
Submitted 21 January, 2021;
originally announced January 2021.
-
An EAGLE view of the missing baryons
Authors:
Toni Tuominen,
Jukka Nevalainen,
Elmo Tempel,
Teet Kuutma,
Nastasha Wijers,
Joop Schaye,
Pekka Heinämäki,
Massimiliano Bonamente,
Punyakoti Ganeshaiah Veena
Abstract:
Context. A significant fraction of the predicted baryons remains undetected in the local universe. We adopted the common assumption that a large fraction of the missing baryons corresponds to the hot (log T(K) = 5.5-7) phase of the Warm Hot Intergalactic Medium (WHIM). We base our missing baryons search on the scenario whereby the WHIM has been heated up via accretion shocks and galactic outflows,…
▽ More
Context. A significant fraction of the predicted baryons remains undetected in the local universe. We adopted the common assumption that a large fraction of the missing baryons corresponds to the hot (log T(K) = 5.5-7) phase of the Warm Hot Intergalactic Medium (WHIM). We base our missing baryons search on the scenario whereby the WHIM has been heated up via accretion shocks and galactic outflows, and is concentrated towards the filaments of the Cosmic Web.
Aims. Our aim is to improve the observational search of the poorly detected hot WHIM.
Methods. We detect the filamentary structure within the EAGLE simulation by applying the Bisous formalism to the galaxy distribution. In addition, we use the MMF/NEXUS+ classification of the large scale environment of the dark matter component in EAGLE. We then study the spatio-thermal distribution of the hot baryons within the extracted filaments.
Results. While the filaments occupy only 5% of the full simulation volume, the diffuse hot intergalactic medium in filaments amounts to 23% $-$ 25% of the total baryon budget, or 79% $-$ 87% of all the hot WHIM. The most optimal filament sample, with a missing baryon mass fraction of 82%, is obtained by selecting Bisous filaments with a high galaxy luminosity density. For these filaments we derived analytic formulae for the radial gas density and temperature profiles, consistent with recent Planck SZ and CMB lensing observations within the central $r$~ 1 Mpc.
Conclusions. Results from EAGLE suggest that the missing baryons are strongly concentrated towards the filament axes. Since the filament finding methods used here are applicable to galaxy surveys, a large fraction of the missing baryons can be localised by focusing the observational efforts on the central 1 Mpc regions of the filaments. Moreover, focusing on high galaxy luminosity density regions will optimise the observational signal.
△ Less
Submitted 16 December, 2020;
originally announced December 2020.
-
Cosmic Ballet III: halo spin evolution in the cosmic web
Authors:
Punyakoti Ganeshaiah Veena,
Marius Cautun,
Rien van de Weygaert,
Elmo Tempel,
Carlos S. Frenk
Abstract:
We explore the evolution of halo spins in the cosmic web using a very large sample of dark matter haloes in the $Λ$CDM Planck-Millennium N-body simulation. We use the NEXUS+ multiscale formalism to identify the hierarchy of filaments and sheets of the cosmic web at several redshifts. We find that at all times the magnitude of halo spins correlates with the web environment, being largest in filamen…
▽ More
We explore the evolution of halo spins in the cosmic web using a very large sample of dark matter haloes in the $Λ$CDM Planck-Millennium N-body simulation. We use the NEXUS+ multiscale formalism to identify the hierarchy of filaments and sheets of the cosmic web at several redshifts. We find that at all times the magnitude of halo spins correlates with the web environment, being largest in filaments, and, for the first time, we show that it also correlates with filament thickness as well as the angle between spin-orientation and the spine of the host filament. For example, massive haloes in thick filaments spin faster than their counterparts in thin filaments, while for low-mass haloes the reverse is true. We also have studied the evolution of alignment between halo spin orientations and the preferential axes of filaments and sheets. The alignment varies with halo mass, with the spins of low-mass haloes being predominantly along the filament spine, while those of high-mass haloes being predominantly perpendicular to the filament spine. On average, for all halo masses, halo spins become more perpendicular to the filament spine at later times. At all redshifts, the spin alignment shows a considerable variation with filament thickness, with the halo mass corresponding to the transition from parallel to perpendicular alignment varying by more than one order of magnitude. The environmental dependence of halo spin magnitude shows little evolution for $z\leq2$ and is likely a consequence of the correlations in the initial conditions or high redshift effects
△ Less
Submitted 18 February, 2021; v1 submitted 20 July, 2020;
originally announced July 2020.
-
Quantifying torque from the Milky Way bar using Gaia DR2
Authors:
Rain Kipper,
Peeter Tenjes,
Taavi Tuvikene,
Punyakoti Ganeshaiah Veena,
Elmo Tempel
Abstract:
We determine the mass of the Milky Way bar and the torque it causes, using Gaia DR2, by applying the orbital arc method. Based on this, we have found that the gravitational acceleration is not directed towards the centre of our Galaxy but a few degrees away from it. We propose that the tangential acceleration component is caused by the bar of the Galaxy. Calculations based on our model suggest tha…
▽ More
We determine the mass of the Milky Way bar and the torque it causes, using Gaia DR2, by applying the orbital arc method. Based on this, we have found that the gravitational acceleration is not directed towards the centre of our Galaxy but a few degrees away from it. We propose that the tangential acceleration component is caused by the bar of the Galaxy. Calculations based on our model suggest that the torque experienced by the region around the Sun is $\approx 2400\, km^2 s^{-2}$ per solar mass. The mass estimate for the bar is $\sim 1.6\pm0.3\times10^{10} M_\odot$. Using greatly improved data from Gaia DR2, we have computed the acceleration field to great accuracy by adapting the oPDF method (Han et al. 2016) locally and used the phase space coordinates of $\sim 4\times10^5$ stars within a distance of 0.5 kpc from the Sun. In the orbital arc method, the first step is to guess an acceleration field and then reconstruct the stellar orbits using this acceleration for all the stars within a specified region. Next, the stars are redistributed along orbits to check if the overall phase space distribution has changed. We repeat this process until we find an acceleration field that results in a new phase space distribution that is the same as the one that we started with; we have then recovered the true underlying acceleration.
△ Less
Submitted 23 April, 2020;
originally announced April 2020.
-
The Cosmic Ballet II: Spin alignment of galaxies and haloes with large-scale filaments in the EAGLE simulation
Authors:
Punyakoti Ganeshaiah Veena,
Marius Cautun,
Elmo Tempel,
Rien van de Weygaert,
Carlos S. Frenk
Abstract:
We investigate the alignment of galaxies and haloes relative to cosmic web filaments using the EAGLE hydrodynamical simulation. We identify filaments by applying the NEXUS+ method to the mass distribution and the Bisous formalism to the galaxy distribution. Both web finders return similar filamentary structures that are well aligned and that contain comparable galaxy populations. EAGLE haloes have…
▽ More
We investigate the alignment of galaxies and haloes relative to cosmic web filaments using the EAGLE hydrodynamical simulation. We identify filaments by applying the NEXUS+ method to the mass distribution and the Bisous formalism to the galaxy distribution. Both web finders return similar filamentary structures that are well aligned and that contain comparable galaxy populations. EAGLE haloes have an identical spin alignment with filaments as their counterparts in dark matter only simulations: a complex mass dependent trend with low mass haloes spinning preferentially parallel to and high mass haloes spinning preferentially perpendicular to filaments. In contrast, galaxy spins do not show such a spin transition and have a propensity for perpendicular alignments at all masses, with the degree of alignment being largest for massive galaxies. This result is valid for both NEXUS+ and Bisous filaments. When splitting by morphology, we find that elliptical galaxies show a stronger orthogonal spin--filament alignment than spiral galaxies of similar mass. The same is true of their haloes, with the host haloes of elliptical galaxies having a larger degree of orthogonal alignment than the host haloes of spirals. Due to the misalignment between galaxy shape and spin, galaxy minor axes are oriented differently with filaments than galaxy spins. We find that the galaxies whose minor axis is perpendicular to a filament are much better aligned with their host haloes. This suggests that many of the same physical processes determine both the galaxy--filament and the galaxy--halo alignments.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
The Cosmic Ballet: spin and shape alignments of haloes in the cosmic web
Authors:
Punyakoti Ganeshaiah Veena,
Marius Cautun,
Rien van de Weygaert,
Elmo Tempel,
Bernard J. T. Jones,
Steven Rieder,
Carlos S. Frenk
Abstract:
We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $Λ$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, wi…
▽ More
We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $Λ$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, with both the major axis of haloes and their peculiar velocity tending to orient along the filament. However, the spin - filament alignment displays a more complex trend changing from preferentially parallel at low masses to preferentially perpendicular at high masses. This "spin flip" occurs at an average mass of $5\times10^{11}~h^{-1}M_\odot$. This mass increases with increasing filament diameter, varying by more than an order of magnitude between the thinnest and thickest filament samples. We also find that the inner parts of haloes have a spin flip mass that is several times smaller than that of the halo as a whole. These results confirm that recent accretion is responsible for the complex behaviour of the halo spin - filament alignment. Low-mass haloes mainly accrete mass along directions perpendicular to their host filament and thus their spins tend to be oriented along the filaments. In contrast, high-mass haloes mainly accrete along their host filaments and have their spins preferentially perpendicular to them. Furthermore, haloes located in thinner filaments are more likely to accrete along their host filaments than haloes of the same mass located in thicker filaments.
△ Less
Submitted 19 September, 2018; v1 submitted 30 April, 2018;
originally announced May 2018.
