Next Article in Journal
Revisiting Vaidya Horizons
Next Article in Special Issue
A Toy Cosmology Using a Hubble-Scale Casimir Effect
Previous Article in Journal
Two-Body Orbit Expansion Due to Time-Dependent Relative Acceleration Rate of the Cosmological Scale Factor
Previous Article in Special Issue
Cosmographic Constraints and Cosmic Fluids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Large Scale Cosmological Anomalies and Inhomogeneous Dark Energy

by
Leandros Perivolaropoulos
Department of Physics, University of Ioannina, Ioannina 45110, Greece
Galaxies 2014, 2(1), 22-61; https://doi.org/10.3390/galaxies2010022
Submission received: 6 November 2013 / Revised: 24 December 2013 / Accepted: 25 December 2013 / Published: 17 January 2014
(This article belongs to the Special Issue Cosmology with Fluid Components)

Abstract

:
A wide range of large scale observations hint towards possible modifications on the standard cosmological model which is based on a homogeneous and isotropic universe with a small cosmological constant and matter. These observations, also known as “cosmic anomalies” include unexpected Cosmic Microwave Background perturbations on large angular scales, large dipolar peculiar velocity flows of galaxies (“bulk flows”), the measurement of inhomogenous values of the fine structure constant on cosmological scales (“alpha dipole”) and other effects. The presence of the observational anomalies could either be a large statistical fluctuation in the context of ΛCDM or it could indicate a non-trivial departure from the cosmological principle on Hubble scales. Such a departure is very much constrained by cosmological observations for matter. For dark energy however there are no significant observational constraints for Hubble scale inhomogeneities. In this brief review I discuss some of the theoretical models that can naturally lead to inhomogeneous dark energy, their observational constraints and their potential to explain the large scale cosmic anomalies.

1. Introduction

The standard cosmological model (ΛCDM) has been established through a wide range of cosmological observations [1] which challenged its validity during the past two decades. The model is based on the following assumptions: cosmological principle (isotropy, homogeneity, general relativity, cold dark matter and baryonic matter), flatness of space, the existence of a cosmological constant and gaussian scale invariant matter perturbations generated during inflation. Based on these assumptions, the model makes well defined predictions which are consistent with the vast majority of cosmological data. Some of these predictions include the following:
  • The Cosmic Microwave Background (CMB) Spectrum: The angular power spectrum of CMB primordial perturbations [2] is in good agreement with the predictions of ΛCDM. However, a few issues related to the orientation and magnitude of low multipole moments (CMB anomalies) constitute remaining puzzles for the standard model [3,4,5,6,7,8,9,10,11].
  • The statistics of CMB perturrbations: These statistics [12] are consistent with the prediction of gaussianity of the standard model [13].
  • Accelerating expansion: Cosmological observations using standard candles [Type Ia Supernovae (SnIa)] [14] and standard rulers (Baryon Acoustic Oscillations) [15] to map the recent accelerating expansion rate of the universe are consistent with the existence of a cosmological constant. No need has appeared for more complicated models based on dynamical dark energy or modified gravity, despite of the continuously improved data. The likelihood of the cosmological constant vs. more complicated homogeneous models has been continuously increasing during the past decade [16].
  • Large Scale Structure: Observations of large scale structure are in good agreement with ΛCDM [17] (basic statistics of galaxies [18], halo power spectrum [19]).
