We have studied the behaviour of stellar streams in the Aquarius fully cosmological N-body simulations of the formation of Milky Way halos. In particular, we have characterised the streams in angle/frequency spaces derived using an approximate but generally well-fitting spherical potential. We have also run several test-particle simulations to understand and guide our interpretation of the different features we see in the Aquarius streams. Our goal is both to establish which deviations of the expected action-angle behaviour of streams exist because of the approximations made on the potential, but also to derive to what degree we can use these coordinates to model streams reliably. We have found that many of the Aquarius streams wrap in angle space along relatively straight lines, and also in frequency space. On the other hand, from our controlled simulations we have been able to establish that deviations from spherical symmetry, the use of incorrect potentials and the inclusion of self-gravity lead to streams in angle space to still be along relatively straight lines but also to depict wiggly behaviour whose amplitude increases as the approximation to the true potential becomes worse. In frequency space streams typically become thicker and somewhat distorted. Therefore, our analysis explains most of the features seen in the approximate angle and frequency spaces for the Aquarius streams with the exception of their somewhat `noisy' and `patchy' morphologies. These are likely due to the interactions with the large number of dark matter subhalos present in the cosmological simulations. Since the measured angle-frequency misalignments of the Aquarius streams can largely be attributed to using the wrong (spherical) potential, determining the mass growth history of these halos will only be feasible once the true potential has been determined robustly.

We study the evolution of streams in a time-dependent spherical gravitational potential. Our goal is to establish what are the imprints of this time evolution on the properties of streams as well as their observability. To this end, we have performed a suite of numerical experiments for a host system that doubles its mass during the integration time and for a variety of initial conditions. In these experiments we found that the most striking imprint is a misalignment of 10 degrees in the angular location of the apocentres of the streams compared to the static case (and to the orbit of the centre of mass), which only becomes apparent for sufficiently long streams. We have also developed an analytic model using action-angle variables which allows us to explain this behaviour and to identify the most important signature of time evolution, namely a difference in the slope defined by the distribution of particles along a stream in frequency and in angle space. Although a difference in slope can arise when the present-day potential is not correctly modelled, this shortcoming can be by-passed because in this case, streams are no longer straight lines in angle space, but depict a wiggly appearance and an implausible energy gradient. The difference in slope due to time-evolution is small, typically $\sim10^{-2}$ and its amplitude depends on the growth rate of the potential, but nonetheless we find that it could be observable if accurate full-space information for nearby long streams is available.

We present a new simple functional form to model the evolution of a spherical mass distribution in a cosmological context. Two parameters control the growth of the system and this is modelled using a redshift dependent exponential for the scale mass and scale radius. In this new model, systems form inside out and the mass of a given shell can be made to never decrease, as generally expected. This feature makes it more suitable for studying the smooth growth of galactic potentials or cosmological halos than other parametrizations often used in the literature. This is further confirmed through a comparison to the growth of dark matter halos in the Aquarius simulations.