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 analyze six X-ray bright elliptical galaxies, observed with Chandra and XMM-Newton, and approximate their gravitational potentials by isothermal spheres ϕ =v 2c logr over a range of radii from ∼0.5 to ∼ 25 kpc. We then compare the circularspeed v c derived from X-ray data with the estimators available from optical data. Inparticular we discuss two simple and robust procedures for evaluating the circularspeed of the galaxy using the observed optical surface brightness and the line-of-sight velocity dispersion proﬁles. The best ﬁtting relation between the circular speedsderived from optical observations of stars and X-ray observations of hot gas is v c,opt ≃η v c,X , where η = 1.10 − 1.15 (depending on the method), suggesting, albeit withlarge statistical and systematic uncertainties, that non-thermal pressure on averagecontributes ∼20-30% of the gas thermal pressure.Key words: Galaxies: Kinematics and Dynamics, X-Rays: Galaxies: Clusters 1 INTRODUCTIONIn spiral galaxies, disk rotation curves oﬀer an accurateand robust way of measuring total gravitational potentialsto distances as large as 10–30 kpc. To a ﬁrst approxi-mation the rotation curves are ﬂat over a broad rangeof radii, suggesting an isothermal (logarithmic) potentialcharacterized by ϕ(r) = v

/pdf/comparison-of-approximately-isothermal-gravitational-2u05s7jlza.pdf

Comparison of approximately isothermal gravitational potentials of elliptical galaxies based on X-ray and optical data

Fabrycky D, 2007, ASTROPHYS J, V669, P1298, DOI 10.1086-521702; Ford EB, 2000, ASTROPHYS J, V535, P385, DOI 10.1086-308815; Gauss C. F., 1866, WERKE, V3, P331; GORIACHEV NN, 1937, METHOD HALPHEN COMPU, P1; Gurkan MA, 2007, MON NOT R ASTRON SOC, V379, P1083, DOI 10.1111-j.1365-2966.2007.11982.x; HAGIHARA Y, 1971, CELESTIAL MECH 1, V2, P199; Halphen G.H., 1888, TRAITE FONCTIONS ELL, VII; HILL GW, 1882, ASTRON PAPERS AM EPH, V1, P317; Hill GW, 1901, AM J MATH, V23, P317, DOI 10.2307-2370133; Holman M, 1997, NATURE, V386, P254, DOI 10.1038-386254a0; Innanen KA, 1997, ASTRON J, V113, P1915, DOI 10.1086-118405; Jacobs V, 2001, ASTROPHYS J, V555, pL25, DOI 10.1086-321728; KOZAI Y, 1962, ASTRON J, V67, P591, DOI 10.1086-108790; LANDAU LD, 2005, CLASSICAL THEORY FIE; LASKAR J, 1986, ASTRON ASTROPHYS, V157, P59; LIBERT AS, 2008, CELEST MECH DYN ASTR, V43, P37; Mazeh T, 1997, ASTROPHYS J, V477, pL103, DOI 10.1086-310536; Murray C. D., 1999, SOLAR SYSTEM DYNAMIC; Murray N, 1999, SCIENCE, V283, P1877, DOI 10.1126-science.283.5409.1877; MUSEN P, 1963, REV GEOPHYS, V1, P85, DOI 10.1029-RG001i001p00085; Peiris HV, 2003, ASTROPHYS J, V599, P237, DOI 10.1086-378638; Rauch K. P., 1996, New Astronomy, V1, DOI 10.1016-S1384-1076(96)00012-7; SAHA P, 1992, ASTRON J, V104, P1633, DOI 10.1086-116347; Sambhus N, 2000, ASTROPHYS J, V542, P143, DOI 10.1086-309497; Sridhar S, 1999, MON NOT R ASTRON SOC, V303, P483, DOI 10.1046-j.1365-8711.1999.02218.x; Touma J, 1998, ASTRON J, V115, P1653, DOI 10.1086-300312; Touma JR, 2002, MON NOT R ASTRON SOC, V333, P583, DOI 10.1046-j.1365-8711.2002.05437.x; TREMAINE S, 2008, ARXIV08090237; Tremaine S, 2001, CELEST MECH DYN ASTR, V79, P231, DOI 10.1023-A:1017540022538; TREMAINE S, 1995, ASTRON J, V110, P628, DOI 10.1086-117548; Wu Y, 2003, ASTROPHYS J, V589, P605, DOI 10.1086-374598

Gauss's method for secular dynamics, softened

We suggest that the irregular structure in Saturn's B ring arises from the formation of shear-free ring particle assemblies of up to ~100 km in radial extent. The characteristic scale of the irregular structure is set by the competition between tidal forces and the yield stress of these assemblies; the required tensile strength of ~105 dyn cm-2 is consistent with the sticking forces observed in laboratory simulations of frosted ice particles. These assemblies could be the nonlinear outcome of a linear instability that occurs in a rotating fluid disk in which the shear stress is a decreasing function of the shear. We show that a simple model of an incompressible, non-Newtonian fluid in shear flow leads to the Cahn-Hilliard equation, which is widely used to model the formation of structure in binary alloys and other systems.

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

On the Origin of Irregular Structure in Saturn's Rings

The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipation. We develop an analytic version of this model for two-planet systems in two limiting cases: the planet mass grows quickly or slowly relative to the characteristic resonant libration time. In both cases the distribution of systems in period ratio develops a characteristic asymmetric peak-trough structure around the resonance, qualitatively similar to that observed in the Kepler sample. We verify this result with numerical integrations of the restricted three-body problem. We show that for the 3:2 resonance, where the observed peak-trough structure is strongest, our simple model is consistent with the observations for a range of mean planet masses 20-100M_{\oplus}. This mass range is higher than expected, by at least a factor of three, from the few Kepler planets with measured masses, but part of this discrepancy could be due to oversimplifications in the dynamical model or uncertainties in the planetary mass-radius relation.

Planets Near Mean-Motion Resonances

It is proposed that the incomplete rings of Neptune consist of a number of short arcs centered on the corotation resonances of a single satellite. The satellite must have a radius of the order of 100 km or more and move on an inclined orbit. Corotation resonances are located at potential maxima. Thus, mechanical energy dissipated by interparticle collisions must be continually replenished to prevent the arcs from spreading. It is shown that each corotation resonance is associated with a nearby Lindblad resonance, which excites the ring particles' orbital eccentricity, thus supplying the energy required to maintain the arc. The ultimate energy reservoir is the satellite's orbital energy. Therefore, interaction with the arcs damps the satellite's orbital inclination. The self-gravity of the arcs limits their contraction and enforces a relation between arc length and mass. The estimated arc masses are so small, of the order of 10^16 g, that the satellite's orbital inclination suffers negligible decay over the age of the solar system. The inferred surface mass densities are comparable to those found in the major rings of Saturn and Uranus.

/pdf/towards-a-theory-for-neptune-s-arc-rings-1kmk5bat9s.pdf

Scott Tremaine

Papers

Comparison of approximately isothermal gravitational potentials of elliptical galaxies based on X-ray and optical data

Gauss's method for secular dynamics, softened

On the Origin of Irregular Structure in Saturn's Rings

Planets Near Mean-Motion Resonances

Towards a theory for Neptune's arc rings