Streams Going Notts: The tidal debris finder comparison project

While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations (Knebe et al., 2011, Onions et al., 2012), this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three non-tracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite (Springel et al., 2008) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the {\em unbound tidal features} associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of $\sim20$. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of $\sim40$. Analysing the results of the substructure finders, we find that the general distribution of {\em substructures} differ significantly from the distribution of bound {\em subhaloes}. Most importantly, both bound and unbound {\em substructures} together constitute $\sim18$ of the host halo mass, which is a factor of $\sim2$ higher than the fraction in self-bound {\em subhaloes}. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history.Comment: 19 pages, 12 figures, accepted for publication in MNRA

Galaxies going MAD: The Galaxy-Finder Comparison Project

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.Comment: 14 interesting pages, 17 beautiful figures, and 2 informative tables accepted for publication in MNRAS (matches published version

Sussing merger trees: a proposed merger tree data format

We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON

Major mergers going Notts: challenges for modern halo finders

Merging haloes with similar masses (i.e. major mergers) pose significant challenges for halo finders. We compare five halo-finding algorithms’ (ahf, hbt, rockstar, subfind, and velociraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z > 1). This raises significant uncertainties for theoretical models that require major mergers for, e.g. galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented

Subhaloes gone Notts: the clustering properties of subhaloes

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and study their differences in the density profile, mass fraction and two-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10 per cent level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5 per cent when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering

“PROCESS” : Systems studies of spherical tokamaks

With a reduced aspect ratio, spherical tokamaks have a number of attractive features for a fusion power plant. This can be studied using systems codes, which allow for the rapid conceptual study of power plants covering everything, from the plasma through to electricity generation. In this paper, we describe models in the systems code PROCESS that have been added specifically for spherical tokamaks. Within PROCESS an alternative relation for the plasma current is included which accounts for the increased ratio of Ip/aB. We have tested this against a series of equilibria created with the free boundary equilibrium code FIESTA, and additionally performed our own fit. We also outline the engineering changes that can be made to the device and describe a water-cooled copper centrepost model. To test our models we recreate the published designs for the Fusion Nuclear Science Facility (FNSF) and a High Temperature Superconducting Pilot Plant (HTS-PP) and find good agreement. We conclude by highlighting the efficiencies needed to produce net electricity from small fusion devices

Solving the puzzle of subhalo spins

Investigating the spin parameter distribution of subhalos in two high-resolution isolated halo simulations, recent work by Onions et al. suggested that typical subhalo spins are consistently lower than the spin distribution found for field halos. To further examine this puzzle, we have analyzed simulations of a cosmological volume with sufficient resolution to resolve a significant subhalo population. We confirm the result of Onions et al. and show that the typical spin of a subhalo decreases with decreasing mass and increasing proximity to the host halo center. We interpret this as the growing influence of tidal stripping in removing the outer layers, and hence the higher angular momentum particles, of the subhalos as they move within the host potential. Investigating the redshift dependence of this effect, we find that the typical subhalo spin is smaller with decreasing redshift. This indicates a temporal evolution, as expected in the tidal stripping scenario

SubHaloes going Notts: The SubHalo-Finder Comparison Project

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.Comment: 16 pages, 7 figures, 2 tables, Accepted for MNRA

Subhaloes gone Notts: subhaloes as tracers of the dark matter halo shape

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way-type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields, we find the former to be more triaxial if the analysis is restricted to massive subhaloes. For three of the four analysed haloes, the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness

Performance study of a 3 x 1 x 1 m(3) dual phase liquid Argon Time Projection Chamber exposed to cosmic rays

This work would not have been possible without the support of the Swiss National Science Foundation, Switzerland; CEA and CNRS/IN2P3, France; KEK and the JSPS program, Japan; Ministerio de Ciencia e Innovacion in Spain under grants FPA2016-77347-C2, SEV-2016-0588 and MdM-2015-0509, Comunidad de Madrid, the CERCA program of the Generalitat de Catalunya and the fellowship (LCF/BQ/DI18/11660043) from "La Caixa" Foundation (ID 100010434); the Programme PNCDI III, CERN-RO, under Contract 2/2020, Romania; the U.S. Department of Energy under Grant No. DE-SC0011686. This project has received funding from the European Union's Horizon 2020 Research and Innovation program under Grant Agreement no. 654168. The authors are also grateful to the French government operated by the National Research Agency (ANR) for the LABEX Enigmass, LABEX Lyon Institute of Origins (ANR-10-LABX-0066) of the Universite de Lyon for its financial support within the program "Investissements d'Avenir" (ANR-11-IDEX-0007).We report the results of the analyses of the cosmic ray data collected with a 4 tonne (3x1x1 m(3)) active mass (volume) Liquid Argon Time-Projection Chamber (TPC) operated in a dual-phase mode. We present a detailed study of the TPC's response, its main detector parameters and performance. The results are important for the understanding and further developments of the dual-phase technology, thanks to the verification of key aspects, such as the extraction of electrons from liquid to gas and their amplification through the entire one square metre readout plain, gain stability, purity and charge sharing between readout views.Swiss National Science Foundation (SNSF)French Atomic Energy CommissionCentre National de la Recherche Scientifique (CNRS)High Energy Accelerator Research Organization (KEK)Ministry of Education, Culture, Sports, Science and Technology, Japan (MEXT)Japan Society for the Promotion of ScienceSpanish Government FPA2016-77347-C2 SEV-2016-0588MdM-2015-0509Comunidad de MadridCERCA program of the Generalitat de CatalunyaLa Caixa Foundation LCF/BQ/DI18/11660043 100010434Programme PNCDI III, RomaniaCERN-RO, Romania 2/2020United States Department of Energy (DOE) SC0011686European Commission 654168Universite de Lyon ANR-10-LABX-0066 ANR-11-IDEX-000