• In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.
• Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.
• ### Cosmology and fundamental physics with the Euclid satellite(1206.1225)

Oct. 22, 2015 hep-ph, gr-qc, astro-ph.CO
Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.
• ### Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos(1503.02689)

March 9, 2015 astro-ph.CO, astro-ph.IM
We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $\epsilon$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\rm fs}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\rm fs}$; filament stability; and to obtain new constraints on the temperature of dark matter.
• ### Milking the spherical cow: on aspherical dynamics in spherical coordinates(1502.07356)

Feb. 25, 2015 astro-ph.CO, astro-ph.GA
Galaxies and the dark matter halos that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealised calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are "maximally stable", i.e. that do not evolve at first order when external potentials (which arise from baryons, large scale tidal fields or infalling substructure) are applied. We show that a spherically-symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealised spherically-symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter halos (formed in collisionless cosmological simulations) nearly attain our maximally-stable limit, but their evolution freezes out before reaching it.
• ### Growing galaxies via superbubble-driven accretion flows(1410.3827)

Oct. 14, 2014 astro-ph.GA
We use a suite a cooling halo simulations to study a new mechanism for rapid accretion of hot halo gas onto star-forming galaxies. Correlated supernovae events create converging 'superbubbles' in the halo gas. Where these collide, the density increases, driving cooling filaments of low metallicity gas that feed the disc. At our current numerical resolution (20 pc) we are only able to resolve the most dramatic events; these could be responsible for the build-up of galaxy discs after the most massive gas-rich mergers have completed (z < 1). As we increase the numerical resolution, we find that the filaments persist for longer, driving continued late-time star formation. This suggests that SNe-driven accretion could act as an efficient mechanism for extracting cold gas from the hot halo, driving late-time star formation in disc galaxies. We show that such filament feeding leads to a peak star formation rate (SFR) of $\sim 3$ M$_{\rm sun}$ yr$^{-1}$, consistent with estimates for the Milky Way. By contrast, direct cooling from the hot halo ('hot-mode' accretion, not present in the simulations that show filament feeding) falls short of the SNe-driven SFR by a factor of 3-4, and is sustained over a shorter time period. The filaments we resolve extend to $\sim$ 50 kpc, reaching column densities of $\sim 10^{18}$ cm$^{-2}$. We show that such structures can plausibly explain the broad dispersion in Mg II absorption seen along sight lines to quasars. Our results suggest a dual role for stellar feedback in galaxy formation, suppressing hot-mode accretion while promoting cold-mode accretion along filaments. This ultimately leads to more star formation, suggesting that the positive feedback effect outweighs the negative. Finally, since the filamentary gas has higher angular momentum than that coming from hot-mode accretion, we show that this leads to the formation of substantially larger gas discs.
• ### Thermal instabilities in cooling galactic coronae: fuelling star formation in galactic discs(1207.3814)

Sept. 4, 2012 astro-ph.CO, astro-ph.GA
We investigate the means by which cold gas can accrete onto Milky Way mass galaxies from a hot corona of gas, using a new smoothed particle hydrodynamics code, 'SPHS'. We find that the 'cold clumps' seen in many classic SPH simulations in the literature are not present in our SPHS simulations. Instead, cold gas condenses from the halo along filaments that form at the intersection of supernovae-driven bubbles from previous phases of star formation. This positive feedback feeds cold gas to the galactic disc directly, fuelling further star formation. The resulting galaxies in the SPH and SPHS simulations differ greatly in their morphology, gas phase diagrams, and stellar content. We show that the classic SPH cold clumps owe to a numerical thermal instability caused by an inability for cold gas to mix in the hot halo. The improved treatment of mixing in SPHS suppresses this instability leading to a dramatically different physical outcome. In our highest resolution SPHS simulation, we find that the cold filaments break up into bound clumps that form stars. The filaments are overdense by a factor of 10-100 compared to the surrounding gas, suggesting that the fragmentation results from a physical non-linear instability driven by the overdensity. This 'fragmenting filament' mode of disc growth has important implications for galaxy formation, in particular the role of star formation in bringing cold gas into disc galaxies.
• ### N-body simulations of gravitational dynamics(1105.1082)

