ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation

Resolving numerically Vlasov-Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm consisting in representing the phase-space sheet with a conforming, self-adaptive simplicial tessellation of which the vertices follow the Lagrangian equations of motion. The algorithm is implemented both in six- and four-dimensional phase-space. Refinement of the tessellation mesh is performed using the bisection method and a local representation of the phase-space sheet at second order relying on additional tracers created when needed at runtime. In order to preserve in the best way the Hamiltonian nature of the system, refinement is anisotropic and constrained by measurements of local Poincar\'e invariants. Resolution of Poisson equation is performed using the fast Fourier method on a regular rectangular grid, similarly to particle in cells codes. To compute the density projected onto this grid, the intersection of the tessellation and the grid is calculated using the method of Franklin and Kankanhalli (1993) generalised to linear order. As preliminary tests of the code, we study in four dimensional phase-space the evolution of an initially small patch in a chaotic potential and the cosmological collapse of a fluctuation composed of two sinusoidal waves. We also perform a "warm" dark matter simulation in six-dimensional phase-space that we use to check the parallel scaling of the code.Comment: Code and illustration movies available at: http://www.vlasix.org/index.php?n=Main.ColDICE - Article submitted to Journal of Computational Physic

Euler complexes and geometry of modular varieties

There is a mysterious connection between the multiple polylogarithms at N-th roots of unity and modular varieties. In this paper we "explain" it in the simplest case of the double logarithm. We introduce an Euler complex data on modular curves. It includes a length two complex on every modular curve. Their second cohomology groups recover the Beilinson-Kato Euler system in K_2 of modular curves. We show that the above connection in the double logarithm case is provided by the specialization at a cusp of the Euler complex data on the modular curve Y_1(N). Furthermore, specializing the Euler complexes at CM points we find new examples of the connection with geometry of modular varieties, this time hyperbolic 3-folds.Comment: Dedicated to Joseph Bernstein for his 60th birthday. The final version. Some corrections were made. To appear in GAFA, special volume dedicated to J. Bernstei