{"article_id": "1009.3123", "article_text": ["for about 20 years the problem of properties of short - term changes of solar activity has been considered extensively .", "many investigators studied the short - term periodicities of the various indices of solar activity .", "several periodicities were detected , but the periodicities about 155 days and from the interval of @xmath3 $ ] days ( @xmath4 $ ] years ) are mentioned most often .", "first of them was discovered by @xcite in the occurence rate of gamma - ray flares detected by the gamma - ray spectrometer aboard the _ solar maximum mission ( smm ) .", "this periodicity was confirmed for other solar flares data and for the same time period @xcite .", "it was also found in proton flares during solar cycles 19 and 20 @xcite , but it was not found in the solar flares data during solar cycles 22 @xcite .", "_    several autors confirmed above results for the daily sunspot area data . @xcite studied the sunspot data from 18741984 .", "she found the 155-day periodicity in data records from 31 years .", "this periodicity is always characteristic for one of the solar hemispheres ( the southern hemisphere for cycles 1215 and the northern hemisphere for cycles 1621 ) .", "moreover , it is only present during epochs of maximum activity ( in episodes of 13 years ) .", "similarinvestigationswerecarriedoutby + @xcite .", "they applied the same power spectrum method as lean , but the daily sunspot area data ( cycles 1221 ) were divided into 10 shorter time series .", "the periodicities were searched for the frequency interval 57115 nhz ( 100200 days ) and for each of 10 time series .", "the authors showed that the periodicity between 150160 days is statistically significant during all cycles from 16 to 21 .", "the considered peaks were remained unaltered after removing the 11-year cycle and applying the power spectrum analysis .", "@xcite used the wavelet technique for the daily sunspot areas between 1874 and 1993 .", "they determined the epochs of appearance of this periodicity and concluded that it presents around the maximum activity period in cycles 16 to 21 .", "moreover , the power of this periodicity started growing at cycle 19 , decreased in cycles 20 and 21 and disappered after cycle 21 .", "similaranalyseswerepresentedby + @xcite , but for sunspot number , solar wind plasma , interplanetary magnetic field and geomagnetic activity index @xmath5 .", "during 1964 - 2000 the sunspot number wavelet power of periods less than one year shows a cyclic evolution with the phase of the solar cycle.the 154-day period is prominent and its strenth is stronger around the 1982 - 1984 interval in almost all solar wind parameters .", "the existence of the 156-day periodicity in sunspot data were confirmed by @xcite .", "they considered the possible relation between the 475-day ( 1.3-year ) and 156-day periodicities .", "the 475-day ( 1.3-year ) periodicity was also detected in variations of the interplanetary magnetic field , geomagnetic activity helioseismic data and in the solar wind speed @xcite .", "@xcite concluded that the region of larger wavelet power shifts from 475-day ( 1.3-year ) period to 620-day ( 1.7-year ) period and then back to 475-day ( 1.3-year ) .", "the periodicities from the interval @xmath6 $ ] days ( @xmath4 $ ] years ) have been considered from 1968 .", "@xcite mentioned a 16.3-month ( 490-day ) periodicity in the sunspot numbers and in the geomagnetic data .", "@xcite analysed the occurrence rate of major flares during solar cycles 19 .", "they found a 18-month ( 540-day ) periodicity in flare rate of the norhern hemisphere .", "@xcite confirmed this result for the @xmath7 flare data for solar cycles 20 and 21 and found a peak in the power spectra near 510540 days .", "@xcite found a 17-month ( 510-day ) periodicity of sunspot groups and their areas from 1969 to 1986 .", "these authors concluded that the length of this period is variable and the reason of this periodicity is still not understood .", "@xcite and + @xcite obtained statistically significant peaks of power at around 158 days for daily sunspot data from 1923 - 1933 ( cycle 16 ) . in this paper the problem of the existence of this periodicity for sunspot data from cycle 16 is considered .", "the daily sunspot areas , the mean sunspot areas per carrington rotation , the monthly sunspot numbers and their fluctuations , which are obtained after removing the 11-year cycle are analysed . in section 2 the properties of the power spectrum methods are described . in section 3 a new approach to the problem of aliases in the power spectrum analysis", "is presented . in section 4 numerical results of the new method of the diagnosis of an echo - effect for sunspot area data are discussed . in section 5 the problem of the existence of the periodicity of about 155 days during the maximum activity period for sunspot data from the whole solar disk and from each solar hemisphere separately is considered .", "to find periodicities in a given time series the power spectrum analysis is applied . in this paper", "two methods are used : the fast fourier transformation algorithm with the hamming window function ( fft ) and the blackman - tukey ( bt ) power spectrum method @xcite .", "the bt method is used for the diagnosis of the reasons of the existence of peaks , which are obtained by the fft method .", "the bt method consists in the smoothing of a cosine transform of an autocorrelation function using a 3-point weighting average .", "such an estimator is consistent and unbiased .", "moreover , the peaks are uncorrelated and their sum is a variance of a considered time series . the main disadvantage of this method is a weak resolution of the periodogram points , particularly for low frequences .", "for example , if the autocorrelation function is evaluated for @xmath8 , then the distribution points in the time domain are : @xmath9 thus , it is obvious that this method should not be used for detecting low frequency periodicities with a fairly good resolution .", "however , because of an application of the autocorrelation function , the bt method can be used to verify a reality of peaks which are computed using a method giving the better resolution ( for example the fft method ) .", "it is valuable to remember that the power spectrum methods should be applied very carefully .", "the difficulties in the interpretation of significant peaks could be caused by at least four effects : a sampling of a continuos function , an echo - effect , a contribution of long - term periodicities and a random noise .", "first effect exists because periodicities , which are shorter than the sampling interval , may mix with longer periodicities . in result , this effect can be reduced by an decrease of the sampling interval between observations .", "the echo - effect occurs when there is a latent harmonic of frequency @xmath10 in the time series , giving a spectral peak at @xmath10 , and also periodic terms of frequency @xmath11 etc .", "this may be detected by the autocorrelation function for time series with a large variance .", "time series often contain long - term periodicities , that influence short - term peaks .", "they could rise periodogram s peaks at lower frequencies .", "however , it is also easy to notice the influence of the long - term periodicities on short - term peaks in the graphs of the autocorrelation functions .", "this effect is observed for the time series of solar activity indexes which are limited by the 11-year cycle .    to find statistically significant periodicities", "it is reasonable to use the autocorrelation function and the power spectrum method with a high resolution . in the case of a stationary time", "series they give similar results .", "moreover , for a stationary time series with the mean zero the fourier transform is equivalent to the cosine transform of an autocorrelation function @xcite .", "thus , after a comparison of a periodogram with an appropriate autocorrelation function one can detect peaks which are in the graph of the first function and do not exist in the graph of the second function .", "the reasons of their existence could be explained by the long - term periodicities and the echo - effect .", "below method enables one to detect these effects .", "( solid line ) and the 95% confidence level basing on thered noise ( dotted line ) .", "the periodogram values are presented on the left axis .", "the lower curve illustrates the autocorrelation function of the same time series ( solid line ) .", "the dotted lines represent two standard errors of the autocorrelation function .", "the dashed horizontal line shows the zero level .", "the autocorrelation values are shown in the right axis . ]     because the statistical tests indicate that the time series is a white noise the confidence level is not marked . ]    . ]", "the method of the diagnosis of an echo - effect in the power spectrum ( de ) consists in an analysis of a periodogram of a given time series computed using the bt method .", "the bt method bases on the cosine transform of the autocorrelation function which creates peaks which are in the periodogram , but not in the autocorrelation function .", "the de method is used for peaks which are computed by the fft method ( with high resolution ) and are statistically significant .", "the time series of sunspot activity indexes with the spacing interval one rotation or one month contain a markov - type persistence , which means a tendency for the successive values of the time series to remember their antecendent values .", "thus , i use a confidence level basing on the red noise of markov @xcite for the choice of the significant peaks of the periodogram computed by the fft method .", "when a time series does not contain the markov - type persistence i apply the fisher test and the kolmogorov - smirnov test at the significance level @xmath12 @xcite to verify a statistically significance of periodograms peaks . the fisher test checks the null hypothesis that the time series is white noise agains the alternative hypothesis that the time series contains an added deterministic periodic component of unspecified frequency . because the fisher test tends to be severe in rejecting peaks as insignificant the kolmogorov - smirnov test is also used .", "the de method analyses raw estimators of the power spectrum .", "they are given as follows    @xmath13    for @xmath14 + where @xmath15 for @xmath16 + @xmath17 is the length of the time series @xmath18 and @xmath19 is the mean value .", "the first term of the estimator @xmath20 is constant .", "the second term takes two values ( depending on odd or even @xmath21 ) which are not significant because @xmath22 for large m. thus , the third term of ( 1 ) should be analysed .", "looking for intervals of @xmath23 for which @xmath24 has the same sign and different signs one can find such parts of the function @xmath25 which create the value @xmath20 .", "let the set of values of the independent variable of the autocorrelation function be called @xmath26 and it can be divided into the sums of disjoint sets : @xmath27 where + @xmath28 + @xmath29 @xmath30 @xmath31 + @xmath32 + @xmath33 @xmath34 @xmath35 @xmath36 @xmath37 @xmath38", "@xmath39 @xmath40    well , the set @xmath41 contains all integer values of @xmath23 from the interval of @xmath42 for which the autocorrelation function and the cosinus function with the period @xmath43 $ ] are positive .", "the index @xmath44 indicates successive parts of the cosinus function for which the cosinuses of successive values of @xmath23 have the same sign .", "however , sometimes the set @xmath41 can be empty .", "for example , for @xmath45 and @xmath46 the set @xmath47 should contain all @xmath48 $ ] for which @xmath49 and @xmath50 , but for such values of @xmath23 the values of @xmath51 are negative .", "thus , the set @xmath47 is empty .    .", "the periodogram values are presented on the left axis .", "the lower curve illustrates the autocorrelation function of the same time series .", "the autocorrelation values are shown in the right axis . ]", "let us take into consideration all sets \\{@xmath52 } , \\{@xmath53 } and \\{@xmath41 } which are not empty . because numberings and power of these sets depend on the form of the autocorrelation function of the given time series , it is impossible to establish them arbitrary .", "thus , the sets of appropriate indexes of the sets \\{@xmath52 } , \\{@xmath53 } and \\{@xmath41 } are called @xmath54 , @xmath55 and @xmath56 respectively . for example", "the set @xmath56 contains all @xmath44 from the set @xmath57 for which the sets @xmath41 are not empty .", "to separate quantitatively in the estimator @xmath20 the positive contributions which are originated by the cases described by the formula ( 5 ) from the cases which are described by the formula ( 3 ) the following indexes are introduced : @xmath58 @xmath59 @xmath60 @xmath61 where @xmath62 @xmath63 @xmath64 taking for the empty sets \\{@xmath53 } and \\{@xmath41 } the indices @xmath65 and @xmath66 equal zero .", "the index @xmath65 describes a percentage of the contribution of the case when @xmath25 and @xmath51 are positive to the positive part of the third term of the sum ( 1 ) .", "the index @xmath66 describes a similar contribution , but for the case when the both @xmath25 and @xmath51 are simultaneously negative .", "thanks to these one can decide which the positive or the negative values of the autocorrelation function have a larger contribution to the positive values of the estimator @xmath20 .", "when the difference @xmath67 is positive , the statement the @xmath21-th peak really exists can not be rejected .", "thus , the following formula should be satisfied : @xmath68    because the @xmath21-th peak could exist as a result of the echo - effect , it is necessary to verify the second condition :    @xmath69\\in c_m.\\ ] ]    .", "the periodogram values are presented on the left axis .", "the lower curve illustrates the autocorrelation function of the same time series ( solid line ) .", "the dotted lines represent two standard errors of the autocorrelation function .", "the dashed horizontal line shows the zero level .", "the autocorrelation values are shown in the right axis . ]    to verify the implication ( 8) firstly it is necessary to evaluate the sets @xmath41 for @xmath70 of the values of @xmath23 for which the autocorrelation function and the cosine function with the period @xmath71 $ ] are positive and the sets @xmath72 of values of @xmath23 for which the autocorrelation function and the cosine function with the period @xmath43 $ ] are negative .", "secondly , a percentage of the contribution of the sum of products of positive values of @xmath25 and @xmath51 to the sum of positive products of the values of @xmath25 and @xmath51 should be evaluated . as a result the indexes @xmath65 for each set", "@xmath41 where @xmath44 is the index from the set @xmath56 are obtained .", "thirdly , from all sets @xmath41 such that @xmath70 the set @xmath73 for which the index @xmath65 is the greatest should be chosen .    the implication ( 8) is true when the set @xmath73 includes the considered period @xmath43 $ ] .", "this means that the greatest contribution of positive values of the autocorrelation function and positive cosines with the period @xmath43 $ ] to the periodogram value @xmath20 is caused by the sum of positive products of @xmath74 for each @xmath75-\\frac{m}{2k},[\\frac{2m}{k}]+\\frac{m}{2k})$ ] .    when the implication ( 8) is false , the peak @xmath20 is mainly created by the sum of positive products of @xmath74 for each @xmath76-\\frac{m}{2k},\\big [ \\frac{2m}{n}\\big ] + \\frac{m}{2k } \\big ) $ ] , where @xmath77 is a multiple or a divisor of @xmath21 .", "it is necessary to add , that the de method should be applied to the periodograms peaks , which probably exist because of the echo - effect .", "it enables one to find such parts of the autocorrelation function , which have the significant contribution to the considered peak .", "the fact , that the conditions ( 7 ) and ( 8) are satisfied , can unambiguously decide about the existence of the considered periodicity in the given time series , but if at least one of them is not satisfied , one can doubt about the existence of the considered periodicity .", "thus , in such cases the sentence the peak can not be treated as true should be used .    using the de method", "it is necessary to remember about the power of the set @xmath78 .", "if @xmath79 is too large , errors of an autocorrelation function estimation appear .", "they are caused by the finite length of the given time series and as a result additional peaks of the periodogram occur . if @xmath79 is too small , there are less peaks because of a low resolution of the periodogram . in applications", "@xmath80 is used . in order to evaluate the value", "@xmath79 the fft method is used .", "the periodograms computed by the bt and the fft method are compared .", "the conformity of them enables one to obtain the value @xmath79 .    .", "the fft periodogram values are presented on the left axis .", "the lower curve illustrates the bt periodogram of the same time series ( solid line and large black circles ) .", "the bt periodogram values are shown in the right axis . ]", "in this paper the sunspot activity data ( august 1923 - october 1933 ) provided by the greenwich photoheliographic results ( gpr ) are analysed .", "firstly , i consider the monthly sunspot number data . to eliminate the 11-year trend from these data ,", "the consecutively smoothed monthly sunspot number @xmath81 is subtracted from the monthly sunspot number @xmath82 where the consecutive mean @xmath83 is given by @xmath84 the values @xmath83 for @xmath85 and @xmath86 are calculated using additional data from last six months of cycle 15 and first six months of cycle 17 .    because of the north - south asymmetry of various solar indices @xcite , the sunspot activity is considered for each solar hemisphere separately .", "analogously to the monthly sunspot numbers , the time series of sunspot areas in the northern and southern hemispheres with the spacing interval @xmath87 rotation are denoted . in order to find periodicities ,", "the following time series are used : + @xmath88  ", "+ @xmath89", "   + @xmath90     + in the lower part of figure [ f1 ] the autocorrelation function of the time series for the northern hemisphere @xmath88 is shown .", "it is easy to notice that the prominent peak falls at 17 rotations interval ( 459 days ) and @xmath25 for @xmath91 $ ] rotations ( [ 81 , 162 ] days ) are significantly negative .", "the periodogram of the time series @xmath88 ( see the upper curve in figures [ f1 ] ) does not show the significant peaks at @xmath92 rotations ( 135 , 162 days ) , but there is the significant peak at @xmath93 ( 243 days ) .", "the peaks at @xmath94 are close to the peaks of the autocorrelation function .", "thus , the result obtained for the periodicity at about @xmath0 days are contradict to the results obtained for the time series of daily sunspot areas @xcite .    for the southern hemisphere ( the lower curve in figure [ f2 ] ) @xmath25 for @xmath95 $ ] rotations ( [ 54 , 189 ] days ) is not positive except @xmath96 ( 135 days ) for which @xmath97 is not statistically significant .", "the upper curve in figures [ f2 ] presents the periodogram of the time series @xmath89 .", "this time series does not contain a markov - type persistence .", "moreover , the kolmogorov - smirnov test and the fisher test do not reject a null hypothesis that the time series is a white noise only .", "this means that the time series do not contain an added deterministic periodic component of unspecified frequency .", "the autocorrelation function of the time series @xmath90 ( the lower curve in figure [ f3 ] ) has only one statistically significant peak for @xmath98 months ( 480 days ) and negative values for @xmath99 $ ] months ( [ 90 , 390 ] days ) .", "however , the periodogram of this time series ( the upper curve in figure [ f3 ] ) has two significant peaks the first at 15.2 and the second at 5.3 months ( 456 , 159 days ) .", "thus , the periodogram contains the significant peak , although the autocorrelation function has the negative value at @xmath100 months .    to explain", "these problems two following time series of daily sunspot areas are considered : + @xmath101  ", "+ @xmath102     + where @xmath103    the values @xmath104 for @xmath105 and @xmath106 are calculated using additional daily data from the solar cycles 15 and 17 .     and the cosine function for @xmath45 ( the period at about 154 days ) .", "the horizontal line ( dotted line ) shows the zero level .", "the vertical dotted lines evaluate the intervals where the sets @xmath107 ( for @xmath108 ) are searched .", "the percentage values show the index @xmath65 for each @xmath41 for the time series @xmath102 ( in parentheses for the time series @xmath101 ) . in the right bottom corner", "the values of @xmath65 for the time series @xmath102 , for @xmath109 are written . ]", "( the 500-day period ) ]    the comparison of the functions @xmath25 of the time series @xmath101 ( the lower curve in figure [ f4 ] ) and @xmath102 ( the lower curve in figure [ f5 ] ) suggests that the positive values of the function @xmath110 of the time series @xmath101 in the interval of @xmath111 $ ] days could be caused by the 11-year cycle .", "this effect is not visible in the case of periodograms of the both time series computed using the fft method ( see the upper curves in figures [ f4 ] and [ f5 ] ) or the bt method ( see the lower curve in figure [ f6 ] ) . moreover , the periodogram of the time series @xmath102 has the significant values at @xmath112 days , but the autocorrelation function is negative at these points .", "@xcite showed that the lomb - scargle periodograms for the both time series ( see @xcite , figures 7 a - c ) have a peak at 158.8 days which stands over the fap level by a significant amount . using the de method the above discrepancies are obvious . to establish the @xmath79 value the periodograms computed by the fft and", "the bt methods are shown in figure [ f6 ] ( the upper and the lower curve respectively ) .", "for @xmath46 and for periods less than 166 days there is a good comformity of the both periodograms ( but for periods greater than 166 days the points of the bt periodogram are not linked because the bt periodogram has much worse resolution than the fft periodogram ( no one know how to do it ) ) . for @xmath46 and @xmath113", "the value of @xmath21 is 13 ( @xmath71=153 $ ] ) .", "the inequality ( 7 ) is satisfied because @xmath114 .", "this means that the value of @xmath115 is mainly created by positive values of the autocorrelation function .", "the implication ( 8) needs an evaluation of the greatest value of the index @xmath65 where @xmath70 , but the solar data contain the most prominent period for @xmath116 days because of the solar rotation .", "thus , although @xmath117 for each @xmath118 , all sets @xmath41 ( see ( 5 ) and ( 6 ) ) without the set @xmath119 ( see ( 4 ) ) , which contains @xmath120 $ ] , are considered . this situation is presented in figure [ f7 ] . in this figure", "two curves @xmath121 and @xmath122 are plotted .", "the vertical dotted lines evaluate the intervals where the sets @xmath107 ( for @xmath123 ) are searched . for such @xmath41 two numbers", "are written : in parentheses the value of @xmath65 for the time series @xmath101 and above it the value of @xmath65 for the time series @xmath102 . to make this figure clear the curves are plotted for the set @xmath124 only .", "( in the right bottom corner information about the values of @xmath65 for the time series @xmath102 , for @xmath109 are written . )", "the implication ( 8) is not true , because @xmath125 for @xmath126 .", "therefore , @xmath43=153\\notin c_6=[423,500]$ ] .", "moreover , the autocorrelation function for @xmath127 $ ] is negative and the set @xmath128 is empty .", "thus , @xmath129 . on the basis of these information one can state , that the periodogram peak at @xmath130 days of the time series @xmath102 exists because of positive @xmath25 , but for @xmath23 from the intervals which do not contain this period .", "looking at the values of @xmath65 of the time series @xmath101 , one can notice that they decrease when @xmath23 increases until @xmath131 .", "this indicates , that when @xmath23 increases , the contribution of the 11-year cycle to the peaks of the periodogram decreases .", "an increase of the value of @xmath65 is for @xmath132 for the both time series , although the contribution of the 11-year cycle for the time series @xmath101 is insignificant .", "thus , this part of the autocorrelation function ( @xmath133 for the time series @xmath102 ) influences the @xmath21-th peak of the periodogram .", "this suggests that the periodicity at about 155 days is a harmonic of the periodicity from the interval of @xmath1 $ ] days .", "( solid line ) and consecutively smoothed sunspot areas of the one rotation time interval @xmath134 ( dotted line ) .", "both indexes are presented on the left axis .", "the lower curve illustrates fluctuations of the sunspot areas @xmath135 .", "the dotted and dashed horizontal lines represent levels zero and @xmath136 respectively .", "the fluctuations are shown on the right axis . ]", "the described reasoning can be carried out for other values of the periodogram .", "for example , the condition ( 8) is not satisfied for @xmath137 ( 250 , 222 , 200 days ) .", "moreover , the autocorrelation function at these points is negative .", "these suggest that there are not a true periodicity in the interval of [ 200 , 250 ] days .", "it is difficult to decide about the existence of the periodicities for @xmath138 ( 333 days ) and @xmath139 ( 286 days ) on the basis of above analysis . the implication ( 8) is not satisfied for @xmath139 and the condition ( 7 ) is not satisfied for @xmath138 , although the function @xmath25 of the time series @xmath102 is significantly positive for @xmath140 .", "the conditions ( 7 ) and ( 8) are satisfied for @xmath141 ( figure [ f8 ] ) and @xmath142 . therefore , it is possible to exist the periodicity from the interval of @xmath1 $ ] days .", "similar results were also obtained by @xcite for daily sunspot numbers and daily sunspot areas .", "she considered the means of three periodograms of these indexes for data from @xmath143 years and found statistically significant peaks from the interval of @xmath1 $ ] ( see @xcite , figure 2 ) .", "@xcite studied sunspot areas from 1876 - 1999 and sunspot numbers from 1749 - 2001 with the help of the wavelet transform .", "they pointed out that the 154 - 158-day period could be the third harmonic of the 1.3-year ( 475-day ) period .", "moreover , the both periods fluctuate considerably with time , being stronger during stronger sunspot cycles .", "therefore , the wavelet analysis suggests a common origin of the both periodicities . this conclusion confirms the de method result which indicates that the periodogram peak at @xmath144 days is an alias of the periodicity from the interval of @xmath1 $ ]", "in order to verify the existence of the periodicity at about 155 days i consider the following time series : + @xmath145     + @xmath146", "   + @xmath147  ", "+ the value @xmath134 is calculated analogously to @xmath83 ( see sect .", "the values @xmath148 and @xmath149 are evaluated from the formula ( 9 ) . in the upper part of figure [ f9 ] the time series of sunspot areas @xmath150 of the one rotation time interval from the whole solar disk and the time series of consecutively smoothed sunspot areas @xmath151", "are showed . in the lower part of figure [ f9 ]", "the time series of sunspot area fluctuations @xmath145 is presented .", "on the basis of these data the maximum activity period of cycle 16 is evaluated .", "it is an interval between two strongest fluctuations e.a .", "@xmath152 $ ] rotations .", "the length of the time interval @xmath153 is 54 rotations .", "if the about @xmath0-day ( 6 solar rotations ) periodicity existed in this time interval and it was characteristic for strong fluctuations from this time interval , 10 local maxima in the set of @xmath154 would be seen .", "then it should be necessary to find such a value of p for which @xmath155 for @xmath156 and the number of the local maxima of these values is 10 .", "as it can be seen in the lower part of figure [ f9 ] this is for the case of @xmath157 ( in this figure the dashed horizontal line is the level of @xmath158 ) .", "figure [ f10 ] presents nine time distances among the successive fluctuation local maxima and the horizontal line represents the 6-rotation periodicity .", "it is immediately apparent that the dispersion of these points is 10 and it is difficult to find even few points which oscillate around the value of 6 .", "such an analysis was carried out for smaller and larger @xmath136 and the results were similar .", "therefore , the fact , that the about @xmath0-day periodicity exists in the time series of sunspot area fluctuations during the maximum activity period is questionable .    .", "the horizontal line represents the 6-rotation ( 162-day ) period . ]    ]    ]    to verify again the existence of the about @xmath0-day periodicity during the maximum activity period in each solar hemisphere separately , the time series @xmath88 and @xmath89 were also cut down to the maximum activity period ( january 1925december 1930 ) .", "the comparison of the autocorrelation functions of these time series with the appriopriate autocorrelation functions of the time series @xmath88 and @xmath89 , which are computed for the whole 11-year cycle ( the lower curves of figures [ f1 ] and [ f2 ] ) , indicates that there are not significant differences between them especially for @xmath23=5 and 6 rotations ( 135 and 162 days ) ) .", "this conclusion is confirmed by the analysis of the time series @xmath146 for the maximum activity period .", "the autocorrelation function ( the lower curve of figure [ f11 ] ) is negative for the interval of [ 57 , 173 ] days , but the resolution of the periodogram is too low to find the significant peak at @xmath159 days .", "the autocorrelation function gives the same result as for daily sunspot area fluctuations from the whole solar disk ( @xmath160 ) ( see also the lower curve of figures [ f5 ] ) . in the case of", "the time series @xmath89 @xmath161 is zero for the fluctuations from the whole solar cycle and it is almost zero ( @xmath162 ) for the fluctuations from the maximum activity period .", "the value @xmath163 is negative .", "similarly to the case of the northern hemisphere the autocorrelation function and the periodogram of southern hemisphere daily sunspot area fluctuations from the maximum activity period @xmath147 are computed ( see figure [ f12 ] ) .", "the autocorrelation function has the statistically significant positive peak in the interval of [ 155 , 165 ] days , but the periodogram has too low resolution to decide about the possible periodicities .", "the correlative analysis indicates that there are positive fluctuations with time distances about @xmath0 days in the maximum activity period .", "the results of the analyses of the time series of sunspot area fluctuations from the maximum activity period are contradict with the conclusions of @xcite .", "she uses the power spectrum analysis only .", "the periodogram of daily sunspot fluctuations contains peaks , which could be harmonics or subharmonics of the true periodicities .", "they could be treated as real periodicities .", "this effect is not visible for sunspot data of the one rotation time interval , but averaging could lose true periodicities .", "this is observed for data from the southern hemisphere .", "there is the about @xmath0-day peak in the autocorrelation function of daily fluctuations , but the correlation for data of the one rotation interval is almost zero or negative at the points @xmath164 and 6 rotations .", "thus , it is reasonable to research both time series together using the correlative and the power spectrum analyses .", "the following results are obtained :    1 .", "a new method of the detection of statistically significant peaks of the periodograms enables one to identify aliases in the periodogram .", "2 .   two effects cause the existence of the peak of the periodogram of the time series of sunspot area fluctuations at about @xmath0 days : the first is caused by the 27-day periodicity , which probably creates the 162-day periodicity ( it is a subharmonic frequency of the 27-day periodicity ) and the second is caused by statistically significant positive values of the autocorrelation function from the intervals of @xmath165 $ ] and @xmath166 $ ] days .", "the existence of the periodicity of about @xmath0 days of the time series of sunspot area fluctuations and sunspot area fluctuations from the northern hemisphere during the maximum activity period is questionable .", "the autocorrelation analysis of the time series of sunspot area fluctuations from the southern hemisphere indicates that the periodicity of about 155 days exists during the maximum activity period .", "i appreciate valuable comments from professor j. jakimiec ."], "abstract_text": ["<S> the short - term periodicities of the daily sunspot area fluctuations from august 1923 to october 1933 are discussed . for these data </S>", "<S> the correlative analysis indicates negative correlation for the periodicity of about @xmath0 days , but the power spectrum analysis indicates a statistically significant peak in this time interval . </S>", "<S> a new method of the diagnosis of an echo - effect in spectrum is proposed and it is stated that the 155-day periodicity is a harmonic of the periodicities from the interval of @xmath1 $ ] days .    the autocorrelation functions for the daily sunspot area fluctuations and for the fluctuations of the one rotation time interval in the northern hemisphere , separately for the whole solar cycle 16 and for the maximum activity period of this cycle do not show differences , especially in the interval of @xmath2 $ ] days . </S>", "<S> it proves against the thesis of the existence of strong positive fluctuations of the about @xmath0-day interval in the maximum activity period of the solar cycle 16 in the northern hemisphere . </S>", "<S> however , a similar analysis for data from the southern hemisphere indicates that there is the periodicity of about @xmath0 days in sunspot area data in the maximum activity period of the cycle 16 only . </S>"], "labels": null, "section_names": ["introduction", "methods of periodicity analysis", "a method of the diagnosis of an echo-effect in the power spectrum", "data analysis", "the periodicity at about 155 days during the maximum activity period", "conclusion", "acknowledgments"], "sections": [["for about 20 years the problem of properties of short - term changes of solar activity has been considered extensively .", "many investigators studied the short - term periodicities of the various indices of solar activity .", "several periodicities were detected , but the periodicities about 155 days and from the interval of @xmath3 $ ] days ( @xmath4 $ ] years ) are mentioned most often .", "first of them was discovered by @xcite in the occurence rate of gamma - ray flares detected by the gamma - ray spectrometer aboard the _ solar maximum mission ( smm ) .", "this periodicity was confirmed for other solar flares data and for the same time period @xcite .", "it was also found in proton flares during solar cycles 19 and 20 @xcite , but it was not found in the solar flares data during solar cycles 22 @xcite .", "_    several autors confirmed above results for the daily sunspot area data . @xcite studied the sunspot data from 18741984 .", "she found the 155-day periodicity in data records from 31 years .", "this periodicity is always characteristic for one of the solar hemispheres ( the southern hemisphere for cycles 1215 and the northern hemisphere for cycles 1621 ) .", "moreover , it is only present during epochs of maximum activity ( in episodes of 13 years ) .", "similarinvestigationswerecarriedoutby + @xcite .", "they applied the same power spectrum method as lean , but the daily sunspot area data ( cycles 1221 ) were divided into 10 shorter time series .", "the periodicities were searched for the frequency interval 57115 nhz ( 100200 days ) and for each of 10 time series .", "the authors showed that the periodicity between 150160 days is statistically significant during all cycles from 16 to 21 .", "the considered peaks were remained unaltered after removing the 11-year cycle and applying the power spectrum analysis .", "@xcite used the wavelet technique for the daily sunspot areas between 1874 and 1993 .", "they determined the epochs of appearance of this periodicity and concluded that it presents around the maximum activity period in cycles 16 to 21 .", "moreover , the power of this periodicity started growing at cycle 19 , decreased in cycles 20 and 21 and disappered after cycle 21 .", "similaranalyseswerepresentedby + @xcite , but for sunspot number , solar wind plasma , interplanetary magnetic field and geomagnetic activity index @xmath5 .", "during 1964 - 2000 the sunspot number wavelet power of periods less than one year shows a cyclic evolution with the phase of the solar cycle.the 154-day period is prominent and its strenth is stronger around the 1982 - 1984 interval in almost all solar wind parameters .", "the existence of the 156-day periodicity in sunspot data were confirmed by @xcite .", "they considered the possible relation between the 475-day ( 1.3-year ) and 156-day periodicities .", "the 475-day ( 1.3-year ) periodicity was also detected in variations of the interplanetary magnetic field , geomagnetic activity helioseismic data and in the solar wind speed @xcite .", "@xcite concluded that the region of larger wavelet power shifts from 475-day ( 1.3-year ) period to 620-day ( 1.7-year ) period and then back to 475-day ( 1.3-year ) .", "the periodicities from the interval @xmath6 $ ] days ( @xmath4 $ ] years ) have been considered from 1968 .", "@xcite mentioned a 16.3-month ( 490-day ) periodicity in the sunspot numbers and in the geomagnetic data .", "@xcite analysed the occurrence rate of major flares during solar cycles 19 .", "they found a 18-month ( 540-day ) periodicity in flare rate of the norhern hemisphere .", "@xcite confirmed this result for the @xmath7 flare data for solar cycles 20 and 21 and found a peak in the power spectra near 510540 days .", "@xcite found a 17-month ( 510-day ) periodicity of sunspot groups and their areas from 1969 to 1986 .", "these authors concluded that the length of this period is variable and the reason of this periodicity is still not understood .", "@xcite and + @xcite obtained statistically significant peaks of power at around 158 days for daily sunspot data from 1923 - 1933 ( cycle 16 ) . in this paper the problem of the existence of this periodicity for sunspot data from cycle 16 is considered .", "the daily sunspot areas , the mean sunspot areas per carrington rotation , the monthly sunspot numbers and their fluctuations , which are obtained after removing the 11-year cycle are analysed . in section 2 the properties of the power spectrum methods are described . in section 3 a new approach to the problem of aliases in the power spectrum analysis", "is presented . in section 4 numerical results of the new method of the diagnosis of an echo - effect for sunspot area data are discussed . in section 5 the problem of the existence of the periodicity of about 155 days during the maximum activity period for sunspot data from the whole solar disk and from each solar hemisphere separately is considered ."], ["to find periodicities in a given time series the power spectrum analysis is applied . in this paper", "two methods are used : the fast fourier transformation algorithm with the hamming window function ( fft ) and the blackman - tukey ( bt ) power spectrum method @xcite .", "the bt method is used for the diagnosis of the reasons of the existence of peaks , which are obtained by the fft method .", "the bt method consists in the smoothing of a cosine transform of an autocorrelation function using a 3-point weighting average .", "such an estimator is consistent and unbiased .", "moreover , the peaks are uncorrelated and their sum is a variance of a considered time series . the main disadvantage of this method is a weak resolution of the periodogram points , particularly for low frequences .", "for example , if the autocorrelation function is evaluated for @xmath8 , then the distribution points in the time domain are : @xmath9 thus , it is obvious that this method should not be used for detecting low frequency periodicities with a fairly good resolution .", "however , because of an application of the autocorrelation function , the bt method can be used to verify a reality of peaks which are computed using a method giving the better resolution ( for example the fft method ) .", "it is valuable to remember that the power spectrum methods should be applied very carefully .", "the difficulties in the interpretation of significant peaks could be caused by at least four effects : a sampling of a continuos function , an echo - effect , a contribution of long - term periodicities and a random noise .", "first effect exists because periodicities , which are shorter than the sampling interval , may mix with longer periodicities . in result , this effect can be reduced by an decrease of the sampling interval between observations .", "the echo - effect occurs when there is a latent harmonic of frequency @xmath10 in the time series , giving a spectral peak at @xmath10 , and also periodic terms of frequency @xmath11 etc .", "this may be detected by the autocorrelation function for time series with a large variance .", "time series often contain long - term periodicities , that influence short - term peaks .", "they could rise periodogram s peaks at lower frequencies .", "however , it is also easy to notice the influence of the long - term periodicities on short - term peaks in the graphs of the autocorrelation functions .", "this effect is observed for the time series of solar activity indexes which are limited by the 11-year cycle .    to find statistically significant periodicities", "it is reasonable to use the autocorrelation function and the power spectrum method with a high resolution . in the case of a stationary time", "series they give similar results .", "moreover , for a stationary time series with the mean zero the fourier transform is equivalent to the cosine transform of an autocorrelation function @xcite .", "thus , after a comparison of a periodogram with an appropriate autocorrelation function one can detect peaks which are in the graph of the first function and do not exist in the graph of the second function .", "the reasons of their existence could be explained by the long - term periodicities and the echo - effect .", "below method enables one to detect these effects .", "( solid line ) and the 95% confidence level basing on thered noise ( dotted line ) .", "the periodogram values are presented on the left axis .", "the lower curve illustrates the autocorrelation function of the same time series ( solid line ) .", "the dotted lines represent two standard errors of the autocorrelation function .", "the dashed horizontal line shows the zero level .", "the autocorrelation values are shown in the right axis . ]     because the statistical tests indicate that the time series is a white noise the confidence level is not marked . ]    . ]"], ["the method of the diagnosis of an echo - effect in the power spectrum ( de ) consists in an analysis of a periodogram of a given time series computed using the bt method .", "the bt method bases on the cosine transform of the autocorrelation function which creates peaks which are in the periodogram , but not in the autocorrelation function .", "the de method is used for peaks which are computed by the fft method ( with high resolution ) and are statistically significant .", "the time series of sunspot activity indexes with the spacing interval one rotation or one month contain a markov - type persistence , which means a tendency for the successive values of the time series to remember their antecendent values .", "thus , i use a confidence level basing on the red noise of markov @xcite for the choice of the significant peaks of the periodogram computed by the fft method .", "when a time series does not contain the markov - type persistence i apply the fisher test and the kolmogorov - smirnov test at the significance level @xmath12 @xcite to verify a statistically significance of periodograms peaks . the fisher test checks the null hypothesis that the time series is white noise agains the alternative hypothesis that the time series contains an added deterministic periodic component of unspecified frequency . because the fisher test tends to be severe in rejecting peaks as insignificant the kolmogorov - smirnov test is also used .", "the de method analyses raw estimators of the power spectrum .", "they are given as follows    @xmath13    for @xmath14 + where @xmath15 for @xmath16 + @xmath17 is the length of the time series @xmath18 and @xmath19 is the mean value .", "the first term of the estimator @xmath20 is constant .", "the second term takes two values ( depending on odd or even @xmath21 ) which are not significant because @xmath22 for large m. thus , the third term of ( 1 ) should be analysed .", "looking for intervals of @xmath23 for which @xmath24 has the same sign and different signs one can find such parts of the function @xmath25 which create the value @xmath20 .", "let the set of values of the independent variable of the autocorrelation function be called @xmath26 and it can be divided into the sums of disjoint sets : @xmath27 where + @xmath28 + @xmath29 @xmath30 @xmath31 + @xmath32 + @xmath33 @xmath34 @xmath35 @xmath36 @xmath37 @xmath38", "@xmath39 @xmath40    well , the set @xmath41 contains all integer values of @xmath23 from the interval of @xmath42 for which the autocorrelation function and the cosinus function with the period @xmath43 $ ] are positive .", "the index @xmath44 indicates successive parts of the cosinus function for which the cosinuses of successive values of @xmath23 have the same sign .", "however , sometimes the set @xmath41 can be empty .", "for example , for @xmath45 and @xmath46 the set @xmath47 should contain all @xmath48 $ ] for which @xmath49 and @xmath50 , but for such values of @xmath23 the values of @xmath51 are negative .", "thus , the set @xmath47 is empty .    .", "the periodogram values are presented on the left axis .", "the lower curve illustrates the autocorrelation function of the same time series .", "the autocorrelation values are shown in the right axis . ]", "let us take into consideration all sets \\{@xmath52 } , \\{@xmath53 } and \\{@xmath41 } which are not empty . because numberings and power of these sets depend on the form of the autocorrelation function of the given time series , it is impossible to establish them arbitrary .", "thus , the sets of appropriate indexes of the sets \\{@xmath52 } , \\{@xmath53 } and \\{@xmath41 } are called @xmath54 , @xmath55 and @xmath56 respectively . for example", "the set @xmath56 contains all @xmath44 from the set @xmath57 for which the sets @xmath41 are not empty .", "to separate quantitatively in the estimator @xmath20 the positive contributions which are originated by the cases described by the formula ( 5 ) from the cases which are described by the formula ( 3 ) the following indexes are introduced : @xmath58 @xmath59 @xmath60 @xmath61 where @xmath62 @xmath63 @xmath64 taking for the empty sets \\{@xmath53 } and \\{@xmath41 } the indices @xmath65 and @xmath66 equal zero .", "the index @xmath65 describes a percentage of the contribution of the case when @xmath25 and @xmath51 are positive to the positive part of the third term of the sum ( 1 ) .", "the index @xmath66 describes a similar contribution , but for the case when the both @xmath25 and @xmath51 are simultaneously negative .", "thanks to these one can decide which the positive or the negative values of the autocorrelation function have a larger contribution to the positive values of the estimator @xmath20 .", "when the difference @xmath67 is positive , the statement the @xmath21-th peak really exists can not be rejected .", "thus , the following formula should be satisfied : @xmath68    because the @xmath21-th peak could exist as a result of the echo - effect , it is necessary to verify the second condition :    @xmath69\\in c_m.\\ ] ]    .", "the periodogram values are presented on the left axis .", "the lower curve illustrates the autocorrelation function of the same time series ( solid line ) .", "the dotted lines represent two standard errors of the autocorrelation function .", "the dashed horizontal line shows the zero level .", "the autocorrelation values are shown in the right axis . ]    to verify the implication ( 8) firstly it is necessary to evaluate the sets @xmath41 for @xmath70 of the values of @xmath23 for which the autocorrelation function and the cosine function with the period @xmath71 $ ] are positive and the sets @xmath72 of values of @xmath23 for which the autocorrelation function and the cosine function with the period @xmath43 $ ] are negative .", "secondly , a percentage of the contribution of the sum of products of positive values of @xmath25 and @xmath51 to the sum of positive products of the values of @xmath25 and @xmath51 should be evaluated . as a result the indexes @xmath65 for each set", "@xmath41 where @xmath44 is the index from the set @xmath56 are obtained .", "thirdly , from all sets @xmath41 such that @xmath70 the set @xmath73 for which the index @xmath65 is the greatest should be chosen .    the implication ( 8) is true when the set @xmath73 includes the considered period @xmath43 $ ] .", "this means that the greatest contribution of positive values of the autocorrelation function and positive cosines with the period @xmath43 $ ] to the periodogram value @xmath20 is caused by the sum of positive products of @xmath74 for each @xmath75-\\frac{m}{2k},[\\frac{2m}{k}]+\\frac{m}{2k})$ ] .    when the implication ( 8) is false , the peak @xmath20 is mainly created by the sum of positive products of @xmath74 for each @xmath76-\\frac{m}{2k},\\big [ \\frac{2m}{n}\\big ] + \\frac{m}{2k } \\big ) $ ] , where @xmath77 is a multiple or a divisor of @xmath21 .", "it is necessary to add , that the de method should be applied to the periodograms peaks , which probably exist because of the echo - effect .", "it enables one to find such parts of the autocorrelation function , which have the significant contribution to the considered peak .", "the fact , that the conditions ( 7 ) and ( 8) are satisfied , can unambiguously decide about the existence of the considered periodicity in the given time series , but if at least one of them is not satisfied , one can doubt about the existence of the considered periodicity .", "thus , in such cases the sentence the peak can not be treated as true should be used .    using the de method", "it is necessary to remember about the power of the set @xmath78 .", "if @xmath79 is too large , errors of an autocorrelation function estimation appear .", "they are caused by the finite length of the given time series and as a result additional peaks of the periodogram occur . if @xmath79 is too small , there are less peaks because of a low resolution of the periodogram . in applications", "@xmath80 is used . in order to evaluate the value", "@xmath79 the fft method is used .", "the periodograms computed by the bt and the fft method are compared .", "the conformity of them enables one to obtain the value @xmath79 .    .", "the fft periodogram values are presented on the left axis .", "the lower curve illustrates the bt periodogram of the same time series ( solid line and large black circles ) .", "the bt periodogram values are shown in the right axis . ]"], ["in this paper the sunspot activity data ( august 1923 - october 1933 ) provided by the greenwich photoheliographic results ( gpr ) are analysed .", "firstly , i consider the monthly sunspot number data . to eliminate the 11-year trend from these data ,", "the consecutively smoothed monthly sunspot number @xmath81 is subtracted from the monthly sunspot number @xmath82 where the consecutive mean @xmath83 is given by @xmath84 the values @xmath83 for @xmath85 and @xmath86 are calculated using additional data from last six months of cycle 15 and first six months of cycle 17 .    because of the north - south asymmetry of various solar indices @xcite , the sunspot activity is considered for each solar hemisphere separately .", "analogously to the monthly sunspot numbers , the time series of sunspot areas in the northern and southern hemispheres with the spacing interval @xmath87 rotation are denoted . in order to find periodicities ,", "the following time series are used : + @xmath88  ", "+ @xmath89", "   + @xmath90     + in the lower part of figure [ f1 ] the autocorrelation function of the time series for the northern hemisphere @xmath88 is shown .", "it is easy to notice that the prominent peak falls at 17 rotations interval ( 459 days ) and @xmath25 for @xmath91 $ ] rotations ( [ 81 , 162 ] days ) are significantly negative .", "the periodogram of the time series @xmath88 ( see the upper curve in figures [ f1 ] ) does not show the significant peaks at @xmath92 rotations ( 135 , 162 days ) , but there is the significant peak at @xmath93 ( 243 days ) .", "the peaks at @xmath94 are close to the peaks of the autocorrelation function .", "thus , the result obtained for the periodicity at about @xmath0 days are contradict to the results obtained for the time series of daily sunspot areas @xcite .    for the southern hemisphere ( the lower curve in figure [ f2 ] ) @xmath25 for @xmath95 $ ] rotations ( [ 54 , 189 ] days ) is not positive except @xmath96 ( 135 days ) for which @xmath97 is not statistically significant .", "the upper curve in figures [ f2 ] presents the periodogram of the time series @xmath89 .", "this time series does not contain a markov - type persistence .", "moreover , the kolmogorov - smirnov test and the fisher test do not reject a null hypothesis that the time series is a white noise only .", "this means that the time series do not contain an added deterministic periodic component of unspecified frequency .", "the autocorrelation function of the time series @xmath90 ( the lower curve in figure [ f3 ] ) has only one statistically significant peak for @xmath98 months ( 480 days ) and negative values for @xmath99 $ ] months ( [ 90 , 390 ] days ) .", "however , the periodogram of this time series ( the upper curve in figure [ f3 ] ) has two significant peaks the first at 15.2 and the second at 5.3 months ( 456 , 159 days ) .", "thus , the periodogram contains the significant peak , although the autocorrelation function has the negative value at @xmath100 months .    to explain", "these problems two following time series of daily sunspot areas are considered : + @xmath101  ", "+ @xmath102     + where @xmath103    the values @xmath104 for @xmath105 and @xmath106 are calculated using additional daily data from the solar cycles 15 and 17 .     and the cosine function for @xmath45 ( the period at about 154 days ) .", "the horizontal line ( dotted line ) shows the zero level .", "the vertical dotted lines evaluate the intervals where the sets @xmath107 ( for @xmath108 ) are searched .", "the percentage values show the index @xmath65 for each @xmath41 for the time series @xmath102 ( in parentheses for the time series @xmath101 ) . in the right bottom corner", "the values of @xmath65 for the time series @xmath102 , for @xmath109 are written . ]", "( the 500-day period ) ]    the comparison of the functions @xmath25 of the time series @xmath101 ( the lower curve in figure [ f4 ] ) and @xmath102 ( the lower curve in figure [ f5 ] ) suggests that the positive values of the function @xmath110 of the time series @xmath101 in the interval of @xmath111 $ ] days could be caused by the 11-year cycle .", "this effect is not visible in the case of periodograms of the both time series computed using the fft method ( see the upper curves in figures [ f4 ] and [ f5 ] ) or the bt method ( see the lower curve in figure [ f6 ] ) . moreover , the periodogram of the time series @xmath102 has the significant values at @xmath112 days , but the autocorrelation function is negative at these points .", "@xcite showed that the lomb - scargle periodograms for the both time series ( see @xcite , figures 7 a - c ) have a peak at 158.8 days which stands over the fap level by a significant amount . using the de method the above discrepancies are obvious . to establish the @xmath79 value the periodograms computed by the fft and", "the bt methods are shown in figure [ f6 ] ( the upper and the lower curve respectively ) .", "for @xmath46 and for periods less than 166 days there is a good comformity of the both periodograms ( but for periods greater than 166 days the points of the bt periodogram are not linked because the bt periodogram has much worse resolution than the fft periodogram ( no one know how to do it ) ) . for @xmath46 and @xmath113", "the value of @xmath21 is 13 ( @xmath71=153 $ ] ) .", "the inequality ( 7 ) is satisfied because @xmath114 .", "this means that the value of @xmath115 is mainly created by positive values of the autocorrelation function .", "the implication ( 8) needs an evaluation of the greatest value of the index @xmath65 where @xmath70 , but the solar data contain the most prominent period for @xmath116 days because of the solar rotation .", "thus , although @xmath117 for each @xmath118 , all sets @xmath41 ( see ( 5 ) and ( 6 ) ) without the set @xmath119 ( see ( 4 ) ) , which contains @xmath120 $ ] , are considered . this situation is presented in figure [ f7 ] . in this figure", "two curves @xmath121 and @xmath122 are plotted .", "the vertical dotted lines evaluate the intervals where the sets @xmath107 ( for @xmath123 ) are searched . for such @xmath41 two numbers", "are written : in parentheses the value of @xmath65 for the time series @xmath101 and above it the value of @xmath65 for the time series @xmath102 . to make this figure clear the curves are plotted for the set @xmath124 only .", "( in the right bottom corner information about the values of @xmath65 for the time series @xmath102 , for @xmath109 are written . )", "the implication ( 8) is not true , because @xmath125 for @xmath126 .", "therefore , @xmath43=153\\notin c_6=[423,500]$ ] .", "moreover , the autocorrelation function for @xmath127 $ ] is negative and the set @xmath128 is empty .", "thus , @xmath129 . on the basis of these information one can state , that the periodogram peak at @xmath130 days of the time series @xmath102 exists because of positive @xmath25 , but for @xmath23 from the intervals which do not contain this period .", "looking at the values of @xmath65 of the time series @xmath101 , one can notice that they decrease when @xmath23 increases until @xmath131 .", "this indicates , that when @xmath23 increases , the contribution of the 11-year cycle to the peaks of the periodogram decreases .", "an increase of the value of @xmath65 is for @xmath132 for the both time series , although the contribution of the 11-year cycle for the time series @xmath101 is insignificant .", "thus , this part of the autocorrelation function ( @xmath133 for the time series @xmath102 ) influences the @xmath21-th peak of the periodogram .", "this suggests that the periodicity at about 155 days is a harmonic of the periodicity from the interval of @xmath1 $ ] days .", "( solid line ) and consecutively smoothed sunspot areas of the one rotation time interval @xmath134 ( dotted line ) .", "both indexes are presented on the left axis .", "the lower curve illustrates fluctuations of the sunspot areas @xmath135 .", "the dotted and dashed horizontal lines represent levels zero and @xmath136 respectively .", "the fluctuations are shown on the right axis . ]", "the described reasoning can be carried out for other values of the periodogram .", "for example , the condition ( 8) is not satisfied for @xmath137 ( 250 , 222 , 200 days ) .", "moreover , the autocorrelation function at these points is negative .", "these suggest that there are not a true periodicity in the interval of [ 200 , 250 ] days .", "it is difficult to decide about the existence of the periodicities for @xmath138 ( 333 days ) and @xmath139 ( 286 days ) on the basis of above analysis . the implication ( 8) is not satisfied for @xmath139 and the condition ( 7 ) is not satisfied for @xmath138 , although the function @xmath25 of the time series @xmath102 is significantly positive for @xmath140 .", "the conditions ( 7 ) and ( 8) are satisfied for @xmath141 ( figure [ f8 ] ) and @xmath142 . therefore , it is possible to exist the periodicity from the interval of @xmath1 $ ] days .", "similar results were also obtained by @xcite for daily sunspot numbers and daily sunspot areas .", "she considered the means of three periodograms of these indexes for data from @xmath143 years and found statistically significant peaks from the interval of @xmath1 $ ] ( see @xcite , figure 2 ) .", "@xcite studied sunspot areas from 1876 - 1999 and sunspot numbers from 1749 - 2001 with the help of the wavelet transform .", "they pointed out that the 154 - 158-day period could be the third harmonic of the 1.3-year ( 475-day ) period .", "moreover , the both periods fluctuate considerably with time , being stronger during stronger sunspot cycles .", "therefore , the wavelet analysis suggests a common origin of the both periodicities . this conclusion confirms the de method result which indicates that the periodogram peak at @xmath144 days is an alias of the periodicity from the interval of @xmath1 $ ]"], ["in order to verify the existence of the periodicity at about 155 days i consider the following time series : + @xmath145     + @xmath146", "   + @xmath147  ", "+ the value @xmath134 is calculated analogously to @xmath83 ( see sect .", "the values @xmath148 and @xmath149 are evaluated from the formula ( 9 ) . in the upper part of figure [ f9 ] the time series of sunspot areas @xmath150 of the one rotation time interval from the whole solar disk and the time series of consecutively smoothed sunspot areas @xmath151", "are showed . in the lower part of figure [ f9 ]", "the time series of sunspot area fluctuations @xmath145 is presented .", "on the basis of these data the maximum activity period of cycle 16 is evaluated .", "it is an interval between two strongest fluctuations e.a .", "@xmath152 $ ] rotations .", "the length of the time interval @xmath153 is 54 rotations .", "if the about @xmath0-day ( 6 solar rotations ) periodicity existed in this time interval and it was characteristic for strong fluctuations from this time interval , 10 local maxima in the set of @xmath154 would be seen .", "then it should be necessary to find such a value of p for which @xmath155 for @xmath156 and the number of the local maxima of these values is 10 .", "as it can be seen in the lower part of figure [ f9 ] this is for the case of @xmath157 ( in this figure the dashed horizontal line is the level of @xmath158 ) .", "figure [ f10 ] presents nine time distances among the successive fluctuation local maxima and the horizontal line represents the 6-rotation periodicity .", "it is immediately apparent that the dispersion of these points is 10 and it is difficult to find even few points which oscillate around the value of 6 .", "such an analysis was carried out for smaller and larger @xmath136 and the results were similar .", "therefore , the fact , that the about @xmath0-day periodicity exists in the time series of sunspot area fluctuations during the maximum activity period is questionable .    .", "the horizontal line represents the 6-rotation ( 162-day ) period . ]    ]    ]    to verify again the existence of the about @xmath0-day periodicity during the maximum activity period in each solar hemisphere separately , the time series @xmath88 and @xmath89 were also cut down to the maximum activity period ( january 1925december 1930 ) .", "the comparison of the autocorrelation functions of these time series with the appriopriate autocorrelation functions of the time series @xmath88 and @xmath89 , which are computed for the whole 11-year cycle ( the lower curves of figures [ f1 ] and [ f2 ] ) , indicates that there are not significant differences between them especially for @xmath23=5 and 6 rotations ( 135 and 162 days ) ) .", "this conclusion is confirmed by the analysis of the time series @xmath146 for the maximum activity period .", "the autocorrelation function ( the lower curve of figure [ f11 ] ) is negative for the interval of [ 57 , 173 ] days , but the resolution of the periodogram is too low to find the significant peak at @xmath159 days .", "the autocorrelation function gives the same result as for daily sunspot area fluctuations from the whole solar disk ( @xmath160 ) ( see also the lower curve of figures [ f5 ] ) . in the case of", "the time series @xmath89 @xmath161 is zero for the fluctuations from the whole solar cycle and it is almost zero ( @xmath162 ) for the fluctuations from the maximum activity period .", "the value @xmath163 is negative .", "similarly to the case of the northern hemisphere the autocorrelation function and the periodogram of southern hemisphere daily sunspot area fluctuations from the maximum activity period @xmath147 are computed ( see figure [ f12 ] ) .", "the autocorrelation function has the statistically significant positive peak in the interval of [ 155 , 165 ] days , but the periodogram has too low resolution to decide about the possible periodicities .", "the correlative analysis indicates that there are positive fluctuations with time distances about @xmath0 days in the maximum activity period .", "the results of the analyses of the time series of sunspot area fluctuations from the maximum activity period are contradict with the conclusions of @xcite .", "she uses the power spectrum analysis only .", "the periodogram of daily sunspot fluctuations contains peaks , which could be harmonics or subharmonics of the true periodicities .", "they could be treated as real periodicities .", "this effect is not visible for sunspot data of the one rotation time interval , but averaging could lose true periodicities .", "this is observed for data from the southern hemisphere .", "there is the about @xmath0-day peak in the autocorrelation function of daily fluctuations , but the correlation for data of the one rotation interval is almost zero or negative at the points @xmath164 and 6 rotations .", "thus , it is reasonable to research both time series together using the correlative and the power spectrum analyses ."], ["the following results are obtained :    1 .", "a new method of the detection of statistically significant peaks of the periodograms enables one to identify aliases in the periodogram .", "2 .   two effects cause the existence of the peak of the periodogram of the time series of sunspot area fluctuations at about @xmath0 days : the first is caused by the 27-day periodicity , which probably creates the 162-day periodicity ( it is a subharmonic frequency of the 27-day periodicity ) and the second is caused by statistically significant positive values of the autocorrelation function from the intervals of @xmath165 $ ] and @xmath166 $ ] days .", "the existence of the periodicity of about @xmath0 days of the time series of sunspot area fluctuations and sunspot area fluctuations from the northern hemisphere during the maximum activity period is questionable .", "the autocorrelation analysis of the time series of sunspot area fluctuations from the southern hemisphere indicates that the periodicity of about 155 days exists during the maximum activity period ."], ["i appreciate valuable comments from professor j. jakimiec ."]]}
{"article_id": "1512.09139", "article_text": ["it is believed that the direct detection of gravitational waves ( gws ) will bring the era of gravitational wave astronomy .", "the interferometer detectors are now under operation and awaiting the first signal of gws  @xcite .", "it is also known that pulsar timing arrays ( ptas ) can be used as a detector for gws @xcite .", "these detectors are used to search for very low frequency ( @xmath0 ) gravitational waves , where the lower limit of the observable frequencies is determined by the inverse of total observation time @xmath1 .", "indeed , the total observation time has a crucial role in ptas , because ptas are most sensitive near the lower edge of observable frequencies @xcite . taking into account its sensitivity ,", "the first direct detection of the gravitational waves might be achieved by ptas .", "the main target of ptas is the stochastic gravitational wave background ( sgwb ) generated by a large number of unresolved sources with the astrophysical origin or the cosmological origin in the early universe .", "the promising sources are super massive black hole binaries  @xcite , cosmic ( super)string  @xcite , and inflation  @xcite .", "previous studies have assumed that the sgwb is isotropic and unpolarized  @xcite .", "these assumptions are reasonable for the primary detection of the sgwb , but the deviation from the isotropy and the polarizations should have rich information of sources of gravitational waves .", "recently , the cross - correlation formalism has been generalized to deal with anisotropy in the sgwb @xcite .", "result of this work enables us to consider arbitrary levels of anisotropy , and a bayesian approach was performed by using this formalism @xcite . on the other hand , for the anisotropy of the sgwb , the cross - correlation formalism has been also developed in the case of interferometer detectors  @xcite .", "as to the polarization , there are works including the ones motivated by the modified gravity  @xcite", ". we can envisage supermassive black hole binaries emit circularly polarized sgwb due to the chern - simons term  @xcite .", "there may also exist cosmological sgwb with circular polarization in the presence of parity violating term in gravity sector  @xcite .    in this paper", ", we investigate the detectability of circular polarization in the sgwb by ptas .", "we characterize sgwb by the so called stokes @xmath2 parameter  @xcite and calculate generalized overlap reduction functions ( orfs ) so that we can probe the circular polarization of the sgwb .", "we also discuss a method to separate the intensity ( @xmath3 mode ) and circular polarization ( @xmath2 mode ) of the sgwb .", "the paper is organized as follows . in section [ sec :", "stokes parameters for a plane gravitational wave ] , we introduce the stokes parameters for monochromatic plane gravitational waves , and clarify the physical meaning of the stokes parameters @xmath3 and @xmath2 . in section [ sec : formulation ] , we formulate the cross - correlation formalism for anisotropic circularly polarized sgwb with ptas .", "the basic framework is essentially a combination of the formalism of @xcite , and the polarization decomposition formula of the sgwb derived in @xcite . in section [ sec : the generalized overlap reduction function for circular polarization ] , we calculate the generalized orfs for the @xmath2 mode .", "the results for @xmath3 mode are consistent with the previous work  @xcite . in section [ sec : separation method ] , we give a method for separation between the @xmath3 mode and @xmath2 mode of the sgwb .", "the final section is devoted to the conclusion . in appendixes , we present analytic results for the generalized overlap reduction functions . in this paper", ", we will use the gravitational units @xmath4 .", "let us consider the stokes parameters for plane waves traveling in the direction @xmath5 , which can be described by @xmath6 \\", ", \\\\ & & h_{xy}(t , z)=h_{yx}(t , z)={\\rm re}[b_{\\times}\\mathrm{e}^{-iw(t - z ) } ] \\ .\\end{aligned}\\ ] ] for an idealized monochromatic plane wave , complex amplitudes @xmath7 and @xmath8 are constants .", "polarization of the plane gws is characterized by the tensor , ( see @xcite and also electromagnetic case @xcite ) @xmath9 where @xmath10 take @xmath11 .", "any @xmath12 hermitian matrix can be expanded by the pauli and the unit matrices with real coefficients .", "hence , the @xmath13 hermitian matrix @xmath14 can be written as @xmath15 where @xmath16 by analogy with electromagnetic cases , @xmath17 and @xmath2 are called stokes parameters . comparing with , we can read off the stokes parameters as @xmath18= b_{+}^{\\ast}b_{\\times}+ b_{\\times}^{\\ast}b_{+},\\\\ v&=&-2{\\rm i m } [ b_{+}^{\\ast}b_{\\times}]=i ( b_{+}^{\\ast}b_{\\times}- b_{\\times}^{\\ast}b_{+}).\\label{stv}\\end{aligned}\\ ] ] apparently , the real parameter @xmath3 is the intensity of gws . in order to reveal the physical meaning of the real parameter @xmath2 , we define the circular polarization bases @xcite @xmath19 from the relation @xmath20 we see @xmath21", "thus , we can rewrite the stokes parameters - as @xmath22 from the above expression , we see that the real parameter @xmath2 characterizes the asymmetry of circular polarization amplitudes .", "the other parameters @xmath23 and @xmath24 have additional information about linear polarizations by analogy with the electromagnetic cases .", "alternatively , we can also define the tensor @xmath25 in circular polarization bases @xmath26 where @xmath27 .", "note that the stokes parameters satisfy a relation @xmath28    next , we consider the transformation of the stokes parameters under rotations around the @xmath5 axis . the rotation around the @xmath5 axis is given by @xmath29 where @xmath30 is the angle of the rotation .", "the gws traveling in the direction @xmath5 @xmath31 transform as @xmath32 where we took the transverse traceless gauge @xmath33 after a short calculation , we obtain @xmath34 using and , the four stokes parameters ( [ sti])-([stv ] ) transform as @xmath35 as you can see , the parameters @xmath23 and @xmath24 depend on the rotation angle @xmath30 .", "this reflects the fact that @xmath23 and @xmath24 parameters characterize linear polarizations .", "note that this transformation is similar to the transformation of electromagnetic case except for the angle @xmath36 and can be rewritten as @xmath37", "in this section , we study anisotropic distribution of sgwb and focus on the detectability of circular polarizations with pulsar timing arrays .", "we combine the analysis of @xcite and that of @xcite . in sec.[subsec : the spectral ] , we derive the power spectral density for anisotropic circularly polarized sgwb @xmath38 .", "then we also derive the dimensionless density parameter @xmath39 which is expressed by the frequency spectrum of intensity @xmath40  @xcite . in sec.[subsec : the signal ] , we extend the generalized orfs to cases with circular polarizations characterized by the parameter @xmath2 . for simplicity ,", "we consider specific anisotropic patterns with @xmath41 expressed by the spherical harmonics @xmath42 .      in the transverse traceless gauge , metric perturbations @xmath43 with a given propagation direction @xmath44", "can be expanded as @xcite @xmath45 where the fourier amplitude satisfies @xmath46 as a consequence of the reality of @xmath43 , @xmath47 , @xmath48 is the frequency of the gws , @xmath49 are spatial indices , @xmath50 label polarizations .", "note that the fourier amplitude @xmath51 satisfies the relation @xmath52 where @xmath53 was defined by .", "the polarized tensors @xmath54 are defined by @xmath55 where @xmath56 and @xmath57 are unit orthogonal vectors perpendicular to @xmath58 .", "the polarization tensors satisfy @xmath59 with polar coordinates , the direction @xmath44 can be represented by @xmath60 and the polarization basis vectors read @xmath61    we assume the fourier amplitudes @xmath62 are random variables , which is stationary and gaussian .", "however , they are not isotropic and unpolarized .", "the ensemble average of fourier amplitudes can be written as @xcite @xmath63 where @xmath64 here , the bracket @xmath65 represents an ensemble average , and @xmath66 is the dirac delta function on the two - sphere .", "the gw power spectral density @xmath38 is a hermitian matrix , and satisfies @xmath67 because of the relation @xmath46 .", "therefore , we have the relations @xmath68 note that the stokes parameters are not exactly the same as the expression of , but they have the relation and characterize the same polarization .", "we further assume that the sgwbs satisfy @xmath69 we also assume the directional dependence of the sgwb is frequency independent @xcite .", "this implies the gw power spectral density is factorized into two parts , one of which depends on the direction while the other depends on the frequency .", "because of the transformations - , the parameters @xmath3 and @xmath2 have spin 0 and the parameters @xmath70 have spin @xmath71  @xcite . to analyze the sgwb on the sky , it is convenient to expand the stokes parameters by spherical harmonics @xmath72", ". however , since @xmath70 parameters have spin @xmath71 , they have to be expanded by the spin - weighted harmonics @xmath73 @xcite .", "thus , we obtain @xmath74 in this paper , we study specific anisotropic patterns with @xmath41 for simplicity .", "therefore , we can neglect @xmath23 and @xmath24 from now on .", "thus , the gw power spectral density becomes @xmath75 where @xmath76 so , we focus on the parameters @xmath3 and @xmath2 . in what follows , we will use the following shorthand notation @xmath77    next , we consider the dimensionless density parameter  @xcite @xmath78 where @xmath79 is the critical density , @xmath80 is the present value of the hubble parameter , @xmath81 is the energy density of gravitational waves , and @xmath82 is the energy density in the frequency range @xmath48 to @xmath83 .", "the bracket @xmath65 represents the ensemble average .", "however , actually , we take a spatial average over the wave lengths @xmath84 of gws or a temporal average over the periods @xmath85 of gws . here", ", we assumed the ergodicity , namely , the ensemble average can be replaced by the temporal average .", "using , , , as well as @xmath46 and @xmath86 , we get @xmath87 then we define @xmath88 hence , the dimensionless quantity @xmath39 in is given by @xmath89 where the spherical harmonics are orthogonal and normalized as @xmath90 using @xmath91 , we obtain @xmath92 without loss of generality , we normalize the monopole moment as @xmath93 so , becomes @xmath94      the time of arrival of radio pulses from the pulsar is affected by gws .", "consider a pulsar with frequency @xmath95 located in the direction @xmath96 . to detect the sgwb ,", "let us consider the redshift of the pulse from a pulsar @xcite @xmath97 where @xmath98 is a frequency detected at the earth and @xmath96 is the direction to the pulsar .", "the unit vector @xmath44 represents the direction of propagation of gravitational plane waves .", "we also defined the difference between the metric perturbations at the pulsar @xmath99 and at the earth @xmath100 as @xmath101 the gravitational plane waves at each point is defined as @xmath102 for the sgwb , the redshift have to be integrated over the direction of propagation of the gravitational waves @xmath44 : @xmath103 we choose a coordinate system @xmath104 and assume that the amplitudes of the metric perturbation at the pulsar and the earth are the same .", "then becomes @xmath105 and therefore , reads @xmath106 where we have defined the pattern functions for pulsars @xmath107 note that our convention for the fourier transformation is @xmath108 therefore , the fourier transformation of can be written as @xmath109    in the actual signals from a pulsar , there exist noises .", "hence , we need to use the correlation analysis .", "we consider the signals from two pulsars @xmath110 where @xmath111 labels the pulsar . here", ", @xmath112 denotes the signal from the pulsar and @xmath113 denotes the noise intrinsic to the measurement .", "we assume the noises are stationary , gaussian and are not correlated between the two pulsars .", "to correlate the signals of two measurements , we define @xmath114 where @xmath1 is the total observation time and @xmath115 is a real filter function which should be optimal to maximize signal - to - noise ratio . in the case of interferometer", ", the optimal filter function falls to zero for large @xmath116 compered to the travel time of the light between the detecters .", "since the signals of two detectors are expected to correlate due to the same effect of the gravitational waves , the optimal filter function should behave this way .", "then , typically one of the detectors is very close to the other compared to the total observation time @xmath1 .", "therefore , the total observation time @xmath1 can be extended to @xmath117 @xcite .", "in contrast , in the case of pta , it is invalid that @xmath1 is very large compered to the travel time of the light between the pulsars .", "nevertheless , we can assume that one of the two @xmath1 can be expanded to @xmath117 , because in situations @xmath118 and @xmath119 it is known that we can ignore the effect of the distance @xmath120 of pulsars .", "in this case , it is clear that any locations of the pulsars are optimal and optimal filter function should behave like as the interferometer case @xcite .    using these assumptions @xmath118 and @xmath119 , we can rewrite as @xmath121 where @xmath122 note that @xmath123 satisfies @xmath124 , because @xmath125 is real .", "moreover , to deal with the unphysical region @xmath126 we require @xmath127 .", "thus , @xmath123 becomes real .", "taking the ensemble average , using @xmath128 , @xmath118 , and assuming the noises in the two measurements are not correlated , we get @xmath129\\ , \\label{s2}\\end{aligned}\\ ] ] where we have defined @xmath130 the functions @xmath131 and @xmath132 are called the generalized orfs , which describe the angular sensitivity of the pulsars for the sgwb . note that , as we already mentioned , we consider the cases of @xmath41 for simplicity", ". then we have assumed @xmath118 and @xmath119 , this assumption implies that approximately becomes @xmath133 due to the rapid oscillation of the phase factor .", "therefore , the distance @xmath120 of the pulsars does not appear in the generalized orfs , and hence the generalized orfs do not depend on the frequency .    as you can see from ,", "the correlation of the two measurements involve both the total intensity and the circular polarization .", "however , the degeneracy can be disentangled by using separation method , which will be discussed in the section [ sec : separation method ] .", "in this section , we consider the generalized orfs for circular polarizations : @xmath134 where we defined @xmath135 in the above , we have used and the fact that the generalized orfs do not depend on frequency . for computation of the generalized orfs for circular polarizations ,", "it is convenient to use the computational frame @xcite defined by @xmath136 where @xmath137 is the angular separation between the two pulsars . using - , , and", ", one can easily show that @xmath138 we therefore get @xmath139 the explicit form of the spherical harmonics reads @xmath140 where @xmath141 is the normalization factor .", "the associated legendre functions are given by @xmath142 and @xmath143 with the legendre functions @xmath144\\ .\\label{pl}\\end{aligned}\\ ] ] using the spherical harmonics , becomes @xmath145 where we have used the fact that the function of @xmath146 is odd parity in the case of @xmath147 and is even parity in the case of @xmath148 .", "note that the generalized orfs for circular polarizations are real functions . in the case of @xmath149 and/or @xmath150 ,", "the integrand in vanishes .", "therefore , we can not detect circular polarizations for these cases .", "this fact for @xmath151 implies that we do not need to consider auto - correlation for a single pulsar .", "this is the reason why we neglected auto - correlation term in .    integrating ,", "we get the following form for @xmath152 : @xmath153 for @xmath154 , we have obtained @xmath155 \\ , \\\\ \\gamma^{v}_{1 - 1}&=&\\gamma^{v}_{11 } \\ , \\end{aligned}\\ ] ] recall that @xmath156 .", "the derivation of this formula for @xmath154 can be found in appendix [ sec : angular integral of the generalized overlap reduction function for dipole circular polarization ] .    for @xmath157 , we derived the following : @xmath158\\ , \\\\ \\gamma^{v}_{2 - 1}&=&\\gamma^{v}_{21}\\ , \\\\", "\\gamma^{v}_{22}&=&-\\frac{\\sqrt{30\\pi}}{6}(1-\\cos\\xi)\\left[2-\\cos\\xi+6\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{v}_{2 - 2}&=&-\\gamma^{v}_{22}\\ , \\end{aligned}\\ ] ] for @xmath159 , the results are @xmath160\\ , \\\\ \\gamma^{v}_{3 - 1}&=&\\gamma^{v}_{31}\\ , \\\\ \\gamma^{v}_{32}&=&\\frac{\\sqrt{210\\pi}}{24}(1-\\cos\\xi)\\left[8 - 5\\cos\\xi-\\cos^2\\xi+24\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{v}_{3 - 2}&=&-\\gamma^{v}_{3 - 2}\\ , \\\\ \\gamma^{v}_{33}&=&-\\frac{\\sqrt{35\\pi}}{16}\\sin\\xi\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\left[11 - 6\\cos\\xi-\\cos^2\\xi+32\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{v}_{3 - 3}&=&\\gamma^{v}_{33}\\ .\\end{aligned}\\ ] ] in fig .", "[ gv ] , we plotted these generalized orfs as a function of the angular separation between the two pulsars @xmath137 .", "it is apparent that considering the @xmath2 mode does not make sense when we only consider the isotropic ( @xmath152 ) orf . on the other hand ,", "when we consider anisotropic ( @xmath161 ) orfs , it is worth taking into account polarizations .", "the polarizations of the sgwb would give us rich information both of super massive black hole binaries and of inflation in the early universe .", "as a function of the angular separation between the two pulsars @xmath137 . in fig .", "[ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( a ) @xmath152     as a function of the angular separation between the two pulsars @xmath137 . in fig . [ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( b ) @xmath154     as a function of the angular separation between the two pulsars @xmath137 . in fig . [ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( c ) @xmath157     as a function of the angular separation between the two pulsars @xmath137 . in fig .", "[ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( d ) @xmath159    using the same procedure described in the above to derive the generalized orfs for circular polarizations , we can also derive the generalized orfs for the intensity @xmath166 where @xmath167 the angular integral in this case was performed in @xcite .", "the results are summarized in appendix [ sec : the generalized overlap reduction function for intensity ] .", "in this section , we separate the @xmath3 mode and @xmath2 mode of the sgwb with correlation analysis @xcite . to this aim ,", "we use four pulsars ( actually we need at least three pulsars ) , and define correlations of @xmath168 @xmath169 where @xmath170 label the pulsars . comparing with , we obtain @xmath171 \\ ,", "\\label{1c12}\\\\ & & c_{34}(f)=\\sum_{lm}^{l=3}\\left[c_{lm}^{i}i(f)\\gamma_{lm,34}^{i}+c_{lm}^{v}v(f)\\gamma_{lm,34}^{v}\\right ] \\ .\\label{1c34}\\end{aligned}\\ ] ] if the @xmath3 mode and @xmath2 mode of the sgwb are dominated by a certain @xmath172 and @xmath173 , and become @xmath174 \\ , \\label{2c12 } \\\\ & & c_{34}(f)=\\left[c _ { l m}^{i}i(f)\\gamma _ { l m,34}^{i}+c _ { l ' m'}^{v}v(f)\\gamma _ { l ' m',34}^{v}\\right ] \\ .\\label{2c34}\\end{aligned}\\ ] ] to separate the intensity and the circular polarization , we take the following linear combinations @xmath175 where we defined coefficients @xmath176 as you can see , @xmath177 contains only @xmath40 , and @xmath178 contains only @xmath179 .", "for the signal @xmath180 , the formulas corresponding to and are given by @xmath181 \\ , \\label{sp}\\end{aligned}\\ ] ] where @xmath182 denotes @xmath3 and @xmath2 .", "we assume @xmath183 and that the noise in the four pulsars are not correlated .", "we also assume that the ensemble average of fourier amplitudes of the noises @xmath184 is of the form @xmath185 where @xmath186 is the noise power spectral density .", "the reality of @xmath187 gives rise to @xmath188 and therefore we obtain @xmath189 . without loss of generality", ", we can assume @xmath190 then we obtain corresponding noises @xmath191 : @xmath192\\ , \\label{np}\\end{aligned}\\ ] ] where @xmath193^{1/2 } \\label{sn12 } \\ , \\quad s_{n,34}(f ) \\equiv [ s_{n,3}(f)s_{n,4}(f)]^{1/2 } \\label{sn34 } \\ .\\end{aligned}\\ ] ] using the inner product @xmath194 \\ , \\end{aligned}\\ ] ] we can rewrite , as @xmath195 therefore , the optimal filter function can be chosen as @xmath196 using , we get optimal signal - to - noise ratio @xmath197^{1/2}\\ .\\label{snr}\\end{aligned}\\ ] ] plugging , , and into , we obtain @xmath198^{1/2}\\ , \\\\", "{ \\rm snr}_{v}&=&\\left[t\\int_{-\\infty}^{\\infty}df\\,\\,\\frac{\\left(c^{v}_{{l}'{m}'}\\right)^{2}v^{2}(f)\\left(\\gamma_{{l}'{m}',34}^{v}\\gamma^{i}_{{l}{m},12}-\\gamma_{{l}'{m}',12}^{v}\\gamma^{i}_{{l}{m},34}\\right)^2}{\\left(\\gamma^{i}_{{l}{m},12}\\right)^2s^{2}_{n,34}(f)+\\left(\\gamma^{i}_{{l}{m},34}\\right)^2s^{2}_{n,12}(f)}\\right]^{1/2}\\ .\\end{aligned}\\ ] ] if we assume all of the noise power spectral densities are the same , becomes @xmath199 thus , the compiled orfs can be defined as @xmath200^{1/2}}\\ , \\\\ \\gamma_{12:34}^{v}&\\equiv&\\frac{\\gamma_{{l}'{m}',34}^{v}\\gamma^{i}_{{l}{m},12}-\\gamma_{{l}'{m}',12}^{v}\\gamma^{i}_{{l}{m},34}}{\\left[\\left(\\gamma^{i}_{{l}{m},12}\\right)^2+\\left(\\gamma^{i}_{{l}{m},34}\\right)^2\\right]^{1/2}}\\ .\\end{aligned}\\ ] ] this compiled orfs @xmath201 and @xmath202 describe the angular sensitivity of the four pulsars for the pure @xmath3 and @xmath2 mode of the sgwb , respectively .", "note that , to do this separation , we must know a priori the coefficients @xmath203 and @xmath204 .", "if we do not assume , the generalized orfs depend on the frequency . in this case", ", it seems difficult to calculate these coefficients .", "we next consider the case that @xmath3 mode and/or @xmath2 mode dominant in two or more @xmath205 . in this case , if we have a priori knowledge of the values of @xmath206 in each of @xmath205 for coefficients", "@xmath203 and @xmath204 , we can separate @xmath3 mode and @xmath2 mode . for example , assume that @xmath3 mode is dominated by @xmath207 , while @xmath2 mode is dominated by @xmath208 , then and become @xmath209\\ , \\label{3c12}\\\\ & & c_{34}(f)=\\left[c^{i}_{00}i(f)\\left(\\gamma_{00,34}^{i}+\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)+c_{11}^{v}v(f)\\gamma_{11,34}^{v}\\right]\\ .\\label{3c34}\\end{aligned}\\ ] ] thus , we can separate @xmath3 mode and @xmath2 mode by using linear combinations @xmath210\\ , \\\\", "d_{v}&\\equiv&a_{v}c_{34}(f)+b_{v}c_{12}(f ) \\nonumber\\\\ & = & c_{11}^{v}v(f)\\left[\\gamma_{11,34}^{v}\\left(\\gamma_{00,12}^{i}+\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,12}^{i}\\right)-\\gamma_{11,12}^{v}\\left(\\gamma_{00,34}^{i}+\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)\\right]\\ , \\end{aligned}\\ ] ] where @xmath211 as in the previous calculations , we can get the compiled orfs @xmath212^{1/2}}\\ , \\label{gi1234}\\\\ \\gamma_{12:34}^{v}&\\equiv&\\frac{\\gamma_{11,34}^{v}\\left(\\gamma_{00,12}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,12}^{i}\\right)-\\gamma_{11,12}^{v}\\left(\\gamma_{00,34}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)}{\\left[\\left(\\gamma_{00,12}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,12}^{i}\\right)^2+\\left(\\gamma_{00,34}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)^2\\right]^{1/2}}\\ .\\label{gv1234}\\end{aligned}\\ ] ]    [ cols=\"^,^ \" , ]     in fig .", "[ cg ] we show some compiled orfs @xmath213 ( left panels ) and @xmath214 ( right panels ) as a function of the two angular separations @xmath137 and @xmath215 for two pulsar pairs , respectively .", "we used the expressions of @xmath2 mode and @xmath3 mode ( see appendix [ sec : the generalized overlap reduction function for intensity ] ) , and we assumed @xmath216 for simplicity . in fig .", "[ cg](a ) and [ cg](b ) , the @xmath3 mode is dominated by @xmath217 and @xmath2 mode is dominated by @xmath218 . in fig .", "[ cg](c ) and [ cg](d ) , the @xmath3 mode is dominated by @xmath219 and @xmath2 mode is dominated by @xmath218 . in fig .", "[ cg](e ) and [ cg](f ) , the @xmath3 mode is dominated by @xmath207 and @xmath2 mode is dominated by @xmath218 . in fig .", "[ cg](e ) and [ cg](f ) , the @xmath3 mode is dominated by @xmath220 and @xmath2 mode is dominated by @xmath218 . by definition , in the case of @xmath221 ,", "the compiled orfs are zero .", "we have studied the detectability of the stochastic gravitational waves with ptas . in most of the previous works ,", "the isotropy of sgwb has been assumed for the analysis .", "recently , however , a stochastic gravitational wave background with anisotropy have been considered .", "the information of the anisotropic pattern of the distribution should contain important information of the sources such as supermassive black hole binaries and the sources in the early universe .", "it is also intriguing to take into account the polarization of sgwb in the pta analysis .", "therefore , we extended the correlation analysis to circularly polarized sgwb and calculated generalized overlap reduction functions for them .", "it turned out that the circular polarization can not be detected for an isotropic background .", "however , when the distribution has anisotropy , we have shown that there is a chance to observe circular polarizations in the sgwb .", "we also discussed how to separate polarized modes from unpolarized modes of gravitational waves .", "if we have a priori knowledge of the abundance ratio for each mode in each of @xmath205 , we can separate @xmath3 mode and @xmath2 mode in general .", "this would be possible if we start from fundamental theory and calculate the spectrum of sgwb .", "in particular , in the case that the signal of lowest @xmath222 is dominant , we performed the separation of @xmath3 mode and @xmath2 mode explicitly .", "this work was supported by grants - in - aid for scientific research ( c ) no.25400251 and \" mext grant - in - aid for scientific research on innovative areas no.26104708 and `` cosmic acceleration''(no.15h05895 ) .", "in this appendix , we perform angular integration of the generalized orf for dipole ( @xmath154 ) circular polarization ( see @xcite ) : @xmath223 where we have defined @xmath224 .", "it is obvious that in the case of @xmath225 , integrand of the generalized orf is zero , because of @xmath226 , then we obtain @xmath227 then , using - , we calculate @xmath228 and we find @xmath229    therefore we only have to consider the dipole generalized orf in the case of @xmath154 , @xmath230 : @xmath231 where @xmath232 first , to calculate @xmath233 , we use contour integral in the complex plane . defining @xmath234 and substituting @xmath235 into , we can rewrite @xmath233 as @xmath236 } \\ , \\end{aligned}\\ ] ] where @xmath237 denotes a unit circle .", "we can factorize the denominator of the integrand and get @xmath238 where @xmath239 hereafter , the upper sign applies when @xmath240 and the lower one applies when @xmath241 .", "note that we only consider the region @xmath242 , so we have used the relation @xmath243 in above expression . in the region", "@xmath244 , @xmath245 is inside the unit circle @xmath237 except for @xmath246 and @xmath247 is outside the unit circle @xmath237 .", "now , we can perform the integral using the residue theorem @xmath248 where @xmath249 the residues inside the unit circle @xmath237 can be evaluated as @xmath250\\right\\ }             = \\frac{i(z_{+}+z_{-})}{2\\sqrt{1-x^2}\\sin\\xi } \\ , \\end{aligned}\\ ] ] @xmath251 thus , we obtain @xmath252 next , we consider @xmath253 defined in", ". using , we can calculate @xmath253 as @xmath254    similarly , we can evaluate @xmath255 given in . to calculate @xmath255 in the complex plane , we again substitute into and obtain @xmath256 we use the residue theorem @xmath257 where @xmath258 the residues inside the unit circle @xmath237 can be calculated as @xmath259\\right\\ }        = \\frac{i(z_{+}^2+z_{-}^2)}{4\\sqrt{1-x^2}\\sin\\xi } \\ , \\end{aligned}\\ ] ] @xmath260 therefore , @xmath255 becomes @xmath261 substituting to , we can calculate @xmath262 : @xmath263    finally , substituting and into , we get the generalized orf for @xmath264 @xmath265\\ .\\end{aligned}\\ ] ]     as a function of the angular separation between the two pulsars @xmath137 .", "[ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) .", "the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( a ) @xmath152     as a function of the angular separation between the two pulsars @xmath137 .", "[ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) .", "the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( b ) @xmath154     as a function of the angular separation between the two pulsars @xmath137 . fig . [ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) .", "the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( c ) @xmath157     as a function of the angular separation between the two pulsars @xmath137 .", "[ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) . the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( d ) @xmath159", "in this appendix , we show orfs for the intensity @xcite .", "the following form for @xmath152 was derived in @xcite , and our expressions are identical to their expressions : @xmath271\\ , \\end{aligned}\\ ] ] for , @xmath154 , we calculated as @xmath272\\ , \\\\ \\gamma^{i}_{11}&=&\\frac{\\sqrt{6\\pi}}{12}\\sin\\xi\\left[1 + 3(1-\\cos\\xi)\\left\\{1+\\frac{4}{1+\\cos\\xi}\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right\\}\\right]\\ , \\\\ \\gamma^{i}_{1 - 1}&=&-\\gamma^{i}_{11}\\ , \\end{aligned}\\ ] ] for @xmath157 , we obtain @xmath273\\ , \\\\ \\gamma^{i}_{21}&=&-\\frac{\\sqrt{30\\pi}}{60}\\sin\\xi\\left[21 - 15\\cos\\xi-5\\cos^2\\xi+60\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{2 - 1}&=&-\\gamma^{i}_{2 - 1}\\ , \\\\ \\gamma^{i}_{22}&=&\\frac{\\sqrt{30\\pi}}{24}(1-\\cos\\xi)\\left[9 - 4\\cos\\xi-\\cos^2\\xi+24\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{2 - 2}&=&\\gamma^{i}_{22}\\ ,", "\\end{aligned}\\ ] ] for @xmath159 , it is straightforward to reach the following @xmath274\\ , \\\\ \\gamma^{i}_{31}&=&\\frac{\\sqrt{21\\pi}}{48}\\sin\\xi(1-\\cos\\xi)\\left[34 + 15\\cos\\xi+5\\cos^2\\xi+\\frac{96}{1+\\cos\\xi}\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{3 - 1}&=&-\\gamma^{i}_{31}\\ , \\\\ \\gamma^{i}_{32}&=&-\\frac{\\sqrt{210\\pi}}{48}(1-\\cos\\xi)\\left[17 - 9\\cos\\xi-3\\cos^2\\xi-\\cos^3\\xi+48\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{3 - 2}&=&\\gamma^{i}_{32}\\ , \\\\ \\gamma^{i}_{33}&=&\\frac{\\sqrt{35\\pi}}{48}\\frac{(1-\\cos\\xi)^2}{\\sin\\xi}\\left[34 - 17\\cos\\xi-4\\cos^2\\xi-\\cos^3\\xi+96\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{3 - 3}&=&-\\gamma^{i}_{33}\\ .\\end{aligned}\\ ] ] these are plotted in fig .", "the generalized orfs of total intensity are different from that of circular polarization in that the value for @xmath149 is non - trivial .", "then the @xmath3 mode orfs for @xmath275 have value even in the case of @xmath151 .", "this implies that we can consider auto - correlation for a single pulsar .", "99 j.  aasi _ et al .", "_ [ ligo scientific collaboration ] , class .", "grav .   * 32 * , 074001 ( 2015 ) doi:10.1088/0264 - 9381/32/7/074001 [ arxiv:1411.4547 [ gr - qc ] ] .", "f.  acernese _ et al . _", "[ virgo collaboration ] , class .", "grav .   * 32 * , no . 2 , 024001 ( 2015 ) doi:10.1088/0264 - 9381/32/2/024001 [ arxiv:1408.3978 [ gr - qc ] ] . k.  somiya [ kagra collaboration ] , class .", "* 29 * , 124007 ( 2012 ) doi:10.1088/0264 - 9381/29/12/124007 [ arxiv:1111.7185 [ gr - qc ] ] .", "s.  l.  detweiler , astrophys .", "j.   * 234 * , 1100 ( 1979 ) .", "doi:10.1086/157593 romani , r. w. ( 1989 ) .", "timing a millisecond pulsar array . timing neutron stars , 113 - 117 .", "l.  lentati _ et al .", "_ , mon .  not .", "soc .   * 453 * , 2576 ( 2015 ) [ arxiv:1504.03692 [ astro-ph.co ] ] .", "z.  arzoumanian _ et al .", "_ [ nanograv collaboration ] , arxiv:1508.03024 [ astro-ph.ga ] . e.  s.  phinney , astro - ph/0108028 . a.  vilenkin , phys .", "rept .   * 121 * , 263 ( 1985 ) .", "s.  kuroyanagi , k.  miyamoto , t.  sekiguchi , k.  takahashi and j.  silk , phys .", "d * 87 * , no . 2 , 023522 ( 2013 ) [ phys .", "d * 87 * , no .", "6 , 069903 ( 2013 ) ] [ arxiv:1210.2829 [ astro-ph.co ] ] .", "m.  maggiore , gr - qc/0008027 .", "l.  p.  grishchuk , phys .", "* 48 * , 1235 ( 2005 ) doi:10.1070/pu2005v048n12abeh005795 [ gr - qc/0504018 ] .", "b.  allen and j.  d.  romano , phys .", "d * 59 * , 102001 ( 1999 ) [ gr - qc/9710117 ] .", "c.  m.  f.  mingarelli , t.  sidery , i.  mandel and a.  vecchio , phys .", "d * 88 * , no . 6 , 062005 ( 2013 ) [ arxiv:1306.5394 [ astro-ph.he ] ] .", "s.  r.  taylor and j.  r.  gair , phys .", "d * 88 * , 084001 ( 2013 ) [ arxiv:1306.5395 [ gr - qc ] ] .", "n.  seto and a.  taruya , phys .", "lett .   * 99 * , 121101 ( 2007 ) doi:10.1103/physrevlett.99.121101 [ arxiv:0707.0535 [ astro - ph ] ] . n.  seto and a.  taruya , phys .  rev .", "d * 77 * , 103001 ( 2008 ) [ arxiv:0801.4185 [ astro - ph ] ] . s.  j.  chamberlin and x.  siemens , phys .", "d * 85 * , 082001 ( 2012 ) doi:10.1103/physrevd.85.082001 [ arxiv:1111.5661 [ astro-ph.he ] ] . j.  r.  gair , j.  d.  romano and s.  r.  taylor , phys .  rev .", "d * 92 * , no .", "10 , 102003 ( 2015 ) doi:10.1103/physrevd.92.102003 [ arxiv:1506.08668 [ gr - qc ] ] .", "r.  jackiw and s.  y.  pi , phys .", "d * 68 * , 104012 ( 2003 ) doi:10.1103/physrevd.68.104012 [ gr - qc/0308071 ] . m.  satoh , s.  kanno and j.  soda , phys .", "d * 77 * , 023526 ( 2008 ) doi:10.1103/physrevd.77.023526 [ arxiv:0706.3585 [ astro - ph ] ] . c.  r.  contaldi , j.  magueijo and l.  smolin , phys .", "* 101 * ( 2008 ) 141101 doi:10.1103/physrevlett.101.141101 [ arxiv:0806.3082 [ astro - ph ] ] .", "t.  takahashi and j.  soda , phys .", "lett .   * 102 * , 231301 ( 2009 ) doi:10.1103/physrevlett.102.231301 [ arxiv:0904.0554 [ hep - th ] ] .", "j.  l.  cook and l.  sorbo , phys .", "d * 85 * , 023534 ( 2012 ) [ phys .", "d * 86 * , 069901 ( 2012 ) ] doi:10.1103/physrevd.86.069901 , 10.1103/physrevd.85.023534 [ arxiv:1109.0022 [ astro-ph.co ] ] .", "i.  obata , t.  miura and j.  soda , phys .", "d * 92 * , no .", "6 , 063516 ( 2015 ) doi:10.1103/physrevd.92.063516 [ arxiv:1412.7620 [ hep - ph ] ] .", "a. p. lightman , w. h. press , r. h. price , and s. a. teukolski , problem book in relativity and gravitation , 2nd ed .", "( princeton university press , 1979 ) .", "m. maggiore , gravitational waves , vol . 1 : theory and experiments ( oxford university press , 2008 )", ". g. b. rybicki , a. p. lightman , radiative processes in astrophysics ( wiley - interscience , 1979 ) .", "l. d. landau and e. m. lifshitz , the classical theory ( pergamon press , 1975 ) . c. misner , k. thorne and j. wheeler , gravitation , ( freeman 1973 ) .", "b.  allen and a.  c.  ottewill , phys .", "d * 56 * , 545 ( 1997 ) doi:10.1103/physrevd.56.545 [ gr - qc/9607068 ] .", "u.  seljak and m.  zaldarriaga , phys .", "lett .   * 78 * , 2054 ( 1997 ) doi:10.1103/physrevlett.78.2054 [ astro - ph/9609169 ] .", "j. n. goldberg et al .", ", journal of mathematical physics * 8 * , 2155 , 1967 . m.  anholm , s.  ballmer , j.  d.  e.  creighton , l.  r.  price and x.  siemens , phys .", "d * 79 * , 084030 ( 2009 ) doi:10.1103/physrevd.79.084030 [ arxiv:0809.0701 [ gr - qc ] ] . l.  g.  book and e.  e.  flanagan , phys .", "d * 83 * , 024024 ( 2011 ) doi:10.1103/physrevd.83.024024 [ arxiv:1009.4192 [ astro-ph.co ] ] . f.  a.  jenet and j.  d.  romano , am .", "j.  phys .", "* 83 * , 635 ( 2015 ) doi:10.1119/1.4916358 [ arxiv:1412.1142 [ gr - qc ] ] .", "r.  w.  hellings and g.  s.  downs , astrophys .", "j.   * 265 * , l39 ( 1983 ) . doi:10.1086/183954"], "abstract_text": ["<S> we study the detectability of circular polarization in a stochastic gravitational wave background from various sources such as supermassive black hole binaries , cosmic strings , and inflation in the early universe with pulsar timing arrays . </S>", "<S> we calculate generalized overlap reduction functions for the circularly polarized stochastic gravitational wave background . </S>", "<S> we find that the circular polarization can not be detected for an isotropic background . however , there is a chance to observe the circular polarization for an anisotropic gravitational wave background . </S>", "<S> we also show how to separate polarized gravitational waves from unpolarized gravitational waves . </S>"], "labels": null, "section_names": ["introduction", "stokes parameters for plane gravitational waves", "formulation", "generalized overlap reduction function for circular polarization", "extraction of circularly polarized components", "conclusion", "generalized overlap reduction function for circular polarization: dipole cases", "generalized overlap reduction function for intensity"], "sections": [["it is believed that the direct detection of gravitational waves ( gws ) will bring the era of gravitational wave astronomy .", "the interferometer detectors are now under operation and awaiting the first signal of gws  @xcite .", "it is also known that pulsar timing arrays ( ptas ) can be used as a detector for gws @xcite .", "these detectors are used to search for very low frequency ( @xmath0 ) gravitational waves , where the lower limit of the observable frequencies is determined by the inverse of total observation time @xmath1 .", "indeed , the total observation time has a crucial role in ptas , because ptas are most sensitive near the lower edge of observable frequencies @xcite . taking into account its sensitivity ,", "the first direct detection of the gravitational waves might be achieved by ptas .", "the main target of ptas is the stochastic gravitational wave background ( sgwb ) generated by a large number of unresolved sources with the astrophysical origin or the cosmological origin in the early universe .", "the promising sources are super massive black hole binaries  @xcite , cosmic ( super)string  @xcite , and inflation  @xcite .", "previous studies have assumed that the sgwb is isotropic and unpolarized  @xcite .", "these assumptions are reasonable for the primary detection of the sgwb , but the deviation from the isotropy and the polarizations should have rich information of sources of gravitational waves .", "recently , the cross - correlation formalism has been generalized to deal with anisotropy in the sgwb @xcite .", "result of this work enables us to consider arbitrary levels of anisotropy , and a bayesian approach was performed by using this formalism @xcite . on the other hand , for the anisotropy of the sgwb , the cross - correlation formalism has been also developed in the case of interferometer detectors  @xcite .", "as to the polarization , there are works including the ones motivated by the modified gravity  @xcite", ". we can envisage supermassive black hole binaries emit circularly polarized sgwb due to the chern - simons term  @xcite .", "there may also exist cosmological sgwb with circular polarization in the presence of parity violating term in gravity sector  @xcite .    in this paper", ", we investigate the detectability of circular polarization in the sgwb by ptas .", "we characterize sgwb by the so called stokes @xmath2 parameter  @xcite and calculate generalized overlap reduction functions ( orfs ) so that we can probe the circular polarization of the sgwb .", "we also discuss a method to separate the intensity ( @xmath3 mode ) and circular polarization ( @xmath2 mode ) of the sgwb .", "the paper is organized as follows . in section [ sec :", "stokes parameters for a plane gravitational wave ] , we introduce the stokes parameters for monochromatic plane gravitational waves , and clarify the physical meaning of the stokes parameters @xmath3 and @xmath2 . in section [ sec : formulation ] , we formulate the cross - correlation formalism for anisotropic circularly polarized sgwb with ptas .", "the basic framework is essentially a combination of the formalism of @xcite , and the polarization decomposition formula of the sgwb derived in @xcite . in section [ sec : the generalized overlap reduction function for circular polarization ] , we calculate the generalized orfs for the @xmath2 mode .", "the results for @xmath3 mode are consistent with the previous work  @xcite . in section [ sec : separation method ] , we give a method for separation between the @xmath3 mode and @xmath2 mode of the sgwb .", "the final section is devoted to the conclusion . in appendixes , we present analytic results for the generalized overlap reduction functions . in this paper", ", we will use the gravitational units @xmath4 ."], ["let us consider the stokes parameters for plane waves traveling in the direction @xmath5 , which can be described by @xmath6 \\", ", \\\\ & & h_{xy}(t , z)=h_{yx}(t , z)={\\rm re}[b_{\\times}\\mathrm{e}^{-iw(t - z ) } ] \\ .\\end{aligned}\\ ] ] for an idealized monochromatic plane wave , complex amplitudes @xmath7 and @xmath8 are constants .", "polarization of the plane gws is characterized by the tensor , ( see @xcite and also electromagnetic case @xcite ) @xmath9 where @xmath10 take @xmath11 .", "any @xmath12 hermitian matrix can be expanded by the pauli and the unit matrices with real coefficients .", "hence , the @xmath13 hermitian matrix @xmath14 can be written as @xmath15 where @xmath16 by analogy with electromagnetic cases , @xmath17 and @xmath2 are called stokes parameters . comparing with , we can read off the stokes parameters as @xmath18= b_{+}^{\\ast}b_{\\times}+ b_{\\times}^{\\ast}b_{+},\\\\ v&=&-2{\\rm i m } [ b_{+}^{\\ast}b_{\\times}]=i ( b_{+}^{\\ast}b_{\\times}- b_{\\times}^{\\ast}b_{+}).\\label{stv}\\end{aligned}\\ ] ] apparently , the real parameter @xmath3 is the intensity of gws . in order to reveal the physical meaning of the real parameter @xmath2 , we define the circular polarization bases @xcite @xmath19 from the relation @xmath20 we see @xmath21", "thus , we can rewrite the stokes parameters - as @xmath22 from the above expression , we see that the real parameter @xmath2 characterizes the asymmetry of circular polarization amplitudes .", "the other parameters @xmath23 and @xmath24 have additional information about linear polarizations by analogy with the electromagnetic cases .", "alternatively , we can also define the tensor @xmath25 in circular polarization bases @xmath26 where @xmath27 .", "note that the stokes parameters satisfy a relation @xmath28    next , we consider the transformation of the stokes parameters under rotations around the @xmath5 axis . the rotation around the @xmath5 axis is given by @xmath29 where @xmath30 is the angle of the rotation .", "the gws traveling in the direction @xmath5 @xmath31 transform as @xmath32 where we took the transverse traceless gauge @xmath33 after a short calculation , we obtain @xmath34 using and , the four stokes parameters ( [ sti])-([stv ] ) transform as @xmath35 as you can see , the parameters @xmath23 and @xmath24 depend on the rotation angle @xmath30 .", "this reflects the fact that @xmath23 and @xmath24 parameters characterize linear polarizations .", "note that this transformation is similar to the transformation of electromagnetic case except for the angle @xmath36 and can be rewritten as @xmath37"], ["in this section , we study anisotropic distribution of sgwb and focus on the detectability of circular polarizations with pulsar timing arrays .", "we combine the analysis of @xcite and that of @xcite . in sec.[subsec : the spectral ] , we derive the power spectral density for anisotropic circularly polarized sgwb @xmath38 .", "then we also derive the dimensionless density parameter @xmath39 which is expressed by the frequency spectrum of intensity @xmath40  @xcite . in sec.[subsec : the signal ] , we extend the generalized orfs to cases with circular polarizations characterized by the parameter @xmath2 . for simplicity ,", "we consider specific anisotropic patterns with @xmath41 expressed by the spherical harmonics @xmath42 .      in the transverse traceless gauge , metric perturbations @xmath43 with a given propagation direction @xmath44", "can be expanded as @xcite @xmath45 where the fourier amplitude satisfies @xmath46 as a consequence of the reality of @xmath43 , @xmath47 , @xmath48 is the frequency of the gws , @xmath49 are spatial indices , @xmath50 label polarizations .", "note that the fourier amplitude @xmath51 satisfies the relation @xmath52 where @xmath53 was defined by .", "the polarized tensors @xmath54 are defined by @xmath55 where @xmath56 and @xmath57 are unit orthogonal vectors perpendicular to @xmath58 .", "the polarization tensors satisfy @xmath59 with polar coordinates , the direction @xmath44 can be represented by @xmath60 and the polarization basis vectors read @xmath61    we assume the fourier amplitudes @xmath62 are random variables , which is stationary and gaussian .", "however , they are not isotropic and unpolarized .", "the ensemble average of fourier amplitudes can be written as @xcite @xmath63 where @xmath64 here , the bracket @xmath65 represents an ensemble average , and @xmath66 is the dirac delta function on the two - sphere .", "the gw power spectral density @xmath38 is a hermitian matrix , and satisfies @xmath67 because of the relation @xmath46 .", "therefore , we have the relations @xmath68 note that the stokes parameters are not exactly the same as the expression of , but they have the relation and characterize the same polarization .", "we further assume that the sgwbs satisfy @xmath69 we also assume the directional dependence of the sgwb is frequency independent @xcite .", "this implies the gw power spectral density is factorized into two parts , one of which depends on the direction while the other depends on the frequency .", "because of the transformations - , the parameters @xmath3 and @xmath2 have spin 0 and the parameters @xmath70 have spin @xmath71  @xcite . to analyze the sgwb on the sky , it is convenient to expand the stokes parameters by spherical harmonics @xmath72", ". however , since @xmath70 parameters have spin @xmath71 , they have to be expanded by the spin - weighted harmonics @xmath73 @xcite .", "thus , we obtain @xmath74 in this paper , we study specific anisotropic patterns with @xmath41 for simplicity .", "therefore , we can neglect @xmath23 and @xmath24 from now on .", "thus , the gw power spectral density becomes @xmath75 where @xmath76 so , we focus on the parameters @xmath3 and @xmath2 . in what follows , we will use the following shorthand notation @xmath77    next , we consider the dimensionless density parameter  @xcite @xmath78 where @xmath79 is the critical density , @xmath80 is the present value of the hubble parameter , @xmath81 is the energy density of gravitational waves , and @xmath82 is the energy density in the frequency range @xmath48 to @xmath83 .", "the bracket @xmath65 represents the ensemble average .", "however , actually , we take a spatial average over the wave lengths @xmath84 of gws or a temporal average over the periods @xmath85 of gws . here", ", we assumed the ergodicity , namely , the ensemble average can be replaced by the temporal average .", "using , , , as well as @xmath46 and @xmath86 , we get @xmath87 then we define @xmath88 hence , the dimensionless quantity @xmath39 in is given by @xmath89 where the spherical harmonics are orthogonal and normalized as @xmath90 using @xmath91 , we obtain @xmath92 without loss of generality , we normalize the monopole moment as @xmath93 so , becomes @xmath94      the time of arrival of radio pulses from the pulsar is affected by gws .", "consider a pulsar with frequency @xmath95 located in the direction @xmath96 . to detect the sgwb ,", "let us consider the redshift of the pulse from a pulsar @xcite @xmath97 where @xmath98 is a frequency detected at the earth and @xmath96 is the direction to the pulsar .", "the unit vector @xmath44 represents the direction of propagation of gravitational plane waves .", "we also defined the difference between the metric perturbations at the pulsar @xmath99 and at the earth @xmath100 as @xmath101 the gravitational plane waves at each point is defined as @xmath102 for the sgwb , the redshift have to be integrated over the direction of propagation of the gravitational waves @xmath44 : @xmath103 we choose a coordinate system @xmath104 and assume that the amplitudes of the metric perturbation at the pulsar and the earth are the same .", "then becomes @xmath105 and therefore , reads @xmath106 where we have defined the pattern functions for pulsars @xmath107 note that our convention for the fourier transformation is @xmath108 therefore , the fourier transformation of can be written as @xmath109    in the actual signals from a pulsar , there exist noises .", "hence , we need to use the correlation analysis .", "we consider the signals from two pulsars @xmath110 where @xmath111 labels the pulsar . here", ", @xmath112 denotes the signal from the pulsar and @xmath113 denotes the noise intrinsic to the measurement .", "we assume the noises are stationary , gaussian and are not correlated between the two pulsars .", "to correlate the signals of two measurements , we define @xmath114 where @xmath1 is the total observation time and @xmath115 is a real filter function which should be optimal to maximize signal - to - noise ratio . in the case of interferometer", ", the optimal filter function falls to zero for large @xmath116 compered to the travel time of the light between the detecters .", "since the signals of two detectors are expected to correlate due to the same effect of the gravitational waves , the optimal filter function should behave this way .", "then , typically one of the detectors is very close to the other compared to the total observation time @xmath1 .", "therefore , the total observation time @xmath1 can be extended to @xmath117 @xcite .", "in contrast , in the case of pta , it is invalid that @xmath1 is very large compered to the travel time of the light between the pulsars .", "nevertheless , we can assume that one of the two @xmath1 can be expanded to @xmath117 , because in situations @xmath118 and @xmath119 it is known that we can ignore the effect of the distance @xmath120 of pulsars .", "in this case , it is clear that any locations of the pulsars are optimal and optimal filter function should behave like as the interferometer case @xcite .    using these assumptions @xmath118 and @xmath119 , we can rewrite as @xmath121 where @xmath122 note that @xmath123 satisfies @xmath124 , because @xmath125 is real .", "moreover , to deal with the unphysical region @xmath126 we require @xmath127 .", "thus , @xmath123 becomes real .", "taking the ensemble average , using @xmath128 , @xmath118 , and assuming the noises in the two measurements are not correlated , we get @xmath129\\ , \\label{s2}\\end{aligned}\\ ] ] where we have defined @xmath130 the functions @xmath131 and @xmath132 are called the generalized orfs , which describe the angular sensitivity of the pulsars for the sgwb . note that , as we already mentioned , we consider the cases of @xmath41 for simplicity", ". then we have assumed @xmath118 and @xmath119 , this assumption implies that approximately becomes @xmath133 due to the rapid oscillation of the phase factor .", "therefore , the distance @xmath120 of the pulsars does not appear in the generalized orfs , and hence the generalized orfs do not depend on the frequency .    as you can see from ,", "the correlation of the two measurements involve both the total intensity and the circular polarization .", "however , the degeneracy can be disentangled by using separation method , which will be discussed in the section [ sec : separation method ] ."], ["in this section , we consider the generalized orfs for circular polarizations : @xmath134 where we defined @xmath135 in the above , we have used and the fact that the generalized orfs do not depend on frequency . for computation of the generalized orfs for circular polarizations ,", "it is convenient to use the computational frame @xcite defined by @xmath136 where @xmath137 is the angular separation between the two pulsars . using - , , and", ", one can easily show that @xmath138 we therefore get @xmath139 the explicit form of the spherical harmonics reads @xmath140 where @xmath141 is the normalization factor .", "the associated legendre functions are given by @xmath142 and @xmath143 with the legendre functions @xmath144\\ .\\label{pl}\\end{aligned}\\ ] ] using the spherical harmonics , becomes @xmath145 where we have used the fact that the function of @xmath146 is odd parity in the case of @xmath147 and is even parity in the case of @xmath148 .", "note that the generalized orfs for circular polarizations are real functions . in the case of @xmath149 and/or @xmath150 ,", "the integrand in vanishes .", "therefore , we can not detect circular polarizations for these cases .", "this fact for @xmath151 implies that we do not need to consider auto - correlation for a single pulsar .", "this is the reason why we neglected auto - correlation term in .    integrating ,", "we get the following form for @xmath152 : @xmath153 for @xmath154 , we have obtained @xmath155 \\ , \\\\ \\gamma^{v}_{1 - 1}&=&\\gamma^{v}_{11 } \\ , \\end{aligned}\\ ] ] recall that @xmath156 .", "the derivation of this formula for @xmath154 can be found in appendix [ sec : angular integral of the generalized overlap reduction function for dipole circular polarization ] .    for @xmath157 , we derived the following : @xmath158\\ , \\\\ \\gamma^{v}_{2 - 1}&=&\\gamma^{v}_{21}\\ , \\\\", "\\gamma^{v}_{22}&=&-\\frac{\\sqrt{30\\pi}}{6}(1-\\cos\\xi)\\left[2-\\cos\\xi+6\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{v}_{2 - 2}&=&-\\gamma^{v}_{22}\\ , \\end{aligned}\\ ] ] for @xmath159 , the results are @xmath160\\ , \\\\ \\gamma^{v}_{3 - 1}&=&\\gamma^{v}_{31}\\ , \\\\ \\gamma^{v}_{32}&=&\\frac{\\sqrt{210\\pi}}{24}(1-\\cos\\xi)\\left[8 - 5\\cos\\xi-\\cos^2\\xi+24\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{v}_{3 - 2}&=&-\\gamma^{v}_{3 - 2}\\ , \\\\ \\gamma^{v}_{33}&=&-\\frac{\\sqrt{35\\pi}}{16}\\sin\\xi\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\left[11 - 6\\cos\\xi-\\cos^2\\xi+32\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{v}_{3 - 3}&=&\\gamma^{v}_{33}\\ .\\end{aligned}\\ ] ] in fig .", "[ gv ] , we plotted these generalized orfs as a function of the angular separation between the two pulsars @xmath137 .", "it is apparent that considering the @xmath2 mode does not make sense when we only consider the isotropic ( @xmath152 ) orf . on the other hand ,", "when we consider anisotropic ( @xmath161 ) orfs , it is worth taking into account polarizations .", "the polarizations of the sgwb would give us rich information both of super massive black hole binaries and of inflation in the early universe .", "as a function of the angular separation between the two pulsars @xmath137 . in fig .", "[ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( a ) @xmath152     as a function of the angular separation between the two pulsars @xmath137 . in fig . [ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( b ) @xmath154     as a function of the angular separation between the two pulsars @xmath137 . in fig . [ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( c ) @xmath157     as a function of the angular separation between the two pulsars @xmath137 . in fig .", "[ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .", "[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .", "[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .", "[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .", "the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title=\"fig:\",width=340 ] ( d ) @xmath159    using the same procedure described in the above to derive the generalized orfs for circular polarizations , we can also derive the generalized orfs for the intensity @xmath166 where @xmath167 the angular integral in this case was performed in @xcite .", "the results are summarized in appendix [ sec : the generalized overlap reduction function for intensity ] ."], ["in this section , we separate the @xmath3 mode and @xmath2 mode of the sgwb with correlation analysis @xcite . to this aim ,", "we use four pulsars ( actually we need at least three pulsars ) , and define correlations of @xmath168 @xmath169 where @xmath170 label the pulsars . comparing with , we obtain @xmath171 \\ ,", "\\label{1c12}\\\\ & & c_{34}(f)=\\sum_{lm}^{l=3}\\left[c_{lm}^{i}i(f)\\gamma_{lm,34}^{i}+c_{lm}^{v}v(f)\\gamma_{lm,34}^{v}\\right ] \\ .\\label{1c34}\\end{aligned}\\ ] ] if the @xmath3 mode and @xmath2 mode of the sgwb are dominated by a certain @xmath172 and @xmath173 , and become @xmath174 \\ , \\label{2c12 } \\\\ & & c_{34}(f)=\\left[c _ { l m}^{i}i(f)\\gamma _ { l m,34}^{i}+c _ { l ' m'}^{v}v(f)\\gamma _ { l ' m',34}^{v}\\right ] \\ .\\label{2c34}\\end{aligned}\\ ] ] to separate the intensity and the circular polarization , we take the following linear combinations @xmath175 where we defined coefficients @xmath176 as you can see , @xmath177 contains only @xmath40 , and @xmath178 contains only @xmath179 .", "for the signal @xmath180 , the formulas corresponding to and are given by @xmath181 \\ , \\label{sp}\\end{aligned}\\ ] ] where @xmath182 denotes @xmath3 and @xmath2 .", "we assume @xmath183 and that the noise in the four pulsars are not correlated .", "we also assume that the ensemble average of fourier amplitudes of the noises @xmath184 is of the form @xmath185 where @xmath186 is the noise power spectral density .", "the reality of @xmath187 gives rise to @xmath188 and therefore we obtain @xmath189 . without loss of generality", ", we can assume @xmath190 then we obtain corresponding noises @xmath191 : @xmath192\\ , \\label{np}\\end{aligned}\\ ] ] where @xmath193^{1/2 } \\label{sn12 } \\ , \\quad s_{n,34}(f ) \\equiv [ s_{n,3}(f)s_{n,4}(f)]^{1/2 } \\label{sn34 } \\ .\\end{aligned}\\ ] ] using the inner product @xmath194 \\ , \\end{aligned}\\ ] ] we can rewrite , as @xmath195 therefore , the optimal filter function can be chosen as @xmath196 using , we get optimal signal - to - noise ratio @xmath197^{1/2}\\ .\\label{snr}\\end{aligned}\\ ] ] plugging , , and into , we obtain @xmath198^{1/2}\\ , \\\\", "{ \\rm snr}_{v}&=&\\left[t\\int_{-\\infty}^{\\infty}df\\,\\,\\frac{\\left(c^{v}_{{l}'{m}'}\\right)^{2}v^{2}(f)\\left(\\gamma_{{l}'{m}',34}^{v}\\gamma^{i}_{{l}{m},12}-\\gamma_{{l}'{m}',12}^{v}\\gamma^{i}_{{l}{m},34}\\right)^2}{\\left(\\gamma^{i}_{{l}{m},12}\\right)^2s^{2}_{n,34}(f)+\\left(\\gamma^{i}_{{l}{m},34}\\right)^2s^{2}_{n,12}(f)}\\right]^{1/2}\\ .\\end{aligned}\\ ] ] if we assume all of the noise power spectral densities are the same , becomes @xmath199 thus , the compiled orfs can be defined as @xmath200^{1/2}}\\ , \\\\ \\gamma_{12:34}^{v}&\\equiv&\\frac{\\gamma_{{l}'{m}',34}^{v}\\gamma^{i}_{{l}{m},12}-\\gamma_{{l}'{m}',12}^{v}\\gamma^{i}_{{l}{m},34}}{\\left[\\left(\\gamma^{i}_{{l}{m},12}\\right)^2+\\left(\\gamma^{i}_{{l}{m},34}\\right)^2\\right]^{1/2}}\\ .\\end{aligned}\\ ] ] this compiled orfs @xmath201 and @xmath202 describe the angular sensitivity of the four pulsars for the pure @xmath3 and @xmath2 mode of the sgwb , respectively .", "note that , to do this separation , we must know a priori the coefficients @xmath203 and @xmath204 .", "if we do not assume , the generalized orfs depend on the frequency . in this case", ", it seems difficult to calculate these coefficients .", "we next consider the case that @xmath3 mode and/or @xmath2 mode dominant in two or more @xmath205 . in this case , if we have a priori knowledge of the values of @xmath206 in each of @xmath205 for coefficients", "@xmath203 and @xmath204 , we can separate @xmath3 mode and @xmath2 mode . for example , assume that @xmath3 mode is dominated by @xmath207 , while @xmath2 mode is dominated by @xmath208 , then and become @xmath209\\ , \\label{3c12}\\\\ & & c_{34}(f)=\\left[c^{i}_{00}i(f)\\left(\\gamma_{00,34}^{i}+\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)+c_{11}^{v}v(f)\\gamma_{11,34}^{v}\\right]\\ .\\label{3c34}\\end{aligned}\\ ] ] thus , we can separate @xmath3 mode and @xmath2 mode by using linear combinations @xmath210\\ , \\\\", "d_{v}&\\equiv&a_{v}c_{34}(f)+b_{v}c_{12}(f ) \\nonumber\\\\ & = & c_{11}^{v}v(f)\\left[\\gamma_{11,34}^{v}\\left(\\gamma_{00,12}^{i}+\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,12}^{i}\\right)-\\gamma_{11,12}^{v}\\left(\\gamma_{00,34}^{i}+\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)\\right]\\ , \\end{aligned}\\ ] ] where @xmath211 as in the previous calculations , we can get the compiled orfs @xmath212^{1/2}}\\ , \\label{gi1234}\\\\ \\gamma_{12:34}^{v}&\\equiv&\\frac{\\gamma_{11,34}^{v}\\left(\\gamma_{00,12}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,12}^{i}\\right)-\\gamma_{11,12}^{v}\\left(\\gamma_{00,34}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)}{\\left[\\left(\\gamma_{00,12}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,12}^{i}\\right)^2+\\left(\\gamma_{00,34}^{i}+\\displaystyle\\frac{c^{i}_{11}}{c^{i}_{00}}\\gamma_{11,34}^{i}\\right)^2\\right]^{1/2}}\\ .\\label{gv1234}\\end{aligned}\\ ] ]    [ cols=\"^,^ \" , ]     in fig .", "[ cg ] we show some compiled orfs @xmath213 ( left panels ) and @xmath214 ( right panels ) as a function of the two angular separations @xmath137 and @xmath215 for two pulsar pairs , respectively .", "we used the expressions of @xmath2 mode and @xmath3 mode ( see appendix [ sec : the generalized overlap reduction function for intensity ] ) , and we assumed @xmath216 for simplicity . in fig .", "[ cg](a ) and [ cg](b ) , the @xmath3 mode is dominated by @xmath217 and @xmath2 mode is dominated by @xmath218 . in fig .", "[ cg](c ) and [ cg](d ) , the @xmath3 mode is dominated by @xmath219 and @xmath2 mode is dominated by @xmath218 . in fig .", "[ cg](e ) and [ cg](f ) , the @xmath3 mode is dominated by @xmath207 and @xmath2 mode is dominated by @xmath218 . in fig .", "[ cg](e ) and [ cg](f ) , the @xmath3 mode is dominated by @xmath220 and @xmath2 mode is dominated by @xmath218 . by definition , in the case of @xmath221 ,", "the compiled orfs are zero ."], ["we have studied the detectability of the stochastic gravitational waves with ptas . in most of the previous works ,", "the isotropy of sgwb has been assumed for the analysis .", "recently , however , a stochastic gravitational wave background with anisotropy have been considered .", "the information of the anisotropic pattern of the distribution should contain important information of the sources such as supermassive black hole binaries and the sources in the early universe .", "it is also intriguing to take into account the polarization of sgwb in the pta analysis .", "therefore , we extended the correlation analysis to circularly polarized sgwb and calculated generalized overlap reduction functions for them .", "it turned out that the circular polarization can not be detected for an isotropic background .", "however , when the distribution has anisotropy , we have shown that there is a chance to observe circular polarizations in the sgwb .", "we also discussed how to separate polarized modes from unpolarized modes of gravitational waves .", "if we have a priori knowledge of the abundance ratio for each mode in each of @xmath205 , we can separate @xmath3 mode and @xmath2 mode in general .", "this would be possible if we start from fundamental theory and calculate the spectrum of sgwb .", "in particular , in the case that the signal of lowest @xmath222 is dominant , we performed the separation of @xmath3 mode and @xmath2 mode explicitly .", "this work was supported by grants - in - aid for scientific research ( c ) no.25400251 and \" mext grant - in - aid for scientific research on innovative areas no.26104708 and `` cosmic acceleration''(no.15h05895 ) ."], ["in this appendix , we perform angular integration of the generalized orf for dipole ( @xmath154 ) circular polarization ( see @xcite ) : @xmath223 where we have defined @xmath224 .", "it is obvious that in the case of @xmath225 , integrand of the generalized orf is zero , because of @xmath226 , then we obtain @xmath227 then , using - , we calculate @xmath228 and we find @xmath229    therefore we only have to consider the dipole generalized orf in the case of @xmath154 , @xmath230 : @xmath231 where @xmath232 first , to calculate @xmath233 , we use contour integral in the complex plane . defining @xmath234 and substituting @xmath235 into , we can rewrite @xmath233 as @xmath236 } \\ , \\end{aligned}\\ ] ] where @xmath237 denotes a unit circle .", "we can factorize the denominator of the integrand and get @xmath238 where @xmath239 hereafter , the upper sign applies when @xmath240 and the lower one applies when @xmath241 .", "note that we only consider the region @xmath242 , so we have used the relation @xmath243 in above expression . in the region", "@xmath244 , @xmath245 is inside the unit circle @xmath237 except for @xmath246 and @xmath247 is outside the unit circle @xmath237 .", "now , we can perform the integral using the residue theorem @xmath248 where @xmath249 the residues inside the unit circle @xmath237 can be evaluated as @xmath250\\right\\ }             = \\frac{i(z_{+}+z_{-})}{2\\sqrt{1-x^2}\\sin\\xi } \\ , \\end{aligned}\\ ] ] @xmath251 thus , we obtain @xmath252 next , we consider @xmath253 defined in", ". using , we can calculate @xmath253 as @xmath254    similarly , we can evaluate @xmath255 given in . to calculate @xmath255 in the complex plane , we again substitute into and obtain @xmath256 we use the residue theorem @xmath257 where @xmath258 the residues inside the unit circle @xmath237 can be calculated as @xmath259\\right\\ }        = \\frac{i(z_{+}^2+z_{-}^2)}{4\\sqrt{1-x^2}\\sin\\xi } \\ , \\end{aligned}\\ ] ] @xmath260 therefore , @xmath255 becomes @xmath261 substituting to , we can calculate @xmath262 : @xmath263    finally , substituting and into , we get the generalized orf for @xmath264 @xmath265\\ .\\end{aligned}\\ ] ]     as a function of the angular separation between the two pulsars @xmath137 .", "[ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) .", "the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( a ) @xmath152     as a function of the angular separation between the two pulsars @xmath137 .", "[ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) .", "the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( b ) @xmath154     as a function of the angular separation between the two pulsars @xmath137 . fig . [ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) .", "the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( c ) @xmath157     as a function of the angular separation between the two pulsars @xmath137 .", "[ gi](a ) shows monopole ( l=0 ) , fig .", "[ gi](b ) shows dipole ( l=1 ) , fig .", "[ gi](c ) shows quadrupole ( l=2 ) and fig .", "[ gi](d ) shows octupole ( l=3 ) . the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title=\"fig:\",width=340 ] ( d ) @xmath159"], ["in this appendix , we show orfs for the intensity @xcite .", "the following form for @xmath152 was derived in @xcite , and our expressions are identical to their expressions : @xmath271\\ , \\end{aligned}\\ ] ] for , @xmath154 , we calculated as @xmath272\\ , \\\\ \\gamma^{i}_{11}&=&\\frac{\\sqrt{6\\pi}}{12}\\sin\\xi\\left[1 + 3(1-\\cos\\xi)\\left\\{1+\\frac{4}{1+\\cos\\xi}\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right\\}\\right]\\ , \\\\ \\gamma^{i}_{1 - 1}&=&-\\gamma^{i}_{11}\\ , \\end{aligned}\\ ] ] for @xmath157 , we obtain @xmath273\\ , \\\\ \\gamma^{i}_{21}&=&-\\frac{\\sqrt{30\\pi}}{60}\\sin\\xi\\left[21 - 15\\cos\\xi-5\\cos^2\\xi+60\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{2 - 1}&=&-\\gamma^{i}_{2 - 1}\\ , \\\\ \\gamma^{i}_{22}&=&\\frac{\\sqrt{30\\pi}}{24}(1-\\cos\\xi)\\left[9 - 4\\cos\\xi-\\cos^2\\xi+24\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{2 - 2}&=&\\gamma^{i}_{22}\\ ,", "\\end{aligned}\\ ] ] for @xmath159 , it is straightforward to reach the following @xmath274\\ , \\\\ \\gamma^{i}_{31}&=&\\frac{\\sqrt{21\\pi}}{48}\\sin\\xi(1-\\cos\\xi)\\left[34 + 15\\cos\\xi+5\\cos^2\\xi+\\frac{96}{1+\\cos\\xi}\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{3 - 1}&=&-\\gamma^{i}_{31}\\ , \\\\ \\gamma^{i}_{32}&=&-\\frac{\\sqrt{210\\pi}}{48}(1-\\cos\\xi)\\left[17 - 9\\cos\\xi-3\\cos^2\\xi-\\cos^3\\xi+48\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{3 - 2}&=&\\gamma^{i}_{32}\\ , \\\\ \\gamma^{i}_{33}&=&\\frac{\\sqrt{35\\pi}}{48}\\frac{(1-\\cos\\xi)^2}{\\sin\\xi}\\left[34 - 17\\cos\\xi-4\\cos^2\\xi-\\cos^3\\xi+96\\left(\\frac{1-\\cos\\xi}{1+\\cos\\xi}\\right)\\log\\left(\\sin\\frac{\\xi}{2}\\right)\\right]\\ , \\\\ \\gamma^{i}_{3 - 3}&=&-\\gamma^{i}_{33}\\ .\\end{aligned}\\ ] ] these are plotted in fig .", "the generalized orfs of total intensity are different from that of circular polarization in that the value for @xmath149 is non - trivial .", "then the @xmath3 mode orfs for @xmath275 have value even in the case of @xmath151 .", "this implies that we can consider auto - correlation for a single pulsar .", "99 j.  aasi _ et al .", "_ [ ligo scientific collaboration ] , class .", "grav .   * 32 * , 074001 ( 2015 ) doi:10.1088/0264 - 9381/32/7/074001 [ arxiv:1411.4547 [ gr - qc ] ] .", "f.  acernese _ et al . _", "[ virgo collaboration ] , class .", "grav .   * 32 * , no . 2 , 024001 ( 2015 ) doi:10.1088/0264 - 9381/32/2/024001 [ arxiv:1408.3978 [ gr - qc ] ] . k.  somiya [ kagra collaboration ] , class .", "* 29 * , 124007 ( 2012 ) doi:10.1088/0264 - 9381/29/12/124007 [ arxiv:1111.7185 [ gr - qc ] ] .", "s.  l.  detweiler , astrophys .", "j.   * 234 * , 1100 ( 1979 ) .", "doi:10.1086/157593 romani , r. w. ( 1989 ) .", "timing a millisecond pulsar array . timing neutron stars , 113 - 117 .", "l.  lentati _ et al .", "_ , mon .  not .", "soc .   * 453 * , 2576 ( 2015 ) [ arxiv:1504.03692 [ astro-ph.co ] ] .", "z.  arzoumanian _ et al .", "_ [ nanograv collaboration ] , arxiv:1508.03024 [ astro-ph.ga ] . e.  s.  phinney , astro - ph/0108028 . a.  vilenkin , phys .", "rept .   * 121 * , 263 ( 1985 ) .", "s.  kuroyanagi , k.  miyamoto , t.  sekiguchi , k.  takahashi and j.  silk , phys .", "d * 87 * , no . 2 , 023522 ( 2013 ) [ phys .", "d * 87 * , no .", "6 , 069903 ( 2013 ) ] [ arxiv:1210.2829 [ astro-ph.co ] ] .", "m.  maggiore , gr - qc/0008027 .", "l.  p.  grishchuk , phys .", "* 48 * , 1235 ( 2005 ) doi:10.1070/pu2005v048n12abeh005795 [ gr - qc/0504018 ] .", "b.  allen and j.  d.  romano , phys .", "d * 59 * , 102001 ( 1999 ) [ gr - qc/9710117 ] .", "c.  m.  f.  mingarelli , t.  sidery , i.  mandel and a.  vecchio , phys .", "d * 88 * , no . 6 , 062005 ( 2013 ) [ arxiv:1306.5394 [ astro-ph.he ] ] .", "s.  r.  taylor and j.  r.  gair , phys .", "d * 88 * , 084001 ( 2013 ) [ arxiv:1306.5395 [ gr - qc ] ] .", "n.  seto and a.  taruya , phys .", "lett .   * 99 * , 121101 ( 2007 ) doi:10.1103/physrevlett.99.121101 [ arxiv:0707.0535 [ astro - ph ] ] . n.  seto and a.  taruya , phys .  rev .", "d * 77 * , 103001 ( 2008 ) [ arxiv:0801.4185 [ astro - ph ] ] . s.  j.  chamberlin and x.  siemens , phys .", "d * 85 * , 082001 ( 2012 ) doi:10.1103/physrevd.85.082001 [ arxiv:1111.5661 [ astro-ph.he ] ] . j.  r.  gair , j.  d.  romano and s.  r.  taylor , phys .  rev .", "d * 92 * , no .", "10 , 102003 ( 2015 ) doi:10.1103/physrevd.92.102003 [ arxiv:1506.08668 [ gr - qc ] ] .", "r.  jackiw and s.  y.  pi , phys .", "d * 68 * , 104012 ( 2003 ) doi:10.1103/physrevd.68.104012 [ gr - qc/0308071 ] . m.  satoh , s.  kanno and j.  soda , phys .", "d * 77 * , 023526 ( 2008 ) doi:10.1103/physrevd.77.023526 [ arxiv:0706.3585 [ astro - ph ] ] . c.  r.  contaldi , j.  magueijo and l.  smolin , phys .", "* 101 * ( 2008 ) 141101 doi:10.1103/physrevlett.101.141101 [ arxiv:0806.3082 [ astro - ph ] ] .", "t.  takahashi and j.  soda , phys .", "lett .   * 102 * , 231301 ( 2009 ) doi:10.1103/physrevlett.102.231301 [ arxiv:0904.0554 [ hep - th ] ] .", "j.  l.  cook and l.  sorbo , phys .", "d * 85 * , 023534 ( 2012 ) [ phys .", "d * 86 * , 069901 ( 2012 ) ] doi:10.1103/physrevd.86.069901 , 10.1103/physrevd.85.023534 [ arxiv:1109.0022 [ astro-ph.co ] ] .", "i.  obata , t.  miura and j.  soda , phys .", "d * 92 * , no .", "6 , 063516 ( 2015 ) doi:10.1103/physrevd.92.063516 [ arxiv:1412.7620 [ hep - ph ] ] .", "a. p. lightman , w. h. press , r. h. price , and s. a. teukolski , problem book in relativity and gravitation , 2nd ed .", "( princeton university press , 1979 ) .", "m. maggiore , gravitational waves , vol . 1 : theory and experiments ( oxford university press , 2008 )", ". g. b. rybicki , a. p. lightman , radiative processes in astrophysics ( wiley - interscience , 1979 ) .", "l. d. landau and e. m. lifshitz , the classical theory ( pergamon press , 1975 ) . c. misner , k. thorne and j. wheeler , gravitation , ( freeman 1973 ) .", "b.  allen and a.  c.  ottewill , phys .", "d * 56 * , 545 ( 1997 ) doi:10.1103/physrevd.56.545 [ gr - qc/9607068 ] .", "u.  seljak and m.  zaldarriaga , phys .", "lett .   * 78 * , 2054 ( 1997 ) doi:10.1103/physrevlett.78.2054 [ astro - ph/9609169 ] .", "j. n. goldberg et al .", ", journal of mathematical physics * 8 * , 2155 , 1967 . m.  anholm , s.  ballmer , j.  d.  e.  creighton , l.  r.  price and x.  siemens , phys .", "d * 79 * , 084030 ( 2009 ) doi:10.1103/physrevd.79.084030 [ arxiv:0809.0701 [ gr - qc ] ] . l.  g.  book and e.  e.  flanagan , phys .", "d * 83 * , 024024 ( 2011 ) doi:10.1103/physrevd.83.024024 [ arxiv:1009.4192 [ astro-ph.co ] ] . f.  a.  jenet and j.  d.  romano , am .", "j.  phys .", "* 83 * , 635 ( 2015 ) doi:10.1119/1.4916358 [ arxiv:1412.1142 [ gr - qc ] ] .", "r.  w.  hellings and g.  s.  downs , astrophys .", "j.   * 265 * , l39 ( 1983 ) . doi:10.1086/183954"]]}
