A long-range force acting only between nonbaryonic particles would be associated with a large violation of the weak equivalence principle. We explore cosmological consequences of this idea, which we label ReBEL (daRk Breaking Equivalence principLe). A high resolution hydrodynamical simulation of the distributions of baryons and dark matter confirms our previous findings that a ReBEL force of comparable strength to gravity on comoving scales of about 1 Mpc/h causes voids between the concentrations of large galaxies to be more nearly empty, suppresses accretion of intergalactic matter onto galaxies at low redshift, and produces an early generation of dense dark matter halos. A preliminary analysis indicates the ReBEL scenario is consistent with the one-dimensional power spectrum of the Lyman-Alpha forest and the three-dimensional galaxy auto-correlation function. Segregation of baryons and DM in galaxies and systems of galaxies is a strong prediction of ReBEL. ReBEL naturally correlates the baryon mass fraction in groups and clusters of galaxies with the system mass, in agreement with recent measurements.

Numerical simulations of the effect of a long-range scalar interaction (LRSI) acting only on nonbaryonic dark matter, with strength comparable to gravity, show patterns of disruption of satellites that can agree with what is seen in the Milky Way. This includes the symmetric Sagittarius stellar stream. The exception presented here to the Kesden and Kamionkowski demonstration that an LRSI tends to produce distinctly asymmetric streams follows if the LRSI is strong enough to separate the stars from the dark matter before tidal disruption of the stellar component, and if stars dominate the mass in the luminous part of the satellite. It requires that the Sgr galaxy now contains little dark matter, which may be consistent with the Sgr stellar velocity dispersion, for in the simulation the dispersion at pericenter exceeds virial. We present other examples of simulations in which a strong LRSI produces satellites with large mass-to-light ratio, as in Draco, or free streams of stars, which might be compared to "orphan" streams.

We address the question of whether or not assembly bias arises in the absence of highly non-linear effects such as tidal stripping of halos near larger mass concentrations. Therefore, we use a simplified dynamical scheme where these effects are not modeled. We choose the punctuated Zel'dovich (PZ) approximation, which prevents orbit mixing by coalescing particles coming within a critical distance of each other. A numerical implementation of this approximation is fast, allowing us to run a large number of simulations to study assembly bias. We measure an assembly bias from 60 PZ simulations, each with 512^3 cold particles in a 128h^-1 Mpc cubic box. The assembly bias estimated from the correlation functions at separations < 5h^-1 Mpc for objects (halos) at z=0 is comparable to that obtained in full N-body simulations. For masses 4x10^11 h^-1 Mo the "oldest" 10% haloes are 3-5 times more correlated than the "youngest" 10%. The bias weakens with increasing mass, also in agreement with full N-body simulations. We find that that halo ages are correlated with the dimensionality of the surrounding linear structures as measured by the parameter (\lambda_1+\lambda_2+\lambda_3)/ (\lambda_1^2+\lambda_2^2+\lambda_3^2)^{1/2} where \lambda_i are proportional to the eigenvalues of the velocity deformation tensor. Our results suggest that assembly bias may already be encoded in the early stages of the evolution.