the existence of gravitational waves is an unambiguous prediction of the theory of general relativity @xcite . yet , despite efforts originating in the early nineteen sixties , gravitational waves have not been detected directly .
nevertheless , observation of the binary pulsar psr 1913 + 16 convinces us that gravitational waves do exist , and that they are correctly described by einstein s theory @xcite .
that gravitational waves have not yet been detected on earth is simply due to their incredible weakness : typical waves would produce in a bulk of matter a strain @xmath0 , where @xmath1 is the extension of the matter , of order @xmath2 @xcite .
needless to say , to measure this effect is a great challenge for experimentalists .
there is reasonable hope that gravitational waves will be detected within the next ten years , thanks to a new generation of detectors which use interferometry to monitor the small displacements induced by the passage of a gravitational wave .
two groups are currently involved in building large - scale interferometers : the american ligo team , and the french - italian virgo team .
the ligo ( laser interferometer gravitational - wave observatory ) project @xcite involves two detectors , to be built in hanford , washington , and in livingston , louisiana .
construction has begun at both sites .
each interferometer has an armlength of approximately 4 km .
ligo should be completed by the turn of the century .
the virgo ( so named after the galaxy cluster ) project @xcite involves a single interferometer , to be built near pisa , italy , with an armlength of approximately 3 km .
virgo should also be completed by the turn of the century .
the basic idea behind interferometric detectors is the following @xcite :    the interferometer is composed of two long ( 4 km for ligo ) vacuum pipes forming the letter * l*. a laser beam is split in two at the corner of the * l * , and is sent into each arm of the interferometer .
each beam then bounces off a mass which is suspended at each end of the * l * ( a mirror has been coated onto each mass ) .
the light is finally recombined at the beam splitter , and its intensity is measured by a photodiode .    when no gravitational wave is present at the interferometer , the length of each arm is so adjusted that when measured by the photodiode , the light s intensity is precisely zero ( the recombined beams are arranged to be precisely out of phase ) . however , when a gravitational wave passes through the interferometer , the armlengths are no longer constant , and the recombined beams no longer precisely out of phase .
more precisely , during the first half of its cycle the gravitational wave increases the length of one arm , and decreases the length of the other . during the second half cycle ,
the first arm is now shorter , and the second arm longer .
the light s intensity therefore oscillates with the gravitational - wave frequency .
the intensity is a measure of @xmath3 , where @xmath1 denotes the interferometer armlength , and @xmath4 the gravitational - wave field .
interferometers are subject to various sources of noise which limit the detector s sensitivity to gravitational waves .
the relative importance of each source depends on the frequency at which the interferometer oscillates @xcite .    at low frequencies ( @xmath5 )
the detector s sensitivity is limited by seismic noise , which is due to the earth s seismic activity . at frequencies larger than 10
hz the seismic noise can be eliminated with sophisticated isolation stacks ; these fail at low frequencies .    at high frequencies ( @xmath6 ) the detector s sensitivity is limited by photon shot noise , which is due to statistical errors in the counting of photons by the photodiode .
this source of noise can be reduced by increasing the laser power , or making use of `` light recycling '' @xcite .    at intermediate frequencies ( @xmath7 between 10 hz and 100 hz )
the noise is dominated by thermal noise , which is due to spurious motions of thermal origin .
for example , the suspended masses are thermally excited and vibrate with their normal - mode frequencies ; this evidently affects the recombined laser beam .
interferometers are therefore broad - band detectors , with good sensitivity in the range @xcite @xmath8 the required sensitivity for full - scale interferometers is approximately @xmath9 at peak sensitivity  a tall order .
[ the subscript @xmath10 stands for `` noise level '' ; we will define @xmath11 precisely below . a plot of @xmath11 , appropriate for an interferometric detector with `` advanced '' sensitivity , is given in fig .  1 . ] for comparison , we may mention that the caltech 40 m prototype has already achieved @xmath12 at peak sensitivity @xmath13 .
it is not implausible that improved technology and a factor of 100 in armlength will permit to reach the desired goal .
the detector noise can be measured when no gravitational wave is present at the interferometer , the typical situation . then the detector output @xmath14 is given by noise alone : @xmath15 where @xmath16 represents the noise .
the noise is a random process @xcite : the function @xmath16 takes purely random values .
consequently , the noise can only be studied using statistical methods . in the presence of a gravitational wave , eq .
( [ 2 ] ) must be replaced by @xmath17 , where @xmath18 is the gravitational - wave field .
the statistical properties of the noise can be determined by careful measurement .
for example , the time average @xmath19 can be constructed .
this mean value can then be subtracted from @xmath16 and , without loss of generality , we can put @xmath20 . also from measurements , the noise s _ autocorrelation function _ @xmath21 can be constructed : @xmath22 @xmath23 gives the mean squared deviation of the noise with respect to the mean value .    in the following
we will assume that the noise is _ stationary _ , in the sense that its statistical properties do not depend on time @xcite .
this means , in particular , that the autocorrelation function does not depend explicitly on the origin of time @xmath24 , but only on the variable @xmath25 , as was expressed in eq .
( [ 4 ] ) .
we shall also assume that the noise satisfies the _ ergodic hypothesis _
, so that time averages can be replaced with ensemble averages @xcite . here , the noise is imagined to be drawn from a representative ensemble , and the probability that it takes a particular realization @xmath16 is given by a specified probability distribution .
the statistical properties of the noise then refer to this ( infinite dimensional ) distribution function .    in full generality , the statistical properties of the noise
can only be summarized by constructing all the higher moments @xmath26 .
if , however , the noise is assumed to be gaussian , in the sense that its probability distribution function is an infinite dimensional gaussian distribution @xcite , then the autocorrelation function contains all the information .
real detector noise is neither strictly stationary nor strictly gaussian .
however , on a timescale of hours , which is long compared with typical gravitational - wave bursts , the noise appears stationary to a good approximation @xcite . and non - gaussian components to the noise can be removed , to a large extent , by cross - correlating detector outputs from two widely separated interferometers @xcite .
it is therefore a satisfactory approximation to take the noise to be stationary and gaussian .    under these assumptions the statistical properties of the detector noise
are fully summarized by @xmath21 .
it is convenient to work instead in the frequency domain , and to define @xcite the noise s _ spectral density _ @xmath27 as @xmath28 the spectral density is defined for @xmath29 only ; as @xmath21 is a real and even function , the negative frequencies only duplicate the information contained in the positive frequencies .    as @xmath16 is dimensionless , the spectral density has dimensions of time . by multiplying @xmath27 with the frequency and taking the square root ( since the spectral density represents the mean squared noise ) , one obtains the _ noise level _
@xmath11 : @xmath30 this gives the equivalent gravitational - wave amplitude which would make the interferometer oscillate at just the noise level @xcite ; this quantity was introduced in the preceding section ( see fig .  1 ) .
coalescing compact binaries , composed of neutron stars and/or black holes , are the most promising source of gravitational waves for interferometric detectors @xcite .
consider the binary pulsar psr 1913 + 16 @xcite .
this system consists of two neutron stars , each of @xmath31 , in orbital motion around each other .
its present orbital period @xmath32 is approximately 8 hours , corresponding to orbital separations of about @xmath33 times the total mass .
( here and throughout we use units such that @xmath34 . )
its present eccentricity is approximately equal to 0.6 .
the binary s orbital period is observed to decay at a rate @xmath35 corresponding precisely to a loss of energy and angular momentum to gravitational waves @xcite . in a timescale of approximately @xmath36 years the orbital period will have decreased to less than a tenth of a second , corresponding to orbital separations smaller than one hundred times the total mass . in this time
, the eccentricity will have been reduced to extremely small values ( by the radiation reaction ) , so that the orbits are practically circular .
the gravitational waves produced then have a frequency larger than 10 hz , and the frequency keeps increasing as the system evolves .
_ during this late stage of orbital evolution , the gravitational waves sweep through the frequency bandwith of interferometric detectors , and thus become visible .
_ by the time the orbital separation becomes as small as a few times the total mass , the neutron stars begin to merge .
the gravitational waves produced during the final merger can not be detected by interferometric detectors , at least in the broad - band configuration described above @xcite : the gravitational - wave frequency is then larger than 1 000 hz , for which the detector noise is large .    of course
, it would be foolish to wait @xmath36 years in order for psr 1913 + 16 to produce gravitational waves with appropriate frequency .
fortunately , interferometric detectors will be sensitive enough to monitor binary coalescences occurring in quite a large volume of the universe , approximately @xmath37 ( corresponding to a radius of 200 mpc @xcite ) .
it has been estimated @xcite that as many as 100 coalescences could occur every year in such a volume ( this includes coalescences of black - hole systems as well ) .
this potentially large event rate is one of the factors that make gravitational waves from coalescing compact binaries especially promising .
the other factor comes from the fact that compact binaries are extremely clean astrophysical systems .
it can be estimated @xcite that tidal interactions between the two stars are completely negligible , up to the point where the objects are about to merge .
in particular , no mass transfer occurs .
the system can therefore be modeled , to extremely good accuracy , as that of two point masses with a limited number of internal properties ( such as mass , spin , and quadrupole moment ) .
the challenge in modeling coalescing compact binaries resides in formulating and solving the equations of motion and wave generation for a general relativistic two - body problem @xcite .
at the crudest level , the gravitational - wave signal corresponding to an inspiraling binary system can be calculated by ( i ) assuming that the orbital motion is newtonian ( with the effects of radiation reaction incorporated ) and ( ii ) using the standard quadrupole formula for wave generation @xcite . as motivated above
, we may also assume that the orbits are circular .
we define @xmath18 to be the gravitational - wave signal .
this is given by a linear combination , appropriate for interferometric detectors , of the two fundamental polarizations , @xmath38 and @xmath39 , of the gravitational - wave field .
at this level of approximation , the signal is given by @xcite @xmath40 as expected , the signal decays as the inverse power of @xmath41 , the distance to the source .    in eq .
( [ 7 ] ) , @xmath42 is a function of all the angles relevant to the problem : position of the source in the sky , orientation of the orbital plane , position and orientation of the detector on earth .
the parameter @xmath43 is called the _ chirp mass _ and represents a particular combination of the masses , given by @xmath44 the waveform depends on the chirp mass only , and not on any other combination of the masses .
the symbol @xmath7 represents the gravitational - wave frequency , which is equal to _ twice _ the orbital frequency . because the system loses energy and angular momentum to gravitational waves , the frequency is not constant , but increases in time according to @xcite @xmath45 as a consequence , eq .
( [ 7 ] ) shows that the amplitude of the signal , which is proportional to @xmath46 , also increases with time .
a signal which increases both in frequency and in amplitude is known as a _ chirp _ , and this is the origin of the term `` chirp mass '' .
finally , the phase function @xmath47 is given by @xmath48 because the frequency is not a constant , the phase accumulates nonlinearly with time .    for a system of two neutron stars ,
the gravitational - wave signal undergoes approximately 16 000 oscillations as it sweeps through the frequency bandwith of an interferometric detector .
the timescale for the frequency sweep is approximately 15 minutes .
the orbital separation ranges from approximately 180 to 10 times the total mass @xmath49 , and the orbital velocity @xmath50 ranges from approximately 0.1 to 0.4 .
this indicates that relativistic corrections must be inserted in eqs .
( [ 7 ] ) and ( [ 9 ] ) in order to obtain a satisfactory degree of accuracy .
how does one go about finding a gravitational - wave signal in a noisy data stream , when typically the signal is not very strong ? and once the signal is found , how does one go about extracting the information it contains ? for signals of precisely known form ,
one goes about this using the technique of _ matched filtering _
signals from inspiraling compact binaries , since they can be calculated with high precision , belong to this class .
the basic idea behind matched filtering is to use our knowledge about the signal in order to go find it in the data stream , after the noisy frequencies have been filtered out .
suppose a signal of known form @xmath51 is present in the data stream . here
, the vector @xmath52 collectively denotes all the parameters characterizing the signal . in the case of inspiraling binaries
, these would be the time of arrival , the initial phase , the distance to the source , the position angles , the chirp mass , and other parameters to be introduced in sec .
the detector output is given by @xmath53 where @xmath16 is the noise , whose statistical properties are fully summarized by the spectral density @xmath27 , as discussed in sec .
the first step in matched filtering @xcite is to pass @xmath54 through a linear filter which removes the noisy frequencies .
the idea here is to use our knowledge about the detector noise contained in @xmath27 in order to discard that part of the detector output for which the detector noise is large .
the output is imagined to be decomposed into fourier modes according to @xmath55 ; the filter suppresses the modes @xmath56 for which @xmath27 is large .    in mathematical terminology ,
a linear filter is a linear operation on a function @xmath54 .
this operation can always be written as @xcite @xmath57 where @xmath58 is the filter function .
the filter which removes the noisy frequencies is the one such that @xmath59 the factor of 2 is conventional .    the next step in matched filtering consists of computing the overlap integral between the filtered output ( [ 13 ] ) and the known signal @xmath60 .
the true value @xmath52 of the source parameters is not known prior to the measurement .
this operation must therefore be repeated for a large number of trial values @xmath61 ; the corresponding signals @xmath60 are known as _
templates_. the overlap integral defines the function @xmath62 given by @xmath63 to obtain the second line we have inserted the fourier decompositions of @xmath4 , @xmath64 , and @xmath65 , and carried out the integrations over time .
the inner product @xmath66 is defined as @xmath67 where @xmath68 and @xmath69 are arbitrary functions of time , with fourier transforms @xmath70 and @xmath71 .
a quantity analogous to @xmath62 can be defined for the noise alone , by replacing @xmath65 to the right of eq .
( [ 15 ] ) by @xmath10 .
operationally , this amounts to filtering the detector output when a gravitational - wave signal is known _ not _ to be present . because the noise is a random function , the integrals @xmath72 are random also . and because the noise is assumed to have zero mean , to take an average over all possible realizations of the noise ( by repeating the measurements many times ) would yield a zero value .
none of these quantities would be especially useful .
we are therefore led to consider the root - mean - square average of these integrals , @xmath73 to go from the first to the second line requires some machinery which will not be presented in this review .
we refer the reader to ref .
@xcite for the missing steps .
the _ signal - to - noise ratio _ is defined to be the ratio of @xmath74 over @xmath75 : @xmath76 it is clear that @xmath77 is a random function , since @xmath54 is itself a random function , being the superposition of signal plus noise .
its expectation value , or average over all possible realizations of the noise , is zero in the absence of signal ( since @xmath20 ) , and @xmath78 in the presence of the signal @xmath51 .
it is important to notice that in eq .
( [ 19 ] ) , the numerator is the overlap integral between the true signal @xmath51 and the templates @xmath60 . in the absence or presence of a signal ,
the variance in the signal - to - noise ratio is precisely equal to unity @xcite , independently of the values of @xmath52 and @xmath61 .
this fully summarizes the statistical properties of the signal - to - noise ratio , since @xmath77 is a gaussian random function .
[ this can be seen from the fact that @xmath54 itself is gaussian , being the superposition of signal plus gaussian noise . ]
it is intuitively clear that choosing a template with @xmath79 will produce the largest possible expectation value of the signal - to - noise ratio .
this statement , known as _
wiener s theorem _
@xcite , can easily be shown to be true by applying schwarz s inequality to the right - hand side of eq .
( [ 19 ] ) .
we have @xmath80 the fact that the signal - to - noise ratio has a variance of unity indicates that a signal can be concluded to be present only if @xmath81 is significantly larger than 1 .
we will come back to this point in the next section .
the true value of the source parameters can therefore be determined by maximizing @xmath82 over all possible values of the trial parameters .
the fact that the signal - to - noise ratio has a variance of unity implies that this determination can only have a limited degree of accuracy .
the statistical errors decrease with increasing @xmath81 ; since this is proportional to the signal s amplitude , a stronger signal gives better accuracy .
maximizing the signal - to - noise ratio is essentially equivalent to maximizing the overlap integral @xmath83 it is easy to see how the choice of the parameters @xmath61 can affect the overlap integral .
consider a toy waveform with three parameters : arrival time , initial phase , and chirp mass .
a mismatch in the arrival times clearly reduces the overlap integral : the signal and the template , taken to be functions of time , might have support in entirely different regions of the time axis , leading to a vanishing overlap . supposing that the arrival times are matched
, a mismatch in the initial phases can also reduce the overlap integral : the signal and the template might be out of phase with each other , leading to an oscillating integrand and a vanishingly small overlap . finally , supposing that both the arrival times and initial phases are matched , a mismatch in the chirp masses would also reduce the overlap integral .
this is because the chirp mass governs the rate at which the signal s frequency changes with time ; cf .
( [ 9 ] ) .
signal and template , starting at the same time with the same phase , might thereafter go out of phase , thereby reducing the overlap .
the first order of business when analysing the output of a gravitational - wave detector is to decide whether or not a signal is present . here we assume that the signal must be of a specific form , corresponding to a coalescing binary system . in this section
we discuss signal detection  how the technique of matched filtering can be employed to find the signal in the noisy data stream . in the next section
we will discuss signal measurement  how matched filtering is used to estimate the value of the source parameters once the signal has been found .    as mentioned in the previous section
, a signal can be concluded to be present if the maximum value of the signal - to - noise ratio , @xmath84 , is significantly larger than unity .
in fact , there exists a threshold value @xmath85 such that a signal is concluded to present , with a certain confidence level , if @xmath86 . to figure out how large this threshold must be is a standard application of the statistical theory of signal detection @xcite , which was developed largely for the purpose of detecting radar signals .
the theory can easily be taken over to the case gravitational - wave signals @xcite . to go into the detail of this theory
would be outside the scope of this review .
we shall simply state that @xmath85 is fixed by selecting a small , acceptable value for the probability that a signal would falsely be concluded to be present  the false alarm probability .
( this is the neyman - pearson criterion @xcite , which is more precisely formulated in terms of the likelihood ratio , the ratio of the probability that a signal is present to the probability that it is absent . )
once @xmath85 is fixed , the level of confidence that a signal is indeed present increases with @xmath86 .
a typical ballpark value for the threshold is @xmath87 .    in the preceding paragraph
it was assumed that the signal - to - noise ratio @xmath88 is computed using template waveforms which are functions of the parameters @xmath52 introduced in sec .
these parameters have direct physical meaning ; they include the chirp mass and other meaningful parameters to be introduced in sec .  12 .
however , since these parameters are _ not _ estimated during the detection stage of the data analysis ( they are estimated only _
after _ a signal has been found ) , there is no particular need to parametrize the templates with @xmath52 .
in fact , it may be desirable , in order to minimize the computational effort , to parametrize the signal in a completely different way .
the new parameters , @xmath89 , would then have no particular physical significance .
what is required is that the new templates , @xmath90 , reproduce the behaviour of the expected gravitational - wave signal . in other words , the templates @xmath51 and
@xmath90 should span the same `` signal space '' , but @xmath90 should do so most efficiently . unlike @xmath51 ,
@xmath90 need not be derived from the field equations of general relativity . with these new _ detection templates _ , the signal - to - noise ratio is defined as @xmath91 where @xmath92 is defined as in eq .
( [ 18 ] ) . as before
, a signal is concluded to be present if @xmath93 .
once the detection templates @xmath90 have been used to conclude that a signal is present in the data stream , they are replaced with the _ measurement templates _ @xmath60 in order to estimate the value of the physical parameters @xmath52 .
the procedure to estimate the source parameters was explained in sec .
as was mentioned , the idea is to maximize over all possible values of the trial parameters the overlap integral @xmath94 we already have discussed the effect of a mismatch in the value of the parameters .
what remains to be discussed is the effect of a mismatch in the functional form of the template with respect to that of the true waveform .
the true waveform is governed by the exact laws of general relativity .
the template , on the other hand , is necessarily constructed using an approximation to the exact laws .
( that approximations must be made in analytic calculations is obvious ; in a numerical treatment the approximation resides in the finite differencing of the field equations . )
this , clearly , must have an effect on our strategy for extracting the information contained in the gravitational - wave signal .
this can be seen simply from the fact that wiener s theorem , as stated in sec .
8 , strictly requires signal and template to have the same functional form ; they are allowed to differ only in the value of their parameters .
suppose that gravitational waves coming from a given source are received without noise , so that the true signal @xmath51 is measured accurately .
suppose also that the true value @xmath52 of the source parameters is known ( god has spoken ) .
then a computation of the overlap integral with @xmath79 but with an approximate template @xmath60 will not yield the maximum possible value for @xmath77 .
this is because the approximation differs from the true signal , both in amplitude and in phase .
since both signal and template undergo a large number of oscillations ( recall that for a system of two neutron stars , this number is approximately 16 000 ) , the overlap integral is most sensitive to phase differences : a slight phase lag causes the integrand to oscillate , thereby severely reducing the signal - to - noise ratio with respect to its maximum possible value .
a gravitational - wave astronomer does nt know before the measurement the true value of the source parameters , and must work with an approximation to the true general - relativistic waveform .
we have seen that the phase lag occurring between the true signal and the approximate template when the parameters are matched reduces the signal - to - noise ratio from its maximum possible value .
it follows that maximizing @xmath77 with approximate templates introduce _ systematic errors _ into the estimation of the source parameters : evaluating the signal - to - noise ratio with @xmath95 will return , for some @xmath96 , a number larger than @xmath97 .
the systematic errors are precisely the value of @xmath96 for which the signal - to - noise ratio is largest
. if the templates are a poor approximation to the true signal , then the systematic errors will be larger than the statistical errors arising because @xmath77 is a random function ( see sec .
we therefore appreciate the need for constructing measurement templates which are as accurate as possible , especially in phase @xcite .
the requirement is that the systematic errors in the estimated parameters must be smaller than the statistical errors .
an estimate for the required degree of accuracy comes from the observation that the overlap integral will be significantly reduced if the template loses phase by as much as one wave cycle with respect to the true signal .
since the total number of wave cycles is approximately 16 000 , we have @xmath98 since the orbital velocity @xmath99 is of order @xmath100 when the gravitational - wave frequency is in the relevant bandwidth , relativistic corrections _ at least _ of order @xmath101 are required to improve the quadrupole - formula expression given in sec .
as we shall see , this is an underestimate .
we have seen that accurate measurement templates are required to make the most of the gravitational - wave signals originating from coalescing compact binaries .
we also have seen that the quadrupole - formula waveform , eq .
( [ 7 ] ) , is not sufficiently accurate ; relativistic corrections at least of order @xmath101 are required .
how does one go about calculating these ?
a possible line of approach is to use a slow - motion approximation to the equations of general relativity .
this is based on the requirement that if @xmath99 is a typical velocity inside the matter source , then @xmath102 we shall call this approximation `` post - newtonian theory '' @xcite .
we point out that for binary systems , post - newtonian theory makes no assumption regarding the relative size of the two masses .
this is to be contrasted with the perturbation approach , discussed in sec .  13 , in which the mass ratio is assumed to be small , but no restriction is put on @xmath99 .    in post - newtonian theory @xcite , one starts by defining fields @xmath103 as @xmath104 where @xmath105 is the inverse of the true metric @xmath106 with determinant @xmath107 , and @xmath108 is the metric of minkowski spacetime .
when the harmonic gauge conditions @xmath109 are imposed , the _ exact _ einstein field equations reduce to @xmath110 here , @xmath111 is the flat - spacetime wave operator , @xmath112 is the stress - energy tensor of the source , and @xmath113 is nonlinear in @xmath103 and represents an effective stress - energy tensor for the gravitational field .    the detailed way in which one solves these equations is quite complicated , and will not be described here .
the essential ideas are these @xcite :    one first integrates the equations in the near zone ( @xmath114 , where @xmath41 is the flat - space radius and @xmath115 the gravitational wavelength ) assuming slow motion , or @xmath116 .
one does this by iterations : the nonlinear terms in eq .
( [ 26 ] ) are first neglected , and the resulting linear equations integrated .
these solutions are then used as input for the next iteration .
this process is continued until the desired degree of accuracy is obtained .
this is the standard post - newtonian approach @xcite .
one next integrates the equations everywhere in the vacuum region outside the source .
this is done once again by iterations , assuming @xmath117 , but assuming nothing about the relative size of @xmath118 with respect to the spatial derivatives .
this is because the vacuum region contains the wave zone , in which the field propagates with the speed of light ; a slow - motion assumption would therefore not do for the field itself .
this is the post - minkowskian approach @xcite .    using the post - minkowskian approach one constructs , by successive approximations , the most general solution to the einstein equations outside the source
this is characterized by two infinite sets of arbitrary multipole moments @xcite , the mass moments @xmath119 and the current moments @xmath120 , were @xmath121 and @xmath122 are the standard spherical - harmonic indices .
( in practice , the fields @xmath103 are equivalently expressed in terms of symmetric - trace - free moments , not spherical - harmonic moments . ) the general solution is then matched to the near - zone solution in the region of common validity , and the multipole moments are thus determined .
finally , one expresses the radiation field  the time - varying , @xmath123 part of the gravitational field  in terms of the derivatives of the mass and current multipole moments @xcite .
this gives the gravitational waveform .
the gravitational - wave luminosity @xmath124 can also be obtained from the radiation field .
to date , the post - newtonian calculation of the waveform has been carried out accurately through second post - newtonian order  @xmath125  beyond the leading - order , quadrupole - formula expressions given in sec .
the complete waveform will not be displayed here .
instead , we will focus solely on the waveform s _
phasing_.    the phasing of the waves can be determined from @xmath126 , the rate of change of the gravitational - wave frequency .
this can be expressed as @xmath127 where @xmath124 is ( minus ) the gravitational - wave luminosity , and @xmath128 relates orbital energy to orbital frequency ( @xmath7 is twice the orbital frequency ; the orbits are assumed to be circular ) .
both these quantities can be expanded in powers of @xmath129 with leading - order terms @xcite @xmath130 where the subscript @xmath75 stands for `` newtonian '' .
we have introduced the reduced mass @xmath131 , the total mass @xmath132 , and the mass ratio @xmath133 as @xmath134 notice that @xmath133 is restricted to the interval @xmath135 , with @xmath136 for @xmath137 .
it is easy to check that eqs .
( [ 27 ] ) and ( [ 29 ] ) reproduce eq .
( [ 9 ] ) above ; the chirp mass can be expressed as @xmath138 .    for simplicity we will , in the following , focus on the quantity @xmath124 .
to second post - newtonian order , the luminosity takes the form @xmath139 . \label{31}\ ] ] here , the terms so and ss are due to spin - orbit and spin - spin interactions , respectively @xcite ; these occur if the masses @xmath140 and @xmath141 are rotating .
let @xmath142 and @xmath143 be the spin angular momentum of each mass , and define the dimensionless quantities @xmath144 , for @xmath145 .
let also @xmath146 be the direction of orbital angular momentum .
then @xcite @xmath147 \hat{l } \cdot \vec{\chi}_a \label{32}\ ] ] is the spin - orbit term , and @xmath148 is the spin - spin term .    in eq .
( [ 31 ] ) , the leading - order term was first calculated in 1963 by peters and mathews @xcite using the standard quadrupole formula .
the first post - newtonian correction , at order @xmath149 , was calculated in 1976 by wagoner and will @xcite .
the @xmath150 term is due to wave propagation effects : as the waves propagate out of the near zone they are scattered by the curvature of spacetime , and this modifies both the amplitude and the phase of the waveform ; this `` tail term '' was first calculated in 1993 by this author @xcite , and then independently by wiseman @xcite and blanchet and schfer @xcite .
the second post - newtonian correction , at order @xmath101 , was calculated in 1994 by blanchet , damour , iyer , will , and wiseman @xcite .
finally , the spin - orbit and spin - spin corrections were calculated in 1993 by kidder , will , and wiseman @xcite .
we see from eq .
( [ 31 ] ) that the post - newtonian corrections bring a number of additional source parameters into the picture .
the waveform no longer depends uniquely upon the chirp mass @xmath43 ; it now depends upon the masses @xmath140 and @xmath141 separately , and upon the spin - orbit and spin - spin parameters ( which stay approximately constant as the system evolves toward coalescence ) .
these new parameters must be included into @xmath52 when the signal is analyzed using matched filtering .
what we have at this point is an expression for the waveform which is accurate to second post - newtonian order .
the question facing us is whether this waveform is sufficiently accurate to be used as measurement templates . in other words ,
are the systematic errors generated by these templates guaranteed to be smaller than the statistical errors ?
evidently , to answer this question is difficult , since we do not have access to the exact waveform in order to make comparisons . in the next section
we will consider a model problem for which the waveform _ can _ be calculated exactly , thereby enabling us to judge the accuracy of the post - newtonian expansion .
we will find , in sec .
15 , that the answer to this question is , most likely , no : the second post - newtonian waveform is not sufficiently accurate .
a different line of attack for solving einstein s equations for a compact binary system is to assume that one of the bodies is very much less massive than the other @xcite .
we therefore demand @xmath151 where @xmath131 is the reduced mass and @xmath132 the total mass . in this limit
@xmath131 is practically equal to the smaller mass @xmath140 , and @xmath132 is practically equal to the larger mass @xmath141 .
in contrast with the post - newtonian approach , we assume _ nothing _ about the size of the velocity @xmath99 .    this approach takes advantage of the fact that when eq .
( [ 34 ] ) is valid , the smaller mass creates only a small perturbation in the gravitational field of the larger mass , which can be taken to be the schwarzschild field .
the total gravitational field can therefore be written as @xmath152 where @xmath153 represents the background schwarzschild metric , and @xmath154 the perturbation .
writing einstein s equations for @xmath106 and linearizing with respect to @xmath154 , one finds that the perturbation must satisfy an inhomogeneous wave equation in the schwarzschild spacetime .
schematically , @xmath155 where @xmath156 is an appropriate curved - spacetime wave operator , and @xmath112 the stress - energy tensor associated with the orbiting mass
.    we will specifically assume that the central body is a schwarzschild black hole of mass @xmath132 .
this assumption is made for simplicity , and removes the need to model the star s interior . as a matter of fact ,
the internal structure of the bodies is irrelevant , except during the last few orbital cycles before coalescence ; this was discussed in sec .
taking advantage of this , we model the orbiting body as a point particle of mass @xmath131 , so that its stress - energy tensor is a dirac distribution with support on the particle s world line . for simplicity , and also because it is physically well motivated ( as explained in sec .  5 ) , we take the world line to be a circular geodesic of the schwarzschild spacetime . with @xmath157 as the usual schwarzschild coordinates , @xmath158 denotes the orbital radius , and @xmath159 is the angular velocity .
we have @xmath160 and the gravitational - wave frequency @xmath7 is given by @xmath161 .
we stress once more that in the perturbation approach , @xmath99 is not required to be small .
the only limitation on @xmath99 comes from the fact that for @xmath162 , circular orbits are no longer stable ; this implies @xmath163 .
black - hole perturbations are conveniently treated with the teukolsky formalism @xcite , in which gravitational perturbations are represented by the complex - valued function @xmath164 , a particular component of the perturbed weyl tensor ; the tensor @xmath154 can be reconstructed from @xmath164 .
the equation satisfied by @xmath164 admits a separation of the variables .
when @xmath164 is expanded in spherical harmonics and decomposed into fourier modes @xmath165 , one obtains an ordinary second - order differential equation  the teukolsky equation  for the radial function @xmath166 . here ,
@xmath121 and @xmath122 are the usual spherical - harmonic indices .
schematically , and omitting the subscript @xmath167 , this equation takes the form @xmath168 where @xmath169 is a second - order differential operator , and @xmath170 the source , constructed from the particle s stress - energy tensor @xcite .    equation ( [ 38 ] ) can be integrated in the standard way by constructing a green s function @xmath171 out of two linearly independent solutions to the homogeneous equation .
these solutions , @xmath172 and @xmath173 , respectively satisfy appropriate boundary conditions at the inner ( @xmath174 ) and outer ( @xmath175 ) boundaries .
schematically , @xmath176 , where @xmath177 ( @xmath178 ) is the lesser ( greater ) of @xmath41 and @xmath179 .
the solution to eq .
( [ 38 ] ) can then be expressed as @xmath180 , and @xmath164 can be reconstructed by summing over all the modes .
finally , the gravitational waveform @xmath4 and the luminosity @xmath124 can be obtained from the asymptotic behaviour of @xmath164 when @xmath181 .    to integrate eq .
( [ 38 ] ) therefore reduces to solving the homogeneous teukolsky equation , @xmath182 , for the functions @xmath172 and @xmath173 .
this , it turns out , is equivalent to integrating the regge - wheeler equation @xcite @xmath183 \biggr\ } x_{\omega\ell}(r ) = 0 , \label{39}\ ] ] for the functions @xmath184 and @xmath185 ; here , @xmath186 .
this comes about because a solution to the homogeneous teukolsky equation can easily be related to a solution to the regge - wheeler equation .
the relation is known as the chandrasekhar transformation @xcite .
the problem of calculating the gravitational waveform for a compact binary system with small mass ratio can therefore be reduced to the simple one of integrating eq .
( [ 39 ] ) .
the relevant dimensionless parameter entering this equation is @xmath187 , where , for circular orbits , @xmath188 @xcite .
we therefore have , using eq .
( [ 37 ] ) , @xmath189 .    for arbitrary values of @xmath187 the regge - wheeler equation must be integrated numerically @xcite .
this must be done separately for each selected value of @xmath99 in the interval @xmath190 ( approximately corresponding , for systems of a few solar masses , to the frequency interval @xmath191 ) .
@xcite    once again we will focus on @xmath124 , the gravitational - wave luminosity . in fig .  2 we present a plot of @xmath192 as a function of @xmath99 .
this is the ratio of the luminosity as calculated exactly ( numerically ) using the perturbation approach , to the quadrupole - formula expression @xmath193 the numerical results are depicted as a solid curve .
apart from negligible numerical errors , these results are _ exact _ to all orders in @xmath99 ; the only assumption used in the calculation concerns the smallness of @xmath194 .
the various broken curves will be described below .
the figure shows that @xmath192 tends toward unity as @xmath195 , and stays within 15% of unity everywhere in the interval @xmath196 .
it is now easy to judge the accuracy of the post - newtonian expansion for @xmath124 in the limit @xmath197 .
one simply evaluates eq .
( [ 31 ] ) in this limit , with @xmath198 , and compare with the numerical results .
this post - newtonian curve is labeled `` 4 '' in the figure , and we see that the numerical results are only imperfectly reproduced .
the consequences of this will be discussed in the next section .
one can in fact do better than this .
the perturbation approach is not only suitable for exact , numerical computations . by combining it with a slow - motion approximation  putting @xmath199 on top of @xmath200
 it also becomes suitable for approximate , analytical computations .