We describe the astrophysical and numerical basis of N-body simulations, both of collisional stellar systems (dense star clusters and galactic centres) and collisionless stellar dynamics (galaxies and large-scale structure). We explain and discuss the state-of-the-art algorithms used for these quite different regimes, attempt to give a fair critique, and point out possible directions of future improvement and development. We briefly touch upon the history of N-body simulations and their most important results.
• ### Resolving Globular Cluster Formation within a Cosmological Context(0908.1254)

Oct. 21, 2009 astro-ph.CO
We place constraints on the formation redshifts for blue globular clusters (BGCs), independent of the details of hydrodynamics and population III star formation. The observed radial distribution of BGCs in the Milky Way Galaxy suggests that they formed in biased dark matter halos at high redshift. As a result, simulations of a ~1 Mpc box up to z~10 must resolve BGC formation in LCDM. We find that most halo stars could be produced from destroyed BGCs and other low-mass clusters that formed at high redshift. We present a proof-of-concept simulation that captures the formation of globular-like star clusters.
• ### Dark Matter Disc Enhanced Neutrino Fluxes from the Sun and Earth(0902.4001)

April 21, 2009 hep-ph, astro-ph.GA, astro-ph.HE
As disc galaxies form in a hierarchical cosmology, massive merging satellites are preferentially dragged towards the disc plane. The material accreted from these satellites forms a dark matter disc that contributes 0.25 - 1.5 times the non-rotating halo density at the solar position. Here, we show the importance of the dark disc for indirect dark matter detection in neutrino telescopes. Previous predictions of the neutrino flux from WIMP annihilation in the Earth and the Sun have assumed that Galactic dark matter is spherically distributed with a Gaussian velocity distribution, the standard halo model. Although the dark disc has a local density comparable to the dark halo, its higher phase space density at low velocities greatly enhances capture rates in the Sun and Earth. For typical dark disc properties, the resulting muon flux from the Earth is increased by three orders of magnitude over the SHM, while for the Sun the increase is an order of magnitude. This significantly increases the sensitivity of neutrino telescopes to fix or constrain parameters in WIMP models. The flux from the Earth is extremely sensitive to the detailed properties of the dark disc, while the flux from the Sun is more robust. The enhancement of the muon flux from the dark disc puts the search for WIMP annihilation in the Earth on the same level as the Sun for WIMP masses < 100 GeV.
• ### The Causes of Halo Shape Changes Induced by Cooling Baryons: Disks Versus Substructures(0707.0737)

June 13, 2008 astro-ph
Cold dark matter cosmogony predicts triaxial dark matter halos, whereas observations find quite round halos. This is most likely due to the condensation of baryons leading to rounder halos. We examine the halo phase space distribution basis for such shape changes. Triaxial halos are supported by box orbits, which pass arbitrarily close to the density center. The decrease in triaxiality caused by baryons is thought to be due to the scattering of these orbits. We test this hypothesis with simulations of disks grown inside triaxial halos. After the disks are grown we check whether the phase space structure has changed by evaporating the disks and comparing the initial and final states. While the halos are substantially rounder when the disk is at full mass, their final shape after the disk is evaporated is not much different from the initial. Likewise, the halo becomes (more) radially anisotropic when the disk is grown, but the final anisotropy is consistent with the initial. Only if the baryons are unreasonably compact or massive does the halo change irreversibly. We show that the character of individual orbits is not generally changed by the growing mass. Thus the central condensation of baryons does not destroy enough box orbits to cause the shape change. Rather, box orbits merely become rounder along with the global potential. However, if angular momentum is transferred to the halo, either via satellites or via bars, a large irreversible change in the halo distribution occurs. The ability of satellites to alter the phase space distribution of the halo is of particular concern to galaxy formation simulations since halo triaxiality can profoundly influence the evolution of disks.
• ### Fundamental differences between SPH and grid methods(astro-ph/0610051)

Oct. 2, 2006 astro-ph
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.