The interplanetary space probe Cassini/Huygens reached Saturn in July 2004 after seven years of cruise phase. The Cosmic Dust Analyzer (CDA) measures the interplanetary, interstellar and planetary dust in our solar system since 1999 and provided unique discoveries. In 1999, CDA detected interstellar dust in the inner solar system followed by the detection of electrical charges of interplanetary dust grains during the cruise phase between Earth and Jupiter. The instrument determined the composition of interplanetary dust and the nanometre sized dust streams originating from Jupiter's moon Io. During the approach to Saturn in 2004, similar streams of submicron grains with speeds in the order of 100 km/s were detected from Saturn's inner and outer ring system and are released to the interplanetary magnetic field. Since 2004 CDA measured more than one million dust impacts characterizing the dust environment of Saturn. The instrument is one of three experiments which discovered the active ice geysers located at the south pole of Saturn's moon Enceladus in 2005. Later, a detailed compositional analysis of the water ice grains in Saturn's E ring system lead to the discovery of large reservoirs of liquid water (oceans) below the icy crust of Enceladus. Finally, the determination of the dust- magnetosphere interaction and the discovery of the extended E ring (at least twice as large as predicted) allowed the definition of a dynamical dust model of Saturn's E ring describing the observed properties. This paper summarizes the discoveries of a ten year story of success based on reliable measurements with the most advanced dust detector flown in space until today. This paper focuses on cruise results and findings achieved at Saturn with a focus on flux and density measurements.

We establish for the first time heuristic correlations between harmonic space phase information and higher order statistics. Using the spherical full-sky maps of the cosmic microwave background as an example we demonstrate that known phase correlations at large spatial scales can gradually be diminished when subtracting a suitable best-fit (Bianchi-) template map of given strength. The weaker phase correlations lead in turn to a vanishing signature of anisotropy when measuring the Minkowski functionals and scaling indices in real-space and comparing them with surrogate maps being free of phase correlations. Those investigations can open a new road to a better understanding of signatures of non-Gaussianities in complex spatial structures by elucidating the meaning of Fourier phase correlations and their influence on higher order statistics.

A novel mechanism of the formation of compact dusty clouds in astrophysical environments is discussed. It is shown that the balance of collective forces operating in space dusty plasmas can result in the effect of dust self-confinement, generating equilibrium spherical clusters. The distribution of dust and plasma density inside such objects and their stability are investigated. Spherical dusty clouds can be formed in a broad range of plasma parameters, suggesting that this process of dust self-organization might be a generic phenomenon occurring in different astrophysical media. We argue that compact dusty clouds can represent condensation seeds for a population of small-scale, cold, gaseous clumps in the diffuse interstellar medium. They could play an important role in regulating its small-scale structure and its thermodynamical evolution.

In order to test for non-Gaussianities with respect to scale-dependencies we use so-called surrogate maps, in which possible phase correlations of the Fourier phases of the original WMAP data and simulations, respectively, are destroyed by applying a shuffling scheme to the maps. A statistical comparison of the original maps with the surrogate maps then allows to test for the existence of higher order correlations in the original maps, also and especially on well-defined Fourier modes. Using Minkowski functionals as an image analysis technique we calculate the deviation between the original data and 500 surrogates for different hemispheres in the sky and find ecliptic hemispherical asymmetries between northern and southern ecliptic sky. We find strong deviations from Gaussianity in the WMAP data when considering the low-l range with l = [2,20]. The analysis technique of the scaling indices leads to a slightly lower deviation. Although the underlying foreground reduction methods of the maps differ from each other, we find similar results for the WMAP 7yr ILC map and the WMAP 7yr (needlet-based) NILC map in the low-l range. Our results point once more to a cosmological nature of the signal. For a higher l range with l = [120,300] the results differ between the two image analysis techniques and between the two maps which makes an intrinsic nature of the signal on this l range less likely. When we decrease the size of the analysed sky regions for the low-l study, we do not find signatures of NG in the northern sky. In the south we find individual spots which show deviations from Gaussianity. In addition, we investigate non-Gaussian CMB simulations that depend on the f_NL parameter of the local type. These simulations cannot account for the detected signatures on the low-l range.

