https://biblio.ugent.be/publication/685594/file/1130605.pdf

Review of Particle Physics

We use high-resolution N-body simulations to study the equilibrium density profiles of dark matter halos in hierarchically clustering universes. We find that all such profiles have the same shape, independent of the halo mass, the initial density fluctuation spectrum, and the values of the cosmological parameters. Spherically averaged equilibrium profiles are well fitted over two decades in radius by a simple formula originally proposed to describe the structure of galaxy clusters in a cold dark matter universe. In any particular cosmology, the two scale parameters of the fit, the halo mass and its characteristic density, are strongly correlated. Low-mass halos are significantly denser than more massive systems, a correlation that reflects the higher collapse redshift of small halos. The characteristic density of an equilibrium halo is proportional to the density of the universe at the time it was assembled. A suitable definition of this assembly time allows the same proportionality constant to be used for all the cosmologies that we have tested. We compare our results with previous work on halo density profiles and show that there is good agreement. We also provide a step-by-step analytic procedure, based on the Press-Schechter formalism, that allows accurate equilibrium profiles to be calculated as a function of mass in any hierarchical model.

https://iopscience.iop.org/article/10.1086/304888/pdf

A Universal Density Profile from Hierarchical Clustering

We discuss the cosmological simulation code GADGET-2, a new massively parallel TreeSPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics (SPH). Our implementation of SPH manifestly conserves energy and entropy in regions free of dissipation, while allowing for fully adaptive smoothing lengths. Gravitational forces are computed with a hierarchical multipole expansion, which can optionally be applied in the form of a TreePM algorithm, where only short-range forces are computed with the ‘tree’ method while long-range forces are determined with Fourier techniques. Time integration is based on a quasi-symplectic scheme where long-range and short-range forces can be integrated with different time-steps. Individual and adaptive short-range time-steps may also be employed. The domain decomposition used in the parallelization algorithm is based on a space-filling curve, resulting in high flexibility and tree force errors that do not depend on the way the domains are cut. The code is efficient in terms of memory consumption and required communication bandwidth. It has been used to compute the first cosmological N-body simulation with more than 10 10 dark matter particles, reaching a homogeneous spatial dynamic range of 10 5 per dimension in a three-dimensional box. It has also been used to carry out very large cosmological SPH simulations that account for radiative cooling and star formation, reaching total particle numbers of more than 250 million. We present the algorithms used by the code and discuss their accuracy and performance using a number of test problems. GADGET-2 is publicly released to the research community. Ke yw ords: methods: numerical ‐ galaxies: interactions ‐ dark matter.

/pdf/the-cosmological-simulation-code-gadget-2-48ba4eyb8k.pdf

The Cosmological simulation code GADGET-2

This paper describes the Seventh Data Release of the Sloan Digital Sky Survey (SDSS), marking the completion of the original goals of the SDSS and the end of the phase known as SDSS-II. It includes 11,663 deg^2 of imaging data, with most of the ~2000 deg^2 increment over the previous data release lying in regions of low Galactic latitude. The catalog contains five-band photometry for 357 million distinct objects. The survey also includes repeat photometry on a 120° long, 2°.5 wide stripe along the celestial equator in the Southern Galactic Cap, with some regions covered by as many as 90 individual imaging runs. We include a co-addition of the best of these data, going roughly 2 mag fainter than the main survey over 250 deg^2. The survey has completed spectroscopy over 9380 deg^2; the spectroscopy is now complete over a large contiguous area of the Northern Galactic Cap, closing the gap that was present in previous data releases. There are over 1.6 million spectra in total, including 930,000 galaxies, 120,000 quasars, and 460,000 stars. The data release includes improved stellar photometry at low Galactic latitude. The astrometry has all been recalibrated with the second version of the USNO CCD Astrograph Catalog, reducing the rms statistical errors at the bright end to 45 milliarcseconds per coordinate. We further quantify a systematic error in bright galaxy photometry due to poor sky determination; this problem is less severe than previously reported for the majority of galaxies. Finally, we describe a series of improvements to the spectroscopic reductions, including better flat fielding and improved wavelength calibration at the blue end, better processing of objects with extremely strong narrow emission lines, and an improved determination of stellar metallicities.