Despite of the above major successes the standard model is challenged by a few puzzling large scale cosmological observations [20,21] which may hint towards required modifications of the model. Some of these challenges of ΛCDM may be summarized as follows:
  • Power Asymmetry of CMB perturbation maps: A hemispherical power asymmetry in the cosmic microwave background (CMB) on various different angular scales (multipole ranges) has been detected [22,23,24,25]. The power in all multipole ranges is consistently found to be significantly higher in the approximate direction towards Galactic longitude and latitude ( l = 237 , b = - 20 ) than in the opposite direction. A more recent study of the WMAP9 data has found a hemispherical directional dependence of ΛCDM cosmological parameters which is maximized at the direction ( l = 227 , b = - 27 ) at the multipole range 2–600 and is statistically significant at the 3 . 4 σ level [26]. A related asymmetry, is the Maximum Temperature Asymmetry (MTA) defined as the maximized temperature difference between opposite pixels in the sky which shows significant alignment with other apparently unrelated asymmetries [27].
  • Large Scale Velocity Flows: Recent studies have indicated the existence of dipole velocity flows on scales of 100 h - 1 Mpc [28,29] with magnitude about 400 km/s using a combination of peculiar velocity surveys, with direction towards l 282 , b 6 . Other studies [30] using the kinematic S-Z effect have found bulk flows on much larger scales (O(1 Gpc)) with magnitude 600–1000 km/s towards a similar direction. These results are inconsistent with the predictions of ΛCDM at a level of 99 % . These studies however have been challenged by other authors which do not confirm these results but find peculiar velocities on these scales consistent with ΛCDM [31]. Even though these studies agree with the direction of the observed flow they disagree on the magnitude and errorbar of the measured velocities.
  • Alignment of low multipoles in the CMB angular power spectrum: The normals to the octopole and quadrupole planes are aligned with the direction of the cosmological dipole at a level inconsistent with Gaussian random, statistically isotropic skies at 99.7% [4]. This inconsistency has been reduced by the recent Planck results to a level of about 98% [22] (the exact level varies slightly depending on the foreground filtering method).
  • Large scale alignment in the QSO optical polarization data: Quasar polarization vectors are not randomly oriented over the sky with a probability often in excess of 99.9%. The alignment effect seems to be prominent along a particular axis in the direction ( l , b ) = ( 267 , 69 ) [32,33,34].
  • Anisotropy in Accelerating Expansion Rate: Recent studies of the accelerating cosmic expansion rate using SnIa as standard candles have indicated that an anisotropic expansion rate fit by a dipole provides a better fit to the data than an isotropic expansion rate at the 2σ level [35] Interestingly, this dark energy dipole is abnormally aligned [35] with the fine structure constant dipole discussed below.
  • Spatial dependence of the value of the fine structure constant α: A spatial cosmic variation of the fine structure constant has been recently identified on redshifts up to z 4 . 2 by analyzing the absorption spectra of quasars. This anisotropy analysis of the fine structure constant α [36,37] is based on a large sample of quasar absorption-line spectra (295 spectra) obtained using UVES (the Ultraviolet and Visual Echelle Spectrograph) on the VLT (Very Large Telescope) in Chile and also previous observations at the Keck Observatory in Hawaii. An apparent variation of α across the sky was found. It was shown to be well fit by an angular dipole model Δ α α = A cos θ + B where θ is the angle with respect to a preferred axis and A , B are the dipole magnitude and an isotropic monopole term. The dipole axis was found to point in the direction ( l , b ) = ( 331 , - 14 ) and the dipole amplitude A was found to be A = ( 0 . 97 ± 0 . 21 ) × 10 - 5 . The statistical significance over an isotropic model was found to be at the 4 . 1 σ level.
  • Large Quasar Group: An elongated structure of quasars with long dimension about 1240 Mpc and mean redshift z ¯ = 1 . 27 has recently been discovered [38]. This structure is a factor of about three larger than the previously known largest structure (Sloan Great Wall ( z ¯ = 0 . 073 and comoving size 450 Mpc [39]) and appears to be inconsistent with the cosmological principle in the context of the standard ΛCDM model in the sense that it is much larger than the scale of homogeneity in the context of ΛCDM (260–370 Mpc [40]).
The above seven large scale puzzles are large scale effects (1 Gpc or larger) and seem to be related to violation of the cosmological principle. In fact preferred cosmological directions may be associated with most of them. As discussed in more detail in Section 3, some of these directions appear to be abnormally close to each other. These directions are shown in Figure 1 [27]. This may be an indication for a common physical origin at least for some of the above anomalies.
These anomalies may be simply large statistical fluctuations that manifest themselves duo to a posteriori selection of particular statistical tests. In view of their diverse nature and their high statistical significance, is not likely that they are all of this nature. In view of the possibility that there is a physical mechanism that leads to their existence, it is of interest to identify classes of physical models that could predict the existence of some or all of the above cosmic anomalies.
Figure 1. Directions in galactic coordinates for the α (blue) and Dark Energy (green) dipoles, for the Dark Flow direction (red) and for the direction of Maximum Temperature Asymmetry (MTA) in the 7 years ILC CMB map degraded to N s i d e = 8 (yellow). The opposite corresponding directions are also shown (from [27]).
Figure 1. Directions in galactic coordinates for the α (blue) and Dark Energy (green) dipoles, for the Dark Flow direction (red) and for the direction of Maximum Temperature Asymmetry (MTA) in the 7 years ILC CMB map degraded to N s i d e = 8 (yellow). The opposite corresponding directions are also shown (from [27]).
Galaxies 02 00022 g001
The origin of the possible violation of the cosmological principle could be either geometric or energetic. A geometric origin would correspond to inherent properties of the cosmic metric (e.g., non-trivial topology on large scales). An energetic origin would correspond to non-trivial properies of the energy momentum tensor of the contents of the universe (e.g., anisotropic dark energy equation of state).
Models based on geometric violation of the cosmological principle include the following:
  • A multiply connected non-trivial cosmic topology [41,42].
  • A rotating universe coupled to an anisotropic scalar field [43].
  • Non-commutative geometry [44].
  • A fundamental anisotropic curvature [45].
Models based on energetic violation of the cosmological principle include:
  • Anisotropic dark energy equation of state [46,47,48] due perhaps to the existence of vector fields [49,50].
  • Dark matter perturbations on a few Gpc scale [51,52]. For example an off center observer in a 1 Gpc void would experience the existence of a preferred cosmological axis through the Lemaître-Tolman-Bondi metric [53,54,55,56,57]. Within this framework, an additional dark energy component is not needed to secure consistency with the cosmological data that indicate accelerating expansion. The basic idea is that the increased expansion rate occurs locally in space rather than at recent cosmological times, a fact that can be achieved by assuming a locally-reduced matter energy density [58,59]. Thus, the observer is placed close to the center of a giant void with dimensions of a few Gpc [60]. Even though this approach is free of dark energy, it is by no means free of fine tuning. Apart from the unnatural assumption of giant-size Gpc voids, which are very unlikely to be produced in any cosmology, these models require the observer to be placed within a very small volume at the center of the void (about 10 - 6 of the total volume of the void). A slightly off-center observer, however, will naturally experience a preferred cosmological direction (towards the center of the void), which may help to resolve some of the observational puzzles of ΛCDM discussed above.
    Such matter perturbations could be induced by statistically anisotropic or non-gaussian primordial perturbations [61,62,63,64]. For example, inflationary perturbations induced by vector fields [65,66,67,68]. Note however that inflationary models with vector fields usually suffer from instabilities due to the existence of ghosts [69,70,71].
  • Large scale electromagnetic fields [72]. For example, the existence of a large scale primordial magnetic field [73,74,75]. Evidence for such a magnetic field has recently been found in CMB maps [76].
  • Dark energy perturbations on scales comparable to the horizon. Even though the sound speed for dark energy is close to unity implying that it can not cluster on scales much smaller than the horizon, it can still produce observable effects due to clustering on horizon scales [77,78]. In addition the possible recent formation of topological defects with Hubble scale core (topological quintessence) could also behave as inhomogeneous dark energy [79].
Models based on matter underdensities (voids) on scales larger than 1 Gpc are severely constrained by observations [54,80,81,82,83] and do not reduce to ΛCDM for any value of their parameters. For radii as large as the horizon they reduce to homogeneous open CDM matter dominated universe and for small radii they reduce to flat CDM models without a cosmological constant. Both of these classes are ruled out by observations. Even for radii of order a few Gpc these models are faced by severe problems that can be summarized as follows:
  • There is no simple physical mechanism to generate such large voids [84].
  • As discussed above, off center observers experience a CMB dipole that is in excess of the observed one. Thus our location is constrained to practically coincide with the center of the void (fine tuning) to within about a fraction of a percent in radious [85,86].
  • The fit for the expansion rate as a function of redshift is worse than ΛCDM requires significant fine tuning of the void profile to obtain a comparable fit with ΛCDM [86].
  • These models predict significant peculiar velocities for distant galaxies within the void. Such velocities could have been detected in CMB maps through the kinematic S-Z effects. The fact that such large radial velocities are not observed imposes the most severe class of constraints in this class of models [82].
In contrast to matter void models, spherically inhomogeneous dark energy models are more natural and reduce to ΛCDM when the inhomogeneity radious exceeds the horizon scale. They constitute one of the two generic generalizations of ΛCDM. Instead of breaking time translation invariance of the cosmological constant (as is the case for quintessence models) inhomogeneous dark energy models break space translation invariance of Λ. This corresponds to addressing the coincidence problem of ΛCDM not as “why now?” question but as a “why here?” question. The motivation for considering this class of inhomogeneous dark energy models may be summarized as follows:
  • It is a new generic generalization of ΛCDM including ΛCDM as a special case.
  • It naturally violates the cosmological principle on large cosmic scales and predicts a preferred axis for off center observers. Thus it has the potential to address at least some of the above discussed cosmic anomalies.
  • There is a well defined physical mechanism that can give rise to this type on dark energy inhomogeneities. This mechanism will be discussed in some detail in Section 3. It is based on applying the principles of topological inflation [87] to the case of late-time acceleration. According to the idea of topological inflation, the false vacuum energy of the core of a topological defect can give rise to accelerating expansion if the core size reaches the Hubble scale when gravity starts dominating the dynamics. Thus, for example, a recently formed global monopole with appropriate scale of symmetry breaking and coupling could naturally produce a Hubble-scale, spherically symmetric, isocurvature dark energy overdensity. By analogy with topological inflation, this mechanism may be called topological quintessence.
In this review we focus on models which can violate the cosmological principle via energy momentum tensor effects (energetic origin violations). Since Hubble scale matter inhomogeneities are severely constrained we consider mainly dark energy Hubble scale inhomogeneities and/or anisotropies which appear to be a more natural source of a mild violation of the cosmological principle.
The structure of this review is the following: In the next section we present a collective analysis of the above discussed cosmic anomalies. After describing in some detail the current observational status of each one of these anomalies we discuss the statistical significance of its incompatibility with large scale homogeneity and isotropy. In Section 3 we describe theoretical models that predict the existence of anisotropic or inhomogeneous dark energy and the potential of these models to explain some of the cosmic anomalies. We discuss two classes of models: Anisotropic dark energy models described by a Bianchi I cosmic metric with anisotropic dark energy equation of state and spherical inhomogeneous dark energy models described by a generalization of the Lemaître-Tolman-Bondi metric. A special case of the later class of models is “topological quintessence” where a well defined physical mechanism generates the spherical dark energy inhomogeneity. Finally in Section 4 we conclude and summarize our results. We also present an outlook discussing remaining open questions and possible future directions of research related to this subject.

2. A Review of Cosmic Anomalies

In this section we discuss in some detail the cosmological observations that hint towards a possible violation of the large scale homogeneity and isotropy of the universe.

2.1. Power Asymmetry of CMB Perturbation Maps

As discussed in the Introduction, the power spectrum of CMB fluctuations as measured by Wilkinson Microwave Anisotropy Probe (WMAP) and Planck is not isotropic. A hemispherical asymmetry has been pointed out in several studies [2,22,26,71,88,89,90] using both the WMAP [2] and the Planck data [22] at several multipole ranges. The persistense of this asymmetry in different multipole ranges along with the good multifrequency component separation (especially in the Planck data) make it unlikely that the asymmetry is due to foregrounds. The observed asymmetry is modeled as a power spectrum dipolar modulation of an otherwise statistically isotropic sky:
Δ T ( n ^ ) = ( 1 + A p ^ · n ^ ) Δ T i s o ( n ^ )
where p ^ and n ^ are unit vectors towards the dipole and the sky directions respectively and Δ T i s o ( n ^ ) is a statistically isotropic temperature fluctuation in the direction n ^ . The best fit dipole is found to have direction ( l , b ) = ( 227 , - 27 ) and amplitude A = 0 . 072 ± 0 . 022 for the multipole range l < 64 ( k < 0.035 Mpc - 1 ) [22]. The derived amplitude is inconsistent with isotropy ( A = 0 ) at the 3 . 4 σ level. If this dipole modulation of the angular power spectrum is not a statistical fluctuation and has a physical cause, then this cause should be such as to change the amplitude of density fluctuations by less than O ( 10 - 3 ) as constrained by the CMB dipole magnitude. This is not easy to achieve for a primordial mechanism since the same mechanism responsible for an O ( 10 - 1 ) (as seen by the value of A) modulation of the primordial density perturbation spectrum should practically leave unaffected the magnitude of these density perturbations. This however would be easier to achieve for a late time mechanism operating after the density perturbations have been created.
On smaller scales ( k 1 Mpc - 1 there is lack of asymmetry in primordial density perturbations as constrained by the SDSS quasar sample [91]. The power spectrum of these quasars is found to be isotropic with A < 0 . 0153 (at 3σ) for a dipole oriented in the direction of the CMB power asymmetry dipole.
Thus, the possible mechanism that leads to the observed CMB power asymmetry at the O ( 10 - 1 ) level should affect the magnitude of the density perturbations at a level less than 10 - 3 and should also not create power asymmetry on scales k 1 Mpc - 1 at a level more than about O ( 10 - 2 ) .
Both of these constraints of this physical mechanism are consistent with a late generation of asymmetry through through large scale dark energy inhomogeneities which affect locally the growth of structures and thus the CMB maps through the ISW effect. On the other hand, it is hard to construct primordial models of inflation that satisfy the above two constraints [92,93].
It has been shown recently [94] that a large part of the asymmetry is due to the ISW effect. In the context of inhomogeneous dark energy, regions of denser dark energy would be associated more pronounced ISW effect. Thus there are two hints that point out to the possible late time origin of the CMB asymmetry:
  • A large part (if not all) of the asymmetry appears to be due to the ISW effect which occurs at late times and is not related to the primordial nature of the CMB perturbations [94].
  • The power spectrum of large scale structure does not show evidence for such an dipole asymmetry on smaller scales. This also hints towards a possible late time origin of the asymmetry.
Dark energy inhomogeneities is clearly a prime candidate that could lead to the generation of late time asymmetry without affecting the spectrum of density perturbations.

2.2. Alignment of Low Multipoles in the CMB Angular Power Spectrum

The alignment of the normals to the octopole and quadrupole planes has been verified by several studies [3,6,7,95,96,97,98]. The majority of these studies conclude that this is not a foreground contamination of the WMAP Internal Linear Combination (ILC) maps but rather is an effect degraded by the foregrounds [99,100] (see however [101] for an alternative point of view). The directions of the quadrupole and the octopole are determined by rotating the coordinate system while determining the multipole coefficients a l m . The coordinate system that maximized the sum | a l l | 2 + | a l - l | 2 is selected and the direction of its z-axis defines the direction of the l multipole moment ( l = 2 for the quadrupole and l = 3 for the octopole). This coordinate system maximizes the power of fluctuations at its x - y plane (equator). The exact direction of these moments varies slightly in accordance with the method used to filter the foregrounds.
According to the recent Plack data [22] and one of the filtering methods used, the octopole plane normal is in the direction ( l , b ) = ( 246 , 66 ) , and the quadrupole plane normal is in the direction ( l , b ) = ( 228 , 60 ) . For comparison the CMB dipole moment is in the direction ( l , b ) = ( 264 , 48 ) [102].
The angle between the quadrupole and the octopole directions is found to in the range 9 13 depending on the foreground cleaning method used [22]. The probability that this occurs by chance in a random simulation of the best-fit ΛCDM model is found to be about 2 % . The alignment however gets significantly degraded by subtracting any large scale specific feature (hot or cold spot) on the CMB map. Thus, it appears the there is no specific CMB feature responsible for the alignment. Instead, the alignment appears to be an unlikely statistical fluctuation of amplitudes and phases which could be justified by the a posteriori selection of the particular statistic [5]. Despite the degraded statistical significance of this anomaly by the Planck data it remains an unlikely feature in the context of ΛCDM hinting towards statistical anisotropy. A well defined theoretical model that could naturally justify its existence would have to naturally produce a combination of large scale features in a single CMB map plane leading to the observed alignment.

2.3. Large Scale Velocity Flows

The kinematic Sunyaev-Zeldovich (kSZ) effect is a spectral distrortion on the CMB photons induced by the coherent velocities of the free electrons in cluster due to the peculiar velocity of the cluster as a whole. A bulk flow of clusters would produce spectral distortions at the locations of clusters whose anisotropy would be well described by a dipole. This method was first applied for the measurement of large scale bulk flows of clusters in [103]. It has the advantage that it is independent of other conventional methods to measure peculiar velocities which require the measurement of distances. In addition it can easily probe scales up to the horizon which is currently impossible with conventional methods. In general however it is not easy to distinguish the small signal of the kSZ dipole from the much larger signal of the intrinsic CMB dipole. Despite of these difficulties the authors of [30,104,105] used a sample of 700 X-Ray selected clusters with redshifts up to z 0 . 2 (scale of about 800 Mpc) and detected a persistent kSZ dipole which they interpreted as evidence for bulk flow of clusters with magnitude 600–1000 km/s in a direction l 296 ± 28 , b 39 ± 14 .
Other studies however have not confirmed this result [106,107,108,109,110] using somewhat different methodologies and datasets. Most of these studies [106,108,110] have obtained a similar dipole but have claimed significantly larger errors which make the derived magnitudes of bulk flows statistically insignificant and consistent with null (and with ΛCDM which predicts much smaller coherent peculiar velocities on these scales). The issue remains controversial [111] and there is currently a debate as to the identification of the proper methodology for the implementation of the kSZ effect to derive the peculiar velocities of clusters.
On smaller scales (50–100 h - 1 Mpc), an application of a Minimal Variance weighting algorithm to a compilation of 4481 peculiar velocity measurements [28] lead to a bulk flow of 407 ± 81 km/s towards l 287 ± 9 , b 8 ± 6 . This result is inconsistent (too large) with respect to ΛCDM normalized on WMAP7 at a level of 98 % while the direction of the flow is consistent with the result of [104]. A more recent analysis [112] using the peculiar velocities of a sample of 245 SnIa has obtained a bulk flow of a lower magnitude ( 248 ± 87 km/s) towards a similar direction l 319 ± 25 , b 7 ± 13 . This magnitude of the bulk flow is marginally consistent with both ΛCDM and with the previous analysis of [28] and its direction is even closer to other anisotropy directions discussed below (dark energy dipole, α dipole and CMB cold spot I [27]). Consistent results with ΛCDM are also obtained by [113] using SnIa peculiar velocities to measure bulk flows.

2.4. Large Scale Alignment in the QSO Optical Polarization Data

The optical polarization properties of 355 quasars with redshifts up to z 2 . 5 and with well defined polarization angles were studied in [33] (see also [32]). The cosmic scales covered are up to about 2 Gpc. It was found that the quasar polarization vectors are not randomly distributed in the sky at a confidence level of 99 . 9 % . Sources of systematics like instrumental and interstellar polarizations were carefully considered and shown not to affect the main conclusions of the analysis. In addition, the polarizations of quasars obtained by different observatories were found to agree within the uncertainties in both polarization degree and angle. The observed large scale angular correlations of quasar polarization in regions sized at least 1 Gpc showed a coherent change of the mean polarization angle with redshift. The effect appeared stronger along an axis close to the CMB dipole even though direction uncertainties are large. If the effect is not due to an unknown systematic or due to a large statistical fluctuation, it may be evidence for departure from the cosmological principle (large scale isotropy).
An interesting physical mechanism [34] that could lead to the observed coherent rotation of quasar polarizations with redshift involves a coupling between the photon and a pseudoscalar (e.g., axion) in a Lagrangian of the form:
L = 1 2 ( μ ϕ ) ( μ ϕ ) - 1 2 m 2 ϕ 2 - 1 4 F μ ν F μ ν - j μ A μ + 1 4 g ϕ F μ ν F ˜ μ ν
where F ˜ μ ν is the dual Maxwell tensor and 1 4 g ϕ F μ ν F ˜ μ ν describes the mixing pseudoscalar axions and photons. In the context of such a Lagrangian, photons with polarization parallel to a weak external magnetic field decay into pseudoscalars leading to a net polarization p. Thus coherent variation of the polarization p over cosmological distances are predicted provided that there is a magnetic field B < 1 n G [114].
If there is a variation of the direction of B with distance then a rotation of the polarization angle with redshift is predicted. In an extended model where the pseudovector ϕ also plays the role of quintessence with induced inhomogeneities a correlation between accelerating expansion and quasar polarization would be anticipated.
Recent data showing a similar polarization alignment of radiowaves [115] on scales up to 1 Gpc however appear to impose constraints on the pseudoscalar coupling mechanism [116]. Similar strong constraints on this mechanism are imposed by the observed lack of circular polarization in the quasar data [117].

2.5. Anisotropy in Accelerating Expansion Rate (Dark Energy Dipole)

Several recent studies [35,118,119,120,121,122,123,124,125,126,127,128,129] have investigated the isotropy of the accelerating cosmological expansion rate at various redshift ranges using samples of SnIa as standard candles. Most studies search for a preferred direction using either a hemisphere comparison method [118,120,121,123,129] or a fit to a dipole angular distribution of the anisotropy [35,122,125,128] (dark energy dipole). There are indications [35] that the dipole fit may be more sensitive in identifying the expansion anisotropy signal than the hemisphere comparison method. The consensus of these studies is that there is an asymmetry axis pointing in the direction ( l 310 , b = - 15 ) which is roughly consistent with the direction of the CMB Cold Spot I [5]. This asymmetry is more prominent for low redshift SnIa ( z 0 . 2 ) perhaps due to the fact that there are significantly more and better quality SnIa data for this redshift range [35]. The magnitude of the asymmetry is marginally consistent with ΛCDM (at about 2σ level). The direction of this asymmetry axis appears to be consistent at 1σ with both the peculiar velocity flows direction and with the alpha dipole direction discussed below (Figure 2).
The approximate coincidence of these four preferred directions (dark flow, dark energy dipole, α dipole and CMB cold spot) may hint towards a possible common physical origin for these anomalies.
Figure 2. The distribution of the Union2 SnIa datapoints in galactic coordinates along with the Dark Energy dipole direction are shown. The datapoints are split in three different redshift bins shown with different shapes. The direction of the α-dipole is also shown with a star. The light blue blob represents the 1-σ error on the Dark Energy dipole direction (from [35]).
Figure 2. The distribution of the Union2 SnIa datapoints in galactic coordinates along with the Dark Energy dipole direction are shown. The datapoints are split in three different redshift bins shown with different shapes. The direction of the α-dipole is also shown with a star. The light blue blob represents the 1-σ error on the Dark Energy dipole direction (from [35]).
Galaxies 02 00022 g002

2.6. Anisotropy in the Values of the Fine Structure Constant α (α Dipole)

The absorption spectra of distant quasars contain two pieces of information: the redshift of the absorber and the value of the fine structure constant α. This information extracted by the Many Multiplet Method has been used to perform an anisotropy analysis of the fine structure constant α [36,37]. The analysis is based on a large sample of quasar absorption-line spectra (295 spectra) obtained using UVES (the Ultraviolet and Visual Echelle Spectrograph) on the VLT (Very Large Telescope) in Chile and also previous observations at the Keck Observatory in Hawaii.
An apparent variation of α across the sky was found. It was shown to be well fit by an angular dipole model Δ α α = A cos θ + B where θ is the angle with respect to a preferred axis and A , B are the dipole magnitude and an isotropic monopole term. The dipole axis was found to point in the direction ( l , b ) = ( 331 , - 14 ) and the dipole amplitude A was found to be A = ( 0 . 97 ± 0 . 21 ) × 10 - 5 . The statistical significance over an isotropic model is at the 4 . 1 σ level. The analysis of [36,37] has received criticism [130,131,132] based mainly on the fact that its quasar sample combines two datasets (Keck and VLT) with different systematic errors which have a small overlapping subset and cover opposite hemispheres on the sky. The axis connecting these two hemispheres has similar direction with the direction of the obtained dipole. The response of the authors of [36,37] was based on the fact that in the equatorial region of the dipole, where both the Keck and VLT samples contribute a number of absorbers, there is no evidence for inconsistency between Keck and VLT.
The angular separation between the Dark Energy Dipole direction and the α dipole is found to be only 11 . 3 ± 11 . 8 . Reference [35] used Monte Carlo simulations to find the probability of obtaining the observed dipole magnitudes with the observed alignment, in the context of an isotropic cosmological model with no correlation between dark energy and fine structure constant α. It was found that this probability is less than one part in 10 6 .
A simple physical model (extended topological quintessence [35,79]) naturally predicts a spherical inhomogeneous distribution for both dark energy density and fine structure constant values. The model is based on the existence of a recently formed giant global monopole with Hubble scale core which also couples non-minimally to electromagnetism. Aligned dipole anisotropies would naturally emerge for an off-centre observer for both the fine structure constant and for dark energy density. This model smoothly reduces to ΛCDM for proper limits of its parameters. Two predictions of this model are (a) a correlation between the existence of strong cosmic electromagnetic fields and the value of α and (b) the existence of a dark flow on Hubble scales due to the repulsive gravity of the global defect core (“Great Repulser”) aligned with the dark energy and α dipoles. The direction of the dark flow is predicted to be towards the spatial region of lower accelerating expansion. Existing data about the dark flow [28,30] are consistent with this prediction. In the next section we present a detailed analysis of this model.

2.7. Large Quasar Groups

According to SDSS DR7 for galaxies with 0 . 22 < z < 0 . 50 [133], isotropy appears to be established on scales larger than 210 Mpc as anticipated from the overall isotropy of the CMB. However, on larger redshifts, isotropy appears to be challenged by recent quasar data.
The Large Quasar Groups (LQG) are large structures of quasars appearing at redshifts z 1 with scales 70–350 Mpc. A new elongated LQG was discovered recently [38] in the DR7QSO catalogue of the Sloan Digital Sky Survey [134]. It has long dimension 1240 Mpc and characteristic size ( V o l u m e ) 1 / 3 about 500 Mpc at a redshift z = 1 . 27 (Huge-LQG). It is located in the direction l = 225 , b = 63 , relatively close (distance of boundaries 140 Mpc) to another fairly large LQG and consists of 73 member quasars. The principal axes of the LQG have lengths of 1240, 640, and 370 Mpc extending to the Gpc-scale. According to [40], the upper limit to the scale of homogeneity for the concordance cosmology is 370 Mpc. Their long dimensions from the inertia tensor of 630 and 780 Mpc are clearly much larger than the anticipated upper limit of homogeneity. It is the currenly known largest structure in the universe and maybe at odds with the cosmological principle [40]. We stress however that other more recent studies have indicated [135,136] that such structures may not be inconsistent with a homogeneous Poisson distribution of quasars. Indeed, such structures seem to be discovered [135,136] by the algorithm of [38] in explicitly homogeneous simulations of a Poisson point process with the same density as the quasar catalogue.

3. Anisotropic/Inhomogeneous Dark Energy Models

In view of the successes of the standard ΛCDM model it is anticipated that a viable theoretical model able to explain the above described anomalies should have the following properties:
  • Deviate from ΛCDM mainly on large cosmological scales.
  • Reduce to ΛCDM for certain values of its parameters.
In addition, it appears reasonable and in fact hinted by some of the above described data that the physical origin of the deviation from isotropy resides at late cosmological times when these large scales are again causally connected (consider, e.g., the possible ISW origin of the CMB power asymmetry [94]). Motivated by these arguments and by the discussion presented in the Introduction we focus on models involving anisotropic and/or inhomogeneous dark energy. In contrast to matter, dark energy can be anisotropic (due to its pressure) without being inhomogeneous and also the possible inhomogeneities of dark energy naturally appear on Hubble scales where the origin of the anomalies seems to reside. In addition, Gpc inhomogeneities (voids) of matter constitute models that do not smoothly reduce to ΛCDM for any value of their parameters. This class of models has been shown to be severely constrained by several cosmological observations and especially the kSZ effect which is predicted to be significantly larger in these models than observations indicate [82]. Thus in what follows we focus on models of anisotropic/inhomogeneosu dark energy.

3.1. Anisotropic Dark Energy Models

A generic generalization of the Friedmann-Lemaître-Robertson-Walker (FLRW) metric that introduces anisotropy but not inhomogeneity is the Bianchi I metric which is of the form:
d s 2 = - d t 2 + a 2 ( t ) d x 2 + b 2 ( t ) d y 2 + c 2 ( t ) d z 2
Such a metric is consistent with an energy momentum tensor of a fluid with homogeneous but anisotropic pressure of the form:
T ν μ = d i a g ( - 1 , w x , w y , w z ) ρ
where ( w x , w y , w z ) are the equation of state componenets corresponding to each component of pressure. The deviation of these components may be parametrised by two parameters [137] γ and δ. Using these parameters, the anisotropic dark energy energy momentum tensor Equation (4) may be expressed as:
T μ ν μ = d i a g ( - 1 , w , w + 3 δ , w + 3 γ ) ρ
The deviations from isotropy may be quantified by defining the average Hubble expansion rate H and the Hubble normalized shear R and S as:
H 1 3 a ˙ a + b ˙ b + c ˙ c
and
R 1 H a ˙ a - b ˙ b , S 1 H a ˙ a - c ˙ c
The generalized Friedmann equation is then obtained as:
H 2 = 8 π G 3 ρ + ρ m 1 - 1 9 R 2 + S 2 - R S
where ρ m is the matter density which in general is allowed to interact with dark energy. Conservation of the energy momentum tensor of anisotropic dark energy implies:
ρ ˙ + 3 1 + w H + δ b ˙ b + γ c ˙ c ρ = - Q H ρ m
where ρ m is the matter density which obeys:
ρ ˙ m + 3 H ρ m = Q H ρ
and a nonzero value of the parameter Q expresses the interaction between dark energy and matter. Such an anisotropic dark energe could emerge for example in the presence of a vector field A μ whose dynamics is described by the action [137]:
S = d 4 x - g R 16 π G - 1 4 F μ ν F μ ν - V
where F μ ν = μ A ν - ν A μ and the potential V is a function of A 2 = A μ A μ [137].
The main observational effect of anisotropic dark energy is its contribution to the CMB anissotropies at the level of low multipole moments. This contribution is introduced through the last scattering redshift which becomes anisotropic in the context of the metric Equation (3). For a photon observed from the a direction p ^ the last scattering redshift z may be obtained using the lightlike geodesic equatios. It is of the form:
1 + z ( p ^ ) = 1 a 1 + p ^ y 2 e y 2 + p ^ z 2 e z 2
where:
e y 2 = a b 2 - 1 , e z 2 = a c 2 - 1
are evaluated at the time of last scattering and the scale factors are all normalized to unity today. Thus, the predicted CMB temperature is of the form:
T ( p ^ ) = T * / ( 1 + z ( p ^ ) )
where T * is the temperature at decoupling. Using the spatial average 4 π T ¯ = d Ω p ^ T ( p ^ ) , the temperature anisotropies is expressed as:
δ T ( p ^ ) T ^ = 1 - T ( p ^ ) T ¯
This anisotropy may be expanded in spherical harmonics to evaluate the contribution to all CMB multipole moments and thus impose constraints on the parameters γ and δ by demanding that these multipoles do not exceed the observed values. Significant fine tuning is required on the anisotropy parameters to secure consistency with the observed low quadrupole value [137,138,139].
Less stringent constraints on the anisotropy parameters may be obtained by considering the isotropy of the accelerating cosmic cosmic expansion rate as determined using the luminosity distances of SnIa. The predicted values of luminosity distances in the context of the metric Equation (3) is obtained [137] at the redshift z in the direction p ^ as:
d L ( z , p ^ ) = ( 1 + z ) t 0 t ( z ) d t p ^ x 2 a 2 + p ^ y 2 b 2 + p ^ z 2 c 2
where z is connected with the scale factors via Equation (12). Using the maximum likelihood method and the observed luminosity distances it is possible to obtain constraints on the anisotropy parameters γ and δ and show that their absolute values are constrained to be less than about 0 . 1 at the 3σ level [137].
An alternative parametrization of the dark energy anisotropy involves the consideration the anisotropic stress σ. This is obtained by expressing the dark energy momentum tensor as [138]:
T μ ν = ρ u μ u ν + p ( g μ ν + u μ u ν ) + Σ μ ν
where where u μ is the four-velocity of the fluid, and Σ μ ν is the anisotropic stress tensor. The anisotropic stress σ is defined as:
( ρ + p ) σ = - ( k ^ i k ^ j - 1 3 δ i j ) Σ i j
The anisotropic stress σ quantifies the deviation of pressure from isotropy. As discussed in [137], only weak constraints can be imposed on σ using recent data.
The introduction of anisotropic dark energy as described above provides the potential for breaking of statistical isotropy and resolution of some of the large scale anomalies discussed in the previous section. In particular, the anisotropic pressure can induce peculiar velocity flows, anisotropy in the SnIa data and significant CMB dipole and quadrupole [137]. The most severe constraints on the level of anisotropy originate from the CMB observed isotropy but these can be softened by considering a time dependent anisotropy which is more significant at late times and consistent with the weaker constraints originating from the SnIa data.

3.2. Inhomogeneous Dark Energy Models

In the context of inhomogeneous cosmologies, the simplest way to introduce a preferred direction in the universe is to consider a Hubble scale spherical dark matter void or dark energy inhomogeneity and an off center observer. In the case of dark matter void, an observer close to the center of the void will detect a larger Hubble parameter for nearby points in space and lower redshift. Thus she may incorrectly interprete it as an accelerating expansion rate since the Hubble expansion will appear higher at lower redshifts. Thus a spatial variation of expansion rate may be misinterpreted as a temporal variation. In the case of a spherical dark energy overdensity the variation of the expansion rate is both spatial and temporal and there is a true accelerating expansion rate amplified by the spatial variation of the dark energy density.
The displacement direction will be naturally preferred and the observer will experience anisotropic accelerated Hubble expansion and Hubble scale velocity flows (Figure 3). In fact the center of the dark matter void or dark energy overdensity will behave as a great repulser inducing velocity flows away from it.
Figure 3. An off-center observer in a spherical inhomogeneity sees a preferred direction of higher accelerating expansion rate.
Figure 3. An off-center observer in a spherical inhomogeneity sees a preferred direction of higher accelerating expansion rate.
Galaxies 02 00022 g003
As a simple toy model we may consider an inhomogeneous fluid of matter and/or dark energy described by the energy-momentum tensor:
T ν μ = diag ( - ρ ( r ) - ρ M ( r , t ) , p r ( r ) , p t ( r ) , p t ( r ) )
The general spherically symmetric inhomogeneous cosmological spacetime is described by the metric [140,141,142]:
d s 2 = - d t 2 + A 2 ( r , t ) d r 2 + B 2 ( r , t ) d θ 2 + sin 2 θ d φ 2
In the special case of a diagonal energy momentum tensor as the one of Equation (19), the above metric takes the simplified form [140,141,142]:
d s 2 = - d t 2 + ( B ( r , t ) ) 2 1 - k ( r ) d r 2 + B 2 ( r , t ) d θ 2 + sin 2 θ d φ 2
where B ( r , t ) is the spacetime dependent scale factor and k ( r ) plays the role of inhomogeneous spherically symmetric spatial curvature. In the limit of a homogeneous spacetime, the LTB metric Equation (21) reduces to the FLRW metric by setting k ( r ) to a constant and B ( r , t ) = r a ( t ) . Using the metric Equation (21) and the energy momentum tensor Equation (19) in the Einstein equations we obtain the generalized Friedman equation:
H 2 ( r , t ) = H 0 2 ( r ) Ω M ( r ) B 0 B 3 + Ω X ( r ) + Ω c ( r ) B 0 B 2
where:
H ( r , t ) B ˙ ( r , t ) B ( r , t )
with B 0 ( r ) B ( r , t 0 ) ( t 0 is the present time), H 0 ( r ) H ( r , t 0 ) , Ω X ( r ) - 8 π G p r ( r ) / 3 H 0 2 ( r ) (which is positive for negative radial pressure) and Ω c ( r ) 1 - Ω X ( r ) - Ω M ( r ) . The profile of the inhomogeneous dark energy pressure, Ω X ( r ) , is in principle arbitrary, and is physically determined by the physical mechanism that induced the inhomogeneous overdensity.
Interesting density profiles of the inhomogeneities (dark matter and dark energy) which parametrize the inside and outside densities ( Ω in and Ω out ) as well as the interface thickness ( Δ r ) are of the form:
Ω M ( r ) = Ω M , out + ( Ω M , in - Ω M , out ) 1 - tanh [ ( r - r 0 ) / 2 Δ r ] 1 + tanh ( r 0 / 2 Δ r )
Ω X ( r ) = Ω X , out + ( Ω X , in - Ω X , out ) 1 - tanh [ ( r - r 0 ) / 2 Δ r ] 1 + tanh ( r 0 / 2 Δ r )
It is straightforward to solve Equation (22) by first determining H 0 ( r ) in units of Gyr - 1 , assuming an age of the universe t 0 = 13 . 7 Gyr (homogeneous Big Bang). We thus use:
H 0 ( r ) = 1 t 0 0 1 d x Ω M ( r ) x - 1 + Ω X ( r ) x 2 + Ω c ( r )
which can be easily derived from Equation (22). With Equation (26) equation and the initial condition B ( r , t 0 ) = r we solve Equation (22) and obtain numerically B ( r , t ) and its derivatives.
For an observer in the center of the spherical inhomogeneity, light travels radially inwards and therefore for these lightlike geodesics d θ = d ϕ = 0 along the geodesic. The remaining two light-like geodesics for radially incoming light rays are of the form [143]:
d t d z = - B ( r , t ) ( 1 + z ) B ˙ ( r , t )
d r d z = c 1 - k ( r ) ( 1 + z ) B ˙ ( r , t )
where c 0 . 3 is the velocity on light in units of Gpc/Gyr used for consistency in order to obtain the luminosity distance in Gpc. Using these lightlike geodesics with initial conditions t ( z = 0 ) = t 0 and r ( z = 0 ) = 0 we find t ( z ) and r ( z ) . The angular diameter and luminosity distances are then obtained as:
d A ( z , r 0 , Ω X , in ) = B ( r ( z ) , t ( z ) )
d L ( z , r 0 , Ω X , in ) = ( 1 + z ) 2 d A ( z )
Using the above form of the luminosity distance, it is straightforward to fit the predicted distance moduli to the observed ones from the Union2 dataset and derive the best-fit values of parameters for the on-center observer. We first consider a dark matter inhomogeneity (void) with no dark energy. Fixing the interface thickness to Δ r = 0 . 35 Mpc, Ω X , in = Ω X , out = 0 (no dark energy), Ω M , out = 1 , it is straightforward to obtain [86] the 1σ and 2σ contours for the parameters Ω M , in (matter underdensity inside the void) and r 0 (size of the void) shown in Figure 4.
Figure 4. The 1σ and 2σ contours for the parameters Ω M , in (matter underdensity inside the void) and r 0 (size of the void) (from [86]).
Figure 4. The 1σ and 2σ contours for the parameters Ω M , in (matter underdensity inside the void) and r 0 (size of the void) (from [86]).
Galaxies 02 00022 g004
This class of Gpc void models had until recently attracted significant attention. They have three main advantages:
  • They are consistent with the observed redshift dependence of the Hubble expansion rate as obtained using standard candles SnIa [143,144] without the requirement of dark energy.
  • They can naturally provide a preferred axis by utilizing a slight displacement of the observer from the center of the spherical void.
  • The coincidence problem is not present since a void can naturally develop at late times along with other structures.
Despite of the above successes, it has been evident that this class of models also suffers from significant problems. As discussed in the Introduction, these problems include the following:
  • In the context of standard cosmology there is no simple mechanism to create Gpc matter voids without conflicting other cosmological data (e.g., CMB fluctuations).
  • If the observer location is not fine tuned close to the center of the spherical void, then a CMB dipole is predicted which is larger than the observed [85,86].
  • There is no simple matter density profile that will provide an equal or better quality of fit compared to ΛCDM [86].
  • This class of models generically predicts large peculiar velocities of free electrons and clusters with respect to the CMB. These velocities could have been detected via inverse Compton scattering of CMB photons onto moving free electrons (kinematic Synyaev-Zeldovich effect). The existing bounds on such velocities effectively rules out this class of models [82,83].
In view of the above problems of matter void models (see also [54]), it is of some interest to consider models where the Gpc void is replaced by a dark energy spherical overdensity with similar size. As discussed in the Introduction, the motivation for this class of inhomogeneous dark energy models may be summarized as follows:
  • The standard cosmological model ΛCDM assumes the existence of homogeneous matter on large scales, the validity of general relativity and the existence of homogeneous dark energy with constant in time energy density (cosmological constant). The generalization of most of these assumptions has been extensively studied in the literature. The main motivation for such generalizations originates at the coincidence problem expressed as a “Why now?” problem: “Why was the energy scale of the cosmological constant tuned to such low values so that it starts dominating at the present time?”. For example general relativity has been generalized to various models of modified gravity, dark energy was allowed to evolve in time as a scalar field (quintessence), matter was allowed to be inhomogeneous on Gpc scales (void models). The least studied generalization is the one based on allowing dark energy to become inhomogeneous on Gpc scales. Such a generalization could in principle address the coincidence problem if the later is expressed as a “Why here?” problem: “Why is the dark energy density in our horizon such that the universe has recently started its accelerating expansion?”. In the same manner that the answer to the “Why now?” question could be time dependence, the answer to the “Why here?” question could be spatial dependence of dark energy density. Therefore, the consideration of inhomogeneous dark energy constitutes a generic generalization of the standard ΛCDM models which reduces back to ΛCDM for a inhomogeneity scale that exceeds the current horizon scale. It is therefore interesting to use cosmological observations to impose constraints on the two basic parameters of this class of models: the scale of the inhomogeneity and the magnitude of the dark energy density inhomogeneity.
  • As in the case of the dark matter void, an off-center observer in a spherically symmetric dark energy inhomogeneity will naturally experience an anisotropy in the accelerating expansion rate. In the case of inhomogeneous dark energy the scale of the inhomogeneity is allowed to be as large as the horizon scale (ΛCDM limit). For large enough inhomogeneity scale, the observer is allowed to be displaced significantly from the center without experiencing a large CMB dipole. Thus, inhomogeneous dark energy models are less subject to fine tuning with respect to the location of the observer [86].
  • Due to negative pressure and speed of sound close to unity, dark energy inhomogeneities are not easy to sustain. However, there is a well defined physical mechanism (topological quintessence [79]) that can lead to sustainable large scale dark energy inhomogeneities, supported by topological considerations. This mechanism can be viewed as a generalization of topological inflation [87] in which early time accelerating expansion (inflation) takes place in the core of a topological defect due to the topologically trapped vacuum energy.
In this brief review I focus on spherically symmetric dark energy inhomogeneities corresponding to global monopoles in the context of topological quintessence. It is straightforward to generalize this assumption by considering other models that are both inhomogeneous and anisotropic. For example the Szekeres quasi-spherical models [145,146], are effectively non-linear FLRW perturbations [147], with metric of the form:
d s 2 = - d t 2 + e 2 A d x 2 + e 2 B ( d y 2 + d z 2 )
where A = A ( t , x , y , z ) , B = B ( t , x , y , z ) . Such models models can be matched to a dust FLRW model across a matter-comoving spherical junction surface [148] and may represent different types of topological defects in the context of topological quintessence (strings or walls).

3.3. Topological Quintessence

Topological quintessence [79] assumes the recent formation of a global topological defect (e.g., a global monopole) with Hubble scale core and with core vacuum energy dansity of the same order as the current density of matter. The dynamics of such a global monopole is determined by the action:
S = d 4 x - g m P l 2 16 π R - 1 2 ( μ Φ a ) 2 - V ( Φ ) + L m
where L m is the Lagrangian density of matter fields, Φ a ( a = 1 , 2 , 3 ) is an O ( 3 ) symmetric scalar field and:
V ( Φ ) = 1 4 λ ( Φ 2 - η 2 ) 2 , Φ Φ a Φ a
The vacuum energy density in the monopole core and the size of the core are determined by the two parameters of the model η (the vacuum expectation value) and λ (the coupling constant). The global monopole field configuration is described by the hedgehog ansatz:
Φ a = Φ ( r , t ) r ^ a Φ ( r , t ) ( sin θ cos φ , sin θ sin φ , cos θ )
with boundary conditions:
Φ ( 0 , t ) = 0 , Φ ( , t ) = η
where a time dependence is allowed, having in mind a cosmological setup of an expanding background.
The general spherically symmetric spacetime around a global monopole may be described by a metric of the form Equation (20). In order to derive the cosmological dynamical equations we need the total energy momentum tensor T μ ν . We assume a cosmological setup at recent cosmological times and therefore T μ ν is dominated by the global monopole vacuum energy and matter, i.e.:
T μ ν = T μ ν ( m o n ) + T μ ν ( m a t )
The energy momentum tensor of the monopole is given by:
T μ ν ( m o n ) = μ Φ a ν Φ a - g μ ν 1 2 ( σ Φ a ) 2 + V ( Φ )
The negative pressure in the monopole core is expected to lead to accelerating expansion when the core size is of the order of the Hubble scale [149]. The monopole energy density is of the form:
ρ m o n = T 00 m o n = Φ ˙ 2 2 + Φ 2 2 A 2 + Φ 2 B 2 r 2 + λ 4 Φ 2 - η 2 2
and has a local maximum at the monopole core where vacuum energy dominates. The energy-momentum tensor of matter is that of a perfect fluid with zero pressure ( p = 0 ):
T μ ν ( m a t ) = ρ m a t u μ u ν
where ρ m a t = ρ m a t ( r , t ) is the matter density and u μ is the velocity 4-vector of the fluid:
u μ = ( 1 1 - v 2 , v A 1 - v 2 , 0 , 0 )
with v = v ( r , t ) the radial matter fluid velocity. An off-center observer in the global monopole will experience anisotropic accelerating cosmological expansion which is maximized in the direction of the monopole center where the vacuum dark energy density is maximum (Figure 5). This is consistent with the mild evidence (2σ) for the existence of a dark energy dipole discussed in Section 2.
Figure 5. An off-center observer in the core of a global monopole will naturally experience an anisotropy in the expansion rate due to the higher vacuum energy density in the direction of the monopole center.
Figure 5. An off-center observer in the core of a global monopole will naturally experience an anisotropy in the expansion rate due to the higher vacuum energy density in the direction of the monopole center.
Galaxies 02 00022 g005
The question that arises is “Can this mechanism also predict the existence of a fine structure constant dipole aligned with the dark energy dipole as seems to be favored by some observations?” The answer to this question is positive. An α dipole aligned with the dark energy dipole can be obtained by considering a dilatonic coupling of the monopole scalar field to electromagnetism. Consider the generalization of the Lagrangian Equation (32) of the form:
S = 1 2 M p 2 R - 1 2 ( μ Φ a ) 2 - V ( Φ ) + 1 4 B ( Φ ) F μ ν 2 + L m - g d 4 x
where M p - 2 = 8 π G is the reduced Planck mass. This action is inspired from the Bekenstein-Sandvik-Barrow-Magueijo class of models [150,151,152] where the dynamics are attributed to a homogeneous scalar field. In our case (extended topological quintessence) the fine structure constant α has obtained dynamics through the inhomogeneous global monopole scalar field. It may be shown that:
α ( Φ ) = e 0 2 4 π B ( Φ ) 2
where e 0 is the bare charge that remains constant throughout the cosmological evolution. By considering a dilatonic coupling B ( Φ ) of the form Consider now a non-minimal coupling of the form:
B ( Φ ) = 1 - ξ Φ 2 η 2
where ξ is constant, we find for small values of Φ / η a spatial variation of α in accordance to:
Δ α α 2 ξ ( Φ 2 - Φ 0 2 ) η 2
where Φ 0 is the field magnitude at the location of the observer. Clearly, α is predicted to vary most in the direction where the field magnitude varies more rapidly. This is the direction towards the center of the monopole which is identical to the direction of the dark energy dipole. Thus, this model naturally predicts the observed alignment between the α dipole and the dark energy dipole.
The observationally required values of the parameters of topological quintessence may be obtained by implementing the requirement that the core size is comparable to the Hubble scale and the vacuum energy density at the monopole core is comparable to the matter density at present. The core scale (where Φ / η < 1 ) of the static global monopole is:
δ λ - 1 / 2 η - 1
while the vacuum energy density in this core region is:
ρ c o r e λ η 4 4
Deep inside the monopole core the expansion is approximately isotropic ( A ( r , t ) B ( r , t ) a i n ( t ) ) and the dynamics of the scale factors are approximated by:
3 H i n ( t ) 2 = 8 π m P l 2 λ η 4 4 + ρ 0 , i n m a t a 0 , i n a i n 3
where a 0 , i n a i n ( t 0 ) is the scale factor at the present time t 0 , ρ 0 , i n m a t is the present density of matter in the monopole core and:
H i n ( t ) a ˙ i n a i n
Similarly, far away from the monopole core ( r δ ), spacetime is also approximately homogeneous and we may set A ( r , t ) B ( r , t ) a o u t ( t ) . In this limit the dynamics of the scale factors are obtained from:
3 H o u t ( t ) 2 = 8 π m P l 2 ρ 0 , o u t m a t a 0 , o u t a o u t 3
corresponding to a flat matter dominated universe. By imposing the physical requirements stated above namely:
δ H 0 - 1 H i n ( t 0 ) - 1
ρ c o r e ρ m a t ( t 0 )
where ρ c o r e is the vacuum energy density in the monopole core, we obtain order of magnitude estimates for the two parameters in the global monopole action:
η O ( 1 ) × M p l
λ 10 - 120 η ¯
where η ¯ η M p l . In what follows we omit the bar on the dimensionless η ¯ .
In an effort to go beyond the above heuristic estimates we consider the full Einstein equations and using the metric Equation (20) and the energy-momentum tensor Equation (36) we obtain the system of Equations [79,149,153]:
- G 0 0 = K 2 2 ( 2 K - 3 K 2 2 ) - 2 B A 2 B - B 2 A 2 B 2 + 2 A B A 3 B - 6 B A 2 B r + 2 A A 3 r - 1 A 2 r 2 + 1 B 2 r 2
= 8 π m P l 2 Φ ˙ 2 2 + Φ 2 2 A 2 + Φ 2 B 2 r 2 + λ 4 ( Φ 2 - η 2 ) 2 + ρ m a t 1 - v 2
1 2 G 01 = K 2 2 + B B + 1 r 3 K 2 2 - K = 4 π m P l 2 Φ ˙ Φ - v 1 - v 2 A ρ m a t
1 2 ( G 1 1 + G 2 2 + G 3 3 - G 0 0 ) = K ˙ - ( K 1 1 ) 2 - 2 ( K 2 2 ) 2
= 8 π m P l 2 Φ ˙ 2 - λ 4 ( Φ 2 - η 2 ) 2 + 1 2 1 + v 2 1 - v 2 ρ m a t
and:
Φ ¨ - K Φ ˙ - Φ A 2 - - A A + 2 B B + 2 r Φ A 2 + 2 Φ B 2 r 2 + λ Φ ( Φ 2 - η 2 ) = 0
where prime (dot) denotes differentiation with respect to radius r (time t) and:
K 1 1 = - A ˙ A , K 2 2 = K 3 3 = - B ˙ B , K = K i i
The conservation of the energy momentum tensor for the matter fluid Equation (39) in the metric Equation (20) leads to the following equations for v and ρ m a t :
v ˙ v = ( v 2 - 1 ) A ˙ A - v A
ρ ˙ m a t ρ m a t = v ˙ v - 2 B ˙ B - v A ( ρ m a t ) ρ m a t - 2 v A B B - 2 v r A
Notice that in the homogeneous FLRW limit ( B ( r , t ) = A ( r , t ) = a ( t ) , Φ ( r , t ) = η ) we obtain the familiar Friedman equations with v = 0 and the usual conservation of matter equation. The system Equations (54)–(57), (59) and (60) consists of six equations with five unknown functions A, B, Φ, v, ρ m a t . Thus, Equation (54) may be considered as a constraint and the system of the other five equations may be solved with initial conditions corresponding to a homogeneous flat matter dominated universe ( A = B = 1 ) with a static global monopole profile f ( r ) ( f ( 0 ) = 0 , f ( ) = 1 ) [79]. The main questions to be addresses in such a cosmological simulation are the following:
  • Does the monopole energy density eventually dominate over matter in the monopole core?
  • Does the possible monopole domination eventually lead to accelerating expansion in the monopole core?
  • Can this cosmological expansion in the core fit the cosmological data?
The answer to all three of the above questions turns out to be positive. The evolution of the densities is shown in Figure 6. As expected, matter develops an underdensity in the core due to the repulsive gravity of the vacuum energy.
Clearly the vacuum energy eventually dominates in the monopole core. The time evolution of the scale factors A and B is shown in Figure 7 for two values of η ¯ . Clearly the core experiences a faster (accelerating) expansion compared to regions outside the monopole core. For η ¯ > 0 . 3 the comoving monopole core size appears to expand in agreement with previous studies of topological inflation.
In Figure 8 we show the evolution of the Hubble expansion rate ratios H A ( r , z ) / H A ( r , 0 ) and H B ( r , z ) / H B ( r , 0 ) at various comoving distances from the monopole center. Superposed is shown the ΛCDM Hubble expansion rate for the best fit values of Ω Λ . We have chosen the present time such that the best fit value at the core center is Ω Λ = 0 . 73 . Clearly the expansion rates at all distances are almost identical with ΛCDM with diminishing choices of Ω Λ . Thus, for a core size larger than a few Gpc we anticipate good agreement of the predictions of the model with SnIa and other cosmological data.
Figure 6. Evolution of the dark energy (left panel) and properly normalized matter (right panel) density profiles.The profiles shown correspond to the initial time t 0 (dotted line), t p / 3 (dot-dashed line), 2 t p / 3 (dashed line), t p (solid line, present time). Notice the matter underdensity that develops in the monopole core while the monopole density appears to slowly collapse to the center in comoving coordinates. The radial unit is the monopole core size δ (from [79]).
Figure 6. Evolution of the dark energy (left panel) and properly normalized matter (right panel) density profiles.The profiles shown correspond to the initial time t 0 (dotted line), t p / 3 (dot-dashed line), 2 t p / 3 (dashed line), t p (solid line, present time). Notice the matter underdensity that develops in the monopole core while the monopole density appears to slowly collapse to the center in comoving coordinates. The radial unit is the monopole core size δ (from [79]).
Galaxies 02 00022 g006
Figure 7. Evolution of the scale factors A ( r , t ) and B ( r , t ) profile for η = 0 . 1 (black lines) and η = 0 . 6 (blue lines) from the time t 0 up to the present time t p . A mesh of dots is superimposed to the curve corresponding to η = 0 . 1 , t = t p to distinguish it from the curve obtained for η = 0 . 6 , t = t p . The profiles shown correspond to t 0 (dotted line), t p / 3 (dot-dashed line), 2 t p / 3 (dashed line), t p (solid line). Higher curves correspond to more recent times. The (rescaled) comoving region of accelerated expansion is slightly smaller in the case of η < 0 . 3 since the monopole core does not expand in this case (from [79]).
Figure 7. Evolution of the scale factors A ( r , t ) and B ( r , t ) profile for η = 0 . 1 (black lines) and η = 0 . 6 (blue lines) from the time t 0 up to the present time t p . A mesh of dots is superimposed to the curve corresponding to η = 0 . 1 , t = t p to distinguish it from the curve obtained for η = 0 . 6 , t = t p . The profiles shown correspond to t 0 (dotted line), t p / 3 (dot-dashed line), 2 t p / 3 (dashed line), t p (solid line). Higher curves correspond to more recent times. The (rescaled) comoving region of accelerated expansion is slightly smaller in the case of η < 0 . 3 since the monopole core does not expand in this case (from [79]).
Galaxies 02 00022 g007
Figure 8. Ratios H A ( r , z ) / H A ( r , 0 ) and H B ( r , z ) / H B ( r , 0 ) for r = 0 (lower solid line), r = 0 . 5 and r = 5 (upper solid lines) along with the corresponding H Λ ( r , z ) / H 0 for Ω Λ = 0 . 73 (thick, dashed line) and for Ω Λ = 0 (dot-dashed line). Also included for comparison H Λ ( r , z ) / H 0 for the best fit values of Ω Λ ( r ) : Ω Λ ( 0 . 5 ) 0 . 27 (left panel) and Ω Λ ( 0 . 5 ) 0 . 54 (right panel). We have set η = 0 . 1 (from [79]).
Figure 8. Ratios H A ( r , z ) / H A ( r , 0 ) and H B ( r , z ) / H B ( r , 0 ) for r = 0 (lower solid line), r = 0 . 5 and r = 5 (upper solid lines) along with the corresponding H Λ ( r , z ) / H 0 for Ω Λ = 0 . 73 (thick, dashed line) and for Ω Λ = 0 (dot-dashed line). Also included for comparison H Λ ( r , z ) / H 0 for the best fit values of Ω Λ ( r ) : Ω Λ ( 0 . 5 ) 0 . 27 (left panel) and Ω Λ ( 0 . 5 ) 0 . 54 (right panel). We have set η = 0 . 1 (from [79]).
Galaxies 02 00022 g008

3.4. Spherical Dark Energy Overdensity

Even though topological quintessence is a well defined physical model based on specific dynamics, the predicted nonlinear evolution of dark energy complicates the detailed comparison with cosmological observation. This comparison, is significantly simplified if we approximate the dark energy inhomogeneity as a static spherically symmetric fluid with negative pressure in a flat background. In particular we may consider the energy momentum tensor:
T ν μ = diag ( - ρ ( r ) - ρ M ( r , t ) , p r ( r ) , p t ( r ) , p t ( r ) )
where ρ M ( r , t ) is the matter density and we have allowed for a general inhomogeneous and spherically symmetric static fluid with:
( T f ) ν μ = diag ( - ρ ( r ) , p r ( r ) , p t ( r ) , p t ( r ) )
(where p t ( r ) is the transverse pressure). This expression is motivated by the energy-momentum of a global monopole, which asymptotically is of the form [154]:
( T mon ) ν μ = diag ( - η / r 2 , - η / r 2 , 0 , 0 )
with η the symmetry-breaking scale related to the formation of the monopole.
Using Einstein equations, it may be shown [86] that the metric corresponding to the diagonal energy momentum tensor Equation (61) is of the LTB form Equation (21). Using Einstein equation we obtain the cosmological Equation (22). In contrast to the matter void here we set Ω c ( r ) = 0 and use the profiles Equation (24) and (25) for the matter and dark energy fluids which are consistent with flatness ( Ω M , out = 1 - Ω X , out = 1 , Figure 9).
Figure 9. The assumed matter and dark energy profiles are consistent with flatness and with matter domination outside the inhomogeneity. Such profiles would develop dynamically in the context of topological quintessence (see Figure 6).
Figure 9. The assumed matter and dark energy profiles are consistent with flatness and with matter domination outside the inhomogeneity. Such profiles would develop dynamically in the context of topological quintessence (see Figure 6).
Galaxies 02 00022 g009
Considering an observer at the center of the inhomogeneity, we may use the geodesics [86,144] to obtain the luminosity distance as described above Equation (29). By comparing the predicted luminosity distance with the Union 2 data [14], we obtain [86] constraints on the model parameters Ω M , in = 1 - Ω X , in and r 0 (the size of the inhomogeneity) using the maximum likelihood method [86]. The allowed range of parameters is significantly larger for the inhomogeneous dark energy model (IDE) and as expected, the model reduces to ΛCDM for a scale of inhomogeneity that is larger than the range of the Union 2 dataset (see Figure 10).
It is straightforward to extend this analysis to the case of an off-center observer [86,155]. In this case the light-like geodesics of the photons coming from the SnIa, depend on two additional parameters: the distance of the observer from the center of the inhomogeneity ( r o b s ) and the incidence angle ξ of the photons with respect to the displacement axis (Figure 11). The angle ξ in turn may be expressed in terms of the direction of the center in galactic coordinates ( l c , b c ) and the direction the supernova that emits the photon ( l , b ) ( ξ = ξ ( l c , b c , l , b ) ). Thus, it is straightforward [86] to obtain an expression of the predicted luminosity distance as d L ( z , r 0 , Ω X , in ; r o b s , ξ ) and use the Union 2 data (with the directions ( l , b ) of each SnIa) with the maximum likelihood method to derive χ 2 ( r 0 , Ω X , in ; r o b s , l c , b c ) .
Figure 10. The 1σ and 2σ parameter contours for the inhomogeneous dark energy (IDE) fluid model obtained using the Union 2 data, compared with the corresponding contours of the usual LTM matter void model. Notice that for an inhomogeneity size larger than the range of the dataset (about 4 Gpc), the model reduces to ΛCDM( Ω m = 2 . 7 , Ω X = 0 . 73 . The transition scale was fixed to Δ r = 0.35 Gpc (from [86]).
Figure 10. The 1σ and 2σ parameter contours for the inhomogeneous dark energy (IDE) fluid model obtained using the Union 2 data, compared with the corresponding contours of the usual LTM matter void model. Notice that for an inhomogeneity size larger than the range of the dataset (about 4 Gpc), the model reduces to ΛCDM( Ω m = 2 . 7 , Ω X = 0 . 73 . The transition scale was fixed to Δ r = 0.35 Gpc (from [86]).
Galaxies 02 00022 g010
Figure 11. The lightlike geodesics for an off-center depend on two additional parameters r o b s and ξ.
Figure 11. The lightlike geodesics for an off-center depend on two additional parameters r o b s and ξ.
Galaxies 02 00022 g011
In Figure 12 we show the dependence of χ 2 on the displacement r o b s of the observer (left panel). Minimization with respect to the direction of the inhomogeneity center has been performed. We have considered four pairs ( r 0 , Ω X , in ) that provide good fits for the on center observer (right panel). Notice that improved fit with respect to ΛCDM (horizontal blue dashed line) is obtained only for large scale inhomogeneities ( r 0 > 3 . 5 Gpc).
Figure 12. The dependence of χ 2 on the displacement r o b s of the observer (left panel). Minimization with respect to the direction of the inhomogeneity center has been performed. We have considered four pairs ( r 0 , Ω X , in ) that provide good fits for the on center observer (right panel). Notice that improved fit with respect to ΛCDM (horizontal blue dashed line) is obtained only for large scale inhomogeneities ( r 0 > 3 . 5 Gpc) (from [86]).
Figure 12. The dependence of χ 2 on the displacement r o b s of the observer (left panel). Minimization with respect to the direction of the inhomogeneity center has been performed. We have considered four pairs ( r 0 , Ω X , in ) that provide good fits for the on center observer (right panel). Notice that improved fit with respect to ΛCDM (horizontal blue dashed line) is obtained only for large scale inhomogeneities ( r 0 > 3 . 5 Gpc) (from [86]).
Galaxies 02 00022 g012
It is straightforward [86] to also use geodesics from the last scattering surface to the present time to derive the CMB photon geodesics for an off center observer. From these we may obtain the last scattering surface redshift z l s as a function of the inhomogeneity parameters and in particular ξ and r o b s (assuming fixed the inhomogeneity parameters r 0 , Ω X , in ). After obtaining z l s ( ξ , r o b s ) we find the angular dependence of the temperature as:
T ( ξ , r o b s ) = T l s 1 + z l s ( ξ , r o b s )
where T l s is the observed temperature of the last scattering surface. This is assumed to be uniform since we are focusing on the additional temperature fluctuations induced by the displacement of the observer. We then convert it to temperature fluctuations from the mean and expand in spherical harmonics to find the multipole coefficients. Comparison with the measured values of the moments (particularly of the dipole) leads to constraints on r o b s . It is straightforward to obtain numerically a 10 ( r o b s ) , the additional CMB fluctuations dipole moment induced by the shift r o b s of the observer. This quantity is shown in Figure 13 for various values of the inhomogeneity scale r 0 . In the same plot we show the value of the measured CMB dipole (dashed line) [156]:
Δ T T ¯ 10 = 3 . 35 × 10 - 3 K 2 . 725 K = ( 1 . 230 ± 0 . 013 ) × 10 - 3 a 10 = 4 π 3 Δ T T ¯ 10 ( 2 . 52 ± 0 . 03 ) × 10 - 3
Figure 13. The predicted CMB dipole moment a 10 ( r o b s ) for Ω X , in = 0 . 69 and various values of the inhomogeneity scale r 0 . The dashed line corresponds to the observed value (whose uncertainty is too small to be shown (from [86]).
Figure 13. The predicted CMB dipole moment a 10 ( r o b s ) for Ω X , in = 0 . 69 and various values of the inhomogeneity scale r 0 . The dashed line corresponds to the observed value (whose uncertainty is too small to be shown (from [86]).
Galaxies 02 00022 g013
Clearly, for r 0 < 5 Gpc the observer is confined to be in a sphere which is a fine tuned small spatial fraction of the dark energy inhomogeneity ( f ( r 0 ) ( r o b s - m a x r 0 ) 3 10 - 6 ) which implies severe fine tuning of the model as in the case of LTB models with matter. This is shown in Figure 14 where f ( r 0 ) is plotted for values of r 0 between 1 and 7 Gpc. As shown in that figure, the fine tuning starts to reduce ( f ( r 0 ) increases) as the size of the inhomogeneity increases beyond 6 Gpc. We anticipate it to disappear as the size of the inhomogeneity reaches ∼14 Gpc (the comoving distance to the last scattering surface).
The predicted CMB maps of additional temperature fluctuations in galactic coordinates corresponding to the full map (a), dipole (b) quadrupole (c) and octopole (d) are shown in Figure 15. Notice the rapid decrease of the magnitude of moments higher than the dipole. Even though the order of magnitude of the dipole and quadrupole is comparable with observations, higher moments have negligible magnitude. This is anticipated since the displacement of the observer is expected to modify the CMB fluctuations only at the largest scales.
Figure 14. The spatial fraction f ( r 0 ) ( r o b s - m a x r 0 ) 3 where the observer needs to be confined in order to be consistent with the value of the observed CMB dipole, for ( r 0 , Ω X , in ) = ( 3 . 37 Gpc , 0 . 69 ) (from [86]).
Figure 14. The spatial fraction f ( r 0 ) ( r o b s - m a x r 0 ) 3 where the observer needs to be confined in order to be consistent with the value of the observed CMB dipole, for ( r 0 , Ω X , in ) = ( 3 . 37 Gpc , 0 . 69 ) (from [86]).
Galaxies 02 00022 g014
Figure 15. The predicted CMB maps of additional temperature fluctuations in galactic coordinates corresponding to the full map (a); dipole (b) quadrupole (c) and octopole (d). We have selected the direction of the preferred axis to coincide with the direction of the observed CMB dipole and used ( r 0 , Ω X , in ) = ( 3 . 37 Gpc , 0 . 69 ) , r o b s = 30 Mpc (from [86]).
Figure 15. The predicted CMB maps of additional temperature fluctuations in galactic coordinates corresponding to the full map (a); dipole (b) quadrupole (c) and octopole (d). We have selected the direction of the preferred axis to coincide with the direction of the observed CMB dipole and used ( r 0 , Ω X , in ) = ( 3 . 37 Gpc , 0 . 69 ) , r o b s = 30 Mpc (from [86]).
Galaxies 02 00022 g015

4. Conclusions

It has become evident during the last few years that early hints for deviation from the cosmological principle and statistical isotropy are being accumulated. This appears to be one of the most likely directions which may lead to new fundamental physics in the coming years. We have reviewed these hints, including the large scale CMB anomalies (CMB asymmetries), the large scale peculiar velocity flows on scales of 100 Mpc and up to a Gpc (dark flow), the preliminary evidence for the spatial variation of the fine structure constant along a preferred axis on Gpc scales, the coherent variation of quasar optical polarization on Gpc scales, the existence of the Large Quasar Structure with a Gpc scale and the mild hints for anisotropic expansion rate. We have pointed out the fact that most of the above cosmological observations are associated with the existence of specific cosmological directions which appear to be abnormally aligned with each other.
The Gpc scale associated with most of these effects indicates that their physical origin has appeared at late cosmological times. This possibility is further amplified by the fact that the CMB power asymmetry is significantly diminished if the ISW effect is subtracted from the CMB maps [94]. It is therefore natural to attempt to associate the physical origin of these effects with large scale dark energy or matter inhomogeneities. Hubble scale matter inhomogeneities in the form of voids are severely constrained by cosmological observations and are practically ruled out. On the other hand, dark energy Hubble scale inhomogeneities are observationally viable and reduce to standard ΛCDM when the inhomogeneity scale exceeds the horizon scale. A simple mechanism that can give rise to a cosmological preferred axis is based on an off-center observer in a spherical dark energy inhomogeneity. We have shown that this mechanism can naturally give rise to the observed alignment of many of the above described cosmic anisotropies.
The repulsive gravity properties of dark energy and its high sound velocity imply that any large scale dark energy inhomogeneity would tend to smooth out. This is a challenge for the above described mechanism. Nontrivial topology however can lead to sustainable Hubble scale dark energy inhomogeneities. If we happen to leave in the core of a Hubble scale recently formed topological defect we would naturally experience Hubble scale dark energy inhomogeneities and the existence of a preferred axis due to our displacement from the defect center. Therefore, topological quintessence constitutes a physical mechanism to produce Hubble scale dark energy inhomogeneities. We have shown that such a mechanism can also give rise to a Dark Energy Dipole, Large Scale Velocity Flows, Fine Structure Constant Dipole and a CMB Temperature Asymmetry aligned with each other. Other effects which may also occur (quasar polarization alignment etc) constitute interesting potential extensions of this research direction.

Acknowledgments

This research has been co-financed by the European Union (European Social Fund–ESF) and Greek national funds through the Operational Program “Education and Lifelong Learning” of the National Strategic Reference Framework (NSRF)–Research Funding Program: ARISTEIA. Investing in the society of knowledge through the European Social Fund.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hamilton, J.-C. What have we learned from observational cosmology? ArXiv E-Prints 2013. arXiv:1304.4446. [Google Scholar]
  2. Komatsu, E.; Smith, K.M.; Dunkley, J.; Bennett, C.L.; Gold, B.; Hinshaw, G.; Jarosik, N.; Larson, D.; Nolta, M.R.; Page, L.; et al. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Cosmological interpretation. Astrophys. J. 2011, 192. [Google Scholar] [CrossRef]
  3. Tegmark, M.; de Oliveira-Costa, A.; Hamilton, A. A high resolution foreground cleaned CMB map from WMAP. Phys. Rev. D 2003, 68, 123523. [Google Scholar] [CrossRef]
  4. Copi, C.J.; Huterer, D.; Schwarz, D.J.; Starkman, G.D. Large angle anomalies in the CMB. Adv. Astron. 2010, 2010, 847541. [Google Scholar] [CrossRef]
  5. Bennett, C.L.; Hill, R.S.; Hinshaw, G.; Larson, D.; Smith, K.M.; Dunkley, J.; Gold, B.; Halpern, M.; Jarosik, N.; Kogut, A.; et al. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Are there cosmic microwave background anomalies? Astrophys. J. 2011, 192. [Google Scholar] [CrossRef]
  6. Schwarz, D.J.; Starkman, G.D.; Huterer, D.; Copi, C.J. Is the low-l microwave background cosmic? Phys. Rev. Lett. 2004, 93, 221301. [Google Scholar] [CrossRef] [PubMed]
  7. Land, K.; Magueijo, J. The axis of evil. Phys. Rev. Lett. 2005, 95, 071301. [Google Scholar] [CrossRef] [PubMed]
  8. Gruppuso, A.; Gorski, K.M. Large scale directional anomalies in the WMAP 5yr ILC map. J. Cosmol. Astropart. Phys. 2010, 2010, 019. [Google Scholar] [CrossRef]
  9. Sarkar, D.; Huterer, D.; Copi, C.J.; Starkman, G.D.; Schwarz, D.J. Missing power vs. low-l alignments in the cosmic microwave background: No correlation in the standard cosmological model. ArXiv E-Prints 2010. arXiv:1004.3784. [Google Scholar]
  10. Hanson, D.; Lewis, A. Estimators for CMB statistical anisotropy. Phys. Rev. D. 2009, 80, 063004. [Google Scholar] [CrossRef]
  11. Copi, C.J.; Huterer, D.; Schwarz, D.J.; Starkman, G.D. The uncorrelated universe: Statistical anisotropy and the vanishing angular correlation function in WMAP years 1–3. Phys. Rev. D 2007, 75, 023507. [Google Scholar] [CrossRef]
  12. Smith, K.M.; Senatore, L.; Zaldarriaga, M. Optimal limits on f N L l o c a l from WMAP 5-year data. J. Cosmol. Astropart. Phys. 2009, 2009, 006. [Google Scholar] [CrossRef] [PubMed]
  13. Ade, P.A.R.; Aghanim, N.; Armitage-Caplan, C.; Arnaud, M.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A.J.; Barreiro, R.B.; et al. Planck 2013 results. I. Overview of products and scientific results. ArXiv E-Prints 2013. arXiv:1303.5062. [Google Scholar]
  14. Amanullah, R.; Lidman, C.; Rubin, D.; Aldering, G.; Astier, P.; Barbary, K.; Burns, M.S.; Conley, A.; Dawson, K.S.; Deustua, S.E.; et al. Spectra and light curves of six type Ia supernovae at 0.511 < z < 1.12 and the Union2 compilation. Astrophys. J. 2010, 716, 712. [Google Scholar] [CrossRef]
  15. Percival, W.J.; Cole, S.; Eisenstein, D.J.; Nichol, R.C.; Peacock, J.A.; Pope, A.C.; Szalay, A.S. Measuring the Baryon Acoustic Oscillation scale using the SDSS and 2dFGRS. Mon. Not. R. Astron. Soc. 2007, 381, 1053–1066. [Google Scholar] [CrossRef] [Green Version]
  16. Bueno Sanchez, J.C.; Nesseris, S.; Perivolaropoulos, L. Comparison of recent SnIa datasets. J. Cosmol. Astropart. Phys. 2009, 2009, 029. [Google Scholar] [CrossRef]
  17. Nesseris, S.; Perivolaropoulos, L. Testing ΛCDM with the growth function δ(a): Current constraints. Phys. Rev. D 2008, 77, 023504. [Google Scholar] [CrossRef]
  18. Trujillo-Gomez, S.; Klypin, A.; Primack, J.; Romanowsky, A.J. LCDM correctly predicts basic statistics of galaxies: Luminosity-velocity relation, baryonic mass-velocity relation, and velocity function. ArXiv E-Prints 2011. arXiv:1005.1289. [Google Scholar]
  19. Reid, B.A.; Percival, W.J.; Eisenstein, D.J.; Verde, L.; Spergel, D.N.; Skibba, R.A.; Bahcall, N.A.; Budavari, T.; Fukugita, M.; Richard Gott, J.; et al. Cosmological constraints from the clustering of the sloan digital sky survey DR7 luminous red galaxies. Mon. Not. R. Astron. Soc. 2010, 404, 60–85. [Google Scholar] [CrossRef] [Green Version]
  20. Perivolaropoulos, L. Six puzzles for LCDM cosmology. ArXiv E-Prints 2008. arXiv:0811.4684. [Google Scholar]
  21. Yang, R.-J.; Zhang, S.N. The age problem in ΛCDM model. Mon. Not. R. Astron. Soc. 2010, 407, 1835–1841. [Google Scholar] [CrossRef]
  22. Ade, P.A.R.; Aghanim, N.; Armitage-Caplan, C.; Arnaud, M.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A.J.; Barreiro, R.B.; et al. Planck 2013 results. XXIII. Isotropy and Statistics of the CMB. ArXiv E-Prints 2013. arXiv:1303.5083. [Google Scholar]
  23. Eriksen, H.K.; Banday, A.J.; Gorski, K.M.; Hansen, F.K.; Lilje, P.B. Hemispherical power asymmetry in the third-year Wilkinson Microwave Anisotropy Probe sky maps. Astrophys. J. 2007, 660, L81. [Google Scholar] [CrossRef]
  24. Paci, F.; Gruppuso, A.; Finelli, F.; Cabella, P.; de Rosa, A.; Mandolesi, N.; Natoli, P. Power asymmetries in the cosmic microwave background temperature and polarization patterns. Mon. Not. R. Astron. Soc. 2010, 407, 399–404. [Google Scholar] [CrossRef] [Green Version]
  25. Hoftuft, J.; Eriksen, H.K.; Banday, A.J.; Gorski, K.M.; Hansen, F.K.; Lilje, P.B. Increasing evidence for hemispherical power asymmetry in the five-year WMAP data. Astrophys. J. 2009, 699, 985. [Google Scholar] [CrossRef]
  26. Axelsson, M.; Fantaye, Y.; Hansen, F.K.; Banday, A.J.; Eriksen, H.K.; Gorski, K.M. Directional dependence of ΛCDM cosmological parameters. Astrophys. J. 2013, 773, L3. [Google Scholar] [CrossRef]
  27. Mariano, A.; Perivolaropoulos, L. CMB maximum temperature asymmetry axis: Alignment with other cosmic asymmetries. Phys. Rev. D 2013, 87, 043511. [Google Scholar] [CrossRef]
  28. Watkins, R.; Feldman, H.A.; Hudson, M.J. Consistently large cosmic flows on scales of 100 Mpc/h: A challenge for the standard LCDM cosmology. Mon. Not. R. Astron. Soc. 2009, 392, 743–756. [Google Scholar] [CrossRef]
  29. Feldman, H.A.; Watkins, R.; Hudson, M.J. Cosmic flows on 100 Mpc/h scales: Standardized minimum variance bulk flow, shear and octupole moments. Mon. Not. R. Astron. Soc. 2010, 407, 2328–2338. [Google Scholar] [CrossRef]
  30. Kashlinsky, A.; Atrio-Barandela, F.; Kocevski, D.; Ebeling, H. A measurement of large-scale peculiar velocities of clusters of galaxies: Results and cosmological implications. Astrophys. J. 2008, 686, L49. [Google Scholar] [CrossRef]
  31. Lavaux, G.; Tully, R.B.; Mohayaee, R.; Colombi, S. Cosmic flow from 2MASS redshift survey: The origin of CMB dipole and implications for LCDM cosmology. Astrophys. J. 2010, 709, 483. [Google Scholar] [CrossRef]
  32. Hutsemekers, D.; Lamy, H. The polarization properties of broad absorption line QSOs: Observational results. ArXiv E-Prints 2001. arXiv:astro-ph/0109528. [Google Scholar]
  33. Hutsemekers, D.; Cabanac, R.; Lamy, H.; Sluse, D. Mapping extreme-scale alignments of quasar polarization vectors. Astron. Astrophys. 2005, 441, 915. [Google Scholar] [CrossRef]
  34. Hutsemekers, D.; Payez, A.; Cabanac, R.; Lamy, H.; Sluse, D.; Borguet, B.; Cudell, J.R. Large-scale alignments of quasar polarization vectors: Evidence at cosmological scales for very light pseudoscalar particles mixing with photons? ArXiv E-Prints 2008. arXiv:0809.3088. [Google Scholar]
  35. Mariano, A.; Perivolaropoulos, L. Is there correlation between fine structure and dark energy cosmic dipoles? Phys. Rev. D 2012, 86, 083517. [Google Scholar] [CrossRef]
  36. Webb, J.K.; King, J.A.; Murphy, M.T.; Flambaum, V.V.; Carswell, R.F.; Bainbridge, M.B. Indications of a spatial variation of the fine structure constant. Phys. Rev. Lett. 2011, 107, 191101. [Google Scholar] [CrossRef] [PubMed]
  37. King, J.A.; Webb, J.K.; Murphy, M.T.; Flambaum, V.V.; Carswell, R.F.; Bainbridge, M.B.; Wilczynska, M.R.; Koch, F.E. Spatial variation in the fine-structure constant—New results from VLT/UVES. Mon. Not. R. Astron. Soc. 2012, 422, 3370–3414. [Google Scholar] [CrossRef]
  38. Clowes, R.G.; Harris, K.A.; Raghunathan, S.; Campusano, L.E.; Soechting, I.K.; Graham, M.J. A structure in the early universe at z ∼ 1.3 that exceeds the homogeneity scale of the R-W concordance cosmology. ArXiv E-Prints 2012. arXiv:1211.6256. [Google Scholar]
  39. Gott, J.R., III; Juric, M.; Schlegel, D.; Hoyle, F.; Vogeley, M.; Tegmark, M.; Bahcall, N.; Brinkmann, J. A map of the universe. Astrophys. J. 2005, 624, 463. [Google Scholar] [CrossRef]
  40. Yadav, J.K.; Bagla, J.S.; Khandai, N. Fractal dimension as a measure of the scale of homogeneity. Mon. Not. R. Astron. Soc. 2010, 405, 2009–2015. [Google Scholar] [CrossRef]
  41. Luminet, J.-P. The shape and topology of the universe. ArXiv E-Prints 2008. arXiv:0802.2236. [Google Scholar]
  42. Bielewicz, P.; Riazuelo, A. The study of topology of the universe using multipole vectors. ArXiv E-Prints 2009. arXiv:0804.2437. [Google Scholar]
  43. Carneiro, S.; Mena Marugan, G.A. Anisotropic cosmologies containing isotropic background radiation. Phys. Rev. D 2001, 64, 083502. [Google Scholar] [CrossRef]
  44. Akofor, E.; Balachandran, A.P.; Jo, S.G.; Joseph, A.; Qureshi, B.A. Direction-dependent CMB power spectrum and statistical anisotropy from noncommutative geometry. J. High Energy Phys. 2008, 2008, 092. [Google Scholar] [CrossRef]
  45. Koivisto, T.S.; Mota, D.F.; Quartin, M.; Zlosnik, T.G. On the possibility of anisotropic curvature in cosmology. Phys. Rev. D 2011, 83, 023509. [Google Scholar] [CrossRef]
  46. Zumalacarregui, M.; Koivisto, T.S.; Mota, D.F.; Ruiz-Lapuente, P. Disformalscalar fields and the dark sector of the universe. J. Cosmol. Astropart. Phys. 2010, 2010, 038. [Google Scholar] [CrossRef]
  47. Koivisto, T.; Mota, D.F. Dark energy anisotropic stress and large scale structure formation. Phys. Rev. D 2006, 73, 083502. [Google Scholar] [CrossRef]
  48. Battye, R.; Moss, A. Anisotropic dark energy and CMB anomalies. Phys. Rev. D 2009, 80, 023531. [Google Scholar] [CrossRef]
  49. Armendariz-Picon, C. Could dark energy be vector-like? J. Cosmol. Astropart. Phys. 2004, 2004, 007. [Google Scholar] [CrossRef]
  50. Esposito-Farese, G.; Pitrou, C.; Uzan, J.-P. Vector theories in cosmology. Phys. Rev. D 2010, 81, 063519. [Google Scholar] [CrossRef]
  51. Rodrigues, D.C. Anisotropic cosmological constant and the CMB quadrupole anomaly. Phys. Rev. D 2008, 77, 023534. [Google Scholar] [CrossRef]
  52. Jimenez, J.B.; Maroto, A.L. Large-scale cosmic flows and moving dark energy model for anisotropization. J. Cosmol. Astropart. Phys. 2009, 2009, 015. [Google Scholar] [CrossRef]
  53. Alexander, S.; Biswas, T.; Notari, A.; Vaid, D. Local void vs. dark energy: Confrontation with wmap and type Ia supernovae. J. Cosmol. Astropart. Phys. 2009, 2009, 025. [Google Scholar] [CrossRef]
  54. Garcia-Bellido, J.; Haugboelle, T. Confronting Lemaître-Tolman-Bondi models with observational cosmology. J. Cosmol. Astropart. Phys. 2008, 2008, 003. [Google Scholar] [CrossRef]
  55. Biswas, T.; Notari, A.; Valkenburg, W. Testing the void against cosmological data: Fitting CMB, BAO, SN and H0. J. Cosmol. Astropart. Phys. 2010, 2010, 030. [Google Scholar] [CrossRef]
  56. Dunsby, P.; Goheer, N.; Osano, B.; Uzan, J.P. How close can an inhomogeneous universe mimic the concordance model? J. Cosmol. Astropart. Phys. 2010, 2010, 017. [Google Scholar] [CrossRef]
  57. Garfinkle, D. The motion of galaxy clusters in inhomogeneous cosmologies. Class. Quantum Gravity 2010, 27, 065002. [Google Scholar] [CrossRef]
  58. Krasinski, A. Inhomogeneous Cosmological Models; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  59. Bolejko, K.; Celerier, M.-N.; Krasinski, A. Inhomogeneous cosmological models: Exact solutions and their applications. Class. Quantum Gravity 2011, 28, 164002. [Google Scholar] [CrossRef]
  60. Alnes, H.; Amarzguioui, M.; Gron, O. An inhomogeneous alternative to dark energy? Phys. Rev. D 2006, 73, 083519. [Google Scholar] [CrossRef]
  61. Armendariz-Picon, C. Creating statistically anisotropic and inhomogeneous perturbations. J. Cosmol. Astropart. Phys. 2007, 2007, 014. [Google Scholar] [CrossRef]
  62. Pullen, A.R.; Kamionkowski, M. Cosmic microwave background statistics for a direction-dependent primordial power spectrum. Phys. Rev. D 2007, 76, 103529. [Google Scholar] [CrossRef]
  63. Ackerman, L.; Carroll, S.M.; Wise, M.B. Imprints of aprimordial preferred direction on the microwave background. Phys. Rev. D 2007, 75, 083502. [Google Scholar] [CrossRef]
  64. Valenzuela-Toledo, C.A. Non-gaussianity and statistical anisotropy in cosmological inflationary models. ArXiv E-Prints 2010. arXiv:1004.5363. [Google Scholar]
  65. Dimopoulos, K.; Karciauskas, M.; Lyth, D.H.; Rodriguez, Y. Statistical anisotropy of the curvature perturbation from vector field perturbations. J. Cosmol. Astropart. Phys. 2009, 2009, 013. [Google Scholar] [CrossRef]
  66. Yokoyama, S.; Soda, J. Primordial statistical anisotropy generated at the end of inflation. J. Cosmol. Astropart. Phys. 2008, 2008, 005. [Google Scholar] [CrossRef]
  67. Golovnev, A.; Vanchurin, V. Cosmological perturbations from vector inflation. Phys. Rev. D 2009, 79, 103524. [Google Scholar] [CrossRef]
  68. Bartolo, N.; Dimastrogiovanni, E.; Matarrese, S.; Riotto, A. Anisotropic bispectrum of curvature perturbations from primordial non-Abelian vector fields. J. Cosmol. Astropart. Phys. 2009, 2009, 015. [Google Scholar] [CrossRef]
  69. Himmetoglu, B.; Contaldi, C.R.; Peloso, M. Instability of anisotropic cosmological solutions supported by vector fields. Phys. Rev. Lett. 2009, 102, 111301. [Google Scholar] [CrossRef] [PubMed]
  70. Himmetoglu, B.; Contaldi, C.R.; Peloso, M. Instability of the ACW model, and problems with massive vectors during inflation. Phys. Rev. D 2009, 79, 063517. [Google Scholar] [CrossRef]
  71. Himmetoglu, B.; Contaldi, C.R.; Peloso, M. Ghost instabilities of cosmological models with vector fields nonminimally coupled to the curvature. Phys. Rev. D 2009, 80, 123530. [Google Scholar] [CrossRef]
  72. Thorsrud, M.; Mota, D.F.; Hervik, S. Cosmology of a scalar field coupled to matter and an isotropy-violating Maxwell field. J. High Energy Phys. 2012, 2012, 066. [Google Scholar] [CrossRef]
  73. Kahniashvili, T.; Lavrelashvili, G.; Ratra, B. CMB temperature anisotropy from broken spatial isotropy due to an homogeneous cosmological magnetic field. Phys. Rev. D 2008, 78, 063012. [Google Scholar] [CrossRef]
  74. Barrow, J.D.; Ferreira, P.G.; Silk, J. Constraints on a primordial magnetic field. Phys. Rev. Lett. 1997, 78, 3610–3613. [Google Scholar] [CrossRef]
  75. Campanelli, L. A model of universe anisotropization. Phys. Rev. D 2009, 80, 063006. [Google Scholar] [CrossRef]
  76. Kim, J.; Naselsky, P. Cosmological Alfvén waves in the recent CMB data, and the observational bound on the primordial vector perturbation. J. Cosmol. Astropart. Phys. 2009, 2009, 041. [Google Scholar] [CrossRef]
  77. Bueno Sanchez, J.C.; Perivolaropoulos, L. Dark energy and matter perturbations in scalar-tensor theories of gravity. J. Phys. Conf. Ser. 2011, 283, 012006. [Google Scholar] [CrossRef]
  78. Bean, R.; Dore, O. Probing dark energy perturbations: The dark energy equation of state and speed of sound as measured by WMAP. Phys. Rev. D 2004, 69, 083503. [Google Scholar] [CrossRef]
  79. Bueno Sanchez, J.C.; Perivolaropoulos, L. Topological quintessence. Phys. Rev. D 2011, 84, 123516. [Google Scholar] [CrossRef]
  80. Marra, V.; Notari, A. Observational constraints on inhomogeneous cosmological models without dark energy. Class. Quantum Gravity 2011, 28, 164004. [Google Scholar] [CrossRef]
  81. Zumalacarregui, M.; Garcia-Bellido, J.; Ruiz-Lapuente, P. Tension in the void: Cosmic rulers strain inhomogeneous cosmologies. J. Cosmol. Astropart. Phys. 2012, 2012, 009. [Google Scholar] [CrossRef] [PubMed]
  82. Caldwell, R.R.; Stebbins, A. A test of the copernican principle. Phys. Rev. Lett. 2008, 100, 191302. [Google Scholar] [CrossRef] [PubMed]
  83. Garcia-Bellido, J.; Haugboelle, T. Looking the void in the eyes—The kinematic Sunyaev–Zeldovich effect in Lemaître-Tolman-Bondi models. J. Cosmol. Astropart. Phys. 2008, 2008, 016. [Google Scholar] [CrossRef]
  84. Ricciardelli, E.; Quilis, V.; Planelles, S. The structure of cosmic voids in a LCDM Universe. ArXiv E-Prints 2013. arXiv:1306.2955. [Google Scholar]
  85. Alnes, H.; Amarzguioui, M. CMB anisotropies seen by an off-center observer in a spherically symmetric inhomogeneous Universe. Phys. Rev. D 2006, 74, 103520. [Google Scholar] [CrossRef]
  86. Grande, J.; Perivolaropoulos, L. Generalized Lemaître-Tolman-Bondi model with inhomogeneous isotropic dark energy: Observational constraints. Phys. Rev. D 2011, 84, 023514. [Google Scholar] [CrossRef]
  87. Vilenkin, A. Topological inflation. Phys. Rev. Lett. 1994, 72, 3137–3140. [Google Scholar] [CrossRef] [PubMed]
  88. Hansen, F.K.; Banday, A.J.; Gorski, K.M.; Eriksen, H.K.; Lilje, P.B. Power asymmetry in cosmic microwave background fluctuations from full sky to sub-degree scales: Is the universe isotropic? Astrophys. J. 2009, 704, 1448. [Google Scholar] [CrossRef]
  89. Eriksen, H.K.; Hansen, F.K.; Banday, A.J.; Gorski, K.M.; Lilje, P.B. Asymmetries in the cosmic microwave background anisotropy field. Astrophys. J. 2004, 605, 14. [Google Scholar] [CrossRef]
  90. Park, C.-G. Non-Gaussian signatures in the temperature fluctuation observed by the Wilkinson Microwave Anisotropy Probe. Mon. Not. R. Astron. Soc. 2004, 349, 313–320. [Google Scholar] [CrossRef]
  91. Hirata, C.M. Constraints on cosmic hemispherical power anomalies from quasars. J. Cosmol. Astropart. Phys. 2009, 2009, 011. [Google Scholar] [CrossRef]
  92. Dai, L.; Jeong, D.; Kamionkowski, M.; Chluba, J. The pesky power asymmetry. Phys. Rev. D 2013, 87, 123005. [Google Scholar] [CrossRef]
  93. Chang, Z.; Wang, S. Inflation and primordial power spectra at anisotropic spacetime inspired by Planck’s constraints on isotropy of CMB. Eur. Phys. J. C 2013, 73, 2516. [Google Scholar] [CrossRef]
  94. Rassat, A.; Starck, J.-L. On Preferred Axes in WMAP cosmic microwave background data after subtraction of the integrated sachs-wolfe effect. ArXiv E-Prints 2013. arXiv:1303.5051. [Google Scholar]
  95. Bielewicz, P.; Gorski, K.M.; Banday, A.J. Low-order multipole maps of CMB anisotropy derived from WMAP. Mon. Not. R. Astron. Soc. 2004, 355, 1283–1302. [Google Scholar] [CrossRef]
  96. Frommert, M.; Ensslin, T.A. The axis of evil—A polarization perspective. Mon. Not. R. Astron. Soc. 2010, 403, 1739–1748. [Google Scholar] [CrossRef]
  97. Land, K.; Magueijo, J. Is the Universe odd? Phys. Rev. D 2005, 72, 101302. [Google Scholar] [CrossRef]
  98. Land, K.; Magueijo, J. Cubic anomalies in WMAP. Mon. Not. R. Astron. Soc. 2005, 357, 994. [Google Scholar] [CrossRef]
  99. Park, C.-G.; Park, C.; Gott, J.R., III. Cleaned 3 year Wilkinson Microwave Anistropy Probe Cosmic Microwave Background map: Magnitude of the quadrupole and alignment of large-scale modes. Astrophys. J. 2007, 660, 959. [Google Scholar] [CrossRef]
  100. De Oliveira-Costa, A.; Tegmark, M. CMB multipole measurements in the presence of foregrounds. Phys. Rev. D 2006, 74, 023005. [Google Scholar] [CrossRef]
  101. Naselsky, P.D.; Verkhodanov, O.V.; Chiang, L.-Y.; Novikov, I.D. Phase analysis of the 1-year WMAP data and its application for the CMB foreground separation. ArXiv E-Prints 2003. arXiv:astro-ph/0310235. [Google Scholar]
  102. Lineweaver, C.H.; Tenorio, L.; Smoot, G.F.; Keegstra, P.; Banday, A.J.; Lubin, P. The dipole observed in the COBE DMR 4 year data. Astrophys. J. 1996, 470, 38. [Google Scholar] [CrossRef]
  103. Kashlinsky, A.; Atrio-Barandela, F. Measuring cosmological bulk flows via the kinematic Sunyaev-Zeldovich effect in the upcoming cosmic microwave background maps. Astrophys. J. 2000, 536, L67. [Google Scholar] [CrossRef] [PubMed]
  104. Kashlinsky, A.; Atrio-Barandela, F.; Ebeling, H.; Edge, A.; Kocevski, D. A new measurement of the bulk flow of X-ray luminous clusters of galaxies. Astrophys. J. 2010, 712, L81. [Google Scholar] [CrossRef]
  105. Atrio-Barandela, F.; Kashlinsky, A.; Ebeling, H.; Kocevski, D. Cosmic Microwave Background filters and the Dark-Flow measurement. ArXiv E-Prints 2013. arXiv:1211.4345. [Google Scholar]
  106. Ade, P.A.R.; Aghanim, N.; Arnaud, M.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Balbi, A.; Banday, A.J.; Barreiro, R.B.; Battaner, E.; et al. Planck intermediate results. XIII. Constraints on peculiar velocities. ArXiv E-Prints 2013. arXiv:1303.5090. [Google Scholar]
  107. Hand, N.; Addison, G.E.; Aubourg, E.; Battaglia, N.; Battistelli, E.S.; Bizyaev, D.; Bond, J.R.; Brewington, H.; Brinkmann, J.; Brown, B.R.; et al. Evidence of galaxy cluster motions with the kinematic Sunyaev-Zel’dovich effect. Phys. Rev. Lett. 2012, 109, 041101. [Google Scholar] [CrossRef] [PubMed]
  108. Keisler, R. The statistical significance of the “dark flow”. Astrophys. J. 2009, 707, L42. [Google Scholar] [CrossRef]
  109. Osborne, S.J.; Mak, D.S.Y.; Church, S.E.; Pierpaoli, E. Measuring the galaxy cluster bulk flow from WMAP data. Astrophys. J. 2011, 737, 98. [Google Scholar] [CrossRef]
  110. Mody, K.; Hajian, A. One thousand and one clusters: Measuring the bulk flow with the planck ESZ and X-ray selected galaxy cluster catalogs. Astrophys. J. 2012, 758, 4. [Google Scholar] [CrossRef]
  111. Atrio-Barandela, F. On the statistical significance of the bulk flow measured by the PLANCK satellite. ArXiv E-Prints 2013. arXiv:1303.6614. [Google Scholar]
  112. Turnbull, S.J.; Hudson, M.J.; Feldman, H.A.; Hicken, M.; Kirshner, R.P.; Watkins, R. Cosmic flows in the nearby universe from Type Ia supernovae. Mon. Not. R. Astron. Soc. 2012, 420, 447–454. [Google Scholar] [CrossRef]
  113. Dai, D.-C.; Kinney, W.H.; Stojkovic, D. Measuring the cosmological bulk flow using the peculiar velocities of supernovae. J. Cosmol. Astropart. Phys. 2011, 2011, 015. [Google Scholar] [CrossRef]
  114. Payez, A.; Cudell, J.R.; Hutsemekers, D. Axions and polarisation of quasars. AIP Conf. Proc. 2008, 1038, 211. [Google Scholar] [CrossRef]
  115. Tiwari, P.; Jain, P. Polarization Alignment in JVAS/CLASS flat spectrum radio surveys. Int. J. Mod. Phys. D 2013, 22, 1350089. [Google Scholar] [CrossRef]
  116. Payez, A. Cornering the axionlike particle explanation of quasar polarisations. Phys. Rev. D 2012, 85, 087701. [Google Scholar] [CrossRef]
  117. Payez, A.; Cudell, J.R.; Hutsemekers, D. Can axionlike particles explain the alignments of the polarisations of light from quasars? Phys. Rev. D 2011, 84, 085029. [Google Scholar] [CrossRef]
  118. Kalus, B.; Schwarz, D.J.; Seikel, M.; Wiegand, A. Constraints on anisotropic cosmic expansion from supernovae. Astron. Astrophys. 2013, 553, A56. [Google Scholar] [CrossRef]
  119. Colin, J.; Mohayaee, R.; Sarkar, S.; Shafieloo, A. Probing the anisotropic local universe and beyond with SNe Ia data. Mon. Not. R. Astron. Soc. 2011, 414, 264–271. [Google Scholar] [CrossRef]
  120. Schwarz, D.J.; Weinhorst, B. (An)isotropy of the Hubble diagram: Comparing hemispheres. Astron. Astrophys. 2007, 474, 717–729. [Google Scholar]
  121. Cai, R.-G.; Tuo, Z.-L. Direction dependence of the deceleration parameter. J. Cosmol. Astropart. Phys. 2012, 2012, 004. [Google Scholar] [CrossRef] [PubMed]
  122. Campanelli, L.; Cea, P.; Fogli, G.L.; Marrone, A. Testing the isotropy of the Universe with type Ia supernovae. Phys. Rev. D 2011, 83, 103503. [Google Scholar] [CrossRef]
  123. Antoniou, I.; Perivolaropoulos, L. Searching for a cosmological preferred axis: Union2 data analysis and comparison with other probes. J. Cosmol. Astropart. Phys. 2010, 2010, 012. [Google Scholar] [CrossRef]
  124. Zhao, W.; Wu, P.X.; Zhang, Y. Anisotropy of cosmic acceleration. Int. J. Mod. Phys. D 2013, 22, 1350060. [Google Scholar] [CrossRef]
  125. Cai, R.-G.; Ma, Y.-Z.; Tang, B.; Tuo, Z.-L. Constraining the anisotropic expansion of the Universe. Phys. Rev. D 2013, 87, 123522. [Google Scholar] [CrossRef]
  126. Blomqvist, M.; Mortsell, E.; Nobili, S. Probing dark energy inhomogeneities with supernovae. J. Cosmol. Astropart. Phys. 2008, 2008, 027. [Google Scholar] [CrossRef]
  127. Blomqvist, M.; Enander, J.; Mortsell, E. Constraining dark energy fluctuations with supernova correlations. J. Cosmol. Astropart. Phys. 2010, 2010, 018. [Google Scholar] [CrossRef]
  128. Cooke, R.; Lynden-Bell, D. Does the Universe accelerate equally in all directions? Mon. Not. R. Astron. Soc. 2010, 401, 1409–1414. [Google Scholar] [CrossRef]
  129. Gupta, S.; Saini, T.D. Direction dependence in supernova data: Constraining isotropy. Mon. Not. R. Astron. Soc. 2010, 407, 651–657. [Google Scholar] [CrossRef]
  130. Cameron, E.; Pettitt, T. On the evidence for cosmic variation of the fine structure constant (I): A Parametric Bayesian model selection analysis of the quasar dataset. ArXiv E-Prints 2013. arXiv:1207.6223. [Google Scholar]
  131. Kanekar, N.; Langston, G.I.; Stocke, J.T.; Carilli, C.L.; Menten, K.L. Constraining fundamental constant evolution with HI and OH lines. Astrophys. J. 2012, 746, L16. [Google Scholar] [CrossRef]
  132. Levshakov, S.A.; Combes, F.; Boone, F.; Agafonova, I.I.; Reimers, D.; Kozlov, M.G. An upper limit to the variation in the fundamental constants at redshift z = 5.2. Astron. Astrophys. 2012, 540, L9. [Google Scholar] [CrossRef]
  133. Marinoni, C.; Bel, J.; Buzzi, A. The scale of cosmic isotropy. J. Cosmol. Astropart. Phys. 2012, 2012, 036. [Google Scholar] [CrossRef]
  134. Schneider, D.P.; Richards, G.T.; Hall, P.B.; Strauss, M.A.; Anderson, S.F.; Boroson, T.A.; Ross, N.P.; Shen, Y.; Brandt, W.N.; Fan, X.; et al. The sloan digital sky survey quasar catalog. V. seventh data release. Astron. J. 2010, 139, 2360. [Google Scholar] [CrossRef]
  135. Nadathur, S. Seeing patterns in noise: Gigaparsec-scale “structures” that do not violate homogeneity. Mon. Not. R. Astron. Soc. 2013, 434, 398–406. [Google Scholar] [CrossRef]
  136. Pilipenko, S.; Malinovsky, A. Testing statistical significance of large quasar groups with sheets model of large scale structure. ArXiv E-Prints 2013. arXiv:1306.3970. [Google Scholar]
  137. Koivisto, T.; Mota, D.F. Accelerating cosmologies with an anisotropic equation of state. Astrophys. J. 2008, 679, 1. [Google Scholar] [CrossRef]
  138. Mota, D.F.; Kristiansen, J.R.; Koivisto, T.; Groeneboom, N.E. Constraining dark energy anisotropic stress. Mon. Not. R. Astron. Soc. 2007, 382, 793–800. [Google Scholar] [CrossRef]
  139. Koivisto, T.; Mota, D.F. Anisotropic dark energy: Dynamics of background and perturbations. J. Cosmol. Astropart. Phys. 2008, 2008, 018. [Google Scholar] [CrossRef]
  140. Lemaître, A.G. The expanding universe. Gen. Relativ. Gravit. 1997, 29, 641–680. [Google Scholar] [CrossRef]
  141. Tolman, R.C. Effect of imhomogeneity on cosmological models. Proc. Natl. Acad. Sci. USA 1934, 20, 169–176. [Google Scholar] [CrossRef] [PubMed]
  142. Bondi, H. Spherically symmetrical models in general relativity. Mon. Not. R. Astron. Soc. 1947, 107, 410–425. [Google Scholar] [CrossRef]
  143. Enqvist, K.; Mattsson, T. The effect of inhomogeneous expansion on the supernova observations. J. Cosmol. Astropart. Phys. 2007, 2007, 019. [Google Scholar] [CrossRef]
  144. Alnes, H.; Amarzguioui, M. The supernova Hubble diagram for off-center observers in a spherically symmetric inhomogeneous Universe. Phys. Rev. D 2007, 75, 023506. [Google Scholar] [CrossRef]
  145. Szekeres, P. Quasispherical Gravitational Collapse. Phys. Rev. D 1975, 12, 2941–2948. [Google Scholar] [CrossRef]
  146. Szekeres, P. A class of inhomogeneous cosmological models. Commun. Math. Phys. 1975, 41, 55–64. [Google Scholar] [CrossRef]
  147. Goode, S.W.; Wainwright, J. Singularities and evolution of the Szekeres cosmological models. Phys. Rev. D 1982, 26, 3315–3326. [Google Scholar] [CrossRef]
  148. Bonnor, W.B. Non-radiative solutions of Einstein’s equations for dust. Commun. Math. Phys. 1976, 51, 191–199. [Google Scholar] [CrossRef]
  149. Cho, I.; Vilenkin, A. Spacetime structure of an inflating global monopole. Phys. Rev. D 1997, 56, 7621–7626. [Google Scholar] [CrossRef]
  150. Bekenstein, J.D. Fine-structure constant: Is it really a constant? Phys. Rev. D 1982, 25, 1527–1539. [Google Scholar] [CrossRef]
  151. Sandvik, H.B.; Barrow, J.D.; Magueijo, J. A simple cosmology with a varying fine structure constant. Phys. Rev. Lett. 2002, 88, 031302. [Google Scholar] [CrossRef] [PubMed]
  152. Avgoustidis, A.; Martins, C.J.A.P.; Monteiro, A.M.R.V.L.; Vielzeuf, P.E.; Luzzi, G. Cosmological effects of scalar-photon couplings: Dark energy and varying-α models. ArXiv E-Prints 2013. arXiv:1305.7031. [Google Scholar]
  153. Sakai, N.; Shinkai, H.-A.; Tachizawa, T.; Maeda, K.-I. Dynamics of topological defects and inflation. Phys. Rev. D 1996, 53, 655–661. [Google Scholar] [CrossRef]
  154. Barriola, M.; Vilenkin, A. Gravitational field of a global monopole. Phys. Rev. Lett. 1989, 63, 341–343. [Google Scholar] [CrossRef] [PubMed]
  155. Balcerzak, A.; Dabrowski, M.P.; Denkiewicz, T. Off-center observers versus supernovae in inhomogeneous pressure universes. ArXiv E-Prints 2013. arXiv:1312.1567. [Google Scholar]
  156. Lineweaver, C.H. The CMB dipole: The most recent measurement and some history. ArXiv E-Prints 1996. arXiv:astro-ph/9609034. [Google Scholar]

Share and Cite

MDPI and ACS Style

Perivolaropoulos, L. Large Scale Cosmological Anomalies and Inhomogeneous Dark Energy. Galaxies 2014, 2, 22-61. https://doi.org/10.3390/galaxies2010022

AMA Style

Perivolaropoulos L. Large Scale Cosmological Anomalies and Inhomogeneous Dark Energy. Galaxies. 2014; 2(1):22-61. https://doi.org/10.3390/galaxies2010022

Chicago/Turabian Style

Perivolaropoulos, Leandros. 2014. "Large Scale Cosmological Anomalies and Inhomogeneous Dark Energy" Galaxies 2, no. 1: 22-61. https://doi.org/10.3390/galaxies2010022

APA Style

Perivolaropoulos, L. (2014). Large Scale Cosmological Anomalies and Inhomogeneous Dark Energy. Galaxies, 2(1), 22-61. https://doi.org/10.3390/galaxies2010022

Article Metrics

Back to TopTop