methods for integrating the regge - wheeler equations analytically in the limit @xmath201 were devised by various authors @xcite .
the basic idea is to proceed by iterations : eq .
( [ 39 ] ) with @xmath202 is solved in terms of spherical bessel functions , and this zeroth order solution is used as input for the first iteration .    using these methods , tagoshi and
sasaki @xcite were able to calculate @xmath124 analytically , accurately through fourth post - newtonian order , two full orders beyond eq .
( [ 31 ] ) .
they derive the rather impressive expression @xmath203 . \label{41}\end{aligned}\ ] ] notice the presence of @xmath204 terms in this expansion , as well as that of the euler number @xmath205 .
notice also that the first four terms reproduce eq .
( [ 31 ] ) in the limit @xmath197 , when @xmath206 .
this , of course , is as it should be .    the broken curves in fig .
2 represent plots of eq .
( [ 41 ] ) truncated to various orders in @xmath99 .
for example , the curve labeled `` 6 '' is a plot of eq .
( [ 41 ] ) with all terms of order @xmath207 and @xmath208 discarded .
we see from the figure that the post - newtonian expansion converges poorly .
( the suspicion , in fact , is that the series is only asymptotic . )
witness in particular the poor quality of the curve `` 5 '' with respect to `` 4 '' , and compare also `` 8 '' to `` 7 '' .
we have seen that in the @xmath197 limit , the post - newtonian expansion for @xmath124 converges poorly .
a similar statement can be made about @xmath128 : the post - newtonian expansion @xmath209 where @xmath210 converges slowly to the exact result @xcite @xmath211 contrary to eq .
( [ 41 ] ) , the expansion ( [ 42 ] ) is actually _ known _ to converge ( for all values of @xmath99 in the interval @xmath196 ) .    in this section
we address the issue as to how much of an obstacle the poor convergence of the post - newtonian expansion poses to the construction of accurate measurement templates @xcite .
we will answer this question for binary systems with small mass ratios , using the results described in the preceding section .
since there is no reason to believe that the convergence of the post - newtonian expansion would be much improved for systems of comparable masses , our results should also apply , at least qualitatively , for such systems .    one way to address
this question is to ask , given a waveform constructed from the exact numerical results , how much signal - to - noise ratio is lost by matched filtering the exact signal with an approximate post - newtonian template ?    from sec .
8 we know that filtering with the exact signal @xmath18 would give @xmath212 , the largest possible value of the signal - to - noise ratio .
on the other hand , filtering with a post - newtonian template @xmath213 , where @xmath10 denotes the order in @xmath99 to which the expansion is taken ( for example , @xmath214 represents a waveform accurate to second post - newtonian order ) , gives the smaller value @xmath215 . from the results of sec .  7 and 8
we obtain @xmath216 the detail of how to compute @xmath217 is presented in ref .
@xcite .
the calculation described here can only be carried out for binary systems with small mass ratios , because only for these do we have access to the exact waveform @xmath18 . nevertheless , in the following we shall let @xmath194 become large , without altering our expressions for the exact and post - newtonian waveforms .
this is done without justification , but reflects the adopted point of view that the quality of the post - newtonian approximation should not be appreciably affected by the finite - mass corrections .    in doing this extrapolation ,
some thought must be given as to the interpretation of the ratio @xmath194 .
in the limit @xmath197 this can be taken to be both the ratio of the individual masses or the ratio of reduced mass to total mass . in the case of comparable masses
, some choice must be made .
we note that as @xmath194 is allowed to grow large , expressions ( [ 40 ] ) and ( [ 43 ] ) for the newtonian quantities @xmath218 and @xmath219 must be replaced by expressions ( [ 29 ] ) , _ which are formally identical_. this shows that when extrapolating to the case of comparable masses , @xmath131 is to be interpreted as the _ reduced mass _ , and @xmath132 as the _ total mass_.    .reduction in signal - to - noise ratio incurred when matched filtering with approximate , post - newtonian templates .
the first column lists the order @xmath10 of the approximation , the second column lists @xmath217 as calculated using the post - newtonian approximation for @xmath128 , and the third column lists @xmath217 as calculated using the exact expression for @xmath128 . [
cols="^,^,^",options="header " , ]     we quote the results for @xmath217 corresponding to a system of two neutron stars , each of 1.4 @xmath220 ; these are displayed in the second column of table 1 .
we see that even at quite a high order in the post - newtonian expansion , only three quarters of the signal - to - noise ratio is reproduced by the post - newtonian template .
this shows that the apparently small discrepancies between the exact and post - newtonian results for @xmath124 and @xmath128 provide a serious obstacle to the construction of accurate measurement templates .
it is interesting to ask how much of the signal - to - noise ratio would be recovered if @xmath128 were kept exact instead of being expressed as a post - newtonian expansion .
the third column of the table displays the values of @xmath217 calculated in this way .
we see that most of the signal - to - noise is recovered : @xmath217 can now be as large as 0.9864 instead of 0.7651 .
why does the exact expression for @xmath128 give such better results ?
it can be established @xcite that this has to do with the following fact : while the exact expression for @xmath128 correctly goes to zero at @xmath221 ( at the innermost circular orbit ) , its post - newtonian analogue fails to do so .
[ for example , the right - hand side of eq .
( [ 42 ] ) , truncated to order @xmath208 , goes to zero at @xmath222 , corresponding to a radius @xmath223 . ]
when corrected for this , the post - newtonian template gives much better results .
we therefore see that the poor convergence of the post - newtonian expansion is a serious obstacle to the construction of accurate measurement templates .
devising ways to extract the information contained in gravitational waves produced during the late inspiral of a compact binary system poses a great challenge to theoretical physicists .
it is not clear that `` simply '' pushing to higher order in post - newtonian theory will be enough to produce sufficiently accurate measurement templates ( such that the systematic errors will be smaller than the statistical errors ) .
there may be a need for theorists to develop alternative ways of dealing with this problem .
it should be stressed that the convergence problem does not arise when constructing _
detection templates_. indeed , as was discussed in sec .
9 , there is no particular need for these templates to be derived from the equations of general relativity . and since the detection templates need not involve any parameters of direct physical significance , the notion of systematic errors does not apply to them .
the only requirement for constructing detection templates is that they should span the appropriate `` signal space '' , and that they should do so the most efficiently .
the poor convergence of the post - newtonian expansion is therefore not an obstacle for _ detecting _ gravitational waves from coalescing compact binaries .
it is only an obstacle for extracting the information that the waves contain . to overcome this obstacle
will undoubtedly be a theorist s challenge for years to come .
this state of affairs is highly interesting from a historical point of view : never before in the history of gravitational physics has experiment demanded such a high degree of sophistication on theoretical calculations .
this work was supported by the natural sciences foundation under grant no .
phy 92 - 22902 and the national aeronautics and space administration under grant no .
nagw 3874 .
this article was completed while visiting the theoretical physics institute of the university of alberta ; the author is most grateful to werner israel for his warm hospitality .
99 for broader reviews on gravitational waves , see k.s .
thorne , in _ 300 years of gravitation _ , edited by s.w . hawking and w. israel ( cambridge university press , cambridge , 1987 ) ; in _ proceedings of the snowmass 95 summer study on particle and nuclear astrophysics and cosmology _ , edited by e.w .
kolb and r. peccei ( world scientific , singapore , 1995 ) .
r.a . hulse and j.h .
taylor , astrophys .
j. * 324 * , 355 ( 1975 ) ; j.h .
taylor , rev .
phys . * 66 * , 711 ( 1994 ) .
a. abramovici _
et al . _ ,
science * 256 * , 325 ( 1992 ) . c. bradaschia _ el al .
instrum . & methods * a289 * , 518 ( 1990 ) .
see , for example , f. reif , _ fundamentals of statistical and thermal physics _
( mcgraw - hill , new - york , 1965 ) .
helstrom , _ statistical theory of signal detection _
( pergamon , oxford , 1968 ) .
schutz , nature ( london ) * 323 * , 310 ( 1986 ) ; class .
quantum grav .
* 6 * , 1761 ( 1989 ) .
other configurations , with narrow - band sensitivity at a chosen frequency , are possible . e.s .
phinney , astrophys .
j. * 380 * , l17 ( 1991 ) ; r. narayan , t. piran , and a. shemi , _ ibid . _ * 379 * , l17 ( 1991 ) ; a.v .
tutukov and l.r .
yungelson , mon . not .
* 260 * , 675 ( 1993 ) . c. kochanek , astrophys .
j. * 398 * , 234 ( 1992 ) ; l. bildsten and c. cutler , _ ibid . _ * 400 * , 175 ( 1992 ) . for an overview ,
see t. damour , in _ 300 years of gravitation _ , edited by s.w .
hawking and w. israel ( cambridge university press , cambridge , 1987 ) .
see , for example , c.w .
misner , k.s .
thorne , and j.a .
wheeler , _ gravitation _ ( freeman , san francisco , 1973 ) , chap .
wainstein and v.d .
zubakov , _ extraction of signals from noise _
( prentice - hall , englewood cliffs , 1962 ) .
the following discussion is based on c. cutler and e.e .
flanagan , phys .
d * 49 * , 2658 ( 1994 ) .
finn , phys .
d * 46 * , 5236 ( 1992 ) .
finn and d.f .
chernoff , phys .
d * 47 * , 2198 ( 1993 ) . c. cutler _
et al . _ ,
70 * , 1984 ( 1993 ) . for an overview , see c.m .
will , in _ relativistic cosmology _ , proceedings of the eighth nishinomiya - yukawa memorial symposium , edited by m. sasaki ( universal academy press , kyoto , 1994 ) . the following discussion is based on the work of blanchet , damour , and iyer , expounded in the following papers : l. blanchet and t. damour , philos .
london a * 320 * , 379 ( 1986 ) ; phys .
d * 37 * , 1410 ( 1988 ) ; ann .
inst . h. poincar ( phys .
thorique ) * 50 * , 377 ( 1989 ) ; l. blanchet , proc .
london a * 409 * , 383 ( 1987 ) ; t. damour and b.r .
iyer , ann .
h. poincar ( phys
. thorique ) * 54 * , 115 ( 1991 ) . k.s .
thorne , rev .
phys . * 52 * , 299 ( 1980 ) .
kidder , c.m .
will , and a.g .
wiseman , phys .
d * 47 * , r4183 ( 1993 ) .
peters and j. mathews , phys . rev . *
131 * , 435 ( 1963 ) .
r.v . wagoner and c.m .
will , astrophys .
j. * 210 * , 764 ( 1976 ) ; * 215 * , 984 ( 1977 ) .
e. poisson , phys .
d * 47 * , 1497 ( 1993 ) .
a.g . wiseman , phys .
d * 48 * , 4757 ( 1993 ) .
l. blanchet and g. shfer , class .
quantum grav .
* 10 * , 2699 ( 1993 ) . l. blanchet , t. damour , b.r .
iyer , c.m .
will , and a.g .
wiseman , phys .
* 74 * , 3515 ( 1995 ) .
the following discussion is based on the work of this author and his collaborators .
see e. poisson and m. sasaki , phys .
d * 51 * , 5753 ( 1995 ) and references therein .
teukolsky , astrophys .
j. * 185 * , 635 ( 1973 ) .
t. regge and j.a .
wheeler , phys .
rev . * 108 * , 1063 ( 1957 ) .
s. chandrasekhar , proc .
london a * 343 * , 289 ( 1975 ) . c. cutler , l.s .
finn , e. poisson , and g.j .
sussman , phys .
d * 47 * , 1511 ( 1993 ) .
m. sasaki , prog .
. phys . * 92 * , 17 ( 1994 ) .
h. tagoshi and m. sasaki , prog .
* 92 * , 745 ( 1994 ) .
this can be derived simply from the equations for circular geodesic motion in schwarzschild .
e. poisson , _ gravitational radiation from a particle in circular orbit around a black hole .
accuracy of the post - newtonian expansion _
d , in press .


