We study the formation of discs via the cooling flow of gas within galactic haloes using smoothed particle hydrodynamic simulations. These simulations resolve mass scales of a few thousand solar ...masses in the gas component for the first time. Thermal instabilities result in the formation of numerous warm clouds that are pressure confined by the hot ambient halo gas. The clouds fall slowly on to the disc through non-spherical accretion from material flowing preferentially down the angular momentum axis. The rotational velocity of the infalling cold gas decreases as a function of height above the disc, closely resembling that of the extra-planar gas recently observed around the spiral galaxy, NGC 891.
Using analytic calculations and N-body simulations we show that in constant density (harmonic) cores, sinking satellites undergo an initial phase of very rapid (super-Chandrasekhar) dynamical ...friction, after which they experience no dynamical friction at all. For density profiles with a central power law profile, ρ∝r−α, the infalling satellite heats the background and causes α to decrease. For α < 0.5 initially, the satellite generates a small central constant density core and stalls as in the α= 0 case. We discuss some astrophysical applications of our results to decaying satellite orbits, galactic bars and mergers of supermassive black hole binaries. In a companion paper we show that a central constant density core can provide a natural solution to the timing problem for Fornax's globular clusters.
The evolution of gravitationally unstable protoplanetary gaseous disks has been studied with the use of three-dimensional smoothed particle hydrodynamics simulations with unprecedented resolution. We ...have considered disks with initial masses and temperature profiles consistent with those inferred for the protosolar nebula and for other protoplanetary disks. We show that long-lasting, self-gravitating protoplanets arise after a few disk orbital periods if cooling is efficient enough to maintain the temperature close to 50 K. The resulting bodies have masses and orbital eccentricities similar to those of detected extrasolar planets.
Remnant planetesimals might have played an important role in reducing the orbital eccentricities of the terrestrial planets after their formation via giant impacts. However, the population and the ...size distribution of remnant planetesimals during and after the giant impact stage are unknown, because simulations of planetary accretion in the runaway growth and giant impact stages have been conducted independently. Here we report results of direct N-body simulations of the formation of terrestrial planets beginning with a compact planetesimal disk. The initial planetesimal disk has a total mass and angular momentum as observed for the terrestrial planets, and we vary the width (0.3 and 0.5 AU) and the number of planetesimals (1000-5000). This initial configuration generally gives rise to three final planets of similar size, and sometimes a fourth small planet forms near the location of Mars. Since a sufficient number of planetesimals remains, even after the giant impact phase, the final orbital eccentricities are as small as those of the Earth and Venus.
We carry out a large set of very high resolution, three-dimensional, smoothed particle hydrodynamics simulations describing the evolution of gravitationally unstable gaseous protoplanetary disks. We ...consider a broad range of initial disk parameters. Disk masses out to 20 AU range from 0.075 to 0.125 M sub(o), roughly consistent with the high end of the mass distribution inferred for disks around T Tauri stars. Minimum outer temperatures range from 30 to 100 K, as expected from studies of the early protosolar nebula and suggested by the modeling of the spectra of protoplanetary disks. The mass of the central star is also varied, although it is usually assumed to be equal to that of the Sun. Overall, the initial disks span minimum Q-parameters between 0.8 and 2, with most models having Q approx 1.4. The disks are evolved assuming either a locally isothermal equation of state or an adiabatic equation of state with varying gamma. Heating by (artificial) viscosity and shocks is included when the adiabatic equation of state is used. When condensations above a specific density threshold appear as a result of gravitational instability in a locally isothermal calculation, the equation of state is switched to adiabatic to account for the increased optical depth. We show that when a disk has a minimum Q-parameter less than 1.4, strong trailing spiral instabilities, typically three- or four-armed modes, form and grow until fragmentation occurs along the arms after about 5 mean disk orbital times. The resulting clumps contract quickly to densities several orders of magnitude higher than the initial disk density, and the densest of them survive even under adiabatic conditions. These clumps are stable to tidal disruption and merge quickly, leaving two to three protoplanets on fairly eccentric orbits (the mean eccentricity being around 0.2) after approx10 super(3) yr. Fragmentation is not strongly dependent on whether the disk starts from a marginally unstable state or gradually achieves it; we show that if the disk is allowed to grow in mass from a very light, very stable state over tens of orbital times, it still fragments at roughly the same mass and temperature as in the standard disk models. We show that the first stages of the instability, until the appearance of the overdensities, can be understood in terms of the maximum unstable Toomre wavelength and the local Jeans length. A high mass and force resolution are needed to correctly resolve both scales and follow the fragmentation process appropriately. Varying disk mass and temperature affects such physical scales and hence the typical masses of the protoplanets that form. Objects smaller than Saturn or a couple of times bigger than Jupiter can both be produced by fragmentation. Their final masses will then depend on the subsequent interactions and mergers with other clumps and on the accretion of disk material. The accretion rate depends on the disk thermodynamics and is negligible with adiabatic conditions. After approx10 super(3) yr the masses range from just below 1M sub(Jup) to more than 7M sub(Jup), well in agreement with those of detected extrasolar planets.
Evolution of the mass function of dark matter haloes Reed, Darren; Gardner, Jeffrey; Quinn, Thomas ...
Monthly notices of the Royal Astronomical Society,
December 2003, Letnik:
346, Številka:
2
Journal Article
Recenzirano
Odprti dostop
We use a high-resolution ΛCDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of 15, in a 50 h−1 Mpc volume ...containing 80 million particles. Our low-redshift results allow us to probe low-σ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth & Tormen mass function provides an excellent match to all of our data except for redshifts of 10 and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 per cent for the 1010–1011 M⊙ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth & Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high-redshift (z≃ 6) QSOs, which are thought to be hosted by haloes at 5σ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50 per cent. We test the sensitivity of our results to force accuracy, starting redshift and halo-finding algorithm.
Polar ring galaxies are flattened stellar systems with an extended ring of gas and stars rotating in a plane almost perpendicular to the central galaxy. We show that their formation can occur ...naturally in a hierarchical universe where most low-mass galaxies are assembled through the accretion of cold gas infalling along megaparsec-scale filamentary structures. Within a large cosmological hydrodynamical simulation, we find a system that closely resembles the classic polar ring galaxy NGC 4650A. How galaxies acquire their gas is a major uncertainty in models of galaxy formation, and recent theoretical work has argued that cold accretion plays a major role. This idea is supported by our numerical simulations and the fact that polar ring galaxies are typically low-mass systems.
Dark matter subhaloes in numerical simulations Reed, Darren; Governato, Fabio; Quinn, Thomas ...
Monthly notices of the Royal Astronomical Society,
June 2005, Letnik:
359, Številka:
4
Journal Article
Recenzirano
Odprti dostop
We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles ...within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 1011 to 1014 M⊙ at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation dn/dv = (1/8)(vcmax/vcmax,host)−4. Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for r ≳ 0.3rvir, where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of P ∝ kn, and demonstrate that a steeper spectral index, n, results in significantly less substructure.
Since the discovery of the "Bullet Cluster," several similar cases have been uncovered that suggest relative velocities well beyond the tail of high speed collisions predicted by the concordance ACDM ...model. However, quantifying such post-merger events with hydrodynamical models requires a wide coverage of possible initial conditions. Here, we show that it is simpler to interpret pre-merger cases, such as A1750, where the gas between the colliding clusters is modestly affected, so that the initial conditions are clear. We analyze publicly available Chandra data confirming a significant increase in the projected X-ray temperature between the two cluster centers in A1750 consistent with our expectations for a merging cluster. We model this system with a self-consistent hydrodynamical simulation of dark matter and gas using the FLASH code. Our simulations reproduce well the X-ray data and the measured redshift difference between the two clusters in the phase before the first core passage viewed at an intermediate projection angle. The deprojected initial relative velocity derived using our model is 1460 km s super(-1), which is considerably higher than the predicted mean impact velocity for simulated massive haloes derived by recent LambdaCDM cosmological simulations, but is within the allowed range. Our simulations demonstrate that such systems can be identified using a multi-wavelength approach and numerical simulations, for which the statistical distribution of relative impact velocities may provide a definitive examination of a broad range of dark matter scenarios.