We demonstrate the feasibility to generate surrogates by Fourier-based methods for an incomplete data set. This is performed for the case of a CMB analysis, where astrophysical foreground emission, mainly present in the Galactic plane, is a major challenge. The shuffling of the Fourier phases for generating surrogates is now enabled by transforming the spherical harmonics into a new set of basis functions that are orthonormal on the cut sky. The results show that non-Gaussianities and hemispherical asymmetries in the CMB as identified in several former investigations, can still be detected even when the complete Galactic plane (|b| < 30{\deg}) is removed. We conclude that the Galactic plane cannot be the dominant source for these anomalies. The results point towards a violation of statistical isotropy.

In the recent years, non-Gaussianity and statistical isotropy of the Cosmic Microwave Background (CMB) was investigated with various statistical measures, first and foremost by means of the measurements of the WMAP satellite. In this Review, we focus on the analyses that were accomplished with a measure of local type, the so-called Scaling Index Method (SIM). The SIM is able to detect structural characteristics of a given data set, and has proven to be highly valuable in CMB analysis. It was used for comparing the data set with simulations as well as surrogates, which are full sky maps generated by randomisation of previously selected features of the original map. During these investigations, strong evidence for non-Gaussianities as well as asymmetries and local features could be detected. In combination with the surrogates approach, the SIM detected the highest significances for non-Gaussianity to date.

We present a model-independent investigation of the WMAP data with respect to scale- dependent non-Gaussianities (NGs) by employing the method of constrained randomization. For generating so-called surrogate maps a shuffling scheme is applied to the Fourier phases of the original data, which allows to test for the presence of higher order correlations (HOCs) on well-defined scales. Using scaling indices as test statistics we find highly significant signatures for non-Gaussianities when considering all scales. We test for NGs in the bands l = [2,20], l = [20,60], l = [60,120] and l = [120,300]. We find highly significant signatures for non-Gaussianities and ecliptic hemispherical asymmetries for l = [2, 20]. We also obtain highly significant deviations from Gaussianity for the band l = [120,300]. The result for the full l-range can be interpreted as a superposition of the signatures found in the bands l = [2, 20] and l = [120, 300]. We find remarkably similar results when analyzing different ILC-like maps. We perform a set of tests to investigate if the detected anomalies can be explained by systematics. While no test can convincingly rule out the intrinsic nature of the anomalies for the low l case, the ILC map making procedure and/or residual noise in the maps can also lead to NGs at small scales. Our investigations prove that there are phase correlations in the WMAP data of the CMB. In the absence of an explanation in terms of Galactic foregrounds or known systematic artefacts, the signatures at low l must so far be taken to be cosmological at high significance. These findings strongly disagree with predictions of isotropic cosmologies with single field slow roll inflation. The task is now to elucidate the origin of the phase correlations and to understand the physical processes leading to these scale-dependent non-Gaussianities - if systematics as cause for them must be ruled out.

We discuss a novel mechanism of dust acceleration which may dominate for particles smaller than $\sim0.1~\mu$m. The acceleration is caused by their direct electrostatic interactions arising from fluctuations of grain charges. The energy source for the acceleration are the irreversible plasma processes occurring on the grain surfaces. We show that this mechanism of charge-fluctuation-induced acceleration likely affects the rate of grain coagulation and shattering of the population of small grains.

Probing Gaussianity represents one of the key questions in modern cosmology, because it allows to discriminate between different models of inflation. We test for large-scale non-Gaussianities in the cosmic microwave background (CMB) in a model-independent way. To this end, so-called first and second order surrogates are generated by first shuffling the Fourier phases belonging to the scales not of interest and then shuffling the remaining phases for the length scales under study. Using scaling indices as test statistics we find highly significant signatures for both non-Gaussianities and asymmetries on large scales for the WMAP data of the CMB. We find remarkably similar results when analyzing different ILC-maps based on the WMAP five and seven year data. Such features being independent from the map-making procedure would disfavor the fundamental principle of isotropy as well as canonical single-field slow-roll inflation - unless there is some undiscovered systematic error in the collection or reduction of the CMB data or yet unknown foreground contributions.