recent works have explored the possibility to show the effects of the gravity onto quantum systems @xcite .
an interesting study concerning the behaviour of a scalar field in the neighbourhood of a black hole @xcite has shown that the entanglement is degraded by the effect of the black hole . in this work ,
i. fuentes - schuller and r.b .
mann have studied a model in which all the ingredients of quantum field theory are present , but the gravitation is only represented as a rindler spacetime corresponding to the uniform surface gravity in the neighbourhood of the event horizon . recently
palmer _ etal _ @xcite have proposed a theory of localized qubits in curved spacetimes . in this paper , we want to reexamine in this framework the problem of qubits around a black hole , especially concerning the entanglement and the quantum teleportation protocol
. moreover we want to analyse the physical meaning of the non - self - adjointness of the localized qubit hamiltonian .
in contrast with the model of fuentes - schuller and mann , quantum field theory is not completely treated , semi - classical approximations in the localized qubit theory induce the lost of the possibility to create and annihilate particles , the lost of the unruh effect , and the lost of the delocalization of the wave packets ( but the non - locality remains in the theory with the entanglement ) .
but we want to treat the complete geometry of the black hole spacetime .
the interest of our approach is the possibility to treat the local effects on a fixed quantum system ( a single qubit boarded in a `` spacecraft '' with definite position and velocity and following a geodesic ) whereas no local information associated with the spacetime geometry is taken into account in the fuentes - schuller mann model since the gravitational field is considered in it as uniform and the qubits are completely delocalized in it .
we want to analyse the entanglement degrading effect with respect to the position and the velocity of the qubits with respect to the black hole . + we consider a qubit realised as the spin of a fermion submitting only to the gravitational field ( no external magnetic field ) and boarded in a `` spacecraft '' following a geodesic around a black hole . by `` gravitational field '' we mean the lorentz connection associated with the spacetime geometry . in order to treat the dynamics of the qubit
, we use the quantum adiabatic approximation because the qubit transport can be considered as slow with respect to its proper quantum time ( the period of its rabi oscillations induced by the gravitational field ) .
section 2 summarizes the localized qubit approach by rewritting it in the language of the fibre bundle theory .
the goal of this reformulation is to provide a description including spacetime geometries and qubit quantum states in a single common geometric structure .
this is achieved by the introduction of the fiber bundles of the quantum adiabatic approximation .
we show that the problem takes place in complex line bundles over a space constituting by the product of the space of lorentz connections by the tangent bundle of the spacetime manifold .
we show that the qubit is submitted to a kind of decoherence process induced by the gravitational field and responsible to the degradation of the entanglement .
section 3 presents in our framework the quantum teleportation protocol with an epr ( einstein podolsky rosen ) qubit pair when one qubit falls to the black hole whereas the other one is comoving with it .
we compute a formula providing the fidelity of the teleportation .
section 4 applies the formalism to two spacetime geometries , firstly to the rindler spacetime used by fuentes - schuller and mann , secondly to a schwarzschild spacetime , where we analyse the fidelity of the quantum teleportation protocol with respect to the geodesic followed by the qubit falling to the black hole .
+ _ throughout this paper , we consider the unit system such that @xmath0 .
_ + _ note about the notations : a fibre bundle of total space @xmath1 and base space @xmath2 is denoted by its projection @xmath3 .
the space of the local sections of a fibre bundle is denoted by @xmath4 . for a map @xmath5
, @xmath6 denotes the map induced by the fibration .
let @xmath7 be a manifold , its tangent bundle is denoted by @xmath8 ( @xmath9 is the tangent space of @xmath7 at the point @xmath10 ) , and its differential @xmath11-form set is denoted by @xmath12 . for a map @xmath13 , @xmath14 and @xmath15 denote the associated tangent and the cotangent maps ( the pusch - forward and the pull - back ) .
@xmath16 and @xmath17 denote the spaces of continuous and differentiable functions of @xmath7 .  @xmath18  between two manifolds denotes a diffeomorphism .
@xmath19 denotes the projection map defined by @xmath20 .
we use the einstein convention concerning the up and down indices repetition .
the greek indices runs in @xmath21 as curved spacetime indices , the capital latin indices runs in @xmath21 as flat minkowski auxiliary spacetime indices , the small latin indices runs in @xmath22 .
@xmath23 denotes the derivation with respect to a proper time @xmath24 . _
in this section we summarize the results ( without details ) of palmer _ etal _
@xcite , we fix the notations , and we embbed the localized qubit theory into the fibre bundle theory ( for an exposition of the fibre bundle theory see for example @xcite ) .
let @xmath7 be an open set of the spacetime manifold endowed with a local coordinates system @xmath25 and a metric tensor @xmath26 ( to simply we will refere to @xmath7 as the spacetime ) .
let @xmath27 be a tetrad field associated with the metric : @xmath28 , where @xmath29 is the minkowski metric .
let @xmath30 be the lorentz connection ( @xmath31 being the christoffel symbols ) .
a dirac field @xmath32 obeys to the dirac - einstein equation : @xmath33 where @xmath34 are the dirac matrices ( in weyl representation ) and @xmath35 is the spinorial covariant derivative defined by @xmath36 where @xmath37 is the @xmath38 representation of @xmath39 ( the coverging group of the lorentz group @xmath40 ) on @xmath41 ( we denote by the same symbol the induced representation of its lie algebra ) , i.e. @xmath42 $ ] .
+ let @xmath43 be the tangent bundle of @xmath7 and @xmath44 be the frame principal @xmath40-bundle of @xmath7 .
let @xmath45 and @xmath46 be the local trivializations of these bundles , which are defined by @xmath47 and @xmath48 ( @xmath49 is the matrix of elements @xmath50 ) .
the tetrad field can be identified as the trivializing local section of @xmath51 : @xmath52 .
let @xmath53 be the principal @xmath39-bundle associated with the local @xmath39 transformations of the spinors ( @xmath54 is an extension of @xmath51 such that @xmath55 ) . let @xmath56 and @xmath57 be the associated vector bundles for the representation @xmath58 and @xmath59 , i.e. @xmath56 is defined by its local trivialization @xmath60 , @xmath61_{g \in sl(2,\mathbb c)}$ ] where @xmath62 is the local trivialization of @xmath63 and @xmath64 is the projection @xmath43 , and @xmath65_{g \in sl(2,\mathbb c)}$ ] ( we have denoted simply the @xmath58-action of @xmath66 on @xmath67 by @xmath68 and the @xmath59-action of @xmath69 on @xmath70 by @xmath71 ) .
+ @xmath72 is a hilbert @xmath73-module endowed with the inner product : @xmath74 where @xmath75 ( @xmath76 , @xmath77 ) .
let @xmath78 be a spacelike hypersurface of @xmath7 and @xmath79 be the set of futur oriented timelike normal vectors to @xmath80 .
the dirac field @xmath32 is a vector of the hilbert space @xmath81 ( the space of the dirac spinor fields living at some time on the instantaneous space @xmath80 ) .
+ let @xmath82 be a geodesic worldline in @xmath7 and @xmath83 be a fiolation of @xmath7 along @xmath82 by spacelike hypersurfaces ( @xmath24 being the proper time along @xmath82 ) . by the wkb ( wentzel kramers brillouin ) approximation associated with the assumption that the compton wavelength is small with respect to the curvature scale @xcite , we have @xmath84 .
this semi - classical approximation suppresses the space delocalisation of the fermion supporting the qubit and the absence of second quantization supresses the particle number ambiguity .
since only the @xmath58-representation is needed to describe a single spin ( qubit ) , we can project onto the space @xmath85 .
we work then with the composite bundle @xcite @xmath86 .
it can be more easy to work with a bundle @xmath87 with structure fibre @xmath88 ( tangent vector model space and spin quantum state space ) .
@xmath89 is defined by the local trivialization @xmath90 , @xmath91 .
we have an action of @xmath39 onto @xmath89 defined by @xmath92 , @xmath93 , where @xmath94 is the group homomorphism associated with the quotient @xmath95 .
+ we can identify the space of local sections of @xmath89 , @xmath96 , to the space of @xmath40-invariant local sections of @xmath97 : @xmath98 .
the restriction of @xmath99 to the invariant sections is important to ensure the physical character of the theory , more precisely to have the following property : @xmath100 i.e. the quantum propabilities are invariant under lorentz transformations .
endowed with the inner product @xmath101 ( with @xmath102 , @xmath103 is the dirac distribution ) , @xmath96 is a hilbert @xmath16-module .
+ at this stage , the localized qubit is described by the spin state @xmath104 which obeys to the schrdinger like equation ( see @xcite ) : @xmath105 where @xmath24 is the proper time along the geodesic followed by the qubit , and @xmath106 ( with @xmath107 and @xmath108 , @xmath109 being the usual pauli matrices ) .
@xmath110 where @xmath111 is the projection of @xmath97 onto @xmath8 and @xmath112 .
the instantaneous qubit hilbert space @xmath113 depends on the four - velocity since the correct inner product induced by the dirac field theory is @xmath114 where the conjugate state is @xmath115 .      in order to interpret the dependence of the instananeous spin hilbert space on the four - velocity , it needs to recall some elementary axioms of the quantum mechanics .
an hilbert space @xmath116 constitutes the space of the states of a quantum system .
but its algebraic dual @xmath117 , i.e. the space of continuous linear functionals of @xmath116 , is the space of the probability amplitudes of the elementary events : @xmath118 is a map from @xmath116 to @xmath119 , such that @xmath120 is the probability of the realization of some measurement event associated with @xmath121 when the quantum system is in the state @xmath70 . by the riesz theorem , we know that @xmath122 , @xmath123 ( up to a renormalization and phase factor ) such that @xmath124 , @xmath125
. we find @xmath126 as an eigenvector of an observable @xmath127 associated with the measure ; @xmath128 corresponds to the event `` the measure of @xmath127 has provided the result @xmath129 '' ( @xmath129 being the eigenvalue associated with the eigenvector @xmath128 ) .
it is then important to note that a `` ket '' @xmath130 characterizes the quantum system as being its state , whereas a `` bra '' @xmath131 characterizes an event for an observer making measures on the quantum system .
+ returning to the localized qubit problem , we want to interpret the difference between the two linear functionals @xmath132 and @xmath133 ( the `` bra '' always denoting in this paper the partial inner product of @xmath134 : @xmath135 , and never @xmath136 ) .
let @xmath137 be a normalized state : @xmath138 where we have introduced the magnetic four - momentum operator @xmath139 , with @xmath140 ; @xmath141 , @xmath142 and @xmath143 .
the formula [ spinchangeref ] is the classical relation between a magnetic four - momentum @xmath144 measured into an inertial frame @xmath145 and the four - momentum @xmath146 measured into its rest frame @xmath147 of four - velocity @xmath148 with respect to @xmath145 ( see for example @xcite ) . for our problem
, @xmath145 is a frame comoving with the black hole and @xmath147 is the frame comoving with the qubit .
+ with this analysis of the normalization with respect to @xmath149 we can say that :    * @xmath132 is the linear functional associated with the ( non - normalized ) probability amplitude to find the spin in the state @xmath70 measured by an observer comoving with the black hole ( @xmath150 ) .
* @xmath151 is the linear functional associated with the ( normalized ) probability amplitude to find the spin in the state @xmath70 measured by an observer comoving with the qubit ( @xmath152 ) .
( with @xmath153 ) .
we consider the schrdinger like equation for the localized qubit : @xmath154 with the hamiltonian @xmath155 ( non - self - adjoint with respect to @xmath156 , @xmath157 ) : @xmath158 @xmath159 ( @xmath160 ) .
we have @xmath161 ( dissipation operator ) and @xmath162 ( hamiltonian of the qubit rotation )
. finally we can write @xmath163 where @xmath164 ( @xmath165 is the complex self - dual lorentz connection ) . from the viewpoint of the qubit ,
the interaction with the gravitational field is similar to a spin submitted to a complexified magnetic field .
let @xmath166 . + to integrate the dynamics involved by the schrdinger like equation , we propose to use the adiabatic approximation for the non - self - adjoint hamiltonians @xcite : @xmath167 where @xmath168 , @xmath169 and @xmath170 are respectively the instantaneous right eigenvectors , left eigenvectors and eigenvalues of @xmath171 : @xmath172 ( the overline denoting the complex conjugation ) , @xmath173 , and @xmath174 are the generators of the non - unitary geometric phases : @xmath175 @xmath176 is the path in @xmath177 defined by @xmath178 for the geodesic @xmath179 followed by the qubit . + a simple calculation shows that @xmath180 @xmath181 @xmath182 the adiabatic approximation is valid if the non - adiabatic coupling is negligible , i.e. @xmath183    let @xmath184 be the space of the @xmath39-connections of the principal bundle @xmath54 .
the eigenvectors can be considered as maps @xmath185 such that @xmath186 , where @xmath187 is the inner product of @xmath7 and @xmath188 is defined by @xmath189 ( @xmath190 constitutes a set of generators of @xmath191 the lie algebra of @xmath39 ) .
the eigenvectors being defined up to an arbitrary normalization and phase factor , they define @xmath119-line bundles @xmath192 with local trivializations @xmath193 with @xmath194 .
@xmath70 obtained by the adiabatic transport formula [ adiabtransp ] is then a local section of @xmath195 over @xmath196 .
note that the left eigenvectors do not define line bundles since their normalization factors are fixed by those of the right eigenvectors .
+ finally the geometric structure in which the qubit transport takes place can be summarized by the following commutative diagram : @xmath197 where @xmath198 .
we can note that @xmath199 is a connection of the principal bundle @xmath54 and @xmath200 are connections of the @xmath201-principal bundles associated with @xmath202 .
we have then three kinds of gauge changes associated with each floor of the composite bundle :    * ground floor : @xmath203 ( diffeomorphism of the spacetime manifold ) , @xmath204 and @xmath205 . *
first floor : @xmath206 ( local lorentz transformation ) , @xmath207 and @xmath208 . * second floor : @xmath209 ( normalization and phase local change ) , @xmath210 .
the different steps of the construction of the localized qubit adiabatic state can be summarized as follows : @xmath211    it can be interesting to note that the holonomy of @xmath199 along @xmath82 ( between @xmath212 and @xmath24 ) is @xmath213 where @xmath214 and @xmath215 denote path and time ordered exponentials ( dyson series ) . @xmath216 which characterize the adiabatic state of the qubit can be viewed as cylindrical functions of the space of lorentz connections , @xmath217 , and @xmath70 as a linear combination of these two cylindrical functions .
it can be interesting to note that ( the topological completion of ) @xmath218 constitutes the kinematical hilbert space of the loop quantum gravity @xcite , this could be indicate a possible connection of the adiabatic localized qubit formalism with a semi - classical limit of the loop quantum gravity .
+ note that the evolution governed by @xmath219 is unitary with respect to @xmath149 ( for a proper observer comoving with the qubit ) ( see @xcite ) : @xmath220 .
but it is not unitary with respect to @xmath221 ( for an observer comoving with the black hole ) .
we begin to examine this point in the next section .
+ _ to simplify the notation , from this point we denote @xmath222 only by @xmath223 . _
the adiabatic transport of a non - self - adjoint two - level quantum system has been extensively studied in the litterature ( see for example @xcite ) .
the interesting effects in the adiabatic transport ( eq . [ adiabtransp ] ) are related to the submanifold @xmath224 of @xmath177 defined by the crossings @xmath225 .
firstly because the validity of the adiabatic approximation ( eq . [ valad ] ) needs to do not approach @xmath224 ( except if @xmath226 ) . secondly because @xmath224 is a kind of hybercone separating the region of @xmath177 for which @xmath227 are real ( and generate only pure phases ) from the region for which @xmath227 are complex ( and generate non - unitary ( for @xmath156 ) evolution modifying the relative weights of the superposition of @xmath228 ) .
since in the self - adjoint case , the geometric phase generator is similar to a magnetic field induced by a magnetic monopole at the eigenvalue crossing point ( see @xcite ) , for the non - unitary case , @xmath224 has been called complex magnetic monopole @xcite ( but @xmath224 is not an isolated point and is associated with exceptional crossings , i.e. @xmath219 is not diagonalizable on @xmath224 )
. let @xmath229 and @xmath230 . @xmath231 . since @xmath232 ,
the complex magnetic monopole is defined by @xmath233 @xmath234 . if the condition @xmath235 is satisfied , outside @xmath224 ( @xmath236 ) , @xmath237 and @xmath238 are just pure phases ; but inside @xmath224 ( @xmath239 ) , @xmath240 and @xmath241 are non - unitary dynamical phases . in this last case , the evolution modifies the weights of the superposition of @xmath228 ( with respect to @xmath221 , i.e. for an observer comoving with the black hole ) .
we will call this effect a * dynamical decoherence * , because the following coherence @xmath242 falls to zero for large @xmath24 ( with @xmath243 , we have neglected the geometric phases and supposed that @xmath244 ) .
we have not considered the effects of the non - unitary geometric phases @xmath245 which can induce a * geometric decoherence * if @xmath246 ( we call it geometric decoherence since the geometric phase depends only on the shape of the followed path @xmath176 and not from the proper time ) .
+ let @xmath247 be the preimage of @xmath224 into @xmath8 .
it is important to note that the complex magnetic monopole for a fixed spacetime geometry @xmath248 is not a submanifold of the spacetime @xmath7 but a submanifold of the tangent bundle @xmath8 .
the complex magnetic monopole around the black hole `` viewed '' by the qubit depends on its four - velocity .
gobally the set of all complex magnetic monopoles is @xmath249 . in some cases
, a class @xmath250 of geodesics can be defined with some first integrals @xmath251 and @xmath252 ( @xmath253 , @xmath254 ) depends only on @xmath251 and @xmath10 . in that case @xmath255
is a submanifold of @xmath7 which is an image ( for the qubits following geodesics of @xmath250 ) of the complex magnetic monopole in the spacetime .      to understand the physical origin of the non - unitarity with respect to @xmath222 ( observer comoving with the black hole ) , consider first a more simple model constituted by a three - level system ( with levels denoted by @xmath256 ) , governed by an hamiltonian @xmath171 and with spontaneous emission from @xmath257 to the `` dark state '' @xmath258 with a rate @xmath259 .
we want to consider the system restricted to @xmath260 ) as a qubit and to forget the dark state @xmath258 .
the system obeys to a master equation @xcite : @xmath261 -\frac{\gamma_-}{2 } \{\sigma_{d0}^+ \sigma_{d0}^- , \rho\ } + \gamma_- \sigma^-_{d0 } \rho \sigma^+_{d0}\ ] ] where @xmath262 is the density matrix of the system , @xmath263 denotes the anti - commutator , @xmath264 and @xmath265 .
the equation can be rewritten as @xmath266 where @xmath267 .
the anti - self - adjoint part of the effective hamiltonian @xmath268 models the lost of population from @xmath257 to the dark state by spontaneous emission , whereas @xmath269 models the gain of population of this dark state .
so , if we forget the dark state in the modelization , the qubit obeys to a schrdinger equation governed by a non - self - adjoint effective hamiltonian @xmath270 .
+ in a curved spacetime there is an ambiguity concerning the particle number . due to the unruh effect
, the vaccum in the rest frame becomes a thermal state in an uniformly accelerated frame @xcite . at the level of the quantum field theory in curved spacetime , the evolution in the black hole frame of the dirac field spontaneously couples the one particle state to the zero particle state ( in the fermionic fock space of the system ) .
but the semi - classical and wkb approximations used in our model forgets this last one ( we want to have one and only one qubit ) .
we have then only two qubit states @xmath271 and @xmath272 ( forming the canonical basis of @xmath134 used in the construction of the different bundles ) , and a dark state : the vacuum @xmath273 . in the same way that for the small example of a three - level system , the qubit is then governed by a non - self - adjoint effective hamiltonian .
we will study this point with more details in a concrete example in section 4 .
let alice and bob be initially at the point @xmath274 of @xmath7 , supposed sufficiently far from the black hole to consider that @xmath7 is flat in the neighbourhood of @xmath274 .
we set @xmath275 and @xmath276 .
bob is supposed comoving with the black hole ( he stays at @xmath277 ) , but alice follows a geodesic going near the event horizon at a point @xmath278 .
alice wants to teleport information when she will be at this point . at the moment @xmath279 when alice leaves bob
, they have an entangled qubit pair in a bell state : @xmath280 @xmath281 and @xmath282 are the initial four - velocities of alice and bob . since alice and bob belongs to two different frames , each one has its proper definition of the qubit states . for bob ,
the linear functionals of finding its qubit in a particular state are @xmath283 and @xmath284 , involving that @xmath285 and @xmath286 .
but for alice , her linear functionals are @xmath287 and @xmath288 since she is not comoving with the black hole .
the qubit states for alice are then defined by @xmath289 ( @xmath290 ) .
it follows that @xmath291 and @xmath292 ( @xmath293 ) .
@xmath294 in the flat region ( @xmath295 ) , we have @xmath296 and @xmath297 .
@xmath298 let @xmath299 the proper time when alice arrives at @xmath278 .
we suppose that the evolution along the geodesic @xmath82 linking @xmath277 to @xmath278 is adiabatic for the alice s qubit ( as eq .
[ adiabtransp ] ) .
we have then for @xmath300 and @xmath301 @xmath302 ( @xmath303 ) .
note that @xmath304 is defined for two proper times , one for alice and one for bob , since their clocks are asynchronous .
the evolution for the bob s qubit is trivial since it is inertial in a flat part of @xmath7 .
@xmath305 where @xmath306 is the alice s four - velocity at @xmath299 .
@xmath307 with @xmath308 alice encodes a quantum information in a qubit @xmath309 ( @xmath310 ) .
the state of the three qubits is then @xmath311 .
alice performs then the operations of the usual teleportation protocol : @xmath312 where @xmath313 and @xmath314 are the cnot and hadamard gates in the alice s frame . after some algebra ,
we find @xmath315 alice performs a measurement of her qubits . to fix the discussion
, we suppose that she finds @xmath316 ( the result can be easily adapted for another result ) .
alice sends to bob by a classical communication chanel what is the operation to perform on his qubit ( in our example , the operation is the identity ) .
bob receives the message at @xmath317 .
the state is then for @xmath318 and @xmath319 : @xmath320 where @xmath321 is the evolution operator for an alice s qubit after @xmath299 .
the fidelity of the quantum teleportation is then @xmath322 the fidelity of the teleportation is clearly degraded by the decoherence induced by the black hole which is encoded in @xmath323 .
+ remark : for a flat spacetime with alice having a constant four - velocity we have @xmath324 , and then @xmath325 ( we refind the efficiency of the usual teleporation protocol ) .
in order to compare with the fuentes - schuller mann model @xcite , we first consider the case of the rindler spacetime defined by the metric @xmath326 which corresponds to a flat spacetime viewed in a noninterial frame uniformly accelerated ( with acceleration parameter @xmath327 ) , or to the surface gravity approximation of a schwarzschid black hole ( @xmath328 and @xmath329 , @xmath330 being the schwarzschild radius ) ( see @xcite ) .
the tetrad fields are @xmath331 and @xmath332 , and the only one non - zero lorentz connection component is @xmath333 .
it follows that @xmath334 and @xmath335 , and then @xmath336 with @xmath337 and @xmath338 @xmath339 .
it is interesting to note that since the qubit moves in the @xmath10-direction , @xmath340 corresponds to a spin parallel to the linear momentum ( so to a positive helicity state ) and @xmath341 to a spin antiparallel to the linear momentum ( a negative helicity state ) ( @xmath342 ) .
the geodesic equations are @xmath343 the first geodesic equation defines the first integral : @xmath344 the second geodesic equation becomes the autonomous equation : @xmath345 the adiabatic transport of a qubit state @xmath346 ( @xmath347 ) is then @xmath348 the dynamical decoherence kills the positive helicity state in favor of the negative helicity state .
we can heuristically understand this fact as follows .
the vacuum of the minkowsky spacetime becomes in the noninertial frame @xmath349 ( see @xcite ) where @xmath350 ( @xmath351 denotes the momentum , @xmath352 the helicity , @xmath353 and @xmath354 denote the two regions separated by the horizon , @xmath355 ) .
it is then associated with a density matrix @xmath356 which is a thermal distribution with temperature @xmath357 ( @xmath358 is the boltzmann constant ) , corresponding to the unruh radiation .
it follows that the fermion is coupled with this thermal bath . for the positive helicity mode ( which is the part of the weyl spinor with positive energy )
, we can write that the density matrix @xmath359 ( for @xmath360 and @xmath273 ) obeys to the master equation ( see for example @xcite ) : @xmath361 where @xmath362 are the fermionic creation / annihilation operators on the positive helicity mode ( @xmath363 and @xmath364 ) , @xmath365 and @xmath366 characterizes the spectral density of the bath . with some assumptions @xmath367 where @xmath368 is a constant ( see @xcite ) .
since @xmath369 is very small ( the unruh temperature @xmath370 is very small ) , the master equation is dominated by the dissipation of @xmath360 in accordance with eq .
[ transprindler ] . by following
the approximation explained section [ origin ] ( by projecting onto @xmath371 , by neglecting the quantum jumps and the @xmath369 terms ) we have @xmath372 we refind @xmath373 by setting @xmath374 .
for the negative helicity mode the problem is quite different since it is associated with the part of the weyl spinor with negative energy .
it follows that the roles of @xmath375 and @xmath376 are inverted in the master equation .
this one is then dominated by the increase of the population of @xmath377 .
we can then postulate an effective hamiltonian creating negative helicity population , as @xmath378 .
+ the geodesic equation ( eq . [ geodesicrindler ] )
has for solution @xmath379 with @xmath380 .
we have then an analytical expression of the dynamical phases : @xmath381 let @xmath382 be the proper time for which the qubit reaches the horizon .
we have @xmath383 .
+ the fidelity of the teleportation protocol for the rindler spacetime is represented fig .
[ fidelity_rindler ] .
average fidelity @xmath384 ( eq . [ fidelity ] ) of the teleportation protocol , for alice following a geodesic of the rindler spacetime and bob being static , with respect to @xmath385 ( first integral of the geodesic ) and @xmath386 the alice s proper time when she realizes its part of the protocol . for each value of @xmath385
, the fidelity is drawn until the proper time when alice reaches the rindler horizon.,width=302 ]    the fidelity of the quantum teleportation falls if alice approaches too close to the rindler horizon .
we consider the metric associated with a static black hole with spherical symmetry : @xmath387 where @xmath370 and @xmath388 are the factors of time dilation and length contraction , the schwarzschild metric being obtained for @xmath389 with @xmath390 the schwarzschild radius .
the tetrad fields are @xmath391 , @xmath392 , @xmath393 and @xmath394 ; and the non - zero components of the lorentz connection are @xmath395 , @xmath396 , @xmath397 and @xmath398 . because of the spherical symmetry , we can restric our attention to the equatorial plane @xmath399 .
the geodesic equations are @xmath400 the first and the last geodesic equations define the first integrals : @xmath401 @xmath97 and @xmath402 being the energy and the angular momentum by mass unit .
we have @xmath403 , @xmath404 , and @xmath405 .
it follows that @xmath406 with @xmath407 .
@xmath408 if @xmath409 ( @xmath410 \text { with } \frac{t'(r_{le } ) r^2_{le}}{t(r_{le})^2 } = \frac{l}{e}\right\}$ ] ) . for the schwarzschild case , the dynamical decoherence disapears for @xmath411 , i.e. if @xmath412 and @xmath413 .
it follows that the qubit is submitted to dynamical decoherence except if it follows a strongly rotating geodesic ( @xmath402 large ) and far from the complex magnetic monopole ( which is a sphere of radius @xmath414 )
. the generators of the geometric phases are @xmath415 ( with @xmath416 , @xmath417 and @xmath418 ) . for the schwarzschild case
, we have @xmath419 .
geometric decoherence is always present except for the radial geodesics ( @xmath420 ) and the circular orbits ( @xmath421 constant ) .
moreover the non - adiabatic coupling is @xmath422 , assuring without any assumption concerning the velocity , the validity of the adiabatic approximation for the radial geodesics ( @xmath420 ) and the circular orbits ( @xmath423 ) .
we consider first the radial geodesics @xmath420 ( @xmath424 ) .
the second geodesic equation ( eq . [ geodesicschwarz ] ) is then reduced to @xmath425 , which has for solution : @xmath426 .
the event horizon is reached at @xmath427 .
the fidelity of the teleportation protocol for this situation is drawn figure [ fidelity_schwarz_radial ] .
average fidelity @xmath428 ( eq . [ fidelity ] ) of the teleportation protocol , for alice following a radial geodesic of the schwarzschild spacetime and bob comoving with the black hole , with respect to @xmath97 ( first integral of the geodesic ) and @xmath386 the alice s propre time when she realizes its part of the protocol .
the fidelity is drawn until the proper time when alice reaches the event horizon.,width=302 ]    as for the rindler spacetime , the fidelity of the quantum teleportation falls if alice approaches too close to the event horizon , because of the decoherence induced by the gravitational field ( all the radial geodesics are inside @xmath429\}$ ] ) .
we consider the circular orbits defined by @xmath430 ( constant ) and @xmath431 .
the second geodesic equation ( eq . [ geodesicschwarz ] ) involves that @xmath432 and the metric ( eq . [ metricschwarz ] ) involves that @xmath433 .
the fidelity of the teleportation protocol for this situation is drawn figure [ fidelity_schwarz_circular ] .
average fidelity @xmath384 ( eq . [ fidelity ] ) of the teleportation protocol , for alice following a circular orbit around the black hole and bob comoving with the black hole , with respect to @xmath434 ( radius of the orbit ) and @xmath386 the alice s propre time when she realizes its part of the protocol .
the red lines indicates the proper times corresponding to the orbital periods .
we start at @xmath435 ( the photon sphere ) since no closed orbit exists under this value.,width=302 ]    for the circular orbits , the effect is essentially caused by the difference of four - velocity between alice and bob ( explaining why the fidelity is almost uniform with respect to @xmath434 and @xmath299 ) .
no decoherence occurs since for all @xmath434 , @xmath237 ( @xmath436\}$ ] , all circular orbits are outside the complex magnetic monopole which is identified with the photon sphere ) .
the adiabatic transport generates a phase difference between @xmath340 and @xmath341 which induces some interferences in the quantum teleportation explaining the small oscillations in the fidelity .
we consider geodesics starting far from the event horizon and almost reaching it by an adiabatic process for the qubit evolution .
since @xmath437 we can suppose that @xmath438 ( @xmath439 being the proper time needed to reach the event horizon and @xmath440 ) .
moreover @xmath441 , @xmath442 , and since @xmath443 , we can compute an evaluation of the fidelity of the teleportation protocol for geodesics almost reaching the event horizon , see figure [ fidelity_schwarz_he ] .
fidelity @xmath444 ( eq . [ fidelity ] ) of the epr teleportation protocol , for alice following a geodesic almost reaching the event horizon and bob comoving with the black hole , with respect to @xmath97 and @xmath402 ( first integrals of the geodesic ) and @xmath445 the angular position of alice when she almost reaches the event horizon ( alice being started from @xmath446 ) , its final radial position @xmath447 being such that @xmath448 .
alice realizes its part of the protocol when she reaches its final point near the event horizon.,title="fig:",width=264]fidelity @xmath444 ( eq . [ fidelity ] ) of the epr teleportation protocol , for alice following a geodesic almost reaching the event horizon and bob comoving with the black hole , with respect to @xmath97 and @xmath402 ( first integrals of the geodesic ) and @xmath445 the angular position of alice when she almost reaches the event horizon ( alice being started from @xmath446 ) , its final radial position @xmath447 being such that @xmath448 .
alice realizes its part of the protocol when she reaches its final point near the event horizon.,title="fig:",width=264 ] + fidelity @xmath444 ( eq . [ fidelity ] ) of the epr teleportation protocol , for alice following a geodesic almost reaching the event horizon and bob comoving with the black hole , with respect to @xmath97 and @xmath402 ( first integrals of the geodesic ) and @xmath445 the angular position of alice when she almost reaches the event horizon ( alice being started from @xmath446 ) , its final radial position @xmath447 being such that @xmath448 .
alice realizes its part of the protocol when she reaches its final point near the event horizon.,title="fig:",width=264]fidelity @xmath444 ( eq . [ fidelity ] ) of the epr teleportation protocol , for alice following a geodesic almost reaching the event horizon and bob comoving with the black hole , with respect to @xmath97 and @xmath402 ( first integrals of the geodesic ) and @xmath445 the angular position of alice when she almost reaches the event horizon ( alice being started from @xmath446 ) , its final radial position @xmath447 being such that @xmath448 .
alice realizes its part of the protocol when she reaches its final point near the event horizon.,title="fig:",width=264 ] + fidelity @xmath444 ( eq . [ fidelity ] ) of the epr teleportation protocol , for alice following a geodesic almost reaching the event horizon and bob comoving with the black hole , with respect to @xmath97 and @xmath402 ( first integrals of the geodesic ) and @xmath445 the angular position of alice when she almost reaches the event horizon ( alice being started from @xmath446 ) , its final radial position @xmath447 being such that @xmath448 .
alice realizes its part of the protocol when she reaches its final point near the event horizon.,title="fig:",width=264]fidelity @xmath444 ( eq . [ fidelity ] ) of the epr teleportation protocol , for alice following a geodesic almost reaching the event horizon and bob comoving with the black hole , with respect to @xmath97 and @xmath402 ( first integrals of the geodesic ) and @xmath445 the angular position of alice when she almost reaches the event horizon ( alice being started from @xmath446 ) , its final radial position @xmath447 being such that @xmath448 .
alice realizes its part of the protocol when she reaches its final point near the event horizon.,title="fig:",width=264 ] +    the fidelity oscillates with the relative angular position of alice when she reaches the event horizon ( for the teleportation of a `` schrdinger cat '' , no oscillation occurs for the teleportation of @xmath257 or @xmath449 ) .
the dependence from @xmath97 and @xmath402 is small except for the small values of these first integrals .
a localized qubit in general relativity is described ( at the adiabatic limit ) by a geometric structure including the description of the quantum states and of the spacetime geometry , i.e. the composite bundle @xmath450 where @xmath184 is the space of lorentz connections . in this bundle , inside a particular submanifold ( the complex magnetic monopole @xmath451 )
decoherence processes appear on the qubit .
this complex magnetic monopole is a sphere surrounding the event horizon for a schwarzschild black hole , with radius decreasing with the increase of the angular momentum of the qubit ( it is infinite for @xmath452 , is equal to @xmath453 ( the photon sphere ) for the circular orbits , and tends to @xmath454 ( the event horizon ) with @xmath455 ) .
we have two different decoherence processes , a dynamical decoherence associated with the non - unitary dynamical phases and a geometric decoherence associated with the non - unitary geometric phases ( and depending only on the shape of the followed path and not from the proper time ) .
the physical origin of these processes is related in a rindler spacetime to the unruh radiation and we can then postulate that in the general case it is related to the hawking radiation ( since the unruh effect can be considered as the near - horizon form of the hawking radiation ) .
we have shown how these decoherence processes degrade the fidelity of the quantum teleportation protocol if alice falls to the event horizon , the adiabatic framework permitting to obtain a simple formula to compute this fidelity with respect to the spacetime position and to the four - velocity of alice when she realizes her part of the protocol .
0 ahmadi f and mehrafarin m 2012 _ j. phys .
a _ * 45 * 485302 .
brushi d e , sabin c , white a , baccetti v , oi d k l and fuentes i 2014 _ new j. phys .
_ * 16 * 053041 .
fuentes - schuller i and mann r b 2005 _ phys .
lett . _ * 95 * 120404 .
palmer m c , takahashi m and westman h f 2012 _ ann .
phys . _ * 327 * 1078 .
nakahara m 1990 _ geometry , topology and physics _
( bristol : iop publishing ) .
viennot d 2010 _ j. math .
phys . _ * 51 * 103501 .
jackson j d 1962 _ classical electrodynamics _
( new york : john wiley & sons ltd ) .
nenciu g and rasche g 1992 _ j. phys . a _ * 25 * 5741 .
rovelli c and vidotto f 2015 _ covariant loop quantum gravity _
( cambridge : cambridge university press ) .
nesterov a i and aceves de la cruz f 2008 _ j. phys . a _ * 41 * 485304 .
mehri - dehnavi h and mostafazadeh a 2008 _ j. math .
phys . _ * 49 * 082105 .
viennot d , jolicard g and killingbeck j p 2008 _ j. phys .
a _ * 41 * 145303 .
dridi g , gurin s , jauslin h r , viennot d and jolicard g 2010 _ phys .
a _ * 82 * 022109 .
breuer h p and petruccione f 2002 _ the theory of open quantum systems _
( new york : oxford university press ) . alsing p m , fuentes - schuller i , mann r b and tessier t e 2006 _ phys .
rev . a _ * 74 * 032326 .
ghosh a , sinha s s and ray d s 2012 _ phys .
e _ * 86 * 011138 .
karmakar a and gangopadhyay g 2010 _ phys .
_ * 85 * 045008 .