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