The Ulysses spacecraft has been orbiting the Sun on a highly inclined ellipse since it encountered Jupiter in February 1992. Since then it made almost three revolutions about the Sun. Here we report on the final three years of data taken by the on-board dust detector. During this time, the dust detector recorded 609 dust impacts of particles with masses 10^-16 g <= m <= 10^-7 g, bringing the mission total to 6719 dust data sets. The impact rate varied from a low value of 0.3 per day at high ecliptic latitudes to 1.5 per day in the inner solar system. The impact direction of the majority of impacts between 2005 and 2007 is compatible with particles of interstellar origin, the rest are most likely interplanetary particles. We compare the interstellar dust measurements from 2005/2006 with the data obtained during earlier periods (1993/1994) and (1999/2000) when Ulysses was traversing the same spatial region at southern ecliptic latitudes but the solar cycle was at a different phase. During these three intervals the impact rate of interstellar grains varied by more than a factor of two. Furthermore, in the two earlier periods the grain impact direction was in agreement with the flow direction of the interstellar helium while in 2005/2006 we observed a shift in the approach direction of the grains by approximately 30 deg away from the ecliptic plane. The reason for this shift remains unclear but may be connected with the configuration of the interplanetary magnetic field during solar maximum. We also find that the dust measurements are in agreement with the interplanetary flux model of Staubach et al. (1997) which was developed to fit a 5-year span of Ulysses data.

The Galileo spacecraft was orbiting Jupiter between Dec 1995 and Sep 2003. The Galileo dust detector monitored the jovian dust environment between about 2 and 370 R_J (jovian radius R_J = 71492 km). We present data from the Galileo dust instrument for the period January 2000 to September 2003. We report on the data of 5389 particles measured between 2000 and the end of the mission in 2003. The majority of the 21250 particles for which the full set of measured impact parameters (impact time, impact direction, charge rise times, charge amplitudes, etc.) was transmitted to Earth were tiny grains (about 10 nm in radius), most of them originating from Jupiter's innermost Galilean moon Io. Their impact rates frequently exceeded 10 min^-1. Surprisingly large impact rates up to 100 min^-1 occurred in Aug/Sep 2000 when Galileo was at about 280 R_J from Jupiter. This peak in dust emission appears to coincide with strong changes in the release of neutral gas from the Io torus. Strong variability in the Io dust flux was measured on timescales of days to weeks, indicating large variations in the dust release from Io or the Io torus or both on such short timescales. Galileo has detected a large number of bigger micron-sized particles mostly in the region between the Galilean moons. A surprisingly large number of such bigger grains was measured in March 2003 within a 4-day interval when Galileo was outside Jupiter's magnetosphere at approximately 350 R_J jovicentric distance. Two passages of Jupiter's gossamer rings in 2002 and 2003 provided the first actual comparison of in-situ dust data from a planetary ring with the results inferred from inverting optical images.

We continue the analysis of non-Gaussianities in the CMB by means of the scaling index method (SIM, Raeth, Schuecker & Banday 2007) by applying this method on the 5-year WMAP data. We compare each of the results with 1000 Monte Carlo simulations mimicing the Gaussian properties of the best fit $\Lambda CDM$-model. Based on the scaling indices, scale-dependent empirical probability distributions, moments of these distributions and $\chi^2$-combinations of them are calculated, obtaining similar results as in the former analysis of the 3-year data: We derive evidence for non-Gaussianity with a probability of up to 97.3% for the mean when regarding the KQ75-masked full sky and summing up over all considered length scales by means of a diagonal $\chi^2$-statistics. Looking at only the northern or southern hemisphere, we obtain up to 98.5% or 96.6%, respectively. For the standard deviation, these results appear as 95.6% for the full sky (99.7% north, 89.4% south) and for a $\chi^2$-combination of both measurements as 97.4% (99.1% north, 95.5% south). By performing an analysis of rotated hemispheres, we detect an obvious asymmetry in the data. In addition to these investigations, we present a method of filling the mask with Gaussian noise to eliminate boundary effects caused by the mask. With the help of this technique, we identify several local features on the map, of which the most significant one turns out to be the well-known cold spot. When excluding all these spots from the analysis, the deviation from Gaussianity increases, which shows that the discovered local anomalies are not the reason of the global detection of non-Gaussianity, but actually were damping the deviations on average. Our analyses per band and per year suggest, however, that it is very unlikely that the detected anomalies are due to foreground effects.

