With the purpose of understanding how time delay interferometry (TDI) combinations can best be used for the characterisation of LISA instrumental noise, we revisit their laser frequency noise cancellation properties. We have developed an algorithm to search for all possible combinations that suppress noise at the same level as the X, Y and Z classical combination. The algorithm calculates delays using symbolic formulas that explicitly include velocities and accelerations of satellites up to the relevant order. In addition, once a combination has been identified, delays are verified by solving numerically the relevant equations using Keplerian orbits and Shapiro delay corrections. We find that the number of combinations that suppress the noise is larger than what was reported in the literature. In particular we find that some combinations that were thought to only partly suppress the noise, in reality do suppress it at the same level of accuracy as the basic X, Y and Z combinations.

LISA Pathfinder (LPF) has been a space-based mission designed to test new technologies that will be required for a gravitational wave observatory in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime (mHz and below), the measurement band of interest for a space-based observatory. The magnetic field can couple to the magnetic susceptibility and remanent magnetic moment from the test masses and disturb them from their geodesic movement. LISA Pathfinder carried on-board a dedicated magnetic measurement subsystem with noise levels of 10 $ \rm nT \ Hz^{-1/2}$ from 1 Hz down to 1 mHz. In this paper we report on the magnetic measurements throughout LISA Pathfinder operations. We characterise the magnetic environment within the spacecraft, study the time evolution of the magnetic field and its stability down to 20 $\mu$Hz, where we measure values around 200 $ \rm nT \ Hz^{-1/2}$ and identify two different frequency regimes, one related to the interplanetary magnetic field and the other to the magnetic field originating inside the spacecraft. Finally, we characterise the non-stationary component of the fluctuations of the magnetic field below the mHz and relate them to the dynamics of the solar wind.

We report on the results of the LISA Pathfinder (LPF) free-fall mode experiment, in which the control force needed to compensate the quasistatic differential force acting on two test masses is applied intermittently as a series of "impulse" forces lasting a few seconds and separated by roughly 350 s periods of true free fall. This represents an alternative to the normal LPF mode of operation in which this balancing force is applied continuously, with the advantage that the acceleration noise during free fall is measured in the absence of the actuation force, thus eliminating associated noise and force calibration errors. The differential acceleration noise measurement presented here with the free-fall mode agrees with noise measured with the continuous actuation scheme, representing an important and independent confirmation of the LPF result. An additional measurement with larger actuation forces also shows that the technique can be used to eliminate actuation noise when this is a dominant factor.

LISA Pathfinder (LPF) was a technology pioneering mission designed to test key technologies required for gravitational wave detection in space. In the low frequency regime (milli-Hertz and below), where space-based gravitational wave observatories will operate, temperature fluctuations play a crucial role since they can couple into the interferometric measurement and the test masses' free-fall accuracy in many ways. A dedicated temperature measurement subsystem, with noise levels in 10$\,\mu$K$\,$Hz$^{-1/2}$ down to $1\,$mHz was part of the diagnostics unit on board LPF. In this paper we report on the temperature measurements throughout mission operations, characterize the thermal environment, estimate transfer functions between different locations and report temperature stability (and its time evolution) at frequencies as low as 10$\,\mu$Hz, where typically values around $1\,$K$\,$Hz$^{-1/2}$ were measured.

Non-recurrent short term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n$^{-1}$ were observed between 2016 February 18 and 2017 July 3 aboard the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5$\times$10$^6$ km from Earth. The energy dependence of three Forbush decreases (FDs) is studied and reported here. A comparison of these observations with others carried out in space down to the energy of a few tens of MeV n$^{-1}$ shows that the same GCR flux parameterization applies to events of different intensity during the main phase. FD observations in L1 with LPF and geomagnetic storm occurrence is also presented. Finally, the characteristics of GCR flux non-recurrent variations (peaks and depressions) of duration $<$ 2 days and their association with interplanetary structures are investigated. It is found that, most likely, plasma compression regions between subsequent corotating high-speed streams cause peaks, while heliospheric current sheet crossing cause the majority of the depressions.

Since the 2017 Nobel Prize in Physics was awarded for the observation of gravitational waves, it is fair to say that the epoch of gravitational wave astronomy (GWs) has begun. However, a number of interesting sources of GWs can only be observed from space. To demonstrate the feasibility of the Laser Interferometer Space Antenna (LISA), a future gravitational wave observatory in space, the LISA Pathfinder satellite was launched on December, 3rd 2015. Measurements of the spurious forces accelerating an otherwise free-falling test mass, and detailed investigations of the individual subsystems needed to achieve the free-fall, have been conducted throughout the mission. This overview article starts with the purpose and aim of the mission, explains satellite hardware and mission operations and ends with a summary of selected important results and an outlook towards LISA. From the LISA Pathfinder experience, we can conclude that the proposed LISA mission is feasible.