/pdf/the-seventh-data-release-of-the-sloan-digital-sky-survey-5cm9fqkek5.pdf

The Seventh Data Release of the Sloan Digital Sky Survey

The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,1603 particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.

/pdf/simulations-of-the-formation-evolution-and-clustering-of-1yahvmv5qk.pdf

Simulations of the formation, evolution and clustering of galaxies and quasars

We investigate the survival of planetesimal discs over Gyr time-scales, using a unified approach that is applicable to all Keplerian discs of solid bodies – dust grains, asteroids, planets, etc. Planetesimal discs can be characterized locally by four parameters: surface density, semimajor axis, planetesimal size and planetesimal radial velocity dispersion. Any planetesimal disc must have survived all dynamical processes, including gravitational instability, dynamical chaos, gravitational scattering, physical collisions, and radiation forces, that would lead to significant evolution over its lifetime. These processes lead to a rich set of constraints that strongly restrict the possible properties of long-lived discs. Within this framework, we also discuss the detection of planetesimal discs using radial velocity measurements, transits, microlensing and the infrared emission from the planetesimals themselves or from dust generated by planetesimal collisions.

/pdf/long-lived-planetesimal-discs-52wn8y8ebp.pdf

Long-lived planetesimal discs

Fitting Selected Random Planetary Systems to Titius–Bode Laws

Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $\lambda \sim 1 \,\mathrm{kpc}$ in typical galactic potentials Such `fuzzy' dark matter (FDM) behaves like cold dark matter (CDM) on much larger scales than the de Broglie wavelength, but may resolve some of the challenges faced by CDM in explaining the properties of galaxies on small scales ($\lesssim 10\,\mathrm{kpc}$) Because of its wave nature, FDM exhibits stochastic density fluctuations on the scale of the de Broglie wavelength that never damp The gravitational field from these fluctuations scatters stars and black holes, causing their orbits to diffuse through phase space We show that this relaxation process can be analyzed quantitatively with the same tools used to analyze classical two-body relaxation in an $N$-body system, and can be described by treating the FDM fluctuations as quasiparticles, with effective mass $\sim 10^7 M_\odot {(1\,\mathrm{kpc}/r)}^2{(10^{-22}\,\mathrm{eV}/m_b)}^3$ in a galaxy with a constant circular speed of $200\,\mathrm{kms}$ This novel relaxation mechanism may stall the inspiral of supermassive black holes or globular clusters due to dynamical friction at radii of a few hundred pc, and can heat and expand the central regions of galaxies These processes can be used to constrain the mass of the light bosons that might comprise FDM

Relaxation in a Fuzzy Dark Matter Halo.

The largest galaxies, and in particular central galaxies in clusters, offer unique insight into understanding the mechanism for the growth of nuclear black holes. We present Hubble Space Telescope kinematics for NGC1399, the central galaxy in Fornax. We find the best-fit model contains a black hole of 5.1 +-0.7 x 10^8 Msun (at a distance of 21.1 Mpc), a factor of over 2 below the correlation of black hole mass and velocity dispersion. We also find a dramatic signature for central tangential anisotropy. The velocity profiles on adjacent sides 0.5" away from the nucleus show strong bimodality, and the central spectrum shows a large drop in the dispersion. Both of these observations point to an orbital distribution that is tangentially biased. The best-fit orbital model suggests a ratio of the tangential to radial internal velocity dispersions of three. This ratio is the largest seen in any galaxy to date and will provide an important measure for the mode by which the central black hole has grown.

/pdf/the-black-hole-mass-and-extreme-orbital-structure-in-ngc1399-30rmusxkch.pdf

Scott Tremaine

Papers

Long-lived planetesimal discs

Fitting Selected Random Planetary Systems to Titius–Bode Laws

Relaxation in a Fuzzy Dark Matter Halo.

The Black Hole Mass and Extreme Orbital Structure in NGC1399

Have interstellar clouds disrupted the oort comet cloud