Future Evolution of Nearby Large-Scale Structure in a Universe Dominated   by a Cosmological Constant
  We simulate the future evolution of the observed inhomogeneities in the local universe assuming that the global expansion rate is dominated by a cosmological constant. We find that within two Hubble times (~ 30 billion years) from the present epoch, large-scale structures will freeze in comoving coordinates and the mass distribution of bound objects will stop evolving. The Local Group will get somewhat closer to the Virgo cluster in comoving coordinates, but will be pulled away from the Virgo in physical coordinates due to the accelerated expansion of the Universe. In the distant future there will only be one massive galaxy within our event horizon, namely the merger product of the Andromeda and the Milky Way galaxies. All galaxies that are not gravitationally bound to the Local Group will recede away from us and eventually exit from our event horizon. More generally, we identify the critical interior overdensity above which a shell of matter around an object will remain bound to it at late times. 
Kentaro Nagamine & Abraham Loeb (Harvard-CfA)