The science operations of the LISA Pathfinder mission has demonstrated the feasibility of sub-femto-g free-fall of macroscopic test masses necessary to build a LISA-like gravitational wave observatory in space. While the main focus of interest, i.e. the optical axis or the $x$-axis, has been extensively studied, it is also of interest to evaluate the stability of the spacecraft with respect to all the other degrees of freedom. The current paper is dedicated to such a study, with a focus set on an exhaustive and quantitative evaluation of the imperfections and dynamical effects that impact the stability with respect to its local geodesic. A model of the complete closed-loop system provides a comprehensive understanding of each part of the in-loop coordinates spectra. As will be presented, this model gives very good agreements with LISA Pathfinder flight data. It allows one to identify the physical noise source at the origin and the physical phenomena underlying the couplings. From this, the performances of the stability of the spacecraft, with respect to its geodesic, are extracted as a function of frequency. Close to $1 mHz$, the stability of the spacecraft on the $X_{SC}$, $Y_{SC}$ and $Z_{SC}$ degrees of freedom is shown to be of the order of $5.0\ 10^{-15} m\ s^{-2}/\sqrt{Hz}$ for X and $4.0 \ 10^{-14} m\ s^{-2}/\sqrt{Hz}$ for Y and Z. For the angular degrees of freedom, the values are of the order $3\ 10^{-12} rad\ s^{-2}/\sqrt{Hz}$ for $\Theta_{SC}$ and $3\ 10^{-13} rad\ s^{-2}/\sqrt{Hz}$ for $H_{SC}$ and $\Phi_{SC}$.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

LISA Pathfinder (LPF) was a European Space Agency mission with the aim to test key technologies for future space-borne gravitational-wave observatories like LISA. The main scientific goal of LPF was to demonstrate measurements of differential acceleration between free-falling test masses at the sub-femto-g level, and to understand the residual acceleration in terms of a physical model of stray forces, and displacement readout noise. A key step toward reaching the LPF goals was the correct calibration of the dynamics of LPF, which was a three-body system composed by two test-masses enclosed in a single spacecraft, and subject to control laws for system stability. In this work, we report on the calibration procedures adopted to calculate the residual differential stray force per unit mass acting on the two test-masses in their nominal positions. The physical parameters of the adopted dynamical model are presented, together with their role on LPF performance. The analysis and results of these experiments show that the dynamics of the system was accurately modeled and the dynamical parameters were stationary throughout the mission. Finally, the impact and importance of calibrating system dynamics for future space-based gravitational wave observatories is discussed.

Galactic cosmic-ray (GCR) energy spectra observed in the inner heliosphere are modulated by the solar activity, the solar polarity and structures of solar and interplanetary origin. A high counting rate particle detector (PD) aboard LISA Pathfinder (LPF), meant for subsystems diagnostics, was devoted to the measurement of galactic cosmic-ray and solar energetic particle integral fluxes above 70 MeV n$^{-1}$ up to 6500 counts s$^{-1}$. PD data were gathered with a sampling time of 15 s. Characteristics and energy-dependence of GCR flux recurrent depressions and of a Forbush decrease dated August 2, 2016 are reported here. The capability of interplanetary missions, carrying PDs for instrument performance purposes, in monitoring the passage of interplanetary coronal mass ejections (ICMEs) is also discussed.

Test mass charging caused by cosmic rays will be a significant source of acceleration noise for space-based gravitational wave detectors like LISA. Operating between December 2015 and July 2017, the technology demonstration mission LISA Pathfinder included a bespoke monitor to help characterise the relationship between test mass charging and the local radiation environment. The radiation monitor made in situ measurements of the cosmic ray flux while also providing information about its energy spectrum. We describe the monitor and present measurements which show a gradual 40% increase in count rate coinciding with the declining phase of the solar cycle. Modulations of up to 10% were also observed with periods of 13 and 26 days that are associated with co-rotating interaction regions and heliospheric current sheet crossings. These variations in the flux above the monitor detection threshold (approximately 70 MeV) are shown to be coherent with measurements made by the IREM monitor on-board the Earth orbiting INTEGRAL spacecraft. Finally we use the measured deposited energy spectra, in combination with a GEANT4 model, to estimate the galactic cosmic ray differential energy spectrum over the course of the mission.