We present a model-independent method to test for scale-dependent non-Gaussianities in combination with scaling indices as test statistics. Therefore, surrogate data sets are generated, in which the power spectrum of the original data is preserved, while the higher order correlations are partly randomised by applying a scale-dependent shuffling procedure to the Fourier phases. We apply this method to the WMAP data of the cosmic microwave background (CMB) and find signatures for non-Gaussianities on large scales. Further tests are required to elucidate the origin of the detected anomalies.

String theoretical arguments led to the hypothesis that the ratio of viscosity to entropy of any physical system has a lower bound. Strongly coupled systems usually have a small viscosity compared to weakly coupled plasmas in which the viscosity is proportional to the mean free path. In the case of a one-component plasma the viscosity as a function of the coupling strength shows a minimum. Here we show that the ratio of viscosity to entropy of a strongly coupled one-component plasma is always above the lower bound predicted by string theory.

The method of constrained randomisation is applied to three-dimensional simulated galaxy distributions. With this technique we generate for a given data set surrogate data sets which have the same linear properties as the original data whereas higher order or nonlinear correlations are not preserved. The analysis of the original and surrogate data sets with measures, which are sensitive to nonlinearities, yields information about the existence of nonlinear correlations in the data. We demonstrate how to generate surrogate data sets from a given point distribution, which have the same linear properties (power spectrum) as well as the same density amplitude distribution. We propose weighted scaling indices as a nonlinear statistical measure to quantify local morphological elements in large scale structure. Using surrogates is is shown that the data sets with the same 2-point correlation functions have slightly different void probability functions and especially a different set of weighted scaling indices. Thus a refined analysis of the large scale structure becomes possible by calculating local scaling properties whereby the method of constrained randomisation yields a vital tool for testing the performance of statistical measures in terms of sensitivity to different topological features and discriminative power.

