the study of the formation , composition , mass and kinematics of galaxy halos are among the most active areas of modern astrophysical research . the formation of galaxy halos are now widely believed to be due to hierarchical structure formation ( e.g. , freeman & bland - hawthorn 2002 ) where galaxies are produced by the combination of a monolithic collapse ( eggen et al .  1962 ) and the accretion of galactic fragments ( searle and zinn 1978 ) . for the halo of the milky way
, all aspects of the halo can probed by studying the nature of the remnants of disrupted dwarf galaxies .
numerous tidal streams and dwarf galaxies have been discovered within the galactic halo in the last 20 years ( e.g. , ibata et al .
1994 ; ibata et al .  2001 ; vivas et al .  2001 ; grillmair 2006 ; newberg et al .  2009 ) .
the most well - studied of these is the sagittarius dwarf galaxy ( sgr , ibata et al .
1994 ) and its associated tidal stream .
the sgr streams have been traced on large scales using blue horizontal branch ( bhb ) stars ( newberg et al . 2003 ) , m - giants ( majewski et al .
2003 ) , and main - sequence turn off ( msto ) stars ( e.g. , belokurov et al .
2006 ; koposov et al .
the sgr system has also been studied using rr lyrae ( rrl , vivas & zinn 2006 ; miceli et al .
2008 ; sesar et al .
2010 ) . other halo streams and structures discovered recently include a virgo stellar stream ( vss ; vivas & zinn 2006 ; vivas et al .
2008 ) , a virgo overdensity ( vod , newberg et al . 2002
; newberg et al .
2007 ) , an overdensity in pisces ( sesar et al .  2007 ; kollmeier et al .
2009 ) and a monoceros stream ( newberg et al . 2002 ; majewski et al .
additionally , a cetus stream has been discovered in the south ( newberg et al .
2009 ; koposov et al .
2012 ) and also evidence for an overdensity dubbed the hercules - aquila cloud ( belokurov et al .
2007 ) .
although the sgr stream system has been the focus of much study , since the structure is extends completely around the galaxy , it is yet to be fully mapped in velocity , metallicity and distance .
furthermore , the cause of the bifurcation in sgr stream stars discovered by belokurov et al .
( 2006 ) is yet to be explained .
additionally , newberg et al .  ( 2003 ) proposed the existence of a stream of stars associated with the sgr dwarf galaxy at galactocentric distances of @xmath3 .
these authors photometrically identified candidate bhb stars in sdss data and found evidence for an overdensity of stars with @xmath4 .
this overdensity was found using sdss data covering @xmath5 , @xmath6 .
newberg et al .  ( 2003 ) attributed this feature to a trailing arm of the sgr stream , yet also noted a possible link to the globular cluster ngc  2419 .
their analysis showed evidence for a feature visible across @xmath7 along the sgr plane .
newberg et al .  ( 2007 ) continued this work and found additional evidence for what they called the sgr trailing tail using bhb candidates selected from sdss data release 5 ( dr5 ; alderman - mccarthy et al .
more recently , ruhland et al .  ( 2011 ) found additional evidence for this overdensity using candidate bhb stars selected from the larger sdss dr7 dataset .
as with newberg et al .
( 2003 , 2007 ) , these authors found a stream of bhb candidates to be located at heliocentric distances of 60 - 80 kpc with @xmath8 significance . within their analysis the feature was seen to span @xmath9 of the sky in the region @xmath10 .
ruhland et al .
( 2011 ) also compared this extended stream with five sets of published numerical simulations of the sgr stream .
they found that none were able to explain the existence of sgr stream stars in the location observed .
current evidence for an outer sgr stream / tidal - tail remains based purely on photometrically selected bhb stars .
as demonstrated by sirko et al .
( 2004 ) and ruhland et al .
( 2011 ) , it is not possible to select a 100% pure set of bhb stars using sdss photometry alone .
depending on the colour cuts used one has to accept varying levels of contamination blue straggler ( bs ) stars which far out number bhb stars at the greatest detection completeness . with strict colour cuts , that remove a large fraction of the bhb stars ,
the level of bs star contamination can be reduced to 30% ( ruhland et al .
however , to retain a large number of bhb stars , a 50% or greater contamination level must be accepted ( brown et al .  2005 ; bell et al .
2010 ) .
in contrast to photometric selection , it is possible to select bhb stars with much less contamination by combining colour selection with high s / n spectra because bhb stars have significantly lower surface gravity than the bs stars ( yanny et al .
the sdss segue - i and segue - ii projects ( yanny et al .
2009 ) undertook spectroscopy of @xmath0 300,000 star with an aim to identifying bhb stars and trace galactic halo structure .
however , while radial velocities are available for stars to @xmath11 ( yanny et al .
2009 ) , the s / n required to accurately measure @xmath12 limits the survey to bhb stars with @xmath13 ( @xmath14 kpc ) .
almost all of the bhb candidates discovered by newberg et al .
( 2003 , 2007 ) and ruhland et al .
( 2011 ) are beyond this limit and do not have measurements of @xmath12 to confirm their nature .
additionally , while the sdss segue surveys cover @xmath15 , the fields cover a patchwork over the northern sky , rather than just the sgr stream region where this feature is located .    like bhb stars ,
rrl stars also exist on horizontal branch ( hb ) . however , unlike bhb stars , rrl exhibit a significant level of intrinsic photometric variability . based on their characteristic variability rrl
can be cleanly separated from other stars .
type - ab rrl stars ( rrab ) have absolute magnitudes of @xmath16 with uncertainties of 6% ( catelan et al .
this makes them excellent distance indicators .
although , significant numbers of repeated observations are required to provide accurate average magnitudes .    to date , a few tens of thousands of rrl are known in dense regions near the galactic bulge and in the magellanic clouds due to microlensing surveys ( soszyski et al .
2009 , pietrukowicz et al .
recently we discovered @xmath0 10,000 rrab over a large fraction of the sky ( drake et al .
2013 ; dr13 ) . however ,
even with the dr13 rrab s most of the galactic halo remains unexplored at heliocentric distances @xmath17 .    in this paper
we outline our search , discovery , and calibration of rrl at distances up to 100 kpc and beyond .
we will compare our discoveries with results from past surveys and undertake an analysis of the distant hb stars uncovered .
the catalina sky survey began in 2004 and uses three telescopes to cover the sky between declination @xmath18 and @xmath19 in order to discover near earth objects ( neos ) and potential hazardous asteroids ( larson et al .
the three telescopes are each considered sub - surveys .
these consist of the catalina schmidt survey ( css ) , the mount lemmon survey ( mls ) and the siding spring survey ( sss ) .
in addition to avoiding high declinations , the galactic plane region is avoided by between 10 and 15 degrees due to reduced source recovery in crowded stellar regions .
all of the survey images are taken in sequences of four observations separated by ten minutes , and all the observations are unfiltered to maximize throughput .
photometry is carried out using the aperture photometry program sextractor ( bertin & arnouts 1996 ) .
all the catalina data is also analyzed for transient sources by the catalina real - time transient survey ( crts , drake et al .
2009 ; djorgovski et al .
2011 ) .    in paper - i ( dr13 ) , we concentrated on analyzing rrab data from the css sub - survey . in this paper
we mainly work with data observed by the mls 1.5 m telescope .
this data predominately covers ecliptic latitudes @xmath20 . for this telescope each image from the @xmath21 catalina ccd camera covers 1.1 @xmath22 on the sky .
observations are taken 21 nights per lunation avoiding bright time .
typical exposures of 30 seconds to reach @xmath23 . in total , the mls source catalog consists of 155 million sources .
as noted in dr13 , the mls data has the same photometric sensitivity as the css data since all observations are taken unfiltered with the same type of ccd camera .
all images are processed with the same software . from dr13 , the colour transformation from catalina photometric system @xmath24 magnitudes to johnson @xmath25
is given by :    @xmath26    the average b - v colour of rrl is about 0.3 mag with stars varying between about 0.1 and 0.5 as they pulsate ( e.g. , clube et al .
1969 ; stepien 1972 ; cacciari et al . 1987
; layden 1997 ; nemec 2004 ) .
as the mls data has no colour information , we adopt the average values in our analysis .
this gives rise to an maximum uncertainty of @xmath27 mags in the @xmath25 magnitudes . combining this with the photometric uncertainty
, we expect a colour - based dispersion of @xmath28 mags in our rrl photometry .      to find rrl candidates we follow the analysis undertaken in dr13 .
that is , we first selected variable sources using a welch - stetson variability index @xmath29 ( welch & stetson 1993 ) .
3.1 million mls sources were selected as variable candidates using @xmath30 .
all sources were processed for periodicity using the lomb - scargle ( ls , lomb 1976 ; scargle 1982 ) periodograms . in total @xmath31 thousand sources
were found to exhibit significant periodicity with a false alarm probability @xmath32
. however , the bulk of these detections are due to the systematic sampling of the data and occur at integer day frequencies .    here , as with dr13 , we are primarily interested in rrab s since they have well known magnitudes and characteristic light curves that are not easily mistaken for other periodic variables . after excluding the period aliases near 0.5 and 1 day , 3087 of the objects found to have periods between 0.34 and 1.4 days .
this period region was deliberately chosen to be broader than the rrab range in order to include sources found by the ls technique at aliases of their true periods .
as with the css sources of dr13 , each of the rrab candidate was run through the afd software ( torreabla et al .
2012 ) to select the best period from among the best 15 given by the ls and aov software ( schwarzenberg - czerny 1989 ) .
briefly , this process involves fourier fitting and iteratively rejecting bad data to determine the best period based on reduced @xmath33 values .
non - rrab sources are rejected using the m - test statistic ( kinemuchi et al ,  2006 ; equation 8) , as well as the fourier fit order .
this initial selection resulted in the detection of 1125 rrab candidates .      upon reviewing the light curves of the rrab candidates we discovered a number of faint rrab in the mls data with @xmath34 .
such rrl were found to be mainly concentrated in the gemini constellation . in figure
[ mlsrrab ] , we plot the phased light curves of eight of these faint stars .
these distant sources are very important for defining distant streams within the galactic halo . from simulations , based on adding artificial rrab light curves to our data
, we discovered that our detection sensitivity was @xmath35 for rrab s with @xmath36 .    to improve our completeness
we investigated the light curves of mls sources with sdss spectroscopic types a0 and f5 that had @xmath37 .
in addition , we used sdss photometry to select stars with the same colour range as the rrab stars that we had discovered .
these stars were matched with mls sources at lower variability and periodicity thresholds ( @xmath38 and @xmath39 ) to select rrab candidates .
to better constrain the variability , we combined the mls photometry with shallower css data .
likewise , we searched the css data for additional faint rrab candidates by selecting objects with rrab colours in sdss photometry . furthermore , since our analysis in dr13 showed that @xmath40 of rrab s were missed because the selection based on ls periods , we reprocessed the css and mls photometry using the aov period - finding software .
we inspected the light curve of each rrl candidate to assess variability and then carried out period finding for all the new variable candidates . in total , approximately six thousand additional sources were searched .
each of the new sources exhibiting some sign of periodicity were searched for improved periods using our fourier - fitting afd process . in a number of cases the rrl
were found to exhibit blazhko phase and amplitude variations ( blazhko 1907 ) .
in total the additional searches yielded an 219 extra mls rrl and 2051 css rrab s .
of these only 425 were previously known . as distant rrl
are important for defining halo streams we retain the 17 clear c - type rrl ( rrc ) that were discovered in mls data beyond 70 kpc .
a number of these have sdss spectra . after removing non - rrl sources
the final set of mls rrl consists of 1,207 stars and is presented in table 1 .
the new rrab s found in css data are given in table 2 .    in figure
[ mlshist ] , we present the magnitude distribution of mls and the css rrl .
the mls histogram shows clear bumps near @xmath41 and @xmath42 .
the first bump is easy to understand as it is clearly due to the sagittarius leading and trailing arms ( as seen in dr13 ) .
the second peak is due to the more distant rrl that are concentrated near the galactic anti - center .
we hereafter call this the gemini stream . in section [ sagstream ] , we will investigate the origin of this feature .
the absolute magnitudes of rrab s are given by catelan & corts ( 2008 ) :    @xmath43}_{\rm zw84}+ 0.948,\ ] ]    where @xmath44}_{\rm zw84}$ ] is the metallicity in the zinn & west ( 1984 ) scale .
the average metallicity for rrab s with sdss spectra found by dr13 was @xmath45= -1.55 $ ] . the comparison given in dr13 , between distant halo rrab s ( @xmath46 kpc ) and the full set of rrab s showed little difference in the distribution or average metalicity . however , css rrab s with sdss metallicities are stars with @xmath47 , corresponding to distance of @xmath48 kpc
. brighter foreground field rrab s are likely to have higher metallicities .
the mls rrab s should have a similar metallicity distribution to the larger sample of css sources discovered in dr13 .
therefore , we adopt the average css metallicity and assume rrab s have magnitude @xmath49 . this value is close to the value of 0.6 that is often adopted ( e.g. , keller et al .
2008 , sesar et al .
the rrl v magnitudes were corrected for extinction using schlegel et al .
( 1998 ) reddening maps .
the dispersion in the metallicity from dr13 is 0.3 dex , which corresponds to a variation of 0.07 magnitudes .
the distances to individual sources are determined using :    @xmath50    here we correct the average rrab @xmath51 magnitudes to static values @xmath52 using values derived from a polynomial fit to the amplitude corrections given by bono et al .
( 1995 ) . combining the uncertainties from the photometric calibration and colour variation of 0.09 , with the variation in metallicity
we derive an overall uncertainty of 0.11 magnitudes , corresponding to a @xmath53 uncertainty in distance . as demonstrated by dr13 , faint rrl have larger uncertainties in their average magnitude .
however , these uncertainties should generally not exceed 0.2 magnitudes ( @xmath54 in distance ) .
the faintest mls rrab in our dataset have @xmath55 corresponding to heliocentric distances of 120 kpc . in figure
[ dist ] , we plot the distribution of the distances of mls rrl and mark the locations of known halo features as well as the gemini stream .
we matched the mls rrab s with past rrl discoveries using the simbad and the international variable star index ( vsx ; watson et al .
2006 ) databases , as well as those presented by dr13 . in total
, there were 605 matches with past surveys of which 316 came from css data ( dr13 ) and 289 were discovered in other surveys . as there is overlap between css and the other previously known rrab s ,
there was css photometry for 540 of the mls rrab s .
of the 289 simbad and vsx rrab s , 278 of the rrab have known periods given by vsx . in figure
[ other ] , we compare the periods derived from mls data with previous determinations .
the figure shows that the agreement between the css and mls periods is much better than with other sources . assuming the differences in periods are normally distributed , the mls - css matches have @xmath56 , while the mls - vsx source matches have @xmath57 . as the level of agreement is good in both cases , we have confidence in the mls rrab periods .
the sekbo survey discovered 2016 candidate rr lyrae in a survey covering 1675 @xmath22 along the ecliptic ( keller et al .
like the mls survey , the sekbo survey covered a band within @xmath58 of the ecliptic .
however , sekbo rrl candidate selection was based on the two colours observed simultaneously by the macho camera and variability was determined from between 2 and 8 epochs of images .
in contrast , the mls rrab candidates are observed unfiltered and have an average of 130 observations ( sampling ranges from 42 to 342 measurements ) .
we matched the mls rrab with the sekbo rrl candidates and found only 103 matches .
however , this is not unexpected since the sampling of mls is highly concentrated within few degrees of the ecliptic with few observations having been taken near the @xmath59 limit .
for example , 70% of the mls rrab s are within @xmath60 of the ecliptic , compared to 27% of the sekbo rrl candidates .    to extend the comparison
, we also matched the sekbo candidates with rrab s from the dr13 sample and the sss rrab s ( torrealba et al .
this yielded an additional 538 rrab matches .
because of overlap between the css , mls and sss rrab catalogs there were 540 unique rrab stars matching sekbo rrl candidates .
clearly this number falls far short of the number of sekbo candidates .
we compared average v magnitudes for the catalina sources with those of the sekbo sources . in figure
[ sekbo ] , we present the v magnitude differences between the measurements .
the average difference between magnitudes is -0.02 magnitudes with @xmath61 magnitudes , corresponding to a @xmath62 uncertainty in source distances .
the internal dispersion for css rrab magnitudes from overlapping fields is @xmath63 mags and thus does not significantly contribute to the observed dispersion .    in order to investigate
the sekbo candidates not found in the catalina rrab catalogs , we matched the remaining 1481 unmatched sekbo objects to the sdss dr8 spectroscopic catalog .
we found 40 matches , of these 35 had measured @xmath12 values and 23 had @xmath64 consistent with rrl ( dr13 ) .
we also matched the sekbo sources with sdss dr8 photometry .
of the 276 sources with sdss photometry , approximately one third had sdss colours beyond the limits observed for dr13 rrab s .
next , we matched the 1481 sekbo sources to photometry from mls , css and sss surveys . among the matches ,
we sub - selected the objects exhibiting significant variability based on the @xmath30 and also objects with ls periodic significance @xmath32 . in total
, we find matches for 90% of the sekbo objects , including 88% of the sources that were not selected as rrab s .
combining the known catalina rrab s and other sources with apparent variability , we find that 74.5% of the 1833 sekbo sources we cover exhibit significant variability . in table 3 , we present the numbers of matches with each of the catalina surveys , as well as the number of objects selected by our variability and periodicity significance thresholds .    for each of the 1314 sekbo rrl candidates with catalina photometry , that were not in our rrab catalogs , we determined a ls period ( regardless of their @xmath29 variability ) .
we inspected both the observed and phase folded light curves of each source and determined a classification based on this photometry .
as there are 2759 catalina light curves matching the 1314 sekbo candidates , the majority of these sources have multiple light curves .
we discovered 551 additional rrl candidates , 140 eclipsing binary candidates ( mainly wuma ) and 140 other types of variable sources ( including qso , cvs , @xmath65 scuti variables and objects where the classification was unclear ) .
this discovery of many non - periodic sources is not unexpected since a number of these are already known , including qsos and cataclysmic variables .
for example , sekbo 106646.2532 ( css080623:140454 - 102702 ) is a cv that was discovered in outburst by crts in 2008 ( kato et al .
2012 ) . additionally , based on sdss dr8 spectroscopy ,
sekbo-105832.627 is a cv , sekbo-117146.1412 is a qso at z=1.02 , and sekbo-096514.1360 is an eclipsing wd - md binary system ( rebassa - mansergas et al .
2010 ) .    in order to obtain the best periods for the clear rrl , and to find rrl that may have been missed during our inspection , we ran every light curve through the afd software ( torrealba et al .
for 307 of the 551 rrl candidates we found good fourier fits and periods consistent with rrab s and rrc s . in table 4 , we present all the new sekbo rrl found using catalina data ( including rrl candidates where the period is uncertain ) .
this table contains 263 rrab , 282 rrc , and 11 rrd candidates .    combining the new rrab s with the 540 from catalina catalogs we find that @xmath66 were missing from our catalogs in good agreement with the simulations presented in dr13 .
the overall sample suggests that 24% of the rrl are rrc stars whereas keller et al .
( 2006 ) estimated that @xmath67 of the sekbo candidates would be rrc s .
however , we note that pietrukowicz et al .
( 2012 ) found @xmath68 rrc among their @xmath69 bulge rrl .
nevertheless , it is possible that some fraction of the rrc candidates we identify may be w  uma variables ( due to the similarity of their light curves ) . the uncertain separation of rrc s and w  uma types is the main reason why rrc s are often excluded as distance indicators .
our results suggest that @xmath70 of the sekbo rrl candidates are likely rrl based on the 1833 objects covered by catalina data .
this is @xmath71 smaller than the @xmath72 non - rrl contamination estimated by prior et al .
( 2009a ) . however , we note that the prior et al .
( 2009a ) sample was based on an average of just nine photometric measurements for @xmath73 ( @xmath53 ) of the sekbo stars .
prior et al .  ( 2009b ) followed 21 additional sekbo rrls overlapping the sgr tidal stream .
more recently , akhter et al .
( 2012 ) followed 137 sekbo candidates and found 57 to be rrls , although they found a high fraction with colours matching the colours of rrl candidates from sdss data .
from catalina photometry we have an average of 121 observations per light curve ( 250 per object ) covering 90% of the sekbo sources .
our analysis suggests that @xmath74 of the sekbo sources are either not variable or have significantly less variability than expected for rrl .
nevertheless , with @xmath75 likely rrl from the sekbo survey , the results shows that large numbers of rrl can be found using very small numbers of observations and colour selection .
the outer halo rrl discovered in this analysis overlap with tidal stream first identified by newberg et al .
( 2003 ) using photometrically selected bhb candidates . as noted above , additional bhb candidates have since been found is sdss dr5 ( newberg et al .  2007 ) and sdss dr7 data ( ruhland et al .
furthermore , ivezic et al .  (
2004 ) discovered many rrl candidates overlapping the sgr stream plane .
the most distant sources coincide with the newberg et al .
( 2003 ) bhb candidates . however , the ivezic et al .
( 2004 ) rrl sample was based on photometric selection combined with limited variability information . of the 1,269 ivezic
( 2004 ) rr lyrae candidates selected in the sdss stripe-82 region , only 483 ( 38% ) were eventually confirmed by sesar et al .
it seems likely the a similar fraction of the sources found in this marginally significant ( @xmath76 ) detection by ivezic et al .
( 2004 ) were indeed rrl .
recently , additional evidence for a distant group of rrl stars was found by sesar et al .
( 2012a ) .
these authors discovered eight rrab between @xmath77 and @xmath78 at heliocentric distances @xmath79 kpc .
this region lies on the edge of the fields where bhb candidates were earlier discovered by newberg et al .
sesar et al .
( 2012a ) selected potential rrab candidates based on variable palomar transient factory ( ptf ) sources with rrl colours based on sdss dr8 photometry . using simulations sesar et al .
( 2012a ) estimated that in their worst - case scenario their search recovered 95% of their simulated sources , with at most one rrab star being missed .
six of the sesar et al .
( 2012a ) rrab s were independently discovered in our searches , while the remaining two were missed due to poor sampling .
for all six known rrab , we find periods matching those given by sesar et al .
( 2012a ) to better than @xmath80 .
in addition to the six ptf rrab , we have discovered seven rrab within the region and distant range noted by sesar et al .
we find 16 more rrab s in this region at heliocentric distances @xmath81 kpc . as we missed two of the ptf rrab s it is possible that there are even more rrab s within this region . in figure
[ ptf ] , we present the locations of the rrab s found within the ptf praesepe fields . among the set of eight rrab s , sesar et al .
( 2012 ) found tentative evidence for two separate groups of rrl based on radial velocities .
one group of four rrab s was noted as having radial velocities in the galactic standard of rest of @xmath82 , while the other had @xmath83 .
although sesar et al .  ( 2012 ) note that there is a 37% chance that these stars are drawn from the same gaussian velocity distribution based on the shapiro & wilk ( 1965 ) sw statistic , they assume that their eight sources are representative of the overall velocity distribution and via numerical simulations find a @xmath84 chance that the sources would exhibit the observed velocities .
these authors also discovered that the metallicities of the two groups were consistent to within uncertainties . in total , we identify 103 rrab s at heliocentric distances beyond 70 kpc .
four of these rrab s appear near @xmath85 , @xmath86 kpc .
these stars are associated with the pisces stream ( sesar et al .  2007 ; watkins et al .
2009 ) . however , we see no evidence for variations in distance . in agreement with the results of kollmeier
et al .  ( 2009 ) .
this supports evidence that this system is not part of a tidal stream itself . of the remaining stars , 82 rrab s lie within the range @xmath87 and @xmath88 .
these objects exhibit a range of distances and thus comprise significant evidence for a stream of stars , suggesting association with a tidally disrupted galaxy .    to better understand the nature of the relation between the distant rrab s and the sgr stream we selected the rrab sources within @xmath89 of the sgr plane ( as defined by majewski et al .
2003 ) . based on koposov et al .
( 2012 ) and dr13 .
most sgr stream stars , including rrab s , lie within these limits .
we combine rrab s from dr13 along with new css rab s and mls rrab s in the range @xmath90 within the sgr streams region .
we transform the rrab distances to galactocentric values and plot the sources in figure [ distg ] . in this plot
the distant stream of rrab s becomes clearer , as does the leading and trailing sgr streams that are observed between @xmath91 and @xmath92 , respectively .    among the set of distant mls rrab s
, the four most distant objects are found within the @xmath93 region @xmath94 , @xmath95 .
these rrab s have distances @xmath96 kpc , their average magnitudes varying by @xmath97 mags and are consistent with no difference at their observed brightness .
this concentrated group is significantly fainter than any of the other mls rrab s , suggesting a possible association between them that varies from the other sources .
however , measurements of radial velocities and metallicities are required before we can rule out association with the other rrab s .    after excluding the group of four very distant rrab s , and four outlier - rrab s near @xmath98kpc with @xmath99 we determine the slope of the outer stream by selecting the rrab s with @xmath100 kpc , in the range @xmath101 .
a simple linear fit gives @xmath102 .
the slope of this line suggests that the distant stream does not meet with the so - called trailing sgr stream lying in the region @xmath103 .    in order to better visualize the direction and slope of the streams , in figure [ polar ] , we present a polar plot of the same rrab sources in the sgr plane coordinate system .
in contrast to figure [ distg ] , here we plot average static star magnitudes for the rrab fainter than @xmath104 ( @xmath105 ) .
clearly a more accurate description of path of this outer stream requires more rrab in the region @xmath106 , which is not covered by the mls photometry .
as noted earlier , many new css rrab were found in this region .
however , as shown in figure [ mlshist ] , and dr13 , rrab s fainter than @xmath107 ( @xmath108 kpc ) are at the limit of css data . discoveries in css data covering this region will be likely biased to the brighter , nearer sources .    to further explore the relationship with the sgr streams , we plot the rrab s in the sgr x - y plane system of majewski et al .
( 2003 ) in figure [ sgrxy ] . here
we mark the proposed path of the gemini stream .
we also plot the law & majewski ( 2010 ) n - body model of the sgr stream system .
the model is reasonable match to the structure of the inner rrab data .
the agreement is increased if the source distances are reduced by @xmath109 . however , as found by ruhland et al .
( 2011 ) , the gemini stream is not explained by this model or other models of the sgr stream system .
nevertheless , it appears that if the path of sgr leading debris was extended , it might better match the location of the sources .
indeed , law & majewski ( 2010 ) found that if the sgr dwarf had been orbiting 2 gyr longer , a much longer leading arm would be present .
nevertheless , their simulations did not produced sources corresponding to the distances and locations of the bhb candidates and rrl .
it is apparent that the sgr trailing arm model is a very poor match to the gemini stream , both in extent and location .      to determine the significance of the gemini stream it is necessary to compare the source density with that expected from halo models .
we adopt the halo density model from sesar et al .
( 2010 ) for our comparison .
the local density of rrl is determined using the nth nearest neighbour method , we find the 8th nearest star to the point where we want to calculate the density and calculate the area that bounds these stars .
we use 8 stars motivated by ivezic et al .
( 2005 ) where they found that this number gives the best results for the underlying density on their improved bayesian method .
we also saw that this number was a good combination of precision and computational efficiency .
we compute the observed number density within our grid ( @xmath110 ) , and density from the sesar et al .
( 2010 ) halo model ( @xmath111 ) . to visualize the over - densities on the sgr plane
, we produce a @xmath112 grid on the mawejeski et al .
( 2003 ) sgr coordinates at @xmath113 ( for the x - y plot ) and @xmath114 ( for the @xmath115-@xmath116 plot ) . in figure
[ sgrgab ] , we plot the resulting density ratios in the x - y plane of the sgr system and in figure [ sgrgabx ] , we plot the densities in the @xmath115-@xmath117 system .
the densities here can be readily contrasted sources in figure [ sgrxy ] .
the main density features are the leading and trailing sgr arms ( marked by a and b , respectively ) . the gemini stream ( marked as c ) , largely has a higher ratio density ratio than the trailing stream system ( marked by b ) .
another possible feature is seen near @xmath118 , @xmath119 kpc ( marked @xmath120 ) .
this has much lower significance than the other features , yet is lactated where models predict the sgr trailing arm should cross the galactic plane and meet with the trailing arm @xmath121 .
since the gemini feature , like much of the sgr stream , is offset from the sgr plane , we have collapsed the source z positions to visualize the stream .    unlike the gemini stream , the leading and trail sgr arms a and b are all ready well known ( law and majewski 2010 , and references therein ) . to quantify the significance of the gemini stream , we estimated the number of stars expected by the model at the position of the stream . for simplicity , we selected the volume that bounds the stream in the sgr coordinate system .
that is @xmath122 , @xmath123 and @xmath124 .
based on the sesar et al .
( 2010 ) model density , within this volume we would expect to find 373 rrl assuming 100% detection efficiency . however , only 106 rrl were found .
assuming pure poisson uncertainties , this suggests the gemini area has @xmath125 under - density . as our detection efficiency is close to @xmath126% based on dr13 this short fall
can not be explained pure by missing rrl .
the detection of no distant rrl across most of the mls survey area strongly suggests that the halo is much less dense than expected by the model .
this result is in agreement with the results of watkins et al .
( 2009 ) and sesar et al .
( 2010 ) based on halo rrl discovered in sdss stripe-82 .
this confirms that the halo density declines more rapidly than suggested by the juric et al .
( 2008 ) halo model .    as a second model comparison ,
we compare the rrl densities with the watkins et al .
( 2009 ) halo model .
their results suggest a break in density occurs @xmath127 kpc . matching the model to the data we find that the watkins et al .
( 2009 ) model requires normalization by a factor of 11 to match our data .
this is in good agreement with the factor of 10 found by sesar et al .
( 2010 ) . in figure [ sgrgabwat ]
, we plot the observed radial density compared with the sesar et al .
( 2010 ) model and the normalized watkins et al .
( 2009 ) model .
the observations are found to be in good agreement with the watkins et al .
( 2009 ) model , although the break in density appears to occur nearer to 50 kpc when averaged over a large area .
however , this is because the densities are enhanced by rrl in the sgr leading and trailing arms .
based on the watkins et al .
( 2009 ) model we expect to find 50 rrl in the region selected above .
the number of gemini stream rrl we find is thus @xmath128 larger than expected from this model .      to investigate differences between the gemini rrab sources and the overall population
, one can infer the oosterhoff type based on the rrab period - amplitude relationship ( smith et al .
2011 , and references therein ) . however , since average colours were used to transform the mls light curves , we need to account for the effect of colour variation of the rrab light curves . in dr13
we found that rrab amplitudes were systematically reduced by 0.15 due to pulsational colour variations within the broad bandpass of css images compared to v - band .
we also found that the css oosterhoff type - i ( ooi ) rrab s exhibit a well - defined amplitude limit that is 0.1 mags higher than the zorotovic et al .
( 2010 ) period - amplitude relationship .
to separate ooi rrab candidates from oosterhoff type - ii s ( ooii s ) , we correct the mls amplitudes to values by adding 0.15 mags .
next we define a set of rrab s that are between 0.1 and 0.25 mags above than the average ooi period - amplitude relationship .
these objects are a mixture of oosterhoff types ( although some may be oosterhoff - intermediate sources ) that we remove from consideration .
we then investigated the variation in mls rrab s oosterhoff types with distance . in figure [ mlsotype ]
we plot the period - amplitude distribution of the mls rrab s .
by dividing these sources into a nearer sample consisting of rrab with @xmath129 and the more distant rrab @xmath130 we see a division in types . in figure [ mlsotype ]
, we also plot the period - amplitude distribution for the gemini rrab s .
we also over - plot the zorotovic et al .
( 2010 ) ooi line and an ooii line that is offset by @xmath131 in @xmath132 to match the css observations .
the nearer rrab exhibit a mixture of oosterhoff types when compared to the zorotovic et al .
( 2010 ) period - amplitude relationships .
this is significant evidence that the rrab s in the range @xmath129 do not come from a single population .
in contrast , the distant mls rrab s all lie near the ooi line , suggesting that they come from a single population that is more metal - rich .    in figure [ mlshbd ]
, we show the locations of gemini stream rrab s after removing the mls rrab s with ambiguous oosterhoff type .
the css rrab s with @xmath129 have also been included , but are not separated by oosterhoff type since their amplitudes and periods are less certain than the mls rrab s .
this figure shows that the gemini rrab s are spread across the width of sgr stream system .
it also shows that the ooii rrab s selected among the nearer rrl set ( @xmath129 ) are distributed across the region .
the nearer ooi group ends around @xmath133 where the stream divides between the nearer and further groups .
analysis of this figure suggests that there may be two or more overlapping populations : one that follows a steep distant gradient and is predominantly ooi , and another that has a shallower gradient and is an ooii population . with respect to this figure
, we once again note that the mls rrl are naturally concentrated toward the ecliptic because of sampling .
so it is not possible to make inferences about changes in density across the sgr stream .
while the sdss photometry has little of the repeated photometry required to unambiguously identify rrab s , it is deeper than both css and mls data ( reaching hb stars to @xmath134 ) .
additionally , sdss data does cover most of the gemini stream observed in the mls data as well as the main region where mls data does not overlap the sgr stream .
therefore , it is possible to use this data to bridge the gap in between the depth of css data and coverage of mls photometry .    unlike previous authors whom searched for bhb stars covering the sgr tidal streams , here we sought to select both bhb and potential rr lyrae stars .
based on the colours of the distant mls rrab s and prior bhb work , we investigated the observed colours of rrl in sdss photometry and selected sdss dr8 stars within @xmath135 , @xmath136 , @xmath137 , sdss object type=6 ( star ) , and @xmath138 .
based on dr13 we know that the bulk of rrab s are located near @xmath139 , @xmath140 and @xmath141 .
however , rrab s with these colours are outnumbered by main sequence turnoff stars ( msto ) by a factor of @xmath142 ( koposov et al .
2012 ) . using maps of sdss source density , we sub - selected stars in the range @xmath143 for @xmath144 .
this selection retains most of the bhb candidates of ruhland et al.(2011 ) as well as a 22% of the distant mls rrab sample and many bs stars .
this colour selection reduces the initial number of sdss point sources in our selection from 1.65 million to 81552 . by selecting stars in the sgr stream region , with @xmath145 ,
the number reduces to 23507 hb candidates . to match the mls photometry we transform the sdss photometry to @xmath25 magnitudes using ivezic et al .
( 2007 ) . as found with newberg et al .
( 2003 , 2007 ) and ruhland et al .
( 2011 ) we found the spatial distribution of sdss hb candidates suffers from a significant crowding and background that is best viewed in source density . in figure [ hess ]
, we provide the hb candidates in the form of the hess ( point - density ) diagram .
this figure shows evidence that the density of hb candidates closely follows the fit to the distant gemini stream of rrab s , in agreement with the distribution of bhb candidates from ruhland et al .
( 2011 ) and others . here
the scaling has been set to match that in figure [ polar ] where the main streams are visible .
in addition to the rrl , the paths expected for bs stars that pass the colour cuts and form a shadow that is @xmath146 mag fainter than the hb stars is also shown .    to obtain another view of mls rrl in the gemini stream in relation to the sdss hb candidates
, we divided the sdss sources into three groups .
these were bright sources , with @xmath147 , intermediate - brightness sources , with @xmath148 , and faint sources , with @xmath149 . in figure
[ mlshb ] , we present the locations of these sources . the bright sources were selected to show the nearby hb stars as well as the bs stars in the monoceros stream ( near the galactic anti - center at @xmath150 limit of the sdss coverage ) .
the intermediate - brightness sources were selected to be indicative of the hb stars in the leading arm of the sgr stream .
the faintest sources were selected to include the hb candidates in the gemini stream . however , this also include bs stars that mirror the distribution of hb stars along the sgr leading arm .
the figure clearly shows the overdensity associated with the sgr streams system .
the separation into two streams is not clear here since there are far fewer bhb stars than the msto stars used by belokurov et al .
( 2006 ) .
we matched the entire mls rrl catalog with the sdss dr8 spectroscopic catalog and found 89 matches
. a much larger sample of rrl with spectra is given in dr13 .
however , here our main purpose was analysis of the distant halo rrl .
we found that 16 of these matches are sources beyond 70 kpc .
as noted earlier , many of the faint rrl candidates were found based on their sdss spectra .
the 16 rrl with sdss spectra include 12 rrab s and four rrc s .
to separate the radial velocities from the velocities due to pulsation , we use the sdss observation times and the fourier fits to derive the phase at which the rrl spectrum was observed .
as noted in dr13 , for sdss spectra , radial velocities are determined by averaging both balmer and metallic lines ( mainly ca lines ) .
we follow dr13 by applying sesar et al .
( 2012b ) velocities corrections for pulsation based corrections derived from both hydrogen and metallic lines .
the average velocity correction for the 12 rrab is @xmath151 , in good agreement with the uncertainty derived from 905 rrab spectra in dr13 ( @xmath152 km / s ) .
we adopt this level of uncertainty for all the sdss rrab spectra .    for the four rrc stars which pulsate in the first overtone mode
the velocity of the pulsation is much smaller than for the rrab s .
based on the rrc s observed by liu & janes ( 1989 ) and jones et al .
( 1988 ) , the amplitude is expected to be approximately @xmath153 . to account for this factor , for these rrl
we increase the observed sdss radial velocity uncertainties by an additional @xmath154 . following law & majewski ( 2010 )
we transform the radial velocities to the galactic standard of rest assuming a solar peculiar motion of ( u , v , w ) = ( 9 , 12 + 220 , 7 ) km / s in the galactic cartesian coordinate system .    in figure
[ mlsvel ] , we plot the radial velocities for the outer halo mls rrab s with spectra , the css rrab s with @xmath155 kpc , and the 10 ptf rrab s found by sesar et al .
( 2012a ) .
we also plot the velocities predicted by law & majewski ( 2010 ) n - body simulations .
the figure shows that two of the rrab s observed near @xmath156 appear to be associated with the leading arm of the sgr stream .
however , for the other gemini stream rrab s , the velocities , like the distances , are not explained by the law & majewski ( 2010 ) n - body simulations . in this
regard , sesar et al .
( 2012a ) noted that the rrab velocity measurements they found suggested the rrab s belonged to two distinct groups .
additionally , sources associated with our group - d of figure [ sgrxy ] ( occurring near @xmath157 ) have radial velocities that appear to vary rapidly with right ascension .
these velocities are also inconsistent with the law & majewski ( 2010 ) model , suggesting that they do not belong to the sgr leading arm .
the association between these feature - d rrab s and the gemini stream rrab s is unclear , although values are similar .      as noted by newberg et al .
( 2003 ) , the outer halo bhb candidates they discovered reside near the unusual globular cluster ( gc ) ngc  2419 .
this system is located at @xmath158 , @xmath159 and distance @xmath160 kpc ( harris 1996 , 2010 edition ) .
the corresponding coordinates in the majewski et al .  ( 2003 ) sgr coordinate system are @xmath161 , @xmath162 , making it well within the limits of the sgr stream system ( koposov et al .
the galactocentric radial velocity is given by newberg et al .
( 2003 ) as @xmath163 , and baumgardt et al .
( 2009 ) find an internal velocity dispersion of @xmath164 .
ngc  2419 is noted as being one of the most metal - poor gcs ( @xmath45 \sim -2.1 $ ] ; mucciarelli et al .
the cluster is notably old , with age 12.3 gyr according to forbes & bridges ( 2010 ) .
the system has the highest luminosity ( @xmath165 mag ) of any gc with a galactocentric distance @xmath166 kpc , apart from the likely sgr dwarf - associated gc , m54 ( cohen et al .
the half - light radius of this cluster is 19 pc making it significantly more extended in the @xmath167 versus @xmath168 plane than other gcs with @xmath166 kpc ( mackey & van den berg 2005 ) .
indeed , the exceptional nature of ngc  2419 relative to outer halo gcs led van den berg & mackey  ( 2004 ) to suggest that the object is the stripped core of a former dwarf spheroidal galaxy ( dsph ) .
based on abundance studies cohen et al .
( 2010 ) and cohen & kirby ( 2012 ) also found that the ngc  2419 appears like no other globular cluster , but rather the core of an accreted dwarf galaxy .    from figures [ distg ] & [ mlsvel ]
we see that the location and radial velocity of ngc  2419 are a relatively good match for the gemini tidal stream .
the average velocities and metallicities of the sesar et al .
( 2012a ) cancer group b ( @xmath169 ) and @xmath45 = -2.1 \pm 0.4~dex$ ] ) and the metallicity of the cancer group b is in reasonable agreement with the values expected for a stream from ngc  2419 .
however , the velocities and metallicities do not provide a strong enough association to link these sources .
another important point to consider in the possible association between the gemini rrl and ngc  2419 is the proximity of the sources .
some of the rrl have @xmath170 near the right ascension of the ngc  2419 ( @xmath158 ) .
thus , if one was to assume that the gemini tidal stream proceeds in the direction of the sgr stream ( as suggested by sdss hb candidates ) , the stream stars would have to be dispersed across the entire @xmath171 of the sgr system between the mls rrab s and ngc  2419 .
sources within this gap are not covered by mls observations because of the coverage limits of the survey .
it is thus remain unclear whether there is any link between the gemini stream and ngc  2419 .
also , as noted earlier , our data shows that the most distant gemini rrab s are ooi type stars while ngc  2419 is well known to be ooii type .
indeed , mucciarelli et al .  ( 2012 ) found that ngc  2419 exhibits very little spread in @xmath45 $ ] ( @xmath172 dex ) suggesting it could not be linked with metal - rich rrl in the gemini structure .
we have performed a periodicity analysis of 3.1 million variable star candidates selected from photometry taken by the mls survey and uncovered 1,207 rrl ( of which 538 are new ) .
comparison of the periods for the @xmath173 previous known rrl shows that the sources are accurately measured .
we have also discovered 2040 new rrab stars in a re - analysis of css photometry .    using catalina surveys photometry we have determined the nature of 90% of the sekbo ( keller et al .
2006 ) rrl candidates and find that 60% are likely to be rrl .
our analysis of the sekbo rrl candidates revealed , selection of a pure set of rrab s with accurate average magnitudes requires many observations .
the importance of repeated observations for characterizing variable star types was also recently , demonstrated by sesar et al .
( 2010 ) , who used rrab light curves to study over - densities that they had earlier been attributed to rrl in sesar et al .
they discovered that a number of the over densities attributed to rrl , based on photometric selection coupled with a small number of observations , were in fact due to intrinsically fainter @xmath65 scuti stars as well as non - variable sources .
overall they found that only 70% of their initial candidates were rrl stars .
similarly , ivezic et al .
( 2005 ) found that although it was possible to completely colour - select a small fraction ( 6% ) of rrl based on sdss photometry , if one wanted to select 60% , a 72% non - rrl contamination rate would result .
similar levels of contamination make it equally difficult to trace halo structures using photometrically selected bhb stars .    in our analysis
we also found a significant group of rrl with average @xmath25 magnitudes @xmath174 . by combining these sources with css rrab s , and photometrically - selected hb candidates ,
we find strong evidence for a tidally disrupted stellar stream that crossing @xmath175 of the sky at galactocentric distances from 70 to 110 kpc .
this result confirms the existence of a stream first noted by newberg et al .
comparison with halo density models shows that the feature is significant and that the halo density declines rapidly beyond 30 - 50 kpc , as previously noted by watkins et al .
( 2009 ) and sesar et al .
however , since these results , like those of watkins et al .
( 2009 ) and sesar et al .
( 2010 ) , are based on a thin slice through the halo , caution has to be taken when interpreting the extent of density regions outside the observed fields .    although the gemini rrl overlap with the sgr stream system , we find that the large distances are inconsistent with existing sgr models .
this result is in agreement with ruhland et al .
( 2011 ) and a recent sample of rrl discovered by sesar et al .
( 2012a ) .
furthermore , we find that the radial velocities of the rrl are inconsistent with simulations of sgr stream . however , we note that models of the sgr streams system as a whole remains poorly constrained by observations . we have investigated the possible relationship between the gemini tidal stream and ngc  2419 as first proposed by newberg et al .
the most recent analyses of ngc  2419 shows significant evidence for it being the nuclear remnant of a disrupted dwarf galaxy ( van den berg & mackey 2004 ; mackey & van den berg 2005 ; cohen et al .
2010 , 2011 ; forbes & bridges 2010 ; cohen & kirby 2012 ) .
although , we find that the distances of many of the rrls and hb candidates are consistent with ngc  2419 , the available velocities and location of the gemini stream are insufficient agreement to link the two structures . furthermore , the most distant of the rrab s discovered appear to be metal rich sources and would thus be inconsistent with stars observed in ngc  2419 . nevertheless , given the location ngc  2419 within the halo , and with 10 degrees of the sgr streams plane , it seems possible that a stream associated with ngc  2419 could join it to the sgr system .
this may in part account for the significant diversity in metallicity observed for varying sgr stellar streams ( law & majewski  2010 ) .
furthermore , as the photometric - selected hb candidates and rrab s cover the sgr stream system and exhibit a distance gradient , it is possible that there is a second galaxy remnant associated with the sgr stream .
such a source might explain the origin of the two intersecting streams of the sgr system that has now been well delineated with msto stars by belokurov et al .
( 2006 ) and kosopov et al .
however , even with 40% uncertainties in the distances to msto stars ( newby et al .
2011 ) , the gemini stream rrab s are twice as distant as expected for msto stream stars ( koposov et al . 2012 ) . alternately
, the gemini stream may originate from the remnant of another disrupted dwarf galaxy that lies beyond the gemini stream stars and is yet to be discovered .
the gemini stream leads into the galactic plane beyond 100kpc .
a highly extincted system in the galactic plane would be very difficult to detect .
future photometric and spectroscopic observations of the hb stars within @xmath58 of ngc  2419 could confirm whether there truly is a tidal stream of rrl associated with ngc  2419 .
for example , the seemingly unique mg and k abundance patterns in ngc  2419 found by mucciarelli et al .
( 2012 ) would chemically tag stars originating from this source , even in the presence of overlapping tidal streams .
if the gemini tidal stream does follow the path expected from the sdss hb candidates , additional deep photometric observations undertaken by projects such as lsst ( abell et al .
2009 ) should reveal numerous additional rrl along this tidal stream .
moreover , if these stars are associated with ngc  2419 , we predict that the rrl will mostly be type oosterhoff ii .
crts and csdr1 are supported by the u.s .
national science foundation under grants ast-0909182 and cns-0540369 .
the css survey is funded by the national aeronautics and space administration under grant no .
nng05gf22 g issued through the science mission directorate near - earth objects observations program .
j. l. p. acknowledges support from nasa through hubble fellowship grant hf-51261.01-a awarded by the stsci , which is operated by aura , inc . for nasa , under contract
nas 5 - 26555 .
support for m.c . and
is provided by the ministry for the economy , development , and tourism s programa inicativa cientfica milenio through grant p07 - 021-f , awarded to the milky way millennium nucleus ; by proyecto basal pfb-06/2007 ; by proyecto fondecyt regular # 1110326 ; and by proyecto anillo act-86 .
sdss - iii is managed by the astrophysical research consortium for the participating institutions of the sdss - iii collaboration funding for sdss - iii has been provided by the alfred p. sloan foundation , the participating institutions , the national science foundation , and the u.s .
department of energy office of science .
the sdss - iii web site is http://www.sdss3.org/.


it has been recognized shortly after the discovery of pulsars that neutron stars have velocities much in excess of those of any other normal stellar populations in our galaxy ( gunn & ostriker 1970 ; minkowski 1970 ) .
however , it is only in the last few years that a significant body of evidence has come into place to support the view that type ii supernovae are nonspherical and neutron stars receive large kick velocities at birth .
lyne and lorimer ( 1994 ) analyzed pulsar velocities in light of new proper - motion measurements ( harrison et al .
1993 ) and up - to - date pulsar distance scale ( taylor & cordes 1993 ) , and concluded that pulsars were born with a mean speed of @xmath9  km  s@xmath7 , much larger than previously thought .
more recent studies of pulsar velocities have adopted more sophisticated satistical methods and included better treatment of selection effects and uncertainties ( lorimer et al .
1997 ; hansen & phinney 1997 ; cordes & chernoff 1997 ) .
they all yielded a result for the pulsar birth velocity in qualitative agreement ( although not in quantitative details ) with that of lyne & lorimer ( 1994 ) . in particular , cordes and chernoff ( 1997 ) found that the three - dimensional space velocities of their sample of 47 pulsars have a bimodal distribution , with characteristic speeds of @xmath10  km  s@xmath7 and @xmath11  km  s@xmath7 ( corresponding to @xmath12 and @xmath13 of the population ) .
they also estimated that pulsars with velocities greater than @xmath14  km  s@xmath7 may be underrepresented owing to selection effects in pulsar surveys [ the uncertainty in the high - velocity end of the pulsar velocity distribution function has also be emphasized by hansen & phinney ( 1997 ) ] .
concrete evidence for the existence of pulsars with velocities of @xmath15  km  s@xmath7 has come from the observation of the guitar nebula pulsar ( b2224 + 65 ) , which produces a bow shock when plowing through the interstellar medium ( cordes , romani & lundgren 1993 ) .
in addition , studies of pulsar - supernova remnant associations have uncovered a number of pulsars having velocities much greater than @xmath14  km  s@xmath7 ( frail et al .  1994 ) , although in many cases the associations are not completely secure .
compelling evidence for supernova asymmetry and pulsar kicks also comes from the detection of geodetic precession in the binary pulsar psr 1913 + 16 ( cordes et al .  1990 ; arzoumanian et al .  1996 ) , and the orbital plane precession in the psr j0045 - 7319/b star binary ( lai et al .
1995 ; kaspi et al .
1996 ) and its fast orbital decay ( which indicates retrograde rotation of the b star with respect to the orbit ; see lai 1996 ) .
these results demonstrate that binary break - up ( as originally suggested by gott , gunn & ostriker 1970 ; see iben & tutukov 1996 ) can not be _ solely _
responsible for the observed pulsar velocities , and that _
natal kicks are required_. in addition , evolutionary studies of neutron star binary population imply the existence of pulsar kicks ( e.g. , deway & cordes 1987 ; fryer & kalogera 1997 ; fryer et al .  1998 ;
see also brandt & podsiadlowski 1995 ) . finally , there are a large number of direct observations of nearby supernovae ( e.g. , cropper et al .  1988 ;
trammell et al .  1993 ; mccray 1993 ; utrobin et al .  1995 ) and supernova remnants ( e.g. , morse , winkler & kirshner 1995 ; aschenbach et al .  1995 ) in radio , optical , and x - ray bands which support the notion that supernova explosions are not spherically symmetric .
the origin of the pulsar velocities is still unknown .
two classes of mechanisms for the _ natal kicks _ have been suggested .
the first class relies on convective instabilities in the collapsed stellar core and within the rebounding shock ( e.g. , burrows & fryxell 1992 ; burrows et al .  1995 ; janka & mller 1994 ,  1996 ; herant et al .
the asymmetries in the matter and temperature distributions associated with the instabilities naturally lead to asymmetric matter ejection and/or asymmetric neutrino emission .
numerical simulations indicate that the local , post - collapse instabilities are not adequate to account for kick velocities higher than a few hundred km  s@xmath7 ( janka & mller 1994 ) .
a variant of this class of models therefore relies on the global asymmetric perturbations seeded in the presupernova cores ( goldreich et al .
1996 ; see also bazan & arnett 1994 ) .
clearly , the magnitude of kick velocity depends on the degree of initial asymmetry in the imploding core ( burrows & hayes 1996 ) . due to various uncertainties in the presupernova stellar models ,
it is not clear at this point whether sufficiently large initial perturbations can be produced in the precollapse core ( lai & goldreich 1998 ) .    in this paper
we focus on the second class of models in which the large pulsar kick velocities arise from asymmetric neutrino emission induced by strong magnetic fields .
since @xmath16 of the neutron star binding energy ( a few times @xmath17  erg ) is released in neutrinos , tapping the neutrino energy would appear to be an efficient means to kick the newly - formed neutron star .
magnetic fields are naturally invoked to break the spherical symmetry in neutrino emission .
but the actual mechanism is unclear .
we first briefly comment on previous works on this and related subjects .
a number of authors have noted that parity violation in weak interactions may lead to asymmetric neutrino emission from proto - neutron stars ( chugai 1984 ; dorofeev et al .  1985 ; vilenkin 1995 ; horowitz & piekarewicz 1997 ) .
however , their studies are largely unsatisfactory for a number of reasons : they either failed to identify the most relevant neutrino emission / interaction processes or the relevant physical conditions in proto - neutron stars , or stopped at estimating the magnetic field effects on neutrino opacities .
chugai ( 1984 ) and vilenkin ( 1995 ) ( see also bezchastnov & haensel 1996 ) considered neutrino - electron scattering and concluded that the effect is extremely small ( e.g. , to obtain @xmath18  km  s@xmath7 would require a magnetic field of at least @xmath8  g ) .
however , neutrino - electron scattering is less important than neutrino - nucleon scattering in determining the characteristics of neutrino transport in proto - neutron stars .
similarly , dorofeev et al .
( 1985 ) considered neutrino emission by urca processes in strong magnetic fields .
but as shown by lai & qian ( 1998 ) , in the bulk interior of the neutron star , the asymmetry associated with neutrino emission is cancelled by the asymmetry associated with neutrino absorption . concerning the parity violation effect in proto - neutron stars ,
horowitz & li ( 1997 ) recently suggested that the asymmetry in neutrino emission may be enhanced due to multiple scatterings of neutrinos by nucleons which are slightly polarized by the magnetic field .
they estimated that a field strength of a few times @xmath19  g is adequate to account for kick velocities of a few hundred km  s@xmath7 .
however , the paper by horowitz & li ( 1997 ) only discusses the idealized situations of scattering media , and ignores various essential physics that is needed in a proper treatment of neutrino transport in proto - neutron stars .
in particular , it does not consider the effect of neutrino absorption which one might suspect to wash out the cumulative effect from multiple scatterings .
our own study of the parity violation effect ( lai & qian 1998 ) was flawed in the treatment of scattering terms in the neutrino transport equation . as shown by arras & lai ( 1998 )
, detailed balance requires that there be no cumulative effect from multiple scatterings in the bulk interior of the proto - neutron star where local thermodynamic equilibrium applies to a good approximation .
enhancement of neutrino emission asymmetry from multiple scatterings obtains only after neutrinos thermally decouple from proto - neutron star matter , and therefore is insignificant .
bisnovatyi - kogan ( 1993 ) attributed pulsar kicks to asymmetric magnetic field distribution in proto - neutron stars .
using the fact that neutron decay rate can be modified by the magnetic field , he inferred that neutrino emissions from opposite sides of the neutron star surface are different .
however , neutron decay is not directly relevant for neutrino emission from a newly - formed neutron star .
roulet ( 1997 ) ( whose paper was posted while this paper was being written ) considered the relevant neutrino absorption processes .
but his study was restricted to calculating the neutrino cross sections as functions of neutrino energy ( corresponding to  2.1 of this paper ) , and therefore was insufficient for addressing the issue of asymmetric neutrino emission .    finally , it is worthwhile to mention a speculative idea on pulsar kicks which relies on nonstandard neutrino physics .
kusenko & segr ( 1996 ) suggested that asymmetric @xmath20 emission could result from the mikheyev - smirnov - wolfenstein flavor transformation between @xmath20 and @xmath0 inside a magnetized proto - neutron star because a magnetic field changes the resonance condition for the flavor transformation .
with a @xmath20 mass of @xmath21  ev , they claimed that magnetic fields of @xmath22  g can give the pulsar a kick velocity of a few hundred km  s@xmath7 .
however , their treatment of neutrino transport was oversimplified : e.g. , they ignored that neutrinos of different energies have different resonance surfaces . furthermore , a simple geometric effect and realistic proto - neutron star conditions easily reduce their estimated kick velocity by an order of magnitude .
most likely , a magnetic field strength of @xmath23  g is needed to produce the observed average pulsar kick velocity via their mechanism ( see qian 1997 ) .
similarly strong magnetic fields are required in variants of the kusenko & segr mechanism , e.g. , those considered by akhmedov et al .
( 1997 ) ( which relies on both the neutrino mass and the neutrino magnetic moment to facilitate the flavor transformation ) .
as should be clear from the brief review in  1.1 , in spite of quite a number of suggestions , there is at present no consensus on the magnitudes of the magnetic field induced asymmetry in neutrino emission from proto - neutron stars and the resulting kick velocities . in this paper
we study the effect of asymmetric magnetic field topology on pulsar kicks .
since the energy levels of @xmath2 and @xmath3 in a magnetic field are quantized , the @xmath0 and @xmath5 absorption opacities near the neutrinosphere and the neutrino - matter decoupling region depend on the local magnetic field strength .
if the _ magnitude _ of magnetic field in the north pole is different from that in the south pole ( the field does not need to be ordered ) , then asymmetric neutrino flux may be generated .
here we ignore the effects of magnetic fields on the proto - neutron star structure through the equation of state  these are secondary effects as far as neutrino transport is concerned .
the organization of this paper is as follows . in  2 we calculate the @xmath0 and @xmath1 absorption cross sections as functions of neutrino energy in strong magnetic fields .
section 3 summarizes the basic features of neutrino transport near the neutron star surface . in
4 we derive and evaluate the `` rosseland mean '' neutrino opacities which are directly related to the local neutrino flux .
we then consider in  5 the change in neutrino flux due to the modified opacities , and estimate the kick velocity resulting from an asymmetric magnetic field topology . in  6 we briefly consider how the emergent neutrino spectra may be modified by the strong magnetic field . in  7 we present our conclusion .
unless noted otherwise , we shall use units in which @xmath24 , @xmath25 , and the boltzmann constant @xmath26 are unity .
in this section we consider how magnetic fields modify the ( total ) neutrino absorption ( @xmath27 and @xmath28 ) cross sections for specific neutrino energies .
we present our calculation of the @xmath0 absorption cross section in detail .
the @xmath1 absorption cross section can be obtained in a similar manner .    to begin with , it is useful to review various energy scales involved in the problem . near the neutrinosphere and the neutrino - matter decoupling region , the matter density is typically @xmath29@xmath19 g@xmath30@xmath31 , and the temperature is about several mev .
the electrons are extremely relativistic and highly degenerate , with fermi energy @xmath32mev , where @xmath33 is the electron fraction ( the number of electrons per nucleon ) and @xmath34 .
the neutrons and protons have fermi energies @xmath35^{2/3}~$]mev and @xmath36mev , respectively , and thus they are essentially nondegenerate .      in strong magnetic fields the transverse motion of the electron is quantized into landau levels .
the electron energy dispersion relation is e_e=^1/2 ,  n=0 ,  1 ,  2 ,  where @xmath37 is the longitudinal momentum ( along the magnetic field ) , and @xmath38 is the landau level index . the dimensionless field strength @xmath39 is defined as b ,  with  b_c = m_e^2c^3e = 4.41410 ^ 13  g. the quantization effect on the proton is extremely small due to the large proton mass , and is neglected throughout this paper .
since the cyclotron radius ( characterizing the size of the landau wavefunction ) @xmath40^{1/2}= [ \hbar/(m_e c)]b^{-1/2}=386\,b^{-1/2}$ ] fm is much greater than the range of the weak interaction [ @xmath41 fm , where @xmath42 gev is the mass of w boson ] , one can use the @xmath43 theory with contact interaction .
locally ( at the point of interaction ) the landau wavefunction resembles a plane wave .
thus we expect that the spin - averaged matrix element is the same as in the case without magnetic fields .
this expectation has been borne out by detailed calculations ( e.g. , fassio - canuto 1969 ; matese & oconnell 1969 ) .
therefore we only need to consider the effect of magnetic field on the electron phase space .
since the neutrino energy ( a few to tens of mev ) of interest is much less than the nucleon mass , we first consider the case where the nucleon mass is taken to be infinity , an approximation adopted in the standard zero - field calculations ( e.g. , tubbs & schramm 1975 ) . the cross section for @xmath44 can be written as a sum over the ( spatial and spin ) states of the electron ( per unit volume ) ^_b(e_)=a_e(1-f_e)(e_+q - e_e ) , where we have defined ag_f^2 ^2_c ( c_v^2 + 3c_a^2 ) , and the other symbols have their usual meanings : @xmath45 is the universal fermi constant , @xmath46 refers to the cabbibo angle , @xmath47 and @xmath48 are the coupling constants , @xmath49 mev is the neutron - proton mass difference , @xmath50 is the @xmath0 energy , and @xmath51 is the pauli blocking factor .
the electron occupation number @xmath52 is f_e=1 + 1 , where @xmath53 is the electron energy , @xmath54 is the electron chemical potential , and @xmath55 is the matter temperature . for zero field ,
the sum becomes @xmath56 , and we have ^_b=0(e_)=a^2(e_+q ) ^1/2(1-f_e ) , with @xmath57 . for @xmath58 , knowing that the degeneracy of the landau level ( per unit area ) is @xmath59 , with @xmath60 for @xmath61 and @xmath62 for @xmath63 , we have ^_b(e_)=a _ -^_n g_nm_e^2b2(1-f_e ) ( e_+q - e_e)=am_e^2b2 ^ 2 _ n=0^n_max g_n e_ep_z , n(1-f_e ) , [ absb]where @xmath57 , and @xmath64^{1/2}$ ] .
the upper limit @xmath65 for the sum is the maximum value of @xmath38 for which @xmath66 is meaningfully defined , i.e. , n_max = int , [ nmax]where int@xmath67 $ ] stands for the integral part of @xmath68 .
henceforth , @xmath65 is similarly defined , and we shall not explicitly write out the lower and upper limits for the sum .
figure 1 depicts @xmath69 for @xmath70 and @xmath71 as functions of @xmath50 for typical conditions ( @xmath72 , @xmath73 mev , and @xmath74 mev ) in a proto - neutron star .
note that ideally , one should compare @xmath69 and @xmath71 at the same matter density . for @xmath75 ( so that @xmath76 ) ,
the matter density is related to @xmath54 via @xmath77 in practice , we have found that for a given @xmath54 , the dependence of @xmath78 on the field strength is rather weak for the range of parameters of interest in this paper .
thus in all of our calculations here and below , we compare the cross sections at a given @xmath54 . a prominent feature of @xmath69 shown in fig .
1 is the presence of `` resonances , '' which occur when @xmath79 ( cf .
eq .  [ [ absb ] ] ) , i.e. , at the neutrino energies for which the final state electron lies in a `` stationary '' landau level ( with zero longitudinal momentum ) .
these resonances reflect the divergence of the density of states of electrons at these stationary landau levels .
the cross section for @xmath80 can be similarly calculated .
one only needs to replace @xmath81 with @xmath82 and @xmath54 with @xmath83 in the expressions for @xmath0 absorption to obtain the desired expressions for @xmath1 absorption .
the numerical results for @xmath1 absorption cross section are shown in fig .  2 for the same physical parameters as used in fig .  1 .
the divergent cross sections at the `` resonances '' shown in figs .  1 and 2 are clearly unphysical .
they are smoothed out by the thermal motion of the neutron . to the leading order , the width of the resonance is @xmath84 , corresponding to the doppler shift of the neutrino energy in the neutron s rest frame .
the recoil energy , of order @xmath85 , is much less effective in smoothing out the resonances , and will be neglected . with this approximation ,
the smoothed cross section can be written as ^_b(e_)v_n^2 dv_nf_nd_n _
b^(e_ ) , [ aver]where @xmath86 and @xmath87 is the velocity distribution function for the neutron . if we replace the factor @xmath51 [ which depends on @xmath88 in @xmath89 with the average value @xmath90 , which can be approximately taken as @xmath51 evaluated at @xmath57 , then the angular integral in eq .
( [ aver ] ) can be carried out analytically . for @xmath91 ,
we have ^_b=0(e_)1-f_ev_n dv_nf_n13 ^ 2e _ ( e_e^2-m_e^2)^3/2|^e_e = e_(1+v_n)+q _ e_e = e_(1-v_n)+q , [ abs0aver]which is reduced to ^_b=0(e_)(e_+q)^21-f_e[abs0aver0]if we neglect the electron mass @xmath92 . for @xmath58 , we find has an integrable singularity at @xmath93 . ]
^_b(e_)1-f_ev_n dv_nf_nm_e^2b2 ^ 2e _ _
n g_n^1/2|^e_e = e_(1+v_n)+q _ e_e = e_(1-v_n)+q .
[ absaver]we could have approximately included the neutron recoil effect by replacing @xmath81 with @xmath94 in eqs .
( [ abs0aver])([absaver ] ) .
but we choose to leave out this effect since as mentioned previously , it is ineffective in smoothing out the resonances in @xmath69 .
we have numerically evaluated eqs .
( [ abs0aver ] ) and ( [ absaver ] ) , and the results are depicted in fig .  1 together with the cross sections before averaging .
similar results for the @xmath1 absorption cross section are shown in fig .
as expected , we find that for @xmath91 , the averaged cross section @xmath95 is almost identical to @xmath71 ( see eq . [ [ abs0aver0 ] ] ) .
however , for @xmath58 , the thermal averaging is essential in obtaining the smoothed physical absorption cross sections .
the oscillatory behaviors of @xmath96 in figs .  1 and 2 can be demonstrated analytically . using the poisson summation formula
( see e.g. , ziman 1979 ) : _ n=0^f(n+12)= _
s=-^(-1)^s_0^f(x)e^2i s xdx , we can rewrite the sum in eq .
( [ absaver ] ) as @xmath97^{1/2}\nonumber\\ & & = { 2\over 3}(e_e^2-m_e^2)^{3/2 } + ( 2m_e^2b)^{3/2}{\rm re}\sum_{s=1}^\infty { 1-i\over 4\pi i s^{3/2}}\exp\left[2\pi i s \left({e_e^2-m_e^2\over 2m_e^2b } \right)\right ] .
\label{summ}\end{aligned}\ ] ] the first term on the right - hand side of eq .
( [ summ ] ) corresponds to the zero - field result , and the second term ( involving the sum of many oscillatory terms ) corresponds to the correction due to the quantization effect .
the sum is clearly convergent , and thus for small @xmath98 , the averaged cross section @xmath96 is reduced to the zero - field value @xmath95 .    for sufficiently large @xmath98 at a given @xmath50 , or , equivalently , for sufficiently small @xmath50 at a given @xmath98
, only the ground state landau level is filled by the electron , i.e. , @xmath99 ( see eq .
[ [ nmax ] ] ) . in this limit , @xmath100 , and
only the first term in eq .
( [ absb ] ) needs to be retained .
we then have ^_b(e_)am_e^2b2 ^ 2 ( e_+q)^1/2(1-f_e ) [ sigbapp]for @xmath0 absorption .
a similar expression for @xmath1 absorption is obtained by replacing @xmath81 with @xmath82 and @xmath54 with @xmath83 in eq .
( [ sigbapp ] ) .
the physical conditions which figs .  1 and 2 are based on ( @xmath74  mev , @xmath72 , and @xmath73  mev , corresponding to a matter density of @xmath101 g  @xmath31 ) are typical of the neutrinosphere of the proto - neutron star . in fig .
3 we plot the @xmath0 and @xmath1 absorption cross sections for a number of slightly different physical parameters , specified by @xmath102 and @xmath55 ( all with @xmath72 ) .
the cross section for elastic neutral - current scattering @xmath103 , @xmath104 ( e.g. , tubbs & schramm 1975 ) , is also shown for comparison .
absorption dominates the transport of @xmath0 for high energies , with the ratio of absorption to scattering cross section per nucleon @xmath105 approaching 5 for @xmath106 . for smaller energies ( @xmath107 ) ,
absorption is severely suppressed by the pauli blocking of the final electron states , and scattering becomes more important . in the case of @xmath1 ,
pauli blocking is negligible , and scattering almost always dominates the transport , with the ratio of scattering to absorption cross section per nucleon @xmath108\simeq 2\,(0.1/y_e)$ ] for @xmath109  mev .    as we can see from figs .  1 - 3 , for a given set of physical parameters of the medium ( @xmath110 or @xmath78 ) , the magnetic field effect on the absorption opacity is more prominent for lower - energy neutrinos than for higher - energy ones .
this is because for a given field strength , as the neutrino energy increases , more and more landau levels become available for the electron in the final state ( see eq .
[ [ nmax ] ] ) . consequently
, the electron phase space and hence the cross section should approach the zero - field values according to the correspondence principle .
we found that typically when the number of landau levels filled by the electron ( @xmath65 ) is greater than 23 , the magnetic field effect on the absorption opacity becomes negligible .
the `` critical ''
neutrino energy is approximately given by e_^(crit)m_e(2bn_max+1)^1/2 6(n_max3)^1/2 ( b10 ^ 15g)^1/2   mev .
for @xmath111 , one can approximate @xmath96 by the zero - field value @xmath95 with good accuracy ( within a few percent ) .
alternatively , we can define a threshold magnetic field , @xmath112 , above which the @xmath0 and @xmath1 absorption opacities are modified ( by more than a few percent ) by the magnetic field : b_thb_c 310 ^ 15(e_10  mev)^2 ( 3n_max )  g , where @xmath113 mev is the average neutrino energy .
in the deeper layers of the neutron star , the typical neutrino energy is higher .
therefore we expect the magnetic field effect on the absorption opacities to be smaller there . in addition , for @xmath0 , absorption becomes unimportant as compared to scattering in the high density region ( although not necessarily in the deepest interior , where @xmath0 can be degenerate ) .
thus , in order to assess the importance of the magnetic field effect on the neutrino emission , we only need to focus on the surface region , the region between the neutrinosphere and the neutrino - matter decoupling layer .
we note that because the absorption opacities in the magnetic field for different neutrino energies oscillate around the the zero - field values , without a detailed calculation of the mean opacity ( appropriately averaged over the neutrino energy spectrum ) one can not even answer the qualitative question such as whether the magnetic field increases or decreases the local neutrino energy flux .    in
36 , we study how the modified absorption cross sections affect the emergent @xmath0 and @xmath1 energy fluxes and spectra from the proto - neutron star .
as discussed in  2.4 , the quantization effect of magnetic fields on neutrino absorption is important in the region only near the neutron star surface ( between the neutrinosphere and the neutrino - matter decoupling layer ) . in this section
we consider neutrino transport in this relatively low density region .    because @xmath0 and @xmath1 with typical energies of 10 mev have comparable absorption and scattering opacities ( see fig .
3 ) , the emergent @xmath0 and @xmath1 energy fluxes and spectra are determined essentially in the same region where these two neutrino species just begin to thermally decouple from the proto - neutron star matter .
the @xmath0 and @xmath1 occupation numbers in this region are then approximately given by @xmath114 + 1 } , \label{occnue}\ ] ] and @xmath115 + 1 } , \label{occnueb}\ ] ] respectively , where @xmath116 , and @xmath117 is the @xmath0 chemical potential . from these occupation numbers ,
we obtain @xmath118 and @xmath119 where for example , @xmath120 and @xmath121 are the number and energy densities , respectively , for the @xmath0
. we can estimate the magnitude of @xmath122 as follows .
the net deleptonization rate of the proto - neutron star is @xmath123 where @xmath124 is the initial electron fraction of the proto - neutron star , @xmath125 is the proto - neutron star mass , @xmath126 is the nucleon mass , @xmath127 is the neutrino diffusion timescale , and @xmath128 is the proto - neutron star radius .
the sum of the @xmath0 and @xmath1 luminosities is @xmath129 from eqs .
( [ nuno])([nulum ] ) , we have @xmath130 where the numerical value is obtained for @xmath131 and @xmath132 km .
therefore , to a good approximation , the @xmath0 and @xmath1 in the decoupling region have fermi - dirac energy distributions with zero chemical potentials .
let @xmath133 be the specific intensity of neutrinos of a given species , where @xmath134 is the spatial position , and @xmath135 specifies the direction of propagation .
the general transport equation takes the form i_e=_e^i_e^(fd ) + ^_e 14d i_e - i_e , [ inu]where @xmath136 , and the subscript `` @xmath137 '' indicates that the relevant physical quantities depend on the neutrino energy @xmath137 .
the opacities @xmath138 and @xmath139 are related to the corresponding cross sections as , e.g. , @xmath140 and @xmath141 for @xmath0 .
note that we have included the effect of stimulated absorption for neutrinos ( e.g. , imshennik & nadezhin 1973 ) by introducing the corrected absorption opacity @xmath142 $ ] ( @xmath55 is the matter temperature ) .
hereafter , we shall suppress the superscript `` @xmath143 '' in @xmath138 , and the notation @xmath144 should be understood as having included the correction due to stimulated absorption .
the quantity @xmath145 in eq .
( [ inu ] ) is the equilibrium neutrino intensity , given by the fermi - dirac distribution i_e^(fd)=e^3c^2h^31(e
/ t)+1 . note that the time derivative term @xmath146 on the left - hand - side of eq .
( [ inu ] ) has been neglected .
this corresponds to an instantaneous redistribution of emission and absorption sources ( i.e. , an instantaneous redistribution of matter temperature )  the timescale for the redistribution is of order the mean free path divided by @xmath25 , much smaller than the neutrino cooling time of the star .
the zeroth order moment of the transport equation can be written as _
e=_e^ c(u_e^(fd)-u_e ) , [ dfnu]where @xmath147 is the ( spectral ) neutrino energy flux , @xmath148 is the energy density , and @xmath149 is the corresponding equilibrium energy density . with the eddington closure relation , the first moment of eq .
( [ inu ] ) becomes _ e =- c3_e^(t)u_e , [ fnu]where @xmath150 .
equations ( [ dfnu ] ) and ( [ fnu ] ) completely specify the neutrino radiation field .
here we are concerned with the total neutrino energy flux emitted from a local surface region of the neutron star .
the magnetic field strength in this local region is assumed to be constant .
since the scale heights for physical quantities such as the matter density and temperature are much smaller than the local radii , we adopt the plane - parallel geometry . integrating eq .
( [ fnu ] ) over the neutrino energy , we find that the total neutrino energy flux ( of a specific species ) @xmath151 is given by f =- c3_rdudr , [ ftot]where @xmath152 is the total neutrino energy density , and @xmath153 is the `` rosseland mean opacity '' defined as = de(1/_e^(t))(d u_e / dr)d e(du_e / dr ) .
[ ross]since neutrino emission and absorption rates are extremely fast even in the regions near the proto - neutron star surface , a radiative equilibrium , @xmath154 , @xmath20 , and @xmath155 is much weaker , the overall radiative equilibrium in the region near the proto - neutron star surface mainly involves matter and @xmath0 and @xmath1 .
since the coupling between matter and @xmath0 is comparable to that between matter and @xmath1 , we further assume that an approximate radiative equilibrium holds separately for @xmath0 and @xmath1 .
] approximately holds , i.e. , _
e^ ( u_e^(fd)-u_e)de0 . > from eq .  ( [ dfnu ] ) , we have @xmath156 , or @xmath157 is approximately constant .
if we introduce the rosseland mean optical depth measured from the stellar surface via @xmath158 , with @xmath157 constant , eq .
( [ ftot ] ) can be integrated to yield cu=3f(_r+23 ) , [ profile]where we have used the boundary condition @xmath159 at the surface ( @xmath160 ) .
equation ( [ profile ] ) resembles the standard eddington profile .
note that outside the decoupling region the neutrino energy density u may not be equal to @xmath161 specified by the matter temperature @xmath55 .
similarly , @xmath162 may differ substantially from @xmath163 ( see  6 ) .
however , we can still approximate @xmath162 by the fermi - dirac function at the neutrino temperature @xmath164 , which is equal to the matter temperature at the decoupling layer .
equation ( [ ross ] ) can then be rewritten as = 307 ^ 4t_^5de1_e^(t ) e^4 ( e / t_)^2 . [ rosst]with the cross sections given in  2 , we can calculate the rosseland mean opacity @xmath165 in magnetic fields using eq .
( [ rosst ] ) .
note that in eq .
( [ rosst ] ) @xmath166 depends on the local matter temperature @xmath55 . only in the regions not far from the decoupling layer , do we have @xmath167 .
figure 4 shows the ratio of the rosseland mean opacity @xmath165 of @xmath0 to the zero - field value @xmath168 as a function of the magnetic field strength for typical conditions near the proto - neutron star surface .
we choose these conditions to be @xmath169  mev , and @xmath170  mev , all with @xmath72 .
figure 5 shows similar results for @xmath1 with similar physical parameters ( @xmath171  mev ; note that the results are insensitive to @xmath54 ) .
we see that deviation of @xmath165 from the zero - field value by more than @xmath172 generally requires field strength @xmath173 .
the behavior of the ratio @xmath174 as a function of @xmath39 is certainly not obvious a priori . since the rosseland mean integral is dominated by neutrinos with the smallest opacities , and since @xmath175 scales linearly with @xmath98 in the regime where only the ground landau level is filled ( see eq .
[ [ sigbapp ] ] ) , we expect that for sufficiently large @xmath98 , the rosseland mean opacity increases linearly with increasing @xmath98 .
this is indeed what is indicated in figs .  4 and 5 .
for `` intermediate '' field strength , the @xmath175 curve for @xmath0 ( see fig
.  1 ) has a significant portion which is below the zero - field value at low neutrino energies . for such intermediate fields
, we would expect @xmath165 to be less than the zero - field value .
this explains the non - monotonic behavior of @xmath174 as a function of @xmath98 shown in fig .
4 . for @xmath1 ,
the non - monotonic behavior is much less pronounced ( see the curve for @xmath176 mev in fig .
5 ) since the energy range in which @xmath177 is smaller .
imagine a proto - neutron star with an asymmetric magnetic field topology , e.g. , the north pole and the south pole have different _ magnitudes _ of magnetic field ( say , by a factor of a few ) .
what is the difference in the neutrino fluxes from these two regions due to the different absorption opacities of @xmath0 and @xmath1 , and what is the resulting kick velocity of the neutron star ?
we recall that the fractional asymmetry @xmath178 in the total neutrino emission ( of all species ) required to generate a kick velocity @xmath179 is @xmath180 \,[e_{\rm tot}/(3\times 10^{53}~{\rm erg})]$ ] , where @xmath181 is the total neutrino energy radiated from the neutron star ( the neutron star mass has been set to @xmath182 ) .
to quantitatively understand how the magnetically - modified opacities change the local neutrino flux , we will have to evaluate these opacities at different depths ( corresponding to different physical conditions ) below the neutron star surface , and to perform complete cooling calculations which take into account these opacities .
however , figs .
4 and 5 clearly indicate that , even near the stellar surface , where one expects the modification to the opacities to be most prominent , and even with the most favorable conditions for the magnetic field to be effective in changing the opacities ( i.e. , low density and low temperature ) , a @xmath172 change in the rosseland mean opacity ( for @xmath0 or @xmath1 ) requires a magnetic field with @xmath183 , or @xmath184
g. in the following , we shall content ourselves with merely estimating the change in the local neutrino flux from the modified opacities .
we first consider the case in which the matter temperature throughout the neutron star is unaffected by the magnetic field .
one might imagine this as a possibility if convection dominates energy transfer between different parts of the star .
under the neutrino - matter decoupling layer the specific neutrino energy density @xmath162 is the same as the fermi - dirac function @xmath163 specified by the matter temperature .
the emergent neutrino flux is then inversely proportional to the rosseland mean optical depth @xmath185 from the decoupling layer to the stellar surface ( see eq .  [ [ profile ] ] ) .
since @xmath185 for @xmath0 or @xmath1 at the decoupling layer is of order a few , the flux @xmath157 is inversely proportional to @xmath153 evaluated near the decoupling region .
if we use the curve for @xmath74  mev and @xmath186  mev in fig .  4 ( which has the most pronounced magnetic correction ) , we find that the fractional change in the @xmath0 flux is @xmath187 at @xmath70 ( if we used the curve for @xmath188 or 20  mev , which more closely represents the physical condition at the decoupling region , the modification @xmath189 would be smaller by a factor of a few or more ) . similarly , using the curve for @xmath74  mev in fig .  5 to maximize the magnetic correction , we find @xmath190 at @xmath70 . since the fluxes of other types of neutrinos are unchanged and @xmath191 depend on the @xmath192-@xmath193 scattering cross sections , which can be affected by the magnetic field .
however , the layers are located much deeper than those of @xmath0 and @xmath5 ( since @xmath192-@xmath193 scattering opacity is small ) , the modification to the @xmath192-@xmath193 scattering opacity due to electron landau levels is therefore smaller than the modification to the absorption opacities of @xmath4 and @xmath5 considered in this paper . ] , and since different neutrino species more or less carry away the same amount of energy , the change in the total neutrino flux is only @xmath194 at @xmath70 . if we now imagine that the north pole of the neutron star has @xmath195 while the south pole has a field strength a factor of a few smaller , then , taking into account the geometric reduction ( @xmath196 ) , the flux asymmetry is about @xmath197 , just enough to give a kick velocity of @xmath21  km  s@xmath7 to the neutron star . with @xmath153
evaluated at more realistic conditions , we expect the resulting @xmath178 to be even smaller .
the assumption that the matter temperature is unaffected by the magnetic field ( and therefore spherically symmetric ) is unlikely to be valid . in reality ,
when a region of the star ( say the north pole ) emits more flux , it also cools faster .
with reduced temperature , this region will tend to emit less flux . in the extreme case
, one would find @xmath198 , with @xmath199 the neutrino energy density in the stellar core , and @xmath185 the rosseland mean optical depth from the core to the surface .
since in the bulk interior of the star , modification to @xmath153 due to the magnetic field is negligible , we would conclude that the change in the neutrino flux is much smaller than what is estimated in the last paragraph .    to summarize , the modifications to the absorption opacities of @xmath0 and @xmath5 due to quantized landau levels give rise to rather small change in neutrino fluxes from the neutron star .
to generate appreciable kick velocities ( @xmath6  km  s@xmath7 ) from this effect requires the difference of field strengths at the two opposite poles to be at least @xmath8  g ( and possibly much larger ) .
we now consider how the oscillatory behaviors ( see figs .
13 ) of the absorption opacities of @xmath0 and @xmath1 may manifest themselves in the radiated neutrino energy spectra . for an
absorption dominated medium , neutrinos with energy @xmath137 decouple from the matter and start free - streaming at a column depth ( in g  @xmath200 ) @xmath201 .
thus the neutrino energy spectrum is given by @xmath202 with @xmath203 evaluated at the depth @xmath204 . in the presence of significant scattering
, decoupling occurs at depth @xmath205 with the effective opacity @xmath206 , whereas the neutrinosphere from which neutrinos start free - streaming is located at @xmath207 . near the decoupling region
we have @xmath208 , and thus from eq .
( [ fnu ] ) the neutrino spectral energy flux is given by f_e~_e^eff_e^(t ) cu_e^(fd)(y~1/_e^eff ) .
the above qualitative argument can be made more rigorous by directly solving eqs .
( [ dfnu ] ) and ( [ fnu ] ) .
consider a plane - parallel atmosphere and define the ( energy - dependent ) optical depth via @xmath209 .
equations ( [ dfnu ] ) and ( [ fnu ] ) can be combined to give u_e=_e^2 [ u_e - u_e^(fd ) ] ,   with  _
e^2=3_e^_e^(t ) . with a constant @xmath210
, the above equation can be solved exactly . for the boundary condition
@xmath211 at @xmath212 , the neutrino energy flux at the surface is given by f_e=_e/31 + 2_e/3_0^d_e_ecu_e^(fd)(_e)(-_e_e ) cu_e^(fd)(_e_e^-1 ) .
[ fnu2 ]    figure 6 shows the spectral energy fluxes of @xmath0 calculated using eq .
( [ fnu2 ] ) for different field strengths ( @xmath70 and @xmath213 ) .
the matter temperature profile is obtained from eq .
( [ profile ] ) assuming @xmath214 . for a total energy flux of @xmath215 ( where @xmath216 is the effective temperature )
, we have @xmath217^{1/4}$ ] .
we choose @xmath218 mev in these examples . for simplicity ,
we have neglected the variations of the opacities as a function of @xmath219 .
thus the temperature at decoupling ( @xmath220 ; see eq .
[ [ fnu2 ] ] ) is given by @xmath221^{1/4}$ ] , where @xmath222 and we take the relevant opacities to have the corresponding values for @xmath72 , @xmath73 mev , and @xmath74 mev
. this simplification should not introduce qualitative changes in the spectrum . as seen from fig .
6 , the spectrum is pinched at the high energies because of the @xmath223-dependence of the opacities .
the magnetic field introduces modulations in the spectrum . for lower magnetic field strengths ,
the modulations occur at the lower energy band . as the magnetic field gets stronger , the modulations shift to the higher energy part of the spectrum .
the distortions and modulations of the @xmath0 and @xmath1 energy spectra by strong magnetic fields may offer an interesting possibility to limit the magnetic field strength in proto - neutron stars through the detection of neutrinos from a galactic supernova .
we note , however , that the @xmath0 spectra shown in fig .  6
are based on very approximate calculations .
they serve only as illustrative examples of the magnetic field effect on the neutrino spectra .
more detailed calculations are needed to obtain realistic spectra which include the magnetic field effect .
in this paper , we have studied the issue of whether asymmetric magnetic field topology in a proto - neutron star can induce asymmetric neutrino emission from the star through the modifications of the neutrino absorption opacities by the magnetic field .
these modifications arise from the quantized landau levels of electrons and positrons produced in strong magnetic fields . by calculating the appropriate mean neutrino opacities in the magnetic field
, we demonstrate that this mechanism is rather inefficient in generating a kick to the neutron star : to obtain appreciable kick velocities ( @xmath6  km  s@xmath7 ) , the difference in the field strengths at the two opposite poles of the star must be at least @xmath8  g.    this work was started while d.l .
held the richard c. tolman postdoctoral fellowship at caltech .
additional support was provided by nasa grant nag 5 - 2756 and nsf grant ast-9417371 , as well as a start - up fund for new faculty at cornell university .
is supported by the david w. morrisroe postdoctoral fellowship at caltech .