This paper is a follow-up of a previous paper about the M82 galaxy and its halo based on Planck observations. As in the case of M82, so also for the M81 galaxy a substantial North-South and East-West temperature asymmetry is found, extending up to galactocentric distances of about $1.5^\circ$. The temperature asymmetry is almost frequency independent and can be interpreted as a Doppler-induced effect related to the M81 halo rotation and/or triggered by the gravitational interaction of the galaxies within the M81 Group. Along with the analogous study of several nearby edge-on spiral galaxies, the CMB temperature asymmetry method thus is shown to act as a direct tool to map the galactic haloes and/or the intergalactic bridges, invisible in other bands or by other methods.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

We report on electrostatic measurements made on board the European Space Agency mission LISA Pathfinder. Detailed measurements of the charge-induced electrostatic forces exerted on free-falling test masses (TMs) inside the capacitive gravitational reference sensor are the first made in a relevant environment for a space-based gravitational wave detector. Employing a combination of charge control and electric-field compensation, we show that the level of charge-induced acceleration noise on a single TM can be maintained at a level close to 1.0 fm/s^2/sqrt(Hz) across the 0.1-100 mHz frequency band that is crucial to an observatory such as LISA. Using dedicated measurements that detect these effects in the differential acceleration between the two test masses, we resolve the stochastic nature of the TM charge build up due to interplanetary cosmic rays and the TM charge-to-force coupling through stray electric fields in the sensor. All our measurements are in good agreement with predictions based on a relatively simple electrostatic model of the LISA Pathfinder instrument.

We used Planck data to study the M33 galaxy and find a substantial temperature asymmetry with respect to its minor axis projected onto the sky plane. This temperature asymmetry correlates well with the HI velocity field at 21 cm, at least within a galactocentric distance of 0.5 degree, and it is found to extend up to about 3 degrees from the galaxy center. We conclude that the revealed effect, that is, the temperature asymmetry and its extension, implies that we detected the differential rotation of the M33 galaxy and of its extended baryonic halo.

Future space borne gravitational wave detectors will require a precise definition of calibration signals to ensure the achievement of their design sensitivity. The careful design of the test signals plays a key role in the correct understanding and characterisation of these instruments. In that sense, methods achieving optimal experiment designs must be considered as complementary to the parameter estimation methods being used to determine the parameters describing the system. The relevance of experiment design is particularly significant for the LISA Pathfinder mission, which will spend most of its operation time performing experiments to characterise key technologies for future space borne gravitational wave observatories. Here we propose a framework to derive the optimal signals ---in terms of minimum parameter uncertainty--- to be injected to these instruments during its calibration phase. We compare our results with an alternative numerical algorithm which achieves an optimal input signal by iteratively improving an initial guess. We show agreement of both approaches when applied to the LISA Pathfinder case.

Planck data towards the galaxy M82 are analyzed in the 70, 100 and 143 GHz bands. A substantial north-south and East-West temperature asymmetry is found, extending up to 1 degree from the galactic center. Being almost frequency-independent, these temperature asymmetries are indicative of a Doppler-induced effect regarding the line-of-sight dynamics on the halo scale, the ejections from the galactic center and, possibly, even the tidal interaction with M81 galaxy. The temperature asymmetry thus acts as a model-independent tool to reveal the bulk dynamics in nearby edge-on spiral galaxies, like the Sunyaev-Zeldovich effect for clusters of galaxies.

Planck data towards the active galaxy Centaurus A are analyzed in the 70, 100 and 143 GHz bands. We find a temperature asymmetry of the northern radio lobe with respect to the southern one that clearly extends at least up to 5 degrees from the Cen A center and diminishes towards the outer regions of the lobes. That transparent parameter - the temperature asymmetry - thus has to carry a principal information, i.e. indication on the line-of-sight bulk motion of the lobes, while the increase of that asymmetry at smaller radii reveals the differential dynamics of the lobes as expected at ejections from the center.

The LISA Pathfinder mission will demonstrate the technology of drag-free test masses for use as inertial references in future space-based gravitational wave detectors. To accomplish this, the Pathfinder spacecraft will perform drag-free flight about a test mass while measuring the acceleration of this primary test mass relative to a second reference test mass. Because the reference test mass is contained within the same spacecraft, it is necessary to apply forces on it to maintain its position and attitude relative to the spacecraft. These forces are a potential source of acceleration noise in the LISA Pathfinder system that are not present in the full LISA configuration. While LISA Pathfinder has been designed to meet it's primary mission requirements in the presence of this noise, recent estimates suggest that the on-orbit performance may be limited by this `suspension noise'. The drift-mode or free-flight experiments provide an opportunity to mitigate this noise source and further characterize the underlying disturbances that are of interest to the designers of LISA-like instruments. This article provides a high-level overview of these experiments and the methods under development to analyze the resulting data.

