1 INTRODUCTION
The position of a face-on spiral galaxy compared to the position of an observer is reflected by the spin patterns of the galaxy, which can seem clockwise (Z) or counterclockwise (S) based on the location of the observer in the universe. Since the spin pattern is merely a matter of the position of the observer, in a sufficiently large population of spiral galaxies, galaxies with Z-wise spin patterns are expected to be identical in their physical properties to galaxies with S-wise spin patterns. However, previous experiments using data from the Sloan Digital Sky Survey (York et al. Reference York2000) show that spiral galaxies with Z-wise spin patterns in SDSS have photometric information that makes them distinguishable from galaxies with S-wise spin patterns (Shamir Reference Shamir2016b, Reference Shamir2017b).
First results (Shamir Reference Shamir2016b) using 13 440 galaxies from Galaxy Zoo 2 (Willett et al. Reference Willett2013) separated manually to Z-wise and S-wise galaxies showed that pattern recognition algorithms are able to predict the spin pattern of a galaxy by its photometric information in ~64% of the cases, significantly higher than mere chance of 50% (P<10−5). Similar results were obtained with 10 281 galaxies annotated automatically (Shamir Reference Shamir2016b). These experiments showed that Z-wise and S-wise spiral galaxies in SDSS are not photometrically indifferent.
Experiments with a larger dataset of SDSS galaxies classified automatically allowed the identification of individual photometric measurements with statistically significant differences between Z-wise and S-wise galaxies. These experiments provided evidence that Z-wise and S-wise galaxies are different by their color (Shamir Reference Shamir2013, Reference Shamir2017c), while exhibiting strong statistical significance for difference in their brightness (Shamir Reference Shamir2017b).
Since the experiments were done with galaxies acquired by one telescope (SDSS) and one image analysis method (Ganalyzer), it is possible that the results are driven by some unknown flaw in the telescope system, photometric pipeline, or automatic image analysis. This paper shows that the results are consistent across different telescopes and different image analysis methods, indicating the possibility that the reported results are not a reflection of a software or hardware flaw. The results also show that the asymmetry between Z-wise and S-wise galaxies changes with the direction of observation, and that opposite hemispheres have inverse asymmetry, suggesting that the asymmetry is aligned with a certain axis.
2 DATA
Several datasets were used to show the asymmetry between galaxies with Z-wise and S-wise spin patterns, taken from two different sky surveys – the Sloan Digital Sky Survey, and the Panoramic Survey Telescope and Rapid Response System (Hodapp et al. Reference Hodapp2004; Tonry et al. Reference Tonry2012).
An initial dataset of 3 · 106 SDSS DR8 galaxies with i magnitude lower than 18 and r Petrosian radius larger than 5.5 arcsec was separated automatically into spiral and elliptical galaxies (Kuminski & Shamir Reference Kuminski and Shamir2016). Then, 740 908 galaxies identified as spiral were separated automatically into Z-wise and S-wise galaxies by applying the Ganalyzer algorithm (Shamir Reference Shamir2011b, Reference Shamir2011a, Reference Shamir2016b).
Ganalyzer works by first transforming the galaxy image to its radial intensity plot (Shamir Reference Shamir2011b), which is a 360×35 image that reflects the brightness variations around the centre of the galaxy. The pixel (x,y) in the radial intensity plot is the median of the 5×5 pixels around (O x + sin(θ) · r, O y − cos(θ) · r) in the original galaxy image, where (O x , O y ) are the coordinates of the galaxy centre, θ is the polar angle (in degrees), and r is the radial distance (in pixels), spanning over 35% of the galaxy radius. That is, the x-coordinate in the radial intensity plot corresponds to the polar angle, and the y-coordinate corresponds to the radial distance from the galaxy centre in the original image (Shamir Reference Shamir2011b).
Then, each line in the radial intensity plot is searched for peaks, and the peaks are grouped by associating each detected peak (x 0, y) in horizontal line y with the nearest peak (x 1, y + 1) in the y+1 horizontal line. The x-coordinate of the two peaks are compared, and if x 1 < x 0 the counter L is incremented, while if x 1 > x 0 the counter R is incremented. If L⩾10 and L⩾3R then the galaxy is determined to have a Z-wise pattern, and if R⩾10 and R⩾3L the galaxy is determined to have a S-wise spin pattern. All other galaxies are assumed to have an undetermined spin pattern and are rejected from the analysis. Detailed description about the Ganalyzer algorithm is available in Shamir (Reference Shamir2011b), Hoehn & Shamir (Reference Hoehn and Shamir2014) and Shamir (Reference Shamir2012).
Applying the algorithm to the 740 908 spiral galaxies (Kuminski & Shamir Reference Kuminski and Shamir2016) provided a dataset of 162 516 galaxies with identified spin patterns (Shamir Reference Shamir2017a, Reference Shamir2017b, Reference Shamir2017c). The spin patterns of the remaining galaxies were not identified (e.g., edge-on galaxies), and were therefore excluded from the experiment. The mean r magnitude of the galaxies in the dataset is ~16.82 (σ ≃0.99), and the mean redshift measured using 10 281 of the galaxies that had spectra is 0.086. The galaxies are in the declination range of (−11.3°, 83.8°). Manual inspection of 400 randomly selected Z-wise galaxies and 400 randomly selected S-wise galaxies showed that 24 galaxies classified as Z-wise and 21 galaxies classified as S-wise did not have identifiable spin patterns.
Figure 1 displays the distribution of the Petrosian radius measured in the r band, the magnitude measured in the r band, and the redshift. Since most galaxies did not have spectra, the redshift is measured among a subset of 10 281 galaxies that had spectroscopic information. The figure shows the distribution of the galaxies with Z-wise patterns, S-wise patterns, and galaxies that were not assigned with an identifiable spin pattern and were excluded from the experiment.
The initial PanSTARRS data includes 3 053 831 celestial objects from data release 1 (Flewelling et al. Reference Flewelling2016)Footnote 1 . That list includes 2 394 452 objects that were identified by PanSTARRS photometric pipeline (Chambers et al. Reference Chambers2016; Tonry et al. Reference Tonry2012; Magnier et al. Reference Magnier2016; Waters et al. Reference Waters2016) as extended sources in all bands. The rest of the galaxies were 659 379 objects not included in the first set, and had r Petrosian radius larger than 5.5 arsec and their PSF magnitude subtracted by their Kron magnitude was larger than 0.05. The brightness of all objects was limited to 19th magnitude in the r band.
The images were downloaded as 120×120 JPG images using the PanSTARRS cutout service in a process that required 62 days to complete. That dataset was reduced by the Ganalyzer algorithm (Shamir Reference Shamir2011b) into 29 013 galaxies annotated by their spin pattern. The remaining galaxies not assigned with an identifiable spin pattern were excluded from the analysis. All galaxies were in the declination range of (−29.99°, 84.94°). Manual inspection of 200 galaxies showed that 11 galaxies did not have identifiable spin patterns. Figure 2 shows the r magnitude and Petrosian radius distribution of the galaxies in the PanSTARRS dataset. Only primary detections (rows with bestDetection=1) were used. For objects with multiple primary detections all detection were averaged. Values such as −9 999 were ignored, as these values are flags and not actual photometric measurements.
In addition to the automatically annotated datasets, 115 349 SDSS galaxies included in the catalog of SDSS spiral galaxies (Kuminski & Shamir Reference Kuminski and Shamir2016) and also had spectroscopic information (Strauss et al. Reference Strauss2002) in SDSS DR8 were annotated manually. Since human classification of spin patterns can be substantially biased (Land et al. Reference Land2008; Lintott et al. Reference Lintott2011), the galaxies were mirrored randomly to balance a possible preference of a certain pattern, and the classification was done by the author in a careful process that required ~450 h of labour. The outcome of the process was a dataset of 20 712 galaxies with Z-wise spin patterns and 20 027 S-wise galaxies. The mean r magnitude of these galaxies is ~16.6 (σ ≃ 0.84), and the mean redshift is ~0.09. Figure 3 shows the distribution of the r magnitude, Petrosian radius, and redshift among the galaxies classified as Z, galaxies classified as S, and galaxies that their spin pattern could not be determined.
The number of galaxies included in both the manually and the automatically classified SDSS datasets is 14 244. Examining the agreement between the two datasets showed that 6 881 galaxies classified manually as Z-wise were also classified automatically as Z-wise, while 276 galaxies classified manually as Z-wise were classified automatically as S-wise. Among the galaxies classified manually as S-wise, 6 775 were classified automatically as S-wise, and 312 were classified by the automatic classifier as Z-wise.
The population of galaxies imaged by SDSS and PanSTARRS is not distributed uniformly in the sky. Table 1 shows the number of galaxies in each of the datasets in each 30° RA range.
3 RESULTS
The galaxy population is distributed such that the right ascension range of (120°, 210°) contains 54 848, 24 281, and 11 508 galaxies in the automatically annotated SDSS dataset, the manually annotated SDSS dataset, and the PanSTARRS dataset, respectively. In the same RA range in the opposite hemisphere (< 30°∨ > 300°), the three datasets contained 51 053, 2 889, and 4 464 galaxies.
Tables 2 and 3 show the mean magnitude of Z-wise and S-wise galaxies in the (120°, 210°) hemisphere and in the opposite (< 30°∨ > 300°) hemispheres, respectively. The tables show the exponential and de Vaucouleurs magnitudes, which are two magnitude models commonly used for galaxies, and are measured by both the SDSS and PanSTARRS photometric pipelines.
While it has been shown that Z-wise and S-wise galaxies imaged by SDSS are photometrically different from each other (Shamir Reference Shamir2013, Reference Shamir2016b, Reference Shamir2017c, Reference Shamir2017b), Tables 2 and 3 show that the asymmetry is consistent across hemispheres, so that the asymmetry between Z-wise and S-wise galaxies in (120°, 210°) is inverse to the asymmetry in the same RA range in the opposite hemisphere (< 30°∨ > 300°). Other magnitude models (PSF magnitude, Kron magnitude, Petrosian magnitude) were tested and showed the same differences, indicating that the asymmetry is not the result of an error in a certain magnitude model.
All three datasets show with strong statistical power that galaxies with S-wise spin pattern in RA range (120°, 210°) are brighter than galaxies with Z-wise spin patterns. The RA range of (< 30°∨ > 300°) shows asymmetry with strong statistical significance observed in the automatically annotated SDSS galaxies. The manually annotated SDSS galaxies do not show statistically significant asymmetry in that RA range, possibly due to the lower number of galaxies, but it is aligned with the automatically annotated SDSS galaxies, showing that Z-wise galaxies in that dataset are brighter in that RA range compared to galaxies with S-wise patterns. The difference in the PanSTARRS galaxies in that RA range is statistically significant when using the one-tailed distribution, given the hypothesis that Z-wise galaxies in that sky region are brighter, as shown by the automatically annotated galaxies. For instance, the chance probability that Z-wise galaxies are brighter in the r exponential magnitude in both the manually annotated SDSS galaxies and the PanSTARRS dataset in the RA range (< 30°∨ > 300°) is ≃ 0.025.
Repeating the experiment after mirroring half of the galaxies randomly provided no statistically significant difference in any of the datasets. The dataset of manually annotated galaxies had spectroscopic information, showing no statistical significance between the mean redshift of the galaxies. The mean redshift of the Z-wise and S-wise galaxies in the RA range of (120°, 210°) is 0.089 ± 0.002 and 0.087 ± 0.002, respectively, and it is 0.082 ± 0.002 and 0.083 ± 0.002 for Z-wise and S-wise galaxies in the opposite hemisphere (< 30°∨ > 300°), showing no statistically significant difference in the redshift.
Figures 4 and 5 show the asymmetry measured using the exponential magnitude and the de Vaucouleurs magnitude, respectively, in the different right ascension ranges. The right ascension ranges of (270°, 300°) and (60°, 90°) in the SDSS manually annotated galaxies were not used due to the low number of galaxies in these ranges, 1 and 20, respectively. The figures show that the asymmetry changes with the right ascension also inside the hemispheres, and that the difference between the brightness of S-wise galaxies and Z-wise galaxies changes with the right ascension.
The tables and graphs show that the asymmetry is strongest at around the RA range of (120°, 210°), and the same RA range in the opposite hemisphere (< 30°∨ > 300°). Figures 6–8 show the histograms of the distribution of the galaxies in different exponential r magnitude ranges in the (120°, 210°) and (< 30°∨ > 300°) RA ranges.
3.1 Spin asymmetry in Galaxy Zoo 1 data
Galaxy Zoo (Lintott et al. Reference Lintott2011) is a citizen science campaign that aimed at classifying SDSS galaxies by their morphology, and also collected manual citizen science annotations of Z-wise and S-wise galaxies. Although the manual annotation of Z-wise and S-wise galaxies was found to be heavily biased (Hayes et al. Reference Hayes, Davis and Silva2017), a subset of galaxies classified by the citizen scientists with a high degree of 95% ‘superclean’ (Land et al. Reference Land2008) can be used to test whether the different magnitude between Z-wise and S-wise galaxies is also consistent when using galaxies annotated manually as part of Galaxy Zoo.
Comparing the Galaxy Zoo galaxies to the dataset described in Section 2 showed that 8 710 galaxies were included in both datasets (4 133 Z-wise galaxies and 4 577 S-wise galaxies). Of these galaxies, 1 199 of the galaxies classified as Z-wise in Galaxy Zoo were also classified as Z-wise in the automatically classified dataset, and 58 galaxies classified as Z-wise were classified automatically as S-wise. Among the galaxies classified by Galaxy Zoo as S-wise, 1 292 were classified automatically as S-wise and 66 as Z-wise.
When comparing to the dataset of manually classified galaxies, 2 045 Z-wise galaxies were also classified as Z-wise in the manually classified galaxies, and 14 galaxies classified as Z-wise were classified as S-wise. Among galaxies classified in the Galaxy Zoo 1 dataset as S-wise, 2 287 were classified as S-wise also in the dataset of manually classified galaxies, and 20 as Z-wise.
Table 4 shows the magnitude of Z-wise and S-wise galaxies measured with 5 923 galaxies (2 804 Z-wise and 3 119 S-wise) within the RA range of (120°, 210°) and had an agreement of more than 95% of the votes.
The table shows that the S-wise galaxies annotated by Galaxy Zoo are also generally brighter than the Z-wise galaxies, which is aligned with the asymmetry observed in the other datasets and shown in Tables 2 and 3. The differences are not statistically significant, with the highest t-test statistical significance of ~0.23 observed in the exponential Z-wise magnitude. It should be noted that the number of Galaxy Zoo galaxies with ‘superclean’ annotations is much lower than the number of galaxies in the other datasets, and therefore statistically significant difference is not expected. The same RA range in the opposite hemisphere had a much lower number of just 687 galaxies with ‘superclean’ agreement between the voters, and the low number does not allow a meaningful analysis.
3.2 Spin asymmetry in Galaxy Zoo 2 data
Galaxy Zoo 2 (Willett et al. Reference Willett2013) provided annotations of the morphology of SDSS galaxies, and a dataset of 13 440 manually annotated Galaxy Zoo 2 galaxies (Shamir Reference Shamir2016a) was used in to show that pattern recognition algorithms can identify between Z-wise and S-wise galaxies with accuracy significantly higher than mere chance (Shamir Reference Shamir2016b). The RA range (120°, 210°) contained 6 974 galaxies. Table 5 shows the magnitude of Z-wise and S-wise Galaxy Zoo 2 galaxies. As the table shows, Galaxy Zoo 2 galaxies show that S-wise galaxies in that RA range are brighter than Z-wise galaxies, and these results are aligned with the observations of the SDSS automatically classified galaxies and PanSTARRS. These results agree with the automatically annotated galaxies, but the statistical power of the difference is not expected to be high due to the much lower number of galaxies. The RA range (< 30°∨ > 300°) in that dataset contained no galaxies.
Comparing the Galaxy Zoo 2 galaxies to the galaxies described in Section 2 showed that 8 710 galaxies were included in both datasets (5 788 Z-wise galaxies and 5 328 S-wise galaxies). Of these galaxies, 1 736 of the galaxies classified as Z-wise in the Galaxy Zoo 2 dataset were also classified as Z-wise in the automatically classified dataset, and 87 galaxies classified as Z-wise were classified automatically as S-wise. Among the galaxies classified by Galaxy Zoo 2 as S-wise, 1 746 were classified automatically as S-wise and 70 as Z-wise.
When comparing to the manually classified galaxies, 2 703 Z-wise galaxies were also classified as Z-wise in the manually classified galaxies, and 41 as S-wise. Among galaxies classified in the Galaxy Zoo 2 dataset as S-wise 2 643 were classified as S-wise, and 44 as Z-wise.
3.3 Redshift dependence
The galaxies in the manually annotated dataset have spectra, allowing to profile the change in the asymmetry in different redshift ranges. Figure 9 shows the asymmetry measured using the exponential and de Vaucouleurs magnitudes in different redshift ranges in the RA range of (120°, 210°). The figure shows that the asymmetry decreases when the redshift gets higher, which can also be attributed to higher inaccuracy of the classification of more distant objects. The asymmetry increases when the redshift is higher than 0.2, but the increase is not statistically significant. Photometric redshift cannot be used as a valid measurement as it depends on the magnitude itself, so magnitude differences between Z-wise and S-wise galaxies would also affect the photometric redshift.
3.4 A possible spin asymmetry axis
The opposite asymmetry in opposite hemispheres indicates that the asymmetry could be aligned with an axis, such that galaxies with spin patterns oriented around that axis are brighter than galaxies with opposite spin patterns. To identify the approximate location of such possible axis, the SDSS automatically and manually annotated galaxies were combined into a single dataset, and the sky was separated into regions with half-width size of 10°, such that no two regions had overlapping parts of the sky. To avoid underpopulated areas, only regions with 3 000 galaxies or more were used. That provided a set of 28 non-overlapping sky regions, and the asymmetry in each sky region was computed.
Then, all possible integer right ascension and declination combinations (α, δ) were tested for cosine dependence by computing the Pearson correlation between vectors D and A, such that Di is the cosine of the angular distances between (α, δ) and the centre of region i, and Ai is the difference between the mean exponential magnitude of the Z-wise galaxies and the mean exponential magnitude of the S-wise galaxies in region i. That was repeated for all five bands, and the correlation of each (α, δ) is the mean correlation measured using all bands. Figure 10 shows the Pearson correlation coefficient measured from all possible integer (α, δ).
The highest Pearson correlation coefficient was ~0.61 observed at (α = 172°, δ = 50°). Repeating the analysis 1 000 times with randomly assigned galaxy spin patterns and selecting the maximum Pearson correlation in each run provided a mean of ~0.28 and standard deviation of ~0.11 (z-score of ~0.41). The correlation coefficient with the z-score of 0.41 from 0.61 is ~0.53, and therefore the 1σ error range of the right ascension is (132°, 224°), and for the declination it is (−26°, 74°). It should be noted that the most probable axis of (α = 172°, δ = 50°) is roughly aligned with the galactic pole (α = 192°, δ = 27°). The right ascension is also close to the right ascension of the CMB dipole (Aghanim et al. Reference Aghanim2014) at 166°, but there is substantial difference in the declination (50° compared to −27°).
The same experiment was repeated with the de Vaucouleurs, and Figure 11 shows the correlation coefficient at each possible integer right ascension and declination combinations. The highest Pearson correlation of ~0.61 is observed at (α = 175°, δ = 50°), which is in close agreement with the exponential magnitude.
4 DISCUSSION
The spin pattern of a spiral galaxy is a crude morphological feature that is sensitive to the perspective of the observer. It is therefore expected that in a sufficiently large population of galaxies the physical properties of galaxies with Z-wise spin patterns will not be different from the physical properties of galaxies with S-wise spin patterns. The observations presented in this paper show the existence of small but significant asymmetry between galaxies with opposite spin patterns. Moreover, the asymmetry in the differences between galaxies with opposite spin patterns is inverse in opposite hemispheres. That is, galaxies with S-wise spin patterns are significantly brighter in the RA range of (120°, 210°), but are significantly dimmer in the same RA range in the opposite hemisphere of (< 30°∨ > 300°). That inverse asymmetry can suggest the existence of an axis around which galaxy spin seem brighter to an Earth-based observer. Clearly, with mean redshift of ~0.09 the region analysed in this study is far larger than the size of a supercluster or any other known structure.
While the observation could be related to incompleteness of current cosmological theories (Kroupa Reference Kroupa2012), the proximity of the most probable axis to the galactic pole might suggest that the asymmetry could also be the effect of solar motion, and the asymmetry can be explained by relativistic beaming (Loeb & Gaudi Reference Loeb and Gaudi2003). With relativistic beaming, stars in a host galaxy that rotate in the same direction of the Sun would seem to an Earth-based observer brighter than stars in a host galaxy that spins in the opposite direction. Since the difference of the bolometric flux is by a fraction of $4\cdot \frac{v}{c}$ (Loeb & Gaudi Reference Loeb and Gaudi2003), assuming $\frac{v}{c}$ of the Sun of ~0.001 and a similar motion of the star in the observed galaxy, the maximum asymmetry can be ~0.8%. Since the asymmetry observed in the statistically significant galaxy population is not greater than 0.6%, and it peaks at around the galactic pole, it is possible that the observed asymmetry is a large-scale observation of relativistic beaming.
The difference in the velocity between the Sun and another Sun-like star in another face-on galaxy rotating in the opposite direction is 2 · 0.001 · c, or z ≃ 0.002. The K corrections for the different filters in that redshift are approximately 0.01, −0.07, −0.02, 0, and −0.03 for the u, g, r, i, and z filters, respectively (Chilingarian & Zolotukhin Reference Chilingarian and Zolotukhin2011). These corrections are far smaller than the observed asymmetry shown in this paper.
ACKNOWLEDGEMENTS
The research was supported in part by NSF grant IIS-1546079. I would like to thank Noah Brosch and Abraham Loeb for the fruitful discussions and kind assistance, as well as the associate editor Stephen Serjeant and the anonymous reviewer for their comments.
The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation.
Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/.
The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.