The dust detector system onboard Galileo records dust impacts in circumjovian space since the spacecraft has been injected into a bound orbit about Jupiter in December 1995. This is the sixth in a series of papers dedicated to presenting Galileo and Ulysses dust data. We present data from the Galileo dust instrument for the period January to December 1996 when the spacecraft completed four orbits about Jupiter (G1, G2, C3 and E4). Data were obtained as high resolution realtime science data or recorded data during a time period of 100 days, or via memory read-outs during the remaining times. Because the data transmission rate of the spacecraft is very low, the complete data set (i. e. all parameters measured by the instrument during impact of a dust particle) for only 2% (5353) of all particles detected could be transmitted to Earth; the other particles were only counted. Together with the data for 2883 particles detected during Galileo's interplanetary cruise and published earlier, complete data of 8236 particles detected by the Galileo dust instrument from 1989 to 1996 are now available. The majority of particles detected are tiny grains (about 10 nm in radius) originating from Jupiter's innermost Galilean moon Io. These grains have been detected throughout the Jovian system and the highest impact rates exceeded $\rm 100 min^{-1}$. A small number of grains has been detected in the close vicinity of the Galilean moons Europa, Ganymede and Callisto which belong to impact-generated dust clouds formed by (mostly submicrometer sized) ejecta from the surfaces of the moons (Kr\"uger et al., Nature, 399, 558, 1999). Impacts of submicrometer to micrometer sized grains have been detected thoughout the Jovian system and especially in the region between the Galilean moons.

The Ulysses spacecraft is orbiting the Sun on a highly inclined ellipse ($ i = 79^{\circ}$, perihelion distance 1.3 AU, aphelion distance 5.4 AU). Between January 1996 and December 1999 the spacecraft was beyond 3 AU from the Sun and crossed the ecliptic plane at aphelion in May 1998. In this four-year period 218 dust impacts were recorded with the dust detector on board. We publish and analyse the complete data set of both raw and reduced data for particles with masses $\rm 10^{-16} g$ to $\rm 10^{-8}$ g. Together with 1477 dust impacts recorded between launch of Ulysses and the end of 1995 published earlier \cite{gruen1995c,krueger1999b}, a data set of 1695 dust impacts detected with the Ulysses sensor between October 1990 and December 1999 is now available. The impact rate measured between 1996 and 1999 was relatively constant with about 0.2 impacts per day. The impact direction of the majority of the impacts is compatible with particles of interstellar origin, the rest are most likely interplanetary particles. The observed impact rate is compared with a model for the flux of interstellar dust particles. The flux of particles several micrometers in size is compared with the measurements of the dust instruments on board Pioneer 10 and Pioneer 11 beyond 3 AU (Humes 1980, JGR, 85, 5841--5852, 1980). Between 3 and 5 AU, Pioneer results predict that Ulysses should have seen five times more ($\rm \sim 10 \mu m$ sized) particles than actually detected.

The gas-to-dust mass ratios found for interstellar dust within the Solar System, versus values determined astronomically for the cloud around the Solar System, suggest that large and small interstellar grains have separate histories, and that large interstellar grains preferentially detected by spacecraft are not formed exclusively by mass exchange with nearby interstellar gas. Observations by the Ulysses and Galileo satellites of the mass spectrum and flux rate of interstellar dust within the heliosphere are combined with information about the density, composition, and relative flow speed and direction of interstellar gas in the cloud surrounding the solar system to derive an in situ value for the gas-to-dust mass ratio, $R_{g/d} = 94^{+46}_{-38}$. Hubble observations of the cloud surrounding the solar system yield a gas-to-dust mass ratio of Rg/d=551+61-251 when B-star reference abundances are assumed. The exclusion of small dust grains from the heliosheath and heliosphere regions are modeled, increasing the discrepancy between interstellar and in situ observations. The shock destruction of interstellar grains is considered, and comparisons are made with interplanetary and presolar dust grains.

Between Jan 1993 and Dec 1995 the Galileo spacecraft traversed interplanetary space between Earth and Jupiter and arrived at Jupiter on 7 Dec 1995. The dust instrument onboard was operating during most of the time. A relatively constant impact rate of interplanetary and interstellar (big) particles of 0.4 impacts per day was detected over the whole three-year time span. In the outer solar system (outside about 2.6 AU) they are mostly of interstellar origin, whereas in the inner solar system they are mostly interplanetary particles. Within about 1.7 AU from Jupiter intense streams of small dust particles were detected with impact rates of up to 20,000 per day whose impact directions are compatible with a Jovian origin. Two different populations of dust particles were detected in the Jovian magnetosphere: small stream particles during Galileo's approach to the planet and big particles concentrated closer to Jupiter between the Galilean satellites. There is strong evidence that the dust stream particles are orders of magnitude smaller in mass and faster than the instrument's calibration, whereas the calibration is valid for the big particles. Because the data transmission rate was very low, the complete data set for only a small fraction (2525) of all detected particles could be transmitted to Earth; the other particles were only counted. Together with the 358 particles published earlier, information about 2883 particles detected by the dust instrument during Galileo's six years' journey to Jupiter is now available.

The Ulysses spacecraft is orbiting the Sun on a highly inclined ellipse ($i = 79^{\circ}$). After its Jupiter flyby in 1992 at a heliocentric distance of 5.4 AU, the spacecraft reapproached the inner solar system, flew over the Sun's south polar region in September 1994, crossed the ecliptic plane at a distance of 1.3 AU in March 1995, and flew over the Sun's north polar region in July 1995. We report on dust impact data obtained with the dust detector onboard Ulysses between January 1993 and December 1995. We publish and analyse the complete data set of 509 recorded impacts of dust particles with masses between $10^{-16}$ g to $10^{-7}$ g. Together with 968 dust impacts from launch until the end of 1992 published earlier (Gr\"un et al., 1995, {\em Planet. Space Sci}, Vol. 43, p. 971-999), information about 1477 particles detected with the Ulysses sensor between October 1990 and December 1995 is now available. The impact rate measured between 1993 and 1995 stayed relatively constant at about 0.4 impacts per day and varied by less than a factor of ten. Most of the impacts recorded outside about 3.5 AU are compatible with particles of interstellar origin. Two populations of interplanetary particles have been recognised: big micrometer-sized particles close to the ecliptic plane and small sub-micrometer-sized particles at high ecliptic latitudes. The observed impact rate is compared with a model for the flux of interstellar dust particles which gives relatively good agreement with the observed impact rate. No change in the instrument's noise characteristics or degradation of the channeltron could be revealed during the three-year period.