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

Near-future surveys promise a dramatic improvement in the number and precision of astrometric, photometric and spectroscopic measurements of stars in the Milky Way's disk. We examine the impact of such surveys on our understanding of the Galaxy by "observing" particle realizations of non-axisymmetric disk distributions orbiting in an axisymmetric halo with appropriate errors and then attempting to recover the underlying potential using a Markov Chain Monte Carlo (MCMC) approach. We demonstrate that the azimuthally averaged gravitational force field in the Galactic plane--and hence, to a lesser extent, the Galactic mass distribution--can be tightly constrained over a large range of radii using a variety of types of surveys so long as the error distribution of the measurements of the parallax, proper motion and radial velocity are well-understood and the disk is surveyed globally. One advantage of our method is that the target stars can be selected non-randomly in real or apparent-magnitude space to ensure just such a global sample without biasing the results. Assuming we can always measure the line-of-sight velocity of a star with at least 1 km/s precision, we demonstrate that the force field can be determined to better than ~1% for Galactocentric radii in the range R=4-20 kpc We conclude that near-future surveys, like SIM Lite, Gaia, and VERA, will provide the first precise mapping of the gravitational force field in the region of the Galactic disk.

Probing the Galactic Potential with Next-Generation Observations of Disk Stars

In galactic nuclei, the gravitational potential is dominated by the central supermassive black hole, so stars follow quasi-Keplerian orbits. These orbits are distorted by gravitational forces from other stars, leading to long-term orbital relaxation. The direct numerical study of these processes is challenging because the fast orbital motion imposed by the central black hole requires very small timesteps. Within the secular approximation of smearing out stars along their underlying Keplerian orbits, a multipole expansion of the pairwise interaction between the stars yields an efficient numerical code to investigate the long-term evolution of their orbital parameters. These new simulations precisely recover the diffusion coefficients of stellar eccentricities obtained through analytical calculations of the secular dynamics. The computational complexity of the present method scales linearly with the total number of stars, so it should prove useful to study the long-term evolution of self-gravitating systems dominated by a central mass.

Secular Dynamics around a Supermassive Black Hole via Multipole Expansion.

We perform a systematic search for sub-parsec binary supermassive black holes (BHs) in normal broad line quasars at z<0.8, using multi-epoch SDSS spectroscopy of the broad Hbeta line. Our working model is that: only one of the two BHs in the binary is active, and dynamically dominates its own broad line region (BLR); the inactive companion BH is orbiting at a distance of a few R_BLR, where R_BLR~0.01-0.1 pc is the BLR size. We search for the expected line-of-sight acceleration of the broad line velocity from binary orbital motion by cross-correlating SDSS spectra from two epochs separated by up to several years in the quasar restframe. Out of ~700 pairs of spectra for which we have good measurements of the velocity shift between two epochs (1-sigma error~40 km/s), we detect 28 systems with significant velocity shifts in broad Hbeta, among which seven are the best candidates for the hypothesized binaries. We use the distribution of the observed accelerations (mostly non-detections) to place constraints on the abundance of such binary systems among the general quasar population. Excess variance in the velocity shift is inferred for observations separated by longer than 0.4 yr (quasar restframe). Attributing all the excess to binary motion would imply that most of the quasars in this sample must be in binaries; that the inactive BH must be on average more massive than the active one; and that the binary separation is at most a few times the size of the BLR. However, if this excess variance is partly or largely due to long-term broad line variability, the requirement of a large population of close binaries is much weakened or even disfavored for massive companions. [Abridged]

Constraining sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. I. The General Quasar Population

We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $\sigma_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{\sigma_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.

/pdf/comparison-of-simple-mass-estimators-for-slowly-rotating-o3dg4vg6j4.pdf

Comparison of simple mass estimators for slowly rotating elliptical galaxies

We describe a technique for solving the combined collisionless Boltzmann and Poisson equations in a discretised, or lattice, phase space. The time and the positions and velocities of `particles' take on integer values, and the forces are rounded to the nearest integer. The equations of motion are symplectic. In the limit of high resolution, the lattice equations become the usual integro-differential equations of stellar dynamics. The technique complements other tools for solving those equations approximately, such as $N$-body simulation, or techniques based on phase-space grids. Equilibria are found in a variety of shapes and sizes. They are true equilibria in the sense that they do not evolve with time, even slowly, unlike existing $N$-body approximations to stellar systems, which are subject to two-body relaxation. They can also be `tailor-made' in the sense that the mass distribution is constrained to be close to some pre-specified function. Their principal limitation is the amount of memory required to store the lattice, which in practice restricts the technique to modeling systems with a high degree of symmetry. We also develop a method for analysing the linear stability of collisionless systems, based on lattice equilibria as an unperturbed model.

Scott Tremaine

Papers

Probing the Galactic Potential with Next-Generation Observations of Disk Stars

Secular Dynamics around a Supermassive Black Hole via Multipole Expansion.

Constraining sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. I. The General Quasar Population

Comparison of simple mass estimators for slowly rotating elliptical galaxies

Lattice stellar dynamics