LISA Pathfinder (LPF), ESA's precursor mission to a gravitational wave observatory, will measure the degree to which two test-masses can be put into free-fall, aiming to demonstrate a residual relative acceleration with a power spectral density (PSD) below 30 fm/s$^2$/Hz$^{1/2}$ around 1 mHz. In LPF data analysis, the measured relative acceleration data series must be fit to other various measured time series data. This fitting is required in different experiments, from system identification of the test mass and satellite dynamics to the subtraction of noise contributions from measured known disturbances. In all cases, the background noise, described by the PSD of the fit residuals, is expected to be coloured, requiring that we perform such fits in the frequency domain. This PSD is unknown {\it a priori}, and a high accuracy estimate of this residual acceleration noise is an essential output of our analysis. In this paper we present a fitting method based on Bayesian parameter estimation with an unknown frequency-dependent background noise. The method uses noise marginalisation in connection with averaged Welch's periodograms to achieve unbiased parameter estimation, together with a consistent, non-parametric estimate of the residual PSD. Additionally, we find that the method is equivalent to some implementations of iteratively re-weighted least-squares fitting. We have tested the method both on simulated data of known PSD, and to analyze differential acceleration from several experiments with the LISA Pathfinder end-to-end mission simulator.

Planck's data acquired during the first 15.4 months of observations towards both the disk and halo of the M31 galaxy are analyzed. We confirm the existence of a temperature asymmetry, previously detected by using the 7-year WMAP data, along the direction of the M31 rotation, therefore indicative of a Doppler-induced effect. The asymmetry extends up to about 10 degrees (about 130 kpc) from the M31 center. We also investigate the recent issue raised in Rubin and Loeb (2014) about the kinetic Sunyaev-Zeldovich effect from the diffuse hot gas in the Local Group, predicted to generate a hot spot of a few degrees size in the CMB maps in the direction of M31, where the free electron optical depth gets the maximum value. We also consider the issue whether in the opposite direction with respect to the M31 galaxy the same effect induces a minimum in temperature in the Planck's maps of the sky. We find that the Planck's data at 100 GHz show an effect even larger than that expected.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We report the results of a deep archive XMM-Newton observation of the Fornax spheroidal galaxy that we analyzed with the aim of fully characterizing the X-ray source population (in most of the cases likely to be background active galactic nuclei) detected towards the target. A cross correlation with the available databases allowed us to find a source that may be associated with a variable star belonging to the galaxy. We also searched for X-ray sources in the vicinity of the Fornax globular clusters GC 3 and GC 4 and found two sources probably associated with the respective clusters. The deep X-ray observation was also suitable for the search of the intermediate-mass black hole (of mass $\simeq 10^{4}$ M$_{\odot}$) expected to be hosted in the center of the galaxy. In the case of Fornax, this search is extremely difficult since the galaxy centroid of gravity is poorly constrained because of the large asymmetry observed in the optical surface brightness. Since we cannot firmly establish the existence of an X-ray counterpart of the putative black hole, we put constraints only on the accretion parameters. In particular, we found that the corresponding upper limit on the accretion efficiency, with respect to the Eddington luminosity, is as low as a few $10^{-5}$.

Baryons constitute about 4% of our universe, but most of them are missing and we do not know where and in what form they are hidden. This constitute the so-called missing baryon problem. A possibility is that part of these baryons are hidden in galactic halos. We show how the 7-year data obtained by the WMAP satellite may be used to trace the halo of the nearby giant spiral galaxy M31. We detect a temperature asymmetry in the M31 halo along the rotation direction up to about 120 kpc. This could be the first detection of a galactic halo in microwaves and may open a new way to probe hidden baryons in these relatively less studied galactic objects using high accuracy CMB measurements.

Data on the cosmic microwave background (CMB) radiation by the Wilkinson Microwave Anisotropy Probe (WMAP) had a profound impact on the understanding of a variety of physical processes in the early phases of the Universe and on the estimation of the cosmological parameters. Here, the 7-year WMAP data are used to trace the disk and the halo of the nearby giant spiral galaxy M31. We analyzed the temperature excess in three WMAP bands (W, V, and Q) by dividing the region of the sky around M31 into several concentric circular areas. We studied the robustness of the detected temperature excess by considering 500 random control fields in the real WMAP maps and simulating 500 sky maps from the best-fitted cosmological parameters. By comparing the obtained temperature contrast profiles with the real ones towards the M31 galaxy, we find that the temperature asymmetry in the M31 disk is fairly robust, while the effect in the halo is weaker. An asymmetry in the mean microwave temperature in the M31 disk along the direction of the M31 rotation is observed with a temperature contrast up to about 130 microK/pixel. We also find a temperature asymmetry in the M31 halo, which is much weaker than for the disk, up to a galactocentric distance of about 10 degrees (120 kpc) with a peak temperature contrast of about 40 microK/pixel. Although the confidence level of the signal is not high, if estimated purely statistically, which could be expected due to the weakness of the effect, the geometrical structure of the temperature asymmetry points towards a definite effect modulated by the rotation of the M31 halo. This result might open a new way to probe these relatively less studied galactic objects using high-accuracy CMB measurements, such as those with the Planck satellite or planned balloon-based experiments, which could prove or disprove our conclusions.

The cosmic microwave background (CMB) radiation is characterized by well-established scales, the 2.7 K temperature of the Planckian spectrum and the $10^{-5}$ amplitude of the temperature anisotropy. These features were instrumental in indicating the hot and equilibrium phases of the early history of the Universe and its large scale isotropy, respectively. We now reveal one more intrinsic scale in CMB properties. We introduce a method developed originally by Kolmogorov, that quantifies a degree of randomness (chaos) in a set of numbers, such as measurements of the CMB temperature in some region. Considering CMB as a composition of random and regular signals, we solve the inverse problem of recovering of their mutual fractions from the temperature sky maps. Deriving the empirical Kolmogorov's function in the Wilkinson Microwave Anisotropy Probe's maps, we obtain the fraction of the random signal to be about 20 per cent, i.e. the cosmological sky is a weakly random one. The paper is dedicated to the memory of Vladimir Arnold (1937-2010).

The Kolmogorov stochasticity parameter is shown to act as a tool to detect point sources in the cosmic microwave background (CMB) radiation temperature maps. Kolmogorov CMB map constructed for the WMAP's 7-year datasets reveals tiny structures which in part coincide with point radio and Fermi/LAT gamma-ray sources. In the first application of this method, we identified several sources not present in the then available 0FGL Fermi catalog. Subsequently they were confirmed in the more recent and more complete 1FGL catalog, thus strengthening the evidence for the power of this methodology.

The power spectrum is obtained for the Kolmogorov stochasticity parameter map for WMAP's cosmic microwave background (CMB) radiation temperature datasets. The interest for CMB Kolmogorov map is that it can carry direct information about voids in the matter distribution, so that the correlations in the distribution of voids have to be reflected in the power spectrum. Although limited by the angular resolution of the WMAP, this analysis shows the possibility of acquiring this crucial information via CMB maps. Even the already obtained behavior, some of which is absent in the simulated maps, can influence the development of views on the void correlations at the large-scale web formation.

The plane-mirror symmetry previously noticed in the Cosmic Microwave Background (CMB) temperature anisotropy maps of Wilkinson Microwave Anisotropy Probe is shown to possess certain anomalous properties. The degree of the randomness determined by the Kolmogorov stochasticity parameter in the both symmetry regions appears to have identical values which, however, essentially differ from the corresponding values for other sky regions. If the mirroring were of cosmological origin, this would imply either additional randomizing properties in those directions of the Universe or their different line-of-sight depth. This analysis also provides a way to test the hypothesis of a link between the nature of dark energy and inhomogeneities.

We continue investigation of the hidden plane-mirror symmetry in the distribution of excursion sets in cosmic microwave background (CMB) temperature anisotropy maps, previously noticed in the three-year data of the Wilkinson microwave anisotropy probe (WMAP), using the WMAP 5 years maps. The symmetry is shown to have higher significance, \chi^2 < 1.7, for low multipoles \ell < 5, while disappearing at larger multipoles, \chi^2 > 3.5 for \ell > 10. The study of the sum and difference maps of temperature inhomogeneity regions, along with simulated maps, confirms its existence.The properties of these mirroring symmetries are compatible with those produced by the Sachs-Wolfe effect in the presence of an anomalously large component of horizon-size density perturbations, independent of one of the spatial coordinates, and/or a slab-like spatial topology of the Universe.

The diffraction effects on gravitational waves propagating through a stellar cluster are analyzed in the relevant approximation of Fresnel diffraction limit. We find that a gravitational wave scintillation effect - similar to the radio source scintillation effect - comes out naturally, implying that the gravitational wave intensity changes in a characteristic way as the observer moves.