Abstract
Recent high-precision meteoritic data infers that Mars finished its accretion rapidly within 10 Myr of the beginning of the Solar System and had an accretion zone that did not entirely ...overlap with the Earth’s. Here we present a detailed study of the accretion zone of planetary embryos from high-resolution simulations of planetesimals in a disk. We found that all simulations with Jupiter and Saturn on their current eccentric orbits (EJS) result in a similar accretion zone between fast-forming Mars and Earth-region embryos. Assuming more circular orbits for Jupiter and Saturn (CJS), on the other hand, has a significantly higher chance of forming Mars with an accretion zone not entirely dominated by Earth and Venus-region embryos; however, CJS in general forms Mars slower than in EJS. By further quantifying the degree of overlap between accretion zones of embryos in different regions with the average overlap coefficient (OVL), we found that the OVL of CJS shows a better match with the OVL from a chondritic isotopic mixing model of Earth and Mars, which indicates that the giant planets are likely to have resided on more circular orbits during gas disk dissipation than they do today, matching their suggested pre-instability orbits. More samples, including those from Mercury and Venus, could potentially confirm this hypothesis.
ABSTRACT
We present new results from a series of ΛCDM simulations of cluster mass haloes resolved with high force and mass resolution. These results are compared with recently published simulations ...from groups using various codes including pkdgrav, art, tpm, grape and gadget. Careful resolution tests show that with 25 million particles within the high‐resolution region we can resolve to about 0.3 per cent of the virial radius and that convergence in radius is proportional to the mean interparticle separation. The density profiles of 26 high‐resolution clusters obtained with the different codes and from different initial conditions agree very well. The average logarithmic slope at one per cent of the virial radius is γ= 1.26 with a scatter of ±0.17. Over the entire resolved regions the density profiles are well fitted by a smooth function that asymptotes to a central cusp ρ∝r−γ, where we find γ= 1.16 ± 0.14 from the mean of the fits to our six highest‐resolution clusters.
The ever increasing size and complexity of data coming from simulations of cosmic structure formation demand equally sophisticated tools for their analysis. During the past decade, the art of object ...finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the 'halo finder comparison project': we investigate in detail the (possible) origin of any deviations across finders. To this extent, we decipher and discuss differences in halo-finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and V
max function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and - to a minor extent - the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing and observables in general.
The dark matter dominated Fornax dwarf spheroidal has five globular clusters orbiting at ∼1 kpc from its centre. In a cuspy cold dark matter halo the globulars would sink to the centre from their ...current positions within a few Gyr, presenting a puzzle as to why they survive undigested at the present epoch. We show that a solution to this timing problem is to adopt a cored dark matter halo. We use numerical simulations and analytic calculations to show that, under these conditions, the sinking time becomes many Hubble times; the globulars effectively stall at the dark matter core radius. We conclude that the Fornax dwarf spheroidal has a shallow inner density profile with a core radius constrained by the observed positions of its globular clusters. If the phase space density of the core is primordial then it implies a warm dark matter particle and gives an upper limit to its mass of ∼0.5 keV, consistent with that required to significantly alleviate the substructure problem.
Terrestrial planets are thought to be the result of a vast number of gravitational interactions and collisions between smaller bodies. We use numerical simulations to show that practically identical ...initial conditions result in a wide array of final planetary configurations. This is a result of the chaotic evolution of trajectories which are highly sensitive to minuscule displacements. We determine that differences between systems evolved from virtually identical initial conditions can be larger than the differences between systems evolved from very different initial conditions. This implies that individual simulations lack predictive power. For example, there is not a reproducible mapping between the initial and final surface density profiles. However, some key global properties can still be extracted if the statistical spread across many simulations is considered. Based on these spreads, we explore the collisional growth and orbital properties of terrestrial planets, which assemble from different initial conditions (we vary the initial planetesimal distribution, planetesimal masses, and giant planet orbits.). Confirming past work, we find that the resulting planetary systems are sculpted by sweeping secular resonances. Configurations with giant planets on eccentric orbits produce fewer and more massive terrestrial planets on tighter orbits than those with giants on circular orbits. This is further enhanced if the initial mass distribution is biased to the inner regions. In all cases, the outer edge of the system is set by the final location of the v6 resonance and we find that the mass distribution peaks at the v5 resonance. Using existing observations, we find that extrasolar systems follow similar trends. Although differences between our numerical modelling and exoplanetary systems remain, we suggest that CoRoT-7, HD 20003 and HD 20781 may host undetected giant planets.
Cold dark matter cosmogony predicts triaxial dark matter halos, whereas observations find quite round halos. This is most likely due to the condensation of baryons leading to rounder halos. We ...examine the halo phase space distribution basis for such shape changes. Triaxial halos are supported by box orbits, which pass arbitrarily close to the density center. The decrease in triaxiality caused by baryons is thought to be due to the scattering of these orbits. We test this hypothesis with simulations of disks grown inside triaxial halos. After the disks are grown we check whether the phase space structure has changed by evaporating the disks and comparing the initial and final states. While the halos are substantially rounder when the disk is at full mass, their final shape after the disk is evaporated is not much different from the initial. Likewise, the halo becomes (more) radially anisotropic when the disk is grown, but the final anisotropy is consistent with the initial. Only if the baryons are unreasonably compact or massive does the halo change irreversibly. We show that the character of individual orbits is not generally changed by the growing mass. Thus, the central condensation of baryons does not destroy enough box orbits to cause the shape change. Rather, box orbits merely become rounder along with the global potential. However, if angular momentum is transferred to the halo, either via satellites or via bars, a large irreversible change in the halo distribution occurs. The ability of satellites to alter the phase space distribution of the halo is of particular concern to galaxy formation simulations since halo triaxiality can profoundly influence the evolution of disks.
Galaxy clusters, the most massive collapsed structures, have been routinely used to determine cosmological parameters. When using clusters for cosmology, the crucial assumption is that they are ...relaxed. However, subarcminute resolution Sunyaev-Zel'dovich (SZ) effect images compared with high-resolution X-ray images of some clusters show significant offsets between the two peaks. We have carried out self-consistent N-body/hydrodynamical simulations of merging galaxy clusters using FLASH to study these offsets quantitatively. We have found that significant displacements result between the SZ and X-ray peaks for large relative velocities for all masses used in our simulations as long as the impact parameters were about 100-250 kpc. Our results suggest that the SZ peak coincides with the peak in the pressure times the line-of-sight characteristic length and not the pressure maximum (as it would for clusters in equilibrium). The peak in the X-ray emission, as expected, coincides with the density maximum of the main cluster. As a consequence, the morphology of the SZ signal, and therefore the offset between the SZ and X-ray peaks, change with viewing angle. As an application, we compare the morphologies of our simulated images to observed SZ and X-ray images and mass surface densities derived from weak-lensing observations of the merging galaxy cluster CL0152-1357, we find that a large relative velocity of 4800 km s super(-1) is necessary to explain the observations. We conclude that an analysis of the morphologies of multi-frequency observations of merging clusters can be used to put meaningful constraints on the initial parameters of the progenitors.
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.