Galaxy filaments as pearl necklaces

Context. Galaxies in the Universe form chains (filaments) that connect groups and clusters of galaxies. The filamentary network includes nearly half of the galaxies and is visually the most striking feature in cosmological maps. Aims. We study the distribution of galaxies along the filamentary network, trying to find specific patterns and regularities. Methods. Galaxy filaments are defined by the Bisous model, a marked point process with interactions. We use the two-point correlation function and the Rayleigh Z-squared statistic to study how galaxies and galaxy groups are distributed along the filaments. Results. We show that galaxies and groups are not uniformly distributed along filaments, but tend to form a regular pattern. The characteristic length of the pattern is around 7 Mpc/h. A slightly smaller characteristic length 4 Mpc/h can also be found, using the Z-squared statistic. Conclusions. We find that galaxy filaments in the Universe are like pearl necklaces, where the pearls are galaxy groups distributed more or less regularly along the filaments. We propose that this well defined characteristic scale could be used to test various cosmological models and to probe environmental effects on the formation and evolution of galaxies.Comment: 8 pages, 9 figures, 1 table, accepted for publication in A&

Sub-grid modeling of pitch-angle diffusion for ion-scale waves in hybrid-Vlasov simulations with Cartesian velocity space

Numerical simulations have grown to play a central role in modern sciences over the years. The ever-improving technology of supercomputers has made large and precise models available. However, this accuracy is often limited by the cost of computational resources. Lowering the simulation's spatial resolution in order to conserve resources can lead to key processes being unresolved. We have shown in a previous study how insufficient spatial resolution of the proton cyclotron instability leads to a misrepresentation of ion dynamics in hybrid-Vlasov simulations. This leads to larger than expected temperature anisotropy and loss-cone shaped velocity distribution functions. In this study, we present a sub-grid numerical model to introduce pitch-angle diffusion in a 3D Cartesian velocity space, at a spatial resolution where the relevant wave-particle interactions were previously not correctly resolved. We show that the method is successfully able to isotropize loss-cone shaped velocity distribution functions, and that this method could be applied to simulations in order to save computational resources and still correctly model wave-particle interactions.Peer reviewe

Electron Signatures of Reconnection in a Global eVlasiator Simulation

Geospace plasma simulations have progressed toward more realistic descriptions of the solar wind-magnetosphere interaction from magnetohydrodynamic to hybrid ion-kinetic, such as the state-of-the-art Vlasiator model. Despite computational advances, electron scales have been out of reach in a global setting. eVlasiator, a novel Vlasiator submodule, shows for the first time how electromagnetic fields driven by global hybrid-ion kinetics influence electrons, resulting in kinetic signatures. We analyze simulated electron distributions associated with reconnection sites and compare them with Magnetospheric Multiscale (MMS) spacecraft observations. Comparison with MMS shows that key electron features, such as reconnection inflows, heated outflows, flat-top distributions, and bidirectional streaming, are in remarkable agreement. Thus, we show that many reconnection-related features can be reproduced despite strongly truncated electron physics and an ion-scale spatial resolution. Ion-scale dynamics and ion-driven magnetic fields are shown to be significantly responsible for the environment that produces electron dynamics observed by spacecraft in near-Earth plasmas.Peer reviewe

Estimating Inner Magnetospheric Radial Diffusion Using a Hybrid-Vlasov Simulation

Radial diffusion coefficients quantify non-adiabatic transport of energetic particles by electromagnetic field fluctuations in planetary radiation belts. Theoretically, radial diffusion occurs for an ensemble of particles that experience irreversible violation of their third adiabatic invariant, which is equivalent to a change in their Roederer L* parameter. Thus, the Roederer L* coordinate is the fundamental quantity from which radial diffusion coefficients can be computed. In this study, we present a methodology to calculate the Lagrangian derivative of L* from global magnetospheric simulations, and test it with an application to Vlasiator, a hybrid-Vlasov model of near-Earth space. We use a Hamiltonian formalism for particles confined to closed drift shells with conserved first and second adiabatic invariants to compute changes in the guiding center drift paths due to electric and magnetic field fluctuations. We investigate the feasibility of this methodology by computing the time derivative of L* for an equatorial ultrarelativistic electron population travelling along four guiding center drift paths in the outer radiation belt during a 5 minute portion of a Vlasiator simulation. Radial diffusion in this simulation is primarily driven by ultralow frequency waves in the Pc3 range (10-45 s period range) that are generated in the foreshock and transmitted through the magnetopause to the outer radiation belt environment. Our results show that an alternative methodology to compute detailed radial diffusion transport is now available and could form the basis for comparison studies between numerical and observational measurements of radial transport in the Earth's radiation belts.Peer reviewe