quantum parameter estimation , originally motivated by optical metrology  @xcite , considers how to estimate parameters associated with the evolution of quantum systems by physical means .
a theory of quantum estimation , derived from the laws of quantum physics and statistics , has been developed and applied to various evolution processes yielding important results for various metrology situations  @xcite .
we consider the single qubit depolarizing channel , which maps a qubit , initially in the state @xmath0 , via @xmath1 } \ ; \hat{i } + \lambda { \hat{\rho}}\label{eq : depol}\ ] ] where @xmath2 the aim of a physical estimation protocol would be to determine @xmath3 as accurately as possible by subjecting a collection of qubits to this evolution .
the depolarizing channel is interesting for various reasons .
first , depolarization is a standard model of certain noise processes and is of general interest for quantum information processing  @xcite .
specific examples appear in nuclear magnetic resonance ( nmr )  @xcite and optical quantum information processing  @xcite .
second , quantum parameter estimation reveals fundamental and quantifiable differences with classical approaches when manifestly quantum resources such as entanglement are used .
for example , in the estimation of the phase associated with a unitary parameter , classical approaches yield an uncertainty that scales as @xmath4 where @xmath5 is the number of times the channel is invoked , whereas using entangled states can result in a lower uncertainty that scales as @xmath6 , thus pointing to a measurable gain provided by entanglement  @xcite .
depolarizing channel parameter estimation has been investigated previously in the context of finding optimal estimation protocols  @xcite where any possible initial state is assumed to be available .
it emerges that the optimal estimation protocol requires qubit pairs , each initially in an entangled pure state , and that states displaying entanglement or correlations amongst more than two qubits will not give any advantages to estimation accuracy  @xcite .
we consider various estimation protocols for the constrained situation where _ when the available initial state of each qubit is mixed .
_ this situation arises in room temperature solution state nmr , where the initial state for any nuclear spin is highly mixed  @xcite .
we ask whether , for mixed initial state situations , protocols involving correlated states can yield gains in estimation accuracy compared to protocols that use product states and whether increasing the number of correlated qubits beyond two is advantageous .
while we will consider these questions for situations involving initial states with all possible ranges of mixedness , we will focus on the special case where the initial state of each qubit is highly mixed .
our approach is motivated by similar work for other types of channels where gains were found when using correlated states  @xcite and is also inspired by various other approaches to investigating parameter estimation with mixed initial states  @xcite .    this article is organized as follows .
section  [ sec : genestimation ] offers a review of quantum estimation theory .
section  [ sec : indest ] applies this to the simplest depolarizing channel parameter channel estimation protocol , establishes a benchmark against which other schemes are to be compared and reviews existing results for parameter estimation using pure initial states .
section  [ sec : indchanneluseest ] considers mixed initial state estimation protocols where there are never any correlations between qubits .
section  [ sec : correst ] , which contains the key general results of our article , considers an estimation protocol that employs correlated states .
section  [ sec : weakpol ] applies these to the special case where the available initial states are highly mixed .
finally section  [ sec : correlation ] briefly discusses the roles of entanglement and quantum discord in the protocol that uses correlated states .
physical parameter estimation requires that suitably prepared quantum systems be subjected to the operation @xmath7 which is dependent on a real parameter @xmath3 that is to be estimated . the generic scheme for doing this , illustrated in fig .
[ fig : generalscheme ] , assumes that @xmath8 copies of the operation , or channel , are available .
these may be invoked on one or more of @xmath5 systems with invocations surrounded by a series of parameter independent unitary transformations .
quantum systems , each represented by a horizontal line .
this procedure uses @xmath8 invocations of the operation @xmath7 and parameter - independent unitary transformations @xmath9 the initial state of each quantum system is @xmath10 and the pre - measurement state of the collection of systems is @xmath11 measurements on individual systems follow @xmath12 [ fig : generalscheme ] ]    each quantum system is prepared in a known initial state @xmath10 .
after the final unitary , a measurement is chosen and performed on each system and the measurement outcomes are recorded .
the final pre - measurement state @xmath13 may depend on the parameter and thus so will the probabilities with which the measurement outcomes occur .
thus information about @xmath3 might be inferred from the measurement outcomes by classical data processing , thus yielding an estimate for the parameter .
the fact that measurement outcomes in quantum physics are usually statistically distributed implies that different runs of the same estimation procedure will yield estimates that fluctuate about a mean value .
the aim of quantum parameter estimation is to tailor the process by choosing the initial state @xmath10 , the parameter independent unitaries @xmath14 and the final measurements so that the estimates will display minimal fluctuations about the actual parameter value .
the accuracy of the estimation is quantified via these statistical fluctuations as follows .
given a set of measurement outcomes , one on each qubit and denoted @xmath15 , an estimate for the parameter @xmath16 is attained via @xmath17 where @xmath18 is an estimator function .
fluctuations in these measurement outcomes from one run of the estimation procedure to another generate fluctuations in @xmath16 .
the extent to which this occurs depends on the construction of the estimator function .
key requirements for an estimation protocol will be that the estimator is unbiased , i.e. @xmath19 , and that it minimizes the variance , @xmath20 where the angle brackets indicate the mean over all possible measurement outcomes . the variance will depend on the choice of estimator .
however , a key result  @xcite is that _ the variance on any estimator is bounded from below via the cramr - rao bound _
@xmath21 where @xmath22 ^ 2\ ;                   \ ;                   \mathrm{d}x_1 \ldots \mathrm{d}x_n .
\label{eq : classicalfisher}\ ] ] is the ( classical ) fisher information associated with the probability distribution for the measurement outcomes @xmath23 .
this is independent of the choice of estimator .
quite generally , there is always an estimator  @xcite , which asymptotically saturates the bound of eq .  .
thus in classical parameter estimation the fisher information quantifies the accuracy of any possible estimate ; the larger @xmath24 , the better the estimate .    in classical estimation , it is easily shown  @xcite that if the same estimation trial is repeated @xmath8 times independently , the the fisher information becomes @xmath25 where @xmath24 is the fisher information for a single trial .
thus the cramr - rao bound implies a reduction in the variance by a factor of @xmath8 , yielding the typical classical reduction in uncertainty of the estimate by a factor of @xmath26    in quantum evolution processes the probability distribution for the measurement outcomes depends on the type of measurement .
thus , given many possible choices of measurement , there will be many possible values for the classical fisher information for the same evolution process .
quantum estimation aims for a measurement which results in the largest classical fisher information .
this is constrained and simplified by the _
quantum cramr - rao bound , _ which states that , for any measurement the resulting classical fisher information satisfies @xmath27 where @xmath28 is the quantum fisher information and is _ independent of the choice of measurement _  @xcite .
the quantum fisher information depends only on @xmath13 and is calculated via  @xcite @xmath29 }      \label{eq : quantumfisher}\ ] ] where the symmetric logarithmic derivative ( sld ) or score operator , @xmath30 is defined implicitly via @xmath31 .
\label{eq : slddefintion}\ ] ] in general the sld is not immediately apparent from @xmath13 ; sometimes it can be computed using relatively simple algebra  @xcite .
however , the sld and quantum fisher information can always be computed whenever an eigenvalue and eigenstate decomposition for @xmath13 can be found .
specifically , if @xmath32 then  @xcite @xmath33 quite generally , there exists a measurement procedure which asymptotically attains the quantum cramr - rao bound  @xcite .
thus the quantum fisher information quantifies the accuracy of the quantum estimation process .
the main task in quantum parameter estimation is to provide a protocol giving a final pre - measurement state @xmath13 with the maximum quantum fisher information ; this approach has been widely adopted elsewhere  @xcite .    before continuing
, we remark that a useful computational alternative to eq .
is  @xcite @xmath34}.              \label{eq : quantumfishertwo}\ ] ] the quantum fisher information will depend on the resources , such as the numbers of qubits , channel invocations , single qubit unitaries or entangling unitaries that are used and various competing estimation protocols can only be compared in terms of their resource use . in this article
we take a common approach , which regards the number of channel invocations @xmath8 in the circuit of fig .
[ fig : generalscheme ] as the only important resource cost  @xcite . in this context , the goal of any quantum estimation protocol is to produce a final state with the maximum quantum fisher information per channel invocation .
the initial state for the entire system can be assumed to be a product of states for the individual systems . if necessary , any of @xmath14 can be used to generate entangled or otherwise correlated states as the procedure unfolds .
we can then consider two broad classes of estimation schemes .
in the first the unitaries @xmath14 never produce any correlation between the individual systems and the entire system evolves through a series of product states .
various classical approaches to estimation , including repetition and averaging , fall into this realm . in the second ,
some of the unitaries @xmath14 produce states where there is some form of quantum correlation , such as entanglement , between the individual systems .
frequently quantum estimation aims to show that using such correlated states results in a parameter estimate with a greater quantum fisher information per channel invocation  @xcite than is possible with non - correlated states ; indirectly this indicates the utility of quantum resources for information processing .      within the framework of fig .
[ fig : generalscheme ] , the typical task is to find the initial state @xmath10 and the unitaries @xmath14 that produce a final state @xmath13 that yields the maximum quantum fisher information per channel invocation .
the convexity of the quantum fisher information guarantees that this will be attained by using a pure initial state  @xcite , and most investigations of depolarizing channel parameter estimation have considered this case  @xcite .
however , it may be that case that the available initial states are not pure .
our aim is to consider various estimation protocols when the available initial states are limited to those that are mixed .
the general initial state of a single qubit is @xmath35 where @xmath36 and @xmath37 are the components of a three dimensional real vector , @xmath38 .
the magnitude of this vector @xmath39 is called the _ polarization or purity _ of the state and it satisfies @xmath40 when @xmath41 the state is pure and as @xmath42 , which we use to quantify the mixedness or purity of the state , decreases the state becomes more mixed .
thus the principal issue of such mixed state quantum estimation is : _ given that the available qubits have initial states each with polarization @xmath42 , how do the accuracies , quantified by the quantum fisher information , of various parameter estimation protocols compare ? _ note that , in this approach , the orientation of the polarization vector @xmath38 is irrelevant , as this can be modified at no cost with a single qubit unitary without changing the polarization @xmath42 .
the simplest parameter estimation scheme , which also forms the baseline against which all other schemes may be compared , is that where the channel is invoked once on a single qubit and is described in fig .
[ fig : sqsc ] . this will be termed the single qubit , single channel ( sqsc ) protocol .    .
a parameter independent unitary transformation @xmath43 precedes the single channel invocation , producing a channel input state , @xmath44 .
a measurement follows the channel invocation .
[ fig : sqsc ] ]    the unitary @xmath43 of fig .
[ fig : sqsc ] only changes the orientation of @xmath38 , leaving the polarization unaltered and thus can be omitted .
then eq .
yields that the state after one channel invocation is @xmath45 .
the quantum fisher information can be computed by first determining the sld using prop .  1 of  @xcite and then substituting into eq .  .
this proposition ( appropriately corrected ) states that if @xmath46 then @xmath47/\operatorname{tr}{{\hat{\rho}_\mathrm{f}}}$ ] while if @xmath48 then @xmath49/\operatorname{tr}{{\hat{\rho}_\mathrm{f } } }                              + \frac{\partial}{\partial \lambda } \ln{(\lvert \alpha/\operatorname{tr}{{\hat{\rho}_\mathrm{f } } } \rvert ) } { \hat{i}}.$ ] straightforward calculation gives that @xmath50 and @xmath51 the only circumstances under which @xmath52 are when @xmath41 _ and _ @xmath53 we henceforth only consider @xmath54 it follows that , _ for a single channel invocation on the state of eq .
_ the sld is @xmath55 substituting into eq .
then gives that _ the quantum fisher information for the sqsc protocol using an initial state with polarization @xmath42 _ is @xmath56 this is the baseline quantum fisher information against which all other estimation schemes are to be compared .
we briefly digress to discuss the previously investigated  @xcite extreme case , where the initial states are pure . if it were possible to modify the initial state polarization , then eq .
, shows that , in accordance with general results  @xcite , the largest quantum fisher information would be attained for a pure input state , i.e.  when @xmath57 giving @xmath58 usually the key issue of interest in quantum parameter estimation is whether a greater quantum fisher information can be attained when entangled multiple qubits are available as input for the channel . for the depolarizing channel
it has been shown that entanglement can yield a greater quantum fisher information .
specifically if the channel is applied to one of two qubits then the optimal quantum fisher information is attained  @xcite when the two qubits are prepared in a maximally entangled pure state , such as @xmath59 .
the resulting quantum fisher information is @xmath60 it is straightforward to show that @xmath61 whenever @xmath62    the remaining question is whether any additional qubits and channel invocations can yield a greater quantum fisher information per channel invocation .
remarkably the result is  @xcite that the optimal depolarizing channel parameter estimation scheme uses pairs of qubits , each in a maximally entangled pure input state .
the channel is applied once to only one of each entangled pair and the resulting quantum fisher information per channel use is that of eq .  .
thus _ when qubits are available in pure initial states then entangling more than two qubits can not yield any advantages for depolarizing channel parameter estimation .
_ however , we stress that this does not necessarily pertain to the situation that we address , in which the available input states are not pure .
the previous discussion motivates the central questions of our work : _ given a collection of qubits , each initially in a state with the same polarization @xmath42 , is it possible to enhance the parameter estimation accuracy by using multiple channel invocations or entanglement or correlations amongst more than two qubits ?
if so to what extent can the accuracy , as measured by the quantum fisher information per channel use , be enhanced ? _
a first strategy for enhancing estimation accuracy is one where the same process is repeated independently , as is done for classical repeat and average schemes .
specifically we consider the protocol of fig .
[ fig : indchanneluse ] , called the _ independent channel use protocol . _ here there are @xmath8 qubits , each in the initial state @xmath10 and subject to a parameter independent single qubit preparatory unitary .
the channel is subsequently invoked exactly once on each qubit . in this protocol
there are never any correlations between the qubits during the entire process ; in this sense it resembles a classical approach .
single qubits , each initially in the state @xmath10 .
qubit @xmath63 is subjected to a single qubit unitary , @xmath64 , which is independent of the channel parameter , producing the input state @xmath65 thereafter a depolarizing channel is applied once to each qubit .
[ fig : indchanneluse ] ]    the quantum fisher information is additive in such cases  @xcite and thus the quantum fisher information for the independent channel use protocol with @xmath8 channel invocations is @xmath66 although this is an improvement over the sqsc protocol , clearly the quantum fisher information _ per channel use _ in the independent channel use protocol is identical to that of the sqsc protocol .
thus _ the independent channel use protocol attains the same estimation accuracy , measured in terms of quantum fisher information per channel invocation , as the sqsc protocol . _    a second strategy , illustrated in fig .
[ fig : seqchanneluse ] , is that where _ only one qubit is available and the channel is invoked sequentially a total of @xmath8 times . _
successive channel invocations are interspersed with parameter - independent unitaries .
is subjected to @xmath8 channel invocations .
each channel invocation is preceded a single qubit unitary , @xmath67 , all of which are independent of the channel parameter .
[ fig : seqchanneluse ] ]    assuming the standard initial state , @xmath10 , the effect of @xmath43 is to rotate the polarization vector , @xmath68 without altering its magnitude .
.   implies that the state after the first channel invocation is @xmath69 . using the same argument , the state after the second channel invocation is @xmath70 where @xmath71 rotates @xmath72 continuing , the state prior to measurement is @xmath73 where @xmath74 a calculation similar to that giving eq
.   yields that the quantum fisher information for @xmath8 sequential channel invocations on a single qubit is @xmath75 thus the quantum fisher information per channel invocation is @xmath76 to compare this to the sqsc protocol we define the _ sequential channel use gain _ as @xmath77 figure  [ fig : gainseqplot ] illustrates the typical gain , plotted here for the @xmath78 case .
this shows that the sequential channel use protocol is advantageous for certain ranges of parameter values and polarizations but disadvantageous for others .    .
[ fig : gainseqplot ] ]    the sequential channel use gain decreases monotonically with @xmath42 and increases monotonically with @xmath3 ( these are proved in appendix  [ app : gainseq ] ) .
a more detailed analysis of this separates the pure initial state case , i.e.  @xmath41 , from mixed initial state cases , i.e.  @xmath79 . for the _ pure initial state case ,
_ it is shown in appendix  [ app : gainseq ] that @xmath80 for the _ mixed initial state case , _ it is also shown appendix  [ app : gainseq ] that @xmath81 in both cases the upper bound is attained when @xmath82    it follows that _ if the available initial state is pure there is no advantage to the sequential channel use protocol .
if the available initial state is mixed then , for certain initial polarizations and the parameter values , there can be an advantage to the sequential channel use protocol .
_    it remains to assess whether for mixed initial states , there is any advantage to increasing the number of channel invocations .
to do so , we compare the quantum fisher information per channel use of eq .   for @xmath83 channel invocations versus that for @xmath8 channel invocations .
an advantage is attained when the ratio of that for @xmath83 to that for @xmath8 exceeds @xmath84 , and straightforward algebra shows that this is equivalent to @xmath85/\lambda^{2m+2}$ ] .
the right hand side of this increases monotonically with @xmath3 and attains a maximum of @xmath84 when @xmath86 .
also it is negative when @xmath87 , implying that when @xmath87 there is no polarization for which an additional channel use will give an advantage .
thus for any polarization there will be some parameter values for which additional invocations in the sequential channel use protocol will be advantageous .
however , the range of possible parameter values for which this occurs , bounded from below by at least @xmath88 , shrinks with increasing channel invocations .
an additional resource available to multiple qubit systems , is the existence of states in which there are correlations between the individual qubits and we now consider whether these can enhance the accuracy of parameter estimation for the depolarizing channel .
we consider the particular correlated state protocol procedure of fig .
[ fig : preparedscheme ] ; this was previously investigated for other channels  @xcite .
we assume that each qubit is initially in the particular state @xmath89 with polarization @xmath42 and that a preparatory unitary @xmath90 , illustrated in fig .
[ fig : uprep ] , is subsequently applied jointly to the qubits .
single qubits , each initially in the state @xmath10 .
these are subjected to a joint preparatory unitary , @xmath90 , producing the channel input state @xmath65 thereafter a depolarizing channel is applied once to each of the uppermost @xmath8 qubits .
the preparatory unitary is illustrated in fig .
[ fig : uprep ] .
[ fig : preparedscheme ] ]    denote the density operator for the system after the preparatory gate by @xmath91     consists of a collection of controlled-@xmath92 operations applied to each distinct pair of qubits followed by a hadamard transformation applied to each qubit .
[ fig : uprep ] ]    it can be shown  @xcite that @xmath93 where @xmath94 and @xmath95                  \\                &                           + i { c_{j(x)}}\ ;                              \biggl [                     { \lvert x \rangle}{\langle n - x \lvert } -                               { \lvert n - x \rangle } { \langle x \lvert }                \biggr ] .
\end{split }   \label{eq : rhoprepped}\ ] ] the counter - diagonal terms of the resulting matrix of eq .
are constructed by first computing the number of zeroes @xmath96 in the binary representation @xmath97 for @xmath98 and then using @xmath99 .
\label{eq : prepcounterdiag }   \end{aligned}\ ] ] the diagonal terms are @xmath100 .
\label{eq : prepdiag }   \end{aligned}\ ] ] as described previously  @xcite , this scheme results in a density operator which has supports consisting of mutually orthogonal two dimensional subspaces ; in terms of matrix representations it can be expressed as a sum of @xmath101 matrices acting on distinct vector spaces .
when the initial states are pure , i.e.  @xmath57 the diagonal and counter - diagonal terms are all zero except @xmath102 and @xmath103 and @xmath104 is the density operator corresponding to the ghz - type state @xmath105    suppose that the channel acts _ once on each of @xmath8 qubits _ corresponding to the least significant ( or rightmost ) digits in the computational basis representation .
the action of the channel on multiple qubits is readily described by noting that eq .
implies that , for a single qubit ,    [ eq : depoloncompbasis ] @xmath106    where @xmath107 and @xmath108 are computational basis states such that @xmath109 and @xmath110 here @xmath111 and @xmath112 . a general state of multiple qubits can be expressed in terms of linear combinations and tensor products of basic operators such as @xmath113 and thus eqs .   can be used to describe the action of any number of channels acting on multiple qubits .
eq   implies the counter - diagonal terms of the matrix representing @xmath13 are attained from those of @xmath104 by multiplication by @xmath114 eq .
implies that the diagonal terms of the matrix representing @xmath13 result from rearrangements of those of @xmath104 multiplied by various factors of @xmath115 and @xmath116 .
the remaining terms are again all @xmath117 .
an additional simplification , demonstrated in appendix  [ app : rhofdiag ] is that @xmath118 thus the final density operator has a form similar to that for @xmath104 , i.e.@xmath119 where @xmath120                  \\                &                           + i { c_{\textrm{f}\ ; x}}\ ;                              \biggl [                     { \lvert x \rangle}{\langle n - x \lvert } -                               { \lvert n - x \rangle } { \langle x \lvert }                \biggr ] .
\end{split }   \label{eq : rhofinal}\ ] ] eq .
yields @xmath121 for the counter - diagonal terms where @xmath96 is the number of zeroes in the bit string for @xmath98 .    determining the diagonal terms
is further simplified by noting that the preparatory unitary is invariant under interchange of any two qubits and that this is also true of the channel operation provided that the channel either acts on each qubit that is interchanged or that it acts on neither of the qubits that are interchanged .
the invariance of the preparatory unitary is readily apparent from its structure , as shown in fig .
[ fig : uprep ] .
more precisely , suppose that @xmath122 is the unitary operation which swaps the qubits labeled @xmath63 and @xmath123 .
then if @xmath124 ( recall that the channel acts on the rightmost @xmath8 bits ) , @xmath125 with an identical result when @xmath126 . thus if the bit string for @xmath98 is @xmath127 , then @xmath128 is invariant under the interchange of any two of @xmath129 and is also invariant under the interchange any two of @xmath130 ( it is not in general invariant under interchanges between one element from one of these sets and another from the other set ) .
for example , if the channel is applied the two least significant qubits out of a total of five qubits then @xmath131 ( or , using a binary representation for the index subscript , @xmath132 ) and @xmath133 ( or @xmath134 ) .
a consequence of this is that @xmath135 only depends on the number of zeroes in the leftmost @xmath136 bits of the binary representation of @xmath98 and also the number of zeroes in the rightmost @xmath8 bits of the binary representation of @xmath98 .
this consideration and a derivation , provided in appendix  [ app : rhofdiagterms ] , yields @xmath137 where @xmath138 is the number of zeroes in the leftmost @xmath136 bits ( i.e.  @xmath139 ) of the binary representation of @xmath98 , @xmath140 is the number of zeroes in the rightmost @xmath8 bits ( i.e.  @xmath141 ) of the binary representation of @xmath98 , @xmath142 and @xmath143 .
clearly @xmath135 depends on @xmath144 and @xmath8 and @xmath3 , via @xmath145 and @xmath146 .
since @xmath3 is implicit throughout this work and @xmath8 is fixed by the number of channel invocations we will find it convenient to express @xmath135 as @xmath147    according to eqs .   and  , @xmath148 is a sum of operators , each acting on a two dimensional space , which have mutually orthogonal supports . thus prop .
2 of  @xcite implies that the quantum fisher information is @xmath149 where @xmath150 is the quantum fisher information associated with @xmath151 a derivation ( see appendix  [ app : qfifirst ] ) based on this and eq .
yields @xmath152 \label{eq : qfidepolone}\end{aligned}\ ] ] where the dot indicates differentiation with respect to @xmath153 although the sum in eq .   is over @xmath98 , the terms indexed as such depend on @xmath154 and @xmath155 and it follows that the quantum fisher information can be expressed as a sum over these . in appendix  [ app : qfifirst ] it is shown that _ if @xmath156 then _
@xmath157 \label{eq : qfidepoltwomlessthann}\end{aligned}\ ] ] where @xmath158 is given by eq .  .
_ if @xmath159
_ then @xmath160 \label{eq : qfidepoltwomequalsn}\end{aligned}\ ] ] where @xmath161 is again given by eq .  .
in both cases @xmath162 is given by eq .  .
these represent the simplest expressions that we have been able to obtain for the quantum fisher information for the correlated state protocol .
the quantum fisher information per channel invocation is @xmath163 we aim to compare the performance of the correlated state protocol to that of the sqsc protocol _ in terms of the quantum fisher information per channel use and given that the initial qubit polarizations are the same .
_ we define the _ correlated state protocol quantum fisher information gain _ as @xmath164 and whenever this exceeds @xmath84 , the correlated state protocol outperforms the sqsc protocol .
the expression for the quantum fisher information is a rational polynomial in @xmath42 and @xmath3 and even with as few as two qubits , this is too complicated to yield insight into the estimation accuracy for the correlated state protocol for all values of polarization and the parameter .
while the same applies to the gain , aspects of its behavior can be illustrated graphically .
first , suppose that there is only one channel invocation ( @xmath165 ) .
[ fig : gainplotssingle ] illustrates the gain for the examples @xmath166 and @xmath167 .     and
@xmath40 [ fig : gainplotssingle ] ]    these suggest that when the channel is invoked once only , the correlated state protocol will be more accurate than the sqsc protocol , assuming the same initial state polarization for both .
however , a general proof of this is elusive ; the methods of  @xcite used to prove a similar fact for the phase - flip channel are not immediately applicable here .
second , suppose that there are multiple channel invocations ( @xmath168 ) .
[ fig : gainplotsmultiple ] illustrates the correlated state gain of eq .   for various numbers of channel invocations on @xmath169 qubits .
all plots are over the range @xmath170 .
[ fig : gainplotsmultiple ] ]    these show that when the channel is invoked multiple times , there are polarizations and parameter values for which the the correlated state protocol will be more accurate than the sqsc protocol but there are also parameter values for which it is inferior .    given that the channel is invoked multiple times , the correlated state protocol can also be compared to the sequential channel protocol that uses the same number of channel invocations .
we define the correlated vs sequential channel gain as @xmath171 where the number of channel invocations is the same for the two protocols
.     qubits .
plots on the left are over the range @xmath172 .
plots on the right provide details for @xmath173 [ fig : gaincorrmultvsseq ] ]    figure  [ fig : gaincorrmultvsseq ] suggests that when the channel is invoked a fixed number of times , the correlated state protocol is more accurate than the sequential protocol , using the same number of channel invocations , over all parameter values and for all polarizations .
in some situations the initial state polarization is very small , i.e.  @xmath174 this is true of room - temperature solution state nmr  @xcite , where typically @xmath175    when  @xmath176 the sqsc protocol quantum fisher information becomes @xmath177 and the sequential channel use protocol quantum fisher information per channel use becomes @xmath178 each to lowest order in @xmath179 in appendix  [ app : smallrqfi ] we show that the correlated state quantum fisher information per channel use becomes @xmath180 to lowest order in @xmath42 .    here
the sequential channel use gain approximates to @xmath181 eq .
shows that the gain exceeds @xmath84 when the depolarizing channel parameter satisfies @xmath182 where @xmath183 is the characteristic parameter cutoff .
the characteristic parameter cutoff , plotted in fig .
[ fig : lambdacutoff ] , is a monotonically increasing function of @xmath8 , with a minimum value of @xmath184 attained as @xmath185 and a maximum value of @xmath84 attained as @xmath186 .
this indicates that for any given parameter value , there will be a maximum number of channel invocations beneath which the sequential channel use protocol offers enhanced parameter estimation accuracy and beyond which it will not do so .
figure  [ fig : lambdacutoff ] shows that as @xmath3 increases , accuracy gains will persist with increasing number of channel invocations .
note that the larger @xmath3 is , the weaker the depolarization is , and for such weak depolarization cases enhanced estimation accuracies can be attained with multiple sequential channel invocations .    .
[ fig : lambdacutoff ] ]    thus , _ for low polarizations , the sequential channel use protocol can yield a greater estimation accuracy per channel invocation than the sqsc protocol .
_ however , this requires that the parameter value exceed the characteristic parameter cutoff and fig .
[ fig : lambdacutoff ] indicates that such advantages will be absent for low enough parameter values and , even when there is an advantage , increasing the number of channel invocations will eventually eliminate this advantage .
the correlated state protocol gain ( versus the sqsc protocol ) approximates to @xmath187 we consider two cases separately .
first , with _ only one channel invocation the correlated state protocol offers an @xmath5-fold gain in accuracy regardless of the parameter value , for low initial polarization states .
_ this is strikingly different from the pure state case where the optimal accuracy is attained  @xcite with a single channel invocation on one of only two qubits ( @xmath166 and @xmath165 ) ; by contrast when the initial state polarization is low and the channel is only invoked on one qubit , then adding more qubits always increases the estimation accuracy .
second , with multiple channel invocations , eq .   indicates that enhanced estimation accuracy will occur if @xmath188 the right hand side forms a new cutoff which is again an increasing function of the number of channel invocations but a decreasing function of the total number of qubits . thus increasing the number of spectator qubits
while leaving the number of channel invocations fixed extends the range of parameter values over which the correlated state protocol is advantageous . furthermore eq .
shows that adding more spectator qubits while keeping the number of channel invocations fixed will increase the quantum fisher information per channel use ; the increase can be arbitrarily large by adding sufficiently many spectator qubits .
_ unlike the pure initial state case , for low polarization initial states , there will always be an advantage to adding additional spectator qubits .
_    separately , when @xmath156 we can assess whether one additional channel invocation yields any advantage .
.   implies that starting with @xmath8 channel invocations and adding one more invocation increases the quantum fisher information if @xmath189 , decreases it if @xmath190 and makes no difference if @xmath191 thus the optimal number of channel invocations in the correlated state protocol is @xmath192 unless @xmath193 is an integer , in which case , @xmath194 examples of this and the associated gains are provided in table  [ tab : optinvocations ] .
ccd @xmath3 & @xmath195 & + @xmath196 & @xmath84 & 1.00 n + @xmath197 & @xmath198 & 1.28 n + @xmath199 & @xmath200 & 2.15 n + @xmath201 & @xmath202 & 3.97 n + @xmath203 & @xmath204 & 18.67
n + @xmath205 & @xmath206 & 37.07 n    when @xmath3 is close enough to @xmath84 , this offers considerable opportunity for improved accuracy . again _ unlike the pure initial state case , for low initial polarization , increasing the number of channel invocations in the correlated state protocol can result in enhanced estimation accuracy , per channel use . _    in many physical situations every qubit is subjected to the depolarizing channel .
when the channel is invoked on all qubits of the correlated state protocol ( @xmath207 ) , the gain becomes @xmath208 and this implies that the correlated state protocol provides enhanced accuracy whenever @xmath209 this lower bound is plotted in fig .
[ fig : lambdacutoff ] .
this is again an increasing function of @xmath8 , approaching @xmath210 as @xmath211 .
the slightly lower cutoff implies that there is somewhat greater scope for enhanced parameter estimation when all qubits are subjected to the channel .
however , for a given parameter value , this cutoff establishes an upper bound on the number of qubits / channel invocations for which the correlated state protocol provides any gain .
we can again ask when adding one more qubit and channel invocation gives any advantage .
equation   with @xmath207 implies that if there are originally @xmath8 channel invocations then adding one more gives a greater quantum fisher information per channel use if @xmath212 where @xmath213 unless @xmath214 is an integer , in which case , @xmath215 examples of this and the associated optimal gains are provided in table  [ tab : optinvocationsall ] .
ccd @xmath3 & @xmath195 & + @xmath216 & @xmath84 & 1.00 + @xmath196 & @xmath217 &
2.16 + @xmath199 & @xmath218 & 15.01 + @xmath201 & @xmath219 & 56.96 + @xmath220 & @xmath221 & 155.03 + @xmath203 & @xmath206 & 1367.00    clearly , for relatively weak depolarizing channels , i.e.  @xmath222 substantial enhancements of estimation accuracy are possible with the correlated state protocol .
finally , the correlated state protocol and sequential channel protocol can be compared in the low initial polarization limit , given the same number of channel invocations . here the correlated vs sequential channel gain as @xmath223 and thus when @xmath224 the correlated state protocol will always offer a better estimation accuracy than the sequential channel protocol provided that the number of channel invocation for each is the same .
the presence of quantum correlations such as entanglement or discord in the correlated state protocol may be assessed for the @xmath166 along the same lines as was done for the pauli channel parameter estimation of  @xcite .
then eqs .  , and yield that with @xmath8 channel invocations , @xmath225 in the computational basis
. the state of the system prior to channel invocation can be obtained from this by setting @xmath53    the presence of entanglement can be assessed by considering partial transpose of @xmath148 and this is @xmath226 the state is separable if @xmath227 is positive  @xcite .
the eigenvalues of this will only be positive if and only if @xmath228 thus the state prior to channel invocation will be separable if and only if @xmath229 the correlated state protocol can give estimation gains as @xmath230 and this clearly corresponds to a state which is separable throughout the protocol . thus the correlated state protocol gains can not always be ascribed to the presence of entanglement .
quantum discord describes correlations between quantum systems in terms of mutual information and measures the extent to which these exceed any possible classical correlations between the systems  @xcite .
it is possible to compute the quantum discord for states with the form  @xcite @xmath231 where @xmath232 are all real and @xmath233 defining    [ eq : evalsforq ] @xmath234    and @xmath235 yields that the quantum discord is  @xcite @xmath236 following the scheme of  @xcite , the final density operator can be brought into the form of eq .   by applying the single qubit unitary @xmath237 to each qubit and
this will not alter the discord .
then @xmath238 and @xmath239 since @xmath40 again the discord for the correlated state protocol prior to channel invocation can be determined by setting @xmath86 and this gives @xmath240 then @xmath241 whenever @xmath242 and @xmath243 when @xmath244 implies that the state prior to channel invocation has non - zero discord whenever the polarization of each qubit is non - zero . in this sense , we can regard the correlated state protocol as using states that have quantum correlations .
we have considered estimation of the depolarizing channel parameter for various situations where the available initial state of each qubit is mixed and we quantified the accuracy of each in terms of the quantum fisher information after channel invocation .
we showed that for certain ranges of initial state polarization and channel parameter , a protocol that uses correlated states prior to channel invocation attains a larger quantum fisher information per channel invocation than either a protocol which uses product states throughout or another protocol which uses a sequence of channel invocations on a single qubit .
for qubits with small initial polarizations we showed that using more than two qubits and more than two channel invocations can be advantageous ; this is in contrast to the situation where qubits in pure initial states are available .
our results suggest that , as with pauli channel parameter estimation , a single invocation of the channel on many correlated qubits gives gains for all parameter values and initial polarizations .
we were only able to prove this for the depolarizing channel in the limit as the initial polarization approaches zero , but we were also unable to generate a counterexample .
this was not central to our results but it does merit further investigation .
the source of the gain can not always be attributed to the presence of entangled or non - separable states prior to channel invocation ; we demonstrated this for the two qubit example .
we also showed that there is non - zero discord present prior to channel invocation in this example
. a full analysis of the evolution of the discord and its possible relationship to gains in the quantum fisher information might provide insight into the source of these accuracy gains ; this is beyond the scope of our article and remains to be investigated .    finally , we chose one particular preparatory unitary transformation for the analysis of the correlated state protocol ; that this yields gains in some cases is interesting enough .
however it is unclear whether this is optimal or whether there exists another parameter independent unitary that offers greater gains .
the behavior of @xmath245 with respect to @xmath42 can be assessed by first setting @xmath246 .
then @xmath247 gives that @xmath248 since @xmath42 increases as @xmath98 increases this proves the fact that @xmath245 is a monotonically decreasing function of @xmath42 .    to assess the behavior of @xmath245 with respect to @xmath3 ,
let @xmath249 and note that @xmath250 then @xmath251 and @xmath252 thus @xmath253 the term @xmath254 is negative since @xmath255 and @xmath256 and thus decreases as @xmath257 increases .
the remaining term @xmath258 increases as @xmath257 increases .
thus @xmath245 decreases as @xmath257 increases .
however , @xmath257 decreases as @xmath3 increases and this proves the fact that @xmath245 is a monotonically increasing function of @xmath3 .    for the pure initial state case , @xmath41 and eq .
gives @xmath259 where @xmath260 .
then @xmath261 gives @xmath262 this is a monotonically decreasing function of @xmath263 and thus attains a maximum of @xmath84 at @xmath264 thus @xmath265 for a pure initial state , where the lower bound stems from the fact that the gain is a ration of positive quantities .
for the mixed initial state case , @xmath245 is a monotonically increasing function of @xmath3 and attains a maximum of @xmath8 when @xmath82 this shows that @xmath266 when @xmath267
eqs .   give that the diagonal terms of @xmath148 only arise from the action of the channel on @xmath268       \label{eq : rhoprepdiag}\ ] ] and that this will be mapped to the diagonal component of the final density operator , denoted @xmath269 and defined via @xmath270 a useful expression for @xmath269 involves the following representation of the channel on any diagonal operator .
eq .   gives @xmath271 thus if the channel acts on the rightmost qubit , labeled @xmath84 , then its action on the diagonal component of the density operator is @xmath272 where @xmath273 and the superscript on @xmath274 indicates action on qubit @xmath84 .
suppose that the channel acts on qubit @xmath84 and then on qubit @xmath198 .
this maps the diagonal component of the prepared density operator as @xmath275                                                                                \nonumber \\                                                                                & &      + q^2        { \hat{x}}_2{\hat{x}}_1 { \hat{\rho}_\textrm{prep diag}}{\hat{x}}_2{\hat{x}}_1                                                                                \end{aligned}\ ] ] where @xmath276 this extends to action of the channel once on each of qubits @xmath277 and gives @xmath278 where the inner sum is over all @xmath279 that satisfy @xmath280 and range over the entire possibilities starting with @xmath281 to @xmath282 if @xmath283 the inner sum is taken to mean @xmath284 .
for example , if @xmath285 and @xmath286 there are three indices @xmath287 and the summation range yields the following possibilities for @xmath288 : @xmath289 this consists of the set of all flips on three out of the four qubits with least significant label .
a more symmetrical alternative is @xmath290 where the sum is over all @xmath291 where no two indices are equal .
we introduce the following notation to render this more compact .
let @xmath292 where , @xmath291 and no two elements in this list are equal . then define @xmath293 it follows that @xmath294 we use this to prove that @xmath295 to do so , note that @xmath296 and let @xmath297 thus , using @xmath298 @xmath299 eq
.   gives    @xmath300    since @xmath301 commutes with @xmath302 .
.   gives @xmath303 this and eq .
let the binary representation of @xmath98 be @xmath305 then @xmath306 and eq .
give    @xmath307    with the same constraints for the inner sum as for eq .   and where @xmath308 .
we shall evaluate the inner sum for any value of @xmath309 , which represents the number of flips on the bits within the binary representation of @xmath98 .
the fact that the diagonal terms of @xmath148 only depend on the number of ( and not their order ) zeroes , @xmath154 , in the @xmath136 most significant bits and also the number of zeroes , @xmath155 , in the @xmath8 least significant bits implies that it is sufficient to consider @xmath135 where the binary representation of @xmath98 is @xmath310 here eqs .   and   give @xmath311 with the same limits for the inner sum .
these limits can be modified by observing that @xmath312 only depends on the number of zeroes and ones in the binary representation of @xmath98 together with the number and arrangement of flips .
suppose that , of the @xmath309 flips , @xmath313 occur on the rightmost cluster of @xmath155 bits .
thus @xmath314 occur on the second to right cluster of @xmath315 bits .
then the bit representation , @xmath316 , contains @xmath154 zeroes from the leftmost cluster , @xmath117 from the next cluster to the right , @xmath317 from the second to right cluster and @xmath318 from the rightmost cluster .
thus the total number of zeroes in this bit string are @xmath319 and @xmath320 the exact location of the flips within each cluster is irrelevant .
thus within the inner sum of eq .
there are @xmath321 ways of attaining @xmath322 , consisting of @xmath323 ways of flipping @xmath313 bits in the rightmost cluster together with @xmath324 ways of flipping @xmath314 in the second to right cluster
. finally not every value of @xmath313 is possible and these are constrained by @xmath325 ( at most @xmath155 flips in the rightmost cluster ) and @xmath326 ( at most @xmath315 flips in the second to right cluster ) . taken together these give @xmath327    combining these yields eq .  .
we shall show that @xmath328\end{aligned}\ ] ] and the result will follow from eq .  .
the expression for @xmath329 comes from the eigenvalue decomposition of @xmath151 the eigenvalues are of @xmath330 @xmath331 and the associated eigenvectors are @xmath332 the fact that the eigenvalues are independent of @xmath3 means that the second term on the right hand side of eq .
is zero and the first term then yields @xmath333 straightforward algebra yields the desired result .
the expression for the quantum fisher information can be rewritten as a sum over @xmath154 and @xmath155 .
the only relevance that each value of @xmath98 has for the summand is that it determines @xmath334 and @xmath335 typically many values of @xmath98 will yield the same combination of values of @xmath154 and @xmath155 . to account for this , note that the sum in eq .
runs from @xmath336 to @xmath337 ( in binary notation ) .
the fact that the leftmost bit is always zero gives two distinct cases .
in the first , where @xmath338 this bit structure reveals that @xmath339 and @xmath340 . for any choice of @xmath154 , there are @xmath341 ways in which the zeroes can appear in the leftmost @xmath136 bits , given that the leftmost single bit must be zero .
there are independently @xmath342 ways in which the zeroes can appear in the rightmost @xmath8 bits .
thus there are @xmath343 values of @xmath98 in the range of the sum which result in exactly the same pair of values @xmath154 and @xmath155 . for each the contribution of
the summand is identical .
this means that the terms within the summand can be expressed in terms of @xmath154 and @xmath155 ( rather than @xmath98 ) and the result sum reduces to that of eq .  .
in the second , where @xmath207 , clearly @xmath344 and @xmath345 a similar counting argument reveals that there are @xmath346 ways in which the zeroes can appear in the rightmost @xmath347 bits , since the rightmost bit must be zero .
a similar argument then results in eq .  .
the correlated state protocol quantum fisher information , eqs .   and  , requires @xmath162 and @xmath348 which in turn requires @xmath349 eqs ,   and   yield , to second order in @xmath42 , @xmath350}{2^{n+1}}.\end{aligned}\ ] ] the latter of these gives @xmath351 where @xmath352 and @xmath353}{2^{n+1}}.\end{aligned}\ ] ] both @xmath354 and @xmath355 are independent of @xmath42 .
then @xmath356 which can be demonstrated indirectly by noting that , when @xmath357 each diagonal term of the final density operator equal @xmath358 however , when @xmath357 @xmath359 and thus @xmath360 .
the channel does not alter this and thus , when @xmath244 , @xmath361 .
it follows that @xmath362 this fact can independently be verified by explicitly evaluating the three summations in @xmath358    note that these imply @xmath363 then eqs  ,   and  , yield that for the terms in eqs .
and   @xmath364 thus if @xmath156 , eqs .
becomes @xmath365 the two summations in eq   give @xmath366 and this yields the result of eq .  .
a similar derivation gives the same result when @xmath367    48ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
+ 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty link:\doibase    10.1103/revmodphys.52.341 [ * * ,   ( ) ] link:\doibase 10.1103/physrevd.23.1693 [ * * , ( ) ] link:\doibase 10.1103/physreva.43.3795 [ * * , ( ) ] @noop _ _  ( ,  ,  ) link:\doibase 10.1103/physrevlett.72.3439 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.79.3865 [ * * ,   ( ) ] http://stacks.iop.org/0305-4470/39/i=26/a=015 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.96.010401 [ * * ,   ( ) ] link:\doibase 10.1142/s0219749909004839 [ * * ,   ( ) ] @noop _ _  ( ,  ,  ) \doibase doi:10.1038/nphoton.2011.35 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.54.3824 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.56.3470 [ * * , ( ) ] @noop _ _  ( , ,  ) link:\doibase 10.1103/physreva.82.012306 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.86.5811 [ * * ,   ( ) ] http://stacks.iop.org/1367-2630/11/i=1/a=013034 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.92.257901 [ * * ,   ( ) ] link:\doibase 10.1002/mop.21210 [ * * ,   ( ) ] link:\doibase 10.1364/josab.25.000668 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.105.180402 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.85.042112 [ * * , ( ) ] link:\doibase 10.1103/physreva.63.042304 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.66.022308 [ * * ,   ( ) ] http://stacks.iop.org/0305-4470/36/i=29/a=314 [ * * ,   ( ) ] http://stacks.iop.org/1751-8121/44/i=20/a=205306 [ * * ,   ( ) ] link:\doibase 10.1103/physrevx.1.021022 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.87.032301 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.72.042327 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.100.100503 [ * * , ( ) ] link:\doibase 10.1140/epjd / e2010 - 00195 - 3 [ * * , ( ) ] link:\doibase    10.1103/physreva.83.063836 [ * * ,   ( ) ] link:\doibase    10.1103/physreva.88.040102 [ * * ,   ( ) ] @noop _ _  ( , ,  ) @noop ( ) ,   in  link:\doibase 10.1007/978 - 4 - 431 - 54493 - 7_4 [ _ _ ] ,  ( ,  )  pp .
http://stacks.iop.org/0305-4470/33/i=24/a=306 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.69.022303 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.70.012317 [ * * ,   ( ) ]
link:\doibase 10.1103/physreva.72.052334 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1109/tit.2008.929940 [ * * , ( ) ] link:\doibase 10.1038/nphys1958 [ * * , ( ) ] link:\doibase    10.1103/physrevlett.98.090501 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.77.1413 [ * * ,   ( ) ] link:\doibase 10.1016/s0375 - 9601(96)00706 - 2 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.88.017901 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.77.042303 [ * * , ( ) ] link:\doibase    10.1103/physreva.84.042313 [ * * ,   ( ) ]


in our planetary system , most objects go around the sun following counterclockwise ( prograde ) orbits as viewed from above the north pole of our star .
however , a large number of comets orbit the sun clockwise ( retrograde ) and a substantial percentage of these objects have highly inclined orbits .
they are believed to have their origin in the oort cloud ( see e.g. duncan 2008 ) , a vast and remote spherical reservoir of cometary material tens of thousands of astronomical units in diameter that completely surrounds the solar system .    out of 3265 objects currently classified by the jet propulsion laboratory ( jpl ) small - body database as comets ,
1918 move in retrograde orbits ( 58.7% ) . for long - period comets
the retrograde fraction is 49.7% , for halley - type comets is 35.7% and for jupiter - family comets is 16.7% .
in contrast , only 53 of the 641502 currently catalogued asteroids orbit the sun clockwise ( 0.008% ) . for asteroids in the outer solar system ( catalogued as centaurs or trans - neptunian objects )
the fraction of retrograde objects is 2.7% ( 48 out of 1804 ) .
therefore , most retrograde asteroids have been found in the region of the giant planets and beyond . at 11.0%
( 30 out of 272 catalogued objects in the jpl small - body database ) , centaurs have the highest retrograde fraction among all the dynamical groups or asteroid families .
centaurs are a transient asteroidal population whose orbits cross those of the outer planets ( see e.g. di sisto and brunini 2007 )
. most centaurs may have their origin in the trans - neptunian belt , some may have come from the oort cloud .
the main source of centaurs remains to be discovered .
they are probably an heterogeneous population with multiple sources ( see e.g. di sisto et al .
2010 ) , perhaps none being dominant ( but see horner and lykawka 2010 ) .
the most intriguing of the known centaurs are those moving in retrograde orbits .
the study of retrograde centaurs , especially those following high inclination trajectories , may represent a rare opportunity to look at objects that have suffered little surface alteration over the age of the solar system ( sheppard 2010 ) .
the three largest known ( in terms of absolute magnitude , @xmath0 ) retrograde centaurs are ( 342842 ) 2008  yb@xmath3 ( @xmath0 = 9.3 mag ) , 2011  mm@xmath4 ( @xmath0 = 9.3 mag ) and 2013  lu@xmath5 ( @xmath0 = 8.1 mag ) ; they all move in near - polar orbits .
minor bodies following eccentric , high - inclination paths spend most of the time well beyond the plane of the solar system and away from the destabilizing influence of the giant planets , although close encounters with these planets are still possible at the nodes .    here , we show that , far from being just dynamical curiosities , the three largest retrograde centaurs may be the key to identify the size threshold for catastrophic collisions in the early solar system . on the other hand ,
they are also submitted to retrograde mean motion resonances that make them dynamically unstable , emphasizing the fact that these objects have not remained in their current orbits for long .
this paper is organized as follows . in sect .
2 , we present a comparative statistical analysis of the properties of known prograde and retrograde centaurs that reveals compelling reasons to single out the largest retrograde objects . the current dynamical status of these large retrograde centaurs is studied in sect .
3 , where the numerical model used in our @xmath6-body simulations is also briefly discussed . in sect .
4 , we consider the implications of our results and in sect .
5 we summarize our conclusions .
following alexandersen et al . ( 2013 ) , let us consider the group of minor bodies whose semi - major axes are between 6 and 34 au .
there are 266 such asteroids currently ( as of 2014 may 5 ) in the jpl small - body database .
this number is slightly less than the total number of known centaurs quoted above . out of them ,
236 ( 88.7% ) follow prograde orbits and 30 ( 11.3% ) follow retrograde orbits .
figure [ pro - ret ] shows the distribution in orbital parameter space and @xmath0 of prograde ( left - hand panels ) and retrograde ( right - hand panels ) centaurs . even if the number of known retrograde centaurs is small , it is highly improbable that the two centaur populations have a common origin .
the distribution in absolute magnitude , which is a proxy for size , clearly shows that retrograde centaurs are unlikely to be the result of gravitational scattering of originally prograde centaurs , unless we assume that smaller objects are more efficiently scattered to become retrograde than larger ones ( see the @xmath7- and @xmath8-panels ) or that smaller retrograde objects are easier to identify than their prograde counterparts .
most prograde orbits have eccentricities in the range 0.20.6 and relatively low inclinations , @xmath920 .
the small but noticeable group of prograde orbits with inclination @xmath1063 in fig .
[ pro - ret ] ( left - hand panel ) is consistent with predictions made by e.g. gallardo ( 2006 ) on the existence of objects submitted to the kozai resonance at that critical inclination
. the relative excess of objects with argument of perihelion close to 0or 180could be due to observational bias because only objects with declinations @xmath1124are observed ( see de la fuente marcos and de la fuente marcos 2014 ) .
the raw distribution in @xmath0 of prograde orbits shows evidence of the so called divot or sudden decrease in numbers for @xmath12 9 mag ( or a diameter of nearly 100 km ) as found by e.g. shankman et al .
this has been interpreted by those authors as the relic signature of the size distribution of the primordial planetesimal population , with the larger objects being the remnants .
a divot in the size distribution of trans - neptunian objects was predicted by fraser ( 2009 ) using the results of collisional evolution calculations .
a lack of objects in the size range 120 km was also anticipated by charnoz and morbidelli ( 2007 ) for minor bodies in the scattered disk and the oort cloud .
this feature had been previously identified and confirmed observationally for the size range 2560 km by bernstein et al .
( 2004 ) , fuentes and holman ( 2008 ) , and fraser and kavelaars ( 2009 ) .
this break in the size distribution may signal the onset of collisional evolution .
schlichting et al .
( 2013 ) consider that , for trans - neptunian objects , the observed size distribution above a radius of @xmath1030 km is primordial .
they also found that there is a significant deficit of bodies with size @xmath1010 km and a strong excess of bodies of nearly 2 km in radius ; this excess leaves a permanent signature in the size distribution of these objects that is still clearly observable after 4.5 gyr of collisional evolution .
fraser et al . (
2014 ) have found that the observed size distributions are difficult to explain within the currently accepted standard model of planetesimal formation .    the distribution in @xmath0 for retrograde orbits in fig .
[ pro - ret ] ( left - hand panel ) exhibits a most unusual behaviour .
there are only three objects ( 10% ) with @xmath1 10 mag : ( 342842 ) 2008  yb@xmath3 ( @xmath0 = 9.3 mag ) , 2011  mm@xmath4 ( @xmath0 = 9.3 mag ) and 2013  lu@xmath5 ( @xmath0 = 8.1 mag ) .
the rest have @xmath2 12 mag ( 73.3% have @xmath2 14 mag ) .
although the number of known retrograde centaurs is small , this result can be seen as a dramatic confirmation of the existence of a break in the size distribution as predicted by e.g. fraser ( 2009 ) .
there is an obvious deficit of objects with @xmath0 in the range 1014 mag although they must be significantly easier to detect , if they do indeed exist , than fainter objects of a size below 10 km .
although the available sample is far from complete and likely biased , it should not be biased in favour of 1-km objects against those in the size range 10100 km .
the scarcity of intermediate - size objects must be real and it may have a collisional , dynamical or even primordial origin ( or a combination of the three ) . in principle , and following fraser ( 2009 ) , schlichting et al .
( 2013 ) or shankman et al .
( 2013 ) , some or all the large retrograde centaurs may represent a remnant population of primordial planetesimals .
in addition , the distribution in @xmath0 of retrograde centaurs suggests that at least two distinct populations of rather different dynamical origin are present . in support of this interpretation , fig .
[ aecen ] ( see also fig .
[ real ] ) shows the distribution in the ( @xmath13 , @xmath7 ) plane of both groups and comets in the same semi - major axis range : the three largest retrograde centaurs clearly stand out ( bottom panel ) .
the brightest ( largest ) objects also appear to be dynamically different .
the distribution of semi - major axes and eccentricities of retrograde objects with @xmath2 12 resembles that of members of the detached or scattered disk populations ( levison and duncan 1997 ; gomes 2011 ) . in terms of the size of their members
, these trans - neptunian populations appear to be of collisional origin and , dynamically , they may have reached their current state after suffering strong gravitational perturbations .
the surprising abundance of retrograde minor bodies at @xmath14 mag may signal a high density of objects with size @xmath15 km as predicted by fraser ( 2009 ) or , more recently , belton ( 2014 ) .    ,
@xmath7 ) plane of prograde ( top panel ) and retrograde ( bottom panel ) centaurs and comets with semi - major axes in the range 634 au .
there are 145 prograde and 20 retrograde comets in that semi - major axis range .
see also fig .
the three largest retrograde centaurs , namely ( 342842 )  2008  yb@xmath3 , 2011  mm@xmath4 and 2013  lu@xmath5 , are clear dynamical outliers ; they have the largest perihelia . ]    , @xmath13 ) plane .
the dark gray area represents the eccentricity / semi - major axis combination with periapsis between the perihelion and aphelion of jupiter , the light gray area shows the equivalent parameter domain if saturn is considered instead of jupiter .
the orange area corresponds to the ( @xmath7 , @xmath13 ) combination with apoapsis between the perihelion and aphelion of uranus and the brown area shows the counterpart for neptune . ]
it may be argued that we unfoundedly invoked the observed decrease in the number of prograde centaurs at magnitude @xmath12 9 to implicitly justify the division into two groups ( @xmath1 10 mag and @xmath2 10 mag ) in fig .
[ aecen ] .
however , such concern is completely unjustified as fig .
[ aecen ] clearly shows .
it is rather obvious from the figure that the three largest retrograde centaurs are clear dynamical outliers because they have the largest perihelia . out of the three largest retrograde centaurs ,
only 342842 is a relatively well - studied object ( sheppard 2010 ; pinilla - alonso et al .
2013 ; bauer et al .
but , what is special or unique about these objects ?
both 342842 and 2013  lu@xmath5 can suffer close encounters with jupiter ; 2011  mm@xmath4 may be considerably more stable .
brasser et al .
( 2012 ) and volk and malhotra ( 2013 ) have found that retrograde or very high inclination centaurs can not have evolved from the trans - neptunian belt due to close encounters with the jovian planets .
instead , they favour a source in the oort cloud for these objects . in the solar system ,
near - polar orbits are almost dynamically frozen in time and , potentially , they could be stable for hundreds of myr .
this is particularly true if the nodes of the orbits of the objects remain relatively far from the paths of the planets .
however , if the objects following these unusual trajectories are also submitted to mean motion resonances , their dynamical stability could be compromised even if they do not undergo close encounters with planetary bodies because they may experience chaotic diffusion .
morais and namouni ( 2013a ) have brought to the attention of the astronomical community that minor bodies following retrograde orbits around the sun may also be trapped in a mean motion resonance with a planet .
the same authors have identified several asteroids currently moving in retrograde resonance with jupiter and saturn ( morais and namouni 2013b ) .
this discovery is of considerable practical interest because of its potential impact on our understanding of the early stages of the evolution of the solar system and also on the dynamics of transient populations like the centaurs .
are the largest retrograde centaurs perhaps submitted to retrograde mean motion resonances with the giant planets ?
the dynamics of resonant retrograde minor bodies has been studied by morais and namouni ( 2013a , 2013b ) . for an object in retrograde resonance with a planet that moves in a prograde orbit
, the resonant argument can be written as @xmath16 , where @xmath17 and @xmath18 are the mean longitudes of the retrograde object and the prograde planet , respectively , @xmath19 is the longitude of perihelion of the retrograde object , @xmath20 , @xmath21 is the longitude of the ascending node and @xmath22 is the argument of perihelion .
the mean longitude of a prograde object is given by @xmath23 , where @xmath24 is the mean anomaly , but the one for a retrograde object is written as @xmath25 .
the numbers @xmath26 , @xmath27 and @xmath28 are integers with @xmath29 . if the two objects are in resonance , @xmath30 , where @xmath31 and @xmath32 are the orbital periods of the retrograde object and the prograde planet , respectively .
if , over a given period of time , the resonant argument , @xmath33 , can take any value in the range 0360 , then it circulates but if it exhibits oscillations ( symmetric or asymmetric ) around a certain value ( usually 0or 180 ) , then we say that it librates .
if the resonant argument librates , then the object is trapped in a @xmath26:@xmath27 retrograde ( or @xmath26:@xmath34 ) mean motion resonance with the planet .
retrograde resonances are weaker than their prograde counterparts ( morais and namouni 2013a , 2013b ) .
0.35truecm    [ cols="<,<,<,<,<",options="header " , ]     [ elements ]    in this work , the orbital evolution of the three largest retrograde centaurs is computed using the hermite integration scheme described by makino ( 1991 ) and implemented by aarseth ( 2003 ) .
the standard version of this @xmath6-body sequential code is publicly available from the ioa web sitesverre / web / pages / nbody.htm ] .
our numerical simulations include the gravitational perturbations of the eight major planets , the moon , the barycentre of the pluto - charon system and the three largest asteroids ( dwarf planet ( 1 ) ceres , ( 2 ) pallas and ( 4 ) vesta ) .
initial positions and velocities are based on the de405 planetary orbital ephemerides ( standish 1998 ) referred to the barycentre of the solar system .
together with the orbital elements , they have been obtained from the jpl on - line solar system data service ( giorgini et al .
1996 ) . in the integrations ,
the relative error in the total energy was always as low as @xmath35 or lower after a simulated time of 1 myr .
the corresponding error in the total angular momentum is several orders of magnitude smaller .
further details on the physical model and the validation of the code can be found in de la fuente marcos and de la fuente marcos ( 2012a ) . besides the calculations performed using the nominal orbits
as provided by the jpl solar system data service , we have followed the evolution of 50 additional control orbits obtained from the nominal ones and the quoted uncertainties ( see table [ elements ] ) .
their orbital elements have been sampled from a six - dimensional gaussian distribution around the nominal ones ( within 3@xmath33 with the exception of 2013  lu@xmath5 ) . in all the figures ,
@xmath36 = 0 coincides with the jd2456800.5 epoch unless explicitly stated .
the dynamical evolution of the objects is followed for @xmath37500 kyr .    , 2011  mm@xmath4 and 2013  lu@xmath5 in the time interval ( -500 , 500 ) kyr .
relatively long episodes of horizontal ( resonant ) oscillations are visible . for 2011  mm@xmath4 ,
no large changes in the orbital elements are observed over the studied time interval .
the dark gray area represents the ( @xmath7 , @xmath13 ) combination with periapsis between the perihelion and aphelion of jupiter , the light gray area shows the equivalent parameter domain if saturn is considered instead of jupiter .
the orange area corresponds to the ( @xmath7 , @xmath13 ) combination with apoapsis between the perihelion and aphelion of uranus and the brown area shows that of neptune .
these results correspond to the nominal orbits in table [ elements ] . ]     but for the ( @xmath7 , @xmath8 ) evolution . ]
figures [ ea ] , [ ei ] and [ all ] evidence that the overall orbital evolution of the three objects is rather different although all of them experience horizontal oscillations in the ( @xmath7 , @xmath13 ) plane as secular resonances modify @xmath7 at constant @xmath13 .
these figures also show that , when @xmath22 is close to 0or 180 , ( 342842 )  2008  yb@xmath3 can experience close encounters with jupiter at perihelion and with uranus at aphelion , 2011  mm@xmath4 can encounter saturn at perihelion and neptune at aphelion , and 2013  lu@xmath5 can find jupiter at perihelion and neptune at aphelion ( see also fig .
[ real ] ) .
encounters can only take place at the nodes of the orbit .
for a retrograde orbit , the distance between the sun and the nodes ( fig .
[ all ] , f - panels ) is given by @xmath38 , where the ` + ' sign denotes the descending node and the ` - ' sign the ascending node . in the following ,
we explore the dynamics of these unusual objects in more detail .
asteroid ( 342842 ) 2008  yb@xmath3 was discovered on 2008 december 18 by r. h. mcnaught observing for the siding spring survey with the 0.5-m uppsala schmidt ( mcnaught et al . 2008 ) . when first observed , it had a @xmath39 magnitude of 18.1 .
it was soon obvious that it follows a near - polar retrograde orbit characterized by @xmath13 = 11.62 au , @xmath7 = 0.44 and @xmath8 = 105(see table [ elements ] ) .
its orbit is well determined with 466 observations spanning a data - arc of 1856 d. it is relatively large with @xmath0 = 9.3 mag or a diameter of 67 km for an albedo of 6.2% , although it is probably larger .
sheppard ( 2010 ) found moderately red colours for this minor body , consistent with those of the neutral group of the centaurs , and concluded that due to its relatively short perihelion , the object may have undergone surface sublimation and interior recrystallization .
pinilla - alonso et al . ( 2013 ) used photometry , spectroscopy and numerical simulations to show that this asteroid has a rotational period of 16.4 hr and , when compared with solar colours , is moderately red with low albedo ( perhaps depleted of ices ) , and it came from the trans - neptunian region .
asteroid 342842 is currently trapped in an accidental 3:10 retrograde resonance with jupiter ( libration around 0 , see fig .
[ yb3 ] , top panel ) .
we call it accidental because the resonant argument circulates with a superimposed short - period oscillation .
this behaviour is similar to the one found for plutino ( 15810 ) 1994 jr@xmath40 , a quasi - satellite of pluto ( de la fuente marcos and de la fuente marcos 2012b ) .
this resonant episode with jupiter started nearly 30 kyr ago and it will cease in about 10 kyr . just after the episode , this minor body will be submitted to a 3:4 retrograde resonance with saturn ( asymmetric libration around 0 , see fig .
[ yb3 ] , bottom panel ) for nearly 10 kyr .
close encounters with jupiter and more distant flybys with saturn trigger the ejection from the retrograde resonances with jupiter and saturn .
all the control orbits exhibit similar behaviour during the time - span plotted in fig .
[ yb3 ] but the dynamical evolution after the close encounters differs significantly between control orbits . in about half the cases , the 3:4 retrograde resonance with saturn is observed prior and after the 3:10 retrograde resonance with jupiter .
the time - scale necessary for two initially infinitesimally close trajectories associated with this object to separate appreciably , or @xmath7-folding time , is @xmath103 kyr .
the object will be ejected towards the trans - neptunian region within the next few hundred kyr .
our results indicate that this object is less dynamically stable than implied in pinilla - alonso et al .
( 2013 ) .    .
]      asteroid 2011  mm@xmath4 was discovered on 2011 june 24 by n. primak , a. schultz , s. watters , j. thiel and t. goggia observing for the pan - starrs 1 project with the 1.8-m ritchey - chretien telescope from haleakala ( woodworth et al .
when discovered , its @xmath41 magnitude was 21.1 .
it follows a path characterized by @xmath13 = 21.21 au , @xmath7 = 0.47 and @xmath8 = 100(see table [ elements ] ) .
its orbit is relatively well determined with 92 observations spanning a data - arc of 396 d. with @xmath0 = 9.3 mag , it is very similar in size to ( 342842 ) 2008  yb@xmath3 although probably larger ( bauer et al . 2013 ) .
asteroid 2011  mm@xmath4 is currently trapped in the 3:10 retrograde resonance with saturn ( asymmetric libration around 90 , see fig .
[ mm4 ] , top panel ) .
this resonant episode started nearly 14 kyr ago and it will end in about 3 kyr . prior to this episode , we observe that the critical angle alternates between circulation and asymmetric libration , indicating that the motion is chaotic . before and after the resonant interaction with saturn , this minor body was submitted to the 5:3 retrograde resonance with neptune ( asymmetric libration around 180 , see fig .
[ mm4 ] , bottom panel ) .
relatively close encounters with saturn and neptune are responsible for the ejection from the resonances . during the entire resonant phase ,
the argument of perihelion remains close to 0 , suggesting a behaviour similar to that of a kozai secular resonance ( kozai 1962 ) although both @xmath7 and @xmath8 remain relatively constant ( see fig .
[ all ] ) . as in the previous case , all the control orbits exhibit similar behaviour during the time - span plotted in fig .
[ mm4 ] . in some cases , the 3:10 retrograde resonance with saturn results in asymmetric libration around 180instead of 90 .
our integrations indicate that this object can remain dynamically stable for several myr although this issue will be studied in more detail later .
its @xmath7-folding time is @xmath1010 kyr though .
it is not submitted to any typical secular resonance ( see fig .
[ rellon ] ) .
the cause of the relative stability of this object in comparison with the others , is connected with the position of its nodes ( see fig . [ all ] , panel f , middle column ) . for a minor body moving in a near - polar orbit , close encounters with major planets
are only possible in the vicinity of the nodes .
when @xmath22 is close to 0or 180 , the perihelion takes place within the ecliptic plane with one node situated at aphelion and the other at perihelion ; however , if @xmath22 is close to 90 or 270 , the perihelion occurs well away from the plane . in the case of 2011  mm@xmath4 ,
when @xmath22 is 0or 180 , encounters with saturn and neptune take place at perihelion and aphelion ; if @xmath22 is 90or 270 , encounters with uranus are possible at both nodes in a manner similar to what is observed for 83982 crantor ( 2002 go@xmath42 ) , a well - studied uranian co - orbital ( de la fuente marcos and de la fuente marcos 2013 ) .
the period of this oscillation of the nodes is nearly 700 kyr . at present time ,
the closest approach to saturn is at 1.08 au but the hill radius of saturn is 0.41 au ; for neptune , the closest approach is at 0.94 au but its hill radius is 0.77 au .
this explains why its motion is rather stable , the planetary perturbations exerted on 2011  mm@xmath4 are currently not strong enough to change the orbit appreciably .    .
]    , of 2011  mm@xmath4 with respect to the giant planets : referred to jupiter ( @xmath43 ) , to saturn ( @xmath44 ) , to uranus ( @xmath45 ) , and to neptune ( @xmath46 ) . in all cases ,
the resonant argument @xmath47 circulates over the entire simulated time interval although it is nearly in secular resonance with jupiter .
these results are for the nominal orbit in table [ elements ] . ]
asteroid 2013  lu@xmath5 was discovered on 2013 june 8 by a. boattini , r. e. hill and j. a. johnson observing with the 1.5-m telescope of the mt .
lemmon survey ( bressi et al . 2013 ) .
it had a @xmath39 magnitude of 21.7 when first observed .
its orbit is very poorly determined with just 27 observations acquired during 7 d and it is mainly included here to encourage follow - up observations .
its orbit is , in principle , the most eccentric of the ones discussed here .
this object is also the largest , with @xmath0 = 8.1 mag or 50150 km .    due to the large uncertainty associated to the orbit of 2013  lu@xmath5 ,
our discussion here is more speculative . assuming an object moving in an orbit like the nominal one ( @xmath13 = 19.25 au , @xmath7 = 0.63 , @xmath8 = 11700 , @xmath21 = 25457 , @xmath22 = 17905 and @xmath24 = 29787 )
, we observe that it is trapped in the 1:1 retrograde resonance with uranus ( asymmetric quasi - satellite behaviour , see fig .
[ lu28 ] , middle panel ) .
the object has remained in the neighbourhood of this resonance for nearly 200 kyr but its dynamical status is rather unstable because it is submitted to additional resonances with neptune and saturn ( see fig .
[ lu28 ] , bottom and top panels ) , 2:1 and 4:11 retrograde , respectively .
our results are consistent with previous studies of the 1:1 resonance at high inclination ( see e.g. brasser et al .
all the control orbits suggest that it is a recent visitor from the trans - neptunian region . in this case , and given its very uncertain parameters , we have restricted our exploration of the orbital domain to orbits very close to the nominal one .
its @xmath7-folding time is @xmath103 kyr .    .
the distribution in absolute magnitude of known retrograde centaurs with semi - major axis in the range 634 au exhibits a clear deficit of objects within 1014 mag and a significant excess at @xmath1014 mag .
this is consistent with theoretical expectations based on collisional models ( fraser 2009 ) .
objects with @xmath1 10 mag form a distinct group in dynamical terms .
the largest known retrograde centaurs move in near - polar , comet - like orbits but no cometary activity has been detected yet on any of them .
the analysis of the short - term orbital evolution of these objects shows that the group is heterogeneous in origin .
all of them are currently trapped in transient retrograde resonances with the jovian planets .
compared with the objects analysed by morais and namouni ( 2013b ) , the duration of the resonant episodes of the objects studied here is longer as they are globally more stable .
asteroid ( 342842 ) 2008  yb@xmath3 will escape towards the trans - neptunian region in the near future and 2013  lu@xmath5 , a recent visitor from the trans - neptunian region , is the first known uranian retrograde co - orbital candidate ( a quasi - satellite ) but further observations are necessary to absolutely confirm that this object is indeed a retrograde companion to uranus .
another interesting feature of the dynamics discussed above is that high - order resonances seem to be preferred by these initial conditions characterized by very high inclinations .
this is , however , not surprising as pointed out by , for example , robutel and laskar ( 2001 ) or gallardo ( 2006 ) . the strength of high - order resonances in general is greater for high - inclination orbits . at high inclination ,
high - order resonances are strong enough to capture objects for relatively long periods of time . in our case ,
the high eccentricities are also a contributing factor to produce stronger dynamical effects ( see e.g. gallardo 2006 ) .
gallardo ( 2006 ) has found that , for a given eccentricity and due to short - period planetary perturbations , high - inclination resonant orbits have a higher probability of survival than low - inclination ones .
this is also consistent with the comparatively weaker stability observed for the retrograde resonances analysed by morais and namouni ( 2013b ) and pointed out above ; the orbits followed by the objects studied in their paper are not almost perpendicular to the ecliptic as in our case and their resonances are low order . in general ,
chaotic diffusion is enhanced for high - inclination , high - eccentricity orbits submitted to high - order resonances .
asteroids 342842 and 2013  lu@xmath5 are not long - term stable but 2011  mm@xmath4 may be considerably more stable . in order to study the long - term stability of 2011  mm@xmath4 we use the regularized mixed variable symplectic ( rmvs ) integrator ( _ swift - rmvs3 _ ) which is part of the swift package ( levison and duncan 1994 ) .
the swift package is available from h. levison web site.hal/swift.html ] the physical model here is similar to the one described in sect .
3 but now the earth and the moon are replaced by the barycentre of the system . the calculations do not consider mercury or any asteroids but they include the barycentre of the pluto - charon system , use a timestep of 50 d and run for @xmath375.0 gyr .
relative errors in the total energy are @xmath48 .
an ejection distance of 10000 au from the sun was used ; any test particle that reached that distance was removed from the calculations . as in the previous simulations , we have followed the evolution of 50 control orbits obtained from the nominal orbital elements of 2011  mm@xmath4 and their quoted uncertainties ( see table [ elements ] ) .
these orbital elements have been sampled within 3@xmath33 from a six - dimensional gaussian distribution around the nominal ones .
results from rmvs and hermite are consistent within the applicable timeframe .
horner et al .
( 2004a , 2004b ) studied the dynamical half - lives of prograde centaurs and found values in the range 0.5432 myr .
the mean half - life of their sample of centaurs was 2.7 myr . for a given orbit , a large number of control orbits , @xmath49 , with parameters close to those of the nominal one were integrated .
these orbits were followed until they resulted in a collision with the sun or the ejection of the test particle from the solar system .
the number of control orbits , @xmath6 , that still remain within the boundaries of the solar system ( in our case , within 10000 au from the sun ) after a time @xmath36 , is given by @xmath50 , where @xmath51 , @xmath17 is the decay constant and @xmath52 is the dynamical half - life of the objects .    using _ swift - rmvs3 _
, we have applied the same approach to 2011  mm@xmath4 in an attempt to study its long - term stability .
asteroid 2011  mm@xmath4 is clearly far more dynamically stable than most centaurs ; its mean dynamical half - life is 23.12 myr .
however , the range in dynamical half - lives forward in time is 5.26310.32 myr and for the past injection into its present orbit is 3.561795.35 myr . in other words ,
this object probably arrived to the region of the giant planets nearly 23 myr ago and certainly not earlier than nearly 2 gyr ago .
therefore , it is not a primordial minor body that has remained around its present location since the formation of the solar system like jovian or martian trojans .
it will be ejected also within 23 myr and not later than 310 myr into the future .    among the control orbits ,
we have found one that is particularly interesting ( @xmath13 = 21.1623 au , @xmath7 = 0.47301 , @xmath8 = 100.446 , @xmath21 = 282.604and @xmath22 = 7.224 , see fig .
[ control ] ) .
it was very long lived into the past ( 1795.35 myr ) and its associated dynamics may signal the presence of a reservoir of asteroids moving in high - inclination prograde orbits that are submitted to a kozai resonance as described by gallardo ( 2006 ) . in fig .
[ control ] , the control orbit begins at a critical inclination of @xmath1063and very high eccentricity .
from there , the eccentricity decreases and the initially prograde orbit turns into a retrograde one . if , somehow , a reservoir of objects moving with orbital inclinations @xmath1063and high eccentricities was created ( and may still survive ) at 200500 au in the remote past , it may be a viable secondary source for retrograde centaurs .
the dynamical evolutionary pathway uncovered by our calculations , even if it has a statistically low probability , shows that the oort cloud is not the only possible source region for retrograde centaurs .
the existence of additional stable reservoirs of rocky and icy bodies at 100300 au can not be simply casted aside .
this interpretation is , however , tentative and should be treated with caution , as there is recent circumstantial evidence on the existence of a massive planet at 250 au from the sun ( trujillo and sheppard 2014 ) . if such planet does indeed exist , the architecture of the entire trans - plutonian region could be very different from the one assumed here , i.e. no truly massive bodies beyond neptune .
a massive perturber in the region 100300 au may affect the dynamical evolution presented in figure [ control ] very significantly .
brasser et al .
( 2012 ) , and volk and malhotra ( 2013 ) argued that retrograde or very high inclination centaurs have their origin in the oort cloud .
fouchard et al . ( 2014 ) also pointed out that the oort cloud , and in particular its inner part , is a potential source of high inclination centaurs .
long - period comets are believed to have their origin in the oort cloud ( see e.g. duncan 2008 ) .
halley - type comets may come from the outer edge of the scattered disc ( levison et al .
2006 ) or the oort cloud ( emelyanenko et al . 2013 ) .
jupiter - family comets appear to have their origin in the scattered disk ( levison and duncan 1997 ) . for the group of minor bodies
whose semi - major axes are between 6 and 34 au , the fraction of retrograde orbits is 11.3% .
there are 165 known comets with semi - major axes in the range 634 au .
out of them , 20 objects or 12.1% move in retrograde orbits .
the fractions of retrograde orbits in both groups are more similar than expected by chance alone . in principle
, this fact argues strongly in favour of a dominant source in the oort cloud for the retrograde centaurs but , again , the possible existence of trans - plutonian planets may change the entire dynamical scenario .
: semi - major axis , @xmath13 , eccentricity , @xmath7 , inclination , @xmath8 , longitude of the perihelion , @xmath21 , and , argument of perihelion , @xmath22 .
the output time - step for this plot is 0.045 myr .
these results correspond to the control orbit defined by @xmath13 = 21.1623 au , @xmath7 = 0.47301 , @xmath8 = 100.446 , @xmath21 = 282.604and @xmath22 = 7.224 . ]
asteroid 2011  mm@xmath4 is a promising candidate to be a representative of a remnant population of primordial planetesimals and it is a good candidate to be submitted to the nodal libration mechanism ( verrier and evans 2008 , 2009 ; farago and laskar 2010 ; doolin and blundell 2011 ) , at least temporarily .
this mechanism , when active , can significantly enhance the stability of near - polar orbits .
minor bodies like 2011  mm@xmath4 may signal the size threshold for catastrophic collisions in the early solar system . returning to the issue of its origin , and as alternative to the oort cloud or the hypothetical high - inclination reservoir pointed out above , we speculate that this and other similar objects may be ancient ejected uranian satellites ( regular or irregular ) . the uranian satellite system is the least massive among those of the jovian planets even if the largest one , titania , is the eighth largest moon in the solar system .
uranus rotates nearly on its side with the axis lying nearly in its orbital plane .
uranus large moons follow regular paths in equatorial - plane orbits .
any object originally bound to uranus and tracing an orbit lying nearly in its equatorial plane will follow an almost polar orbit when ejected .
in this research , we have used a comparative statistical analysis of the properties of known prograde and retrograde centaurs to reveal compelling reasons to single out the largest retrograde objects .
the dynamical evolution of these objects has been subsequently studied using @xmath6-body simulations .
the main conclusions of our work can be summarized as follows :    * the distribution in absolute magnitude of known retrograde centaurs with semi - major axis in the range 634 au exhibits a clear deficit of objects in the range 10 - 14 mag and a significant excess at @xmath1014 mag . * the largest retrograde centaurs move in almost polar , very eccentric paths with their nodal points currently located near perihelion and aphelion . within the retrograde centaur population ,
they are clear outliers both in terms of dynamics and size . *
the largest retrograde centaurs are currently trapped in transient retrograde , mostly high - order , resonances with the jovian planets .
these resonances are strong enough to induce substantial chaotic diffusion and make their orbits dynamically unstable . *
asteroid ( 342842 ) 2008  yb@xmath3 appears to be significantly less dynamically stable than usually assumed . *
asteroid 2011  mm@xmath4 may be a visitor from the oort cloud but an origin in a relatively stable closer reservoir can not be ruled out . *
asteroid 2013  lu@xmath5 is the first known uranian retrograde co - orbital candidate ( a quasi - satellite ) .
* this population may have its origin in the oort cloud but alternative or concurrent sources could be ancient ejected uranian satellites or a hypothetical reservoir of minor bodies moving in high - inclination , high - eccentricity orbits at 200500 au and submitted to the kozai resonance .
our dynamical analysis has been restricted to just three objects and their associated orbital parameter space . however ,
a number of distinctive properties that should be common to any retrograde centaurs have been identified .
in particular , the role of high - order resonances at high inclination as anticipated by gallardo ( 2006 ) seems to be determinant on the short - term evolution of this population .
we would like to thank s.  j.  aarseth , h.  f.  levison and m.  j.  duncan for providing the codes used in this research , and the referee for his / her quick and helpful report .
this work was partially supported by the spanish ` comunidad de madrid ' under grant cam s2009/esp-1496 .
we thank m. j. fernndez - figueroa , m. rego fernndez and the department of astrophysics of the universidad complutense de madrid ( ucm ) for providing computing facilities .
most of the calculations and part of the data analysis were completed on the ` servidor central de clculo ' of the ucm and we thank s. cano alsa for his help during this stage . in preparation of this paper , we made use of the nasa astrophysics data system , the astro - ph e - print server and the mpc data server .
aarseth , s.j . : gravitational n - body simulations .
cambridge university press , cambridge , p.  27
( 2003 ) alexandersen , m. , gladman , b. , greenstreet , s. , kavelaars , j.j . , petit , j .-
m . , gwyn , s. : science  * 341 * , 994 ( 2013 ) bauer , j.m .
, et al . :  * 773 * , 22 ( 2013 ) belton , m.j.s . :  * 231 * , 168 ( 2014 ) bernstein , g.m . , trilling , d.e . ,
allen , r.l . , brown , m.e . ,
holman , m. , malhotra , r. :  * 128 * , 1364 ( 2004 ) brasser , r. , heggie , d.c . , mikkola , s. : celest .
* 88 * , 123 ( 2004 ) brasser , r. , schwamb , m.e . ,
lykawka , p.s . ,
gomes , r.s .
:  * 420 * , 3396 ( 2012 ) bressi , t.h .
, et al . :
mpec 2013-l69 ( 2013 ) charnoz , s. , morbidelli , a. :  * 188 * , 468 ( 2007 ) de la fuente marcos , c. , de la fuente marcos , r. :  * 427 * , 728 ( 2012a ) de la fuente marcos , c. , de la fuente marcos , r. :  * 427 * , l85 ( 2012b ) de la fuente marcos , c. , de la fuente marcos , r. :  * 551 * , a114 ( 2013 ) de la fuente marcos , c. , de la fuente marcos , r. :  * 441 * , 2280 ( 2014 ) di sisto , r.p . ,
brunini , a. :  * 190 * , 224 ( 2007 ) di sisto , r.p . ,
brunini , a. , de ela , g.c .
:  * 519 * , a112 ( 2010 ) doolin , s. , blundell , k. m. :  * 418 * , 2656 ( 2011 ) duncan , m. j. :  * 138 * , 109 ( 2008 ) emelyanenko , v.v . , asher , d.j . , bailey , m.e .
: earth , moon , and planets  * 110 * , 105 ( 2013 ) farago , f. , laskar , j. :  * 401 * , 1189 ( 2010 ) fouchard , m. , rickman , h. , froeschl , c. , valsecchi , g. b. :  * 231 * , 99 ( 2014 ) fraser , w.c . :  * 706 * , 119 ( 2009 ) fraser , w.c .
, kavelaars , j.j .
:  * 137 * , 72 ( 2009 ) fraser , w.c . ,
brown , m.e . ,
morbidelli , a. , parker , a. , batygin , k. :  * 782 * , 100 ( 2014 ) fuentes , c.i . , holman , m.j . :  * 136 * , 83 ( 2008 ) gallardo , t. :  * 181 * , 205 ( 2006 ) giorgini , j.d .
, et  al .
:  * 28 * , 1158 ( 1996 ) gomes , r.s . :  * 215 * , 661 ( 2011 ) horner , j. , lykawka , p.s . :  * 402 * , 13 ( 2010 ) horner , j. , evans , n.w . ,
bailey , m.e .
:  * 354 * , 798 ( 2004a ) horner , j. , evans , n.w . ,
bailey , m.e . :  * 355 * , 321 ( 2004b ) kozai , y. :  * 67 * , 591 ( 1962 ) levison , h. , duncan , m.j . :
* 108 * , 18 ( 1994 ) levison , h.f . ,
duncan , m.j . :  * 127 * , 13 ( 1997 ) levison , h.f . ,
duncan , m.j . ,
dones , l. , gladman , b.j .
:  * 184 * , 619 ( 2006 ) makino , j. :  * 369 * , 200 ( 1991 ) mcnaught , r.h . , et al . :
mpec 2008-y38 ( 2008 ) morais , m.h.m . ,
namouni , f. : celest .
* 117 * , 405 ( 2013a ) morais , m.h.m . ,
namouni , f. :  * 436 * , l30 ( 2013b ) pinilla - alonso , n. , et al . ;
* 550 * , a13 ( 2013 ) robutel , p. , laskar , j. :  * 152 * , 4 ( 2001 ) schlichting , h.e . ,
fuentes , c.i . , trilling , d.e .
:  * 146 * , 36 ( 2013 ) shankman , c. , gladman , b.j . , kaib , n. , kavelaars , j.j . , petit , j.m . :  * 764 * , l2 ( 2013 ) sheppard , s.s . :  * 139 * , 1394 ( 2010 ) standish , e.m . :
jpl planetary and lunar ephemerides , de405/le405 , interoffice memo .
312.f-98 - 048 , jet propulsion laboratory , pasadena , ca , usa ( 1998 ) trujillo , c.a . ,
sheppard , s.s .
:  * 507 * , 471 ( 2014 ) verrier , p.e . ,
evans , n.w . :  * 390 * , 1377 ( 2008 ) verrier , p.e . ,
evans , n.w . :  * 394 * , 1721 ( 2009 ) volk , k. , malhotra , r. :  * 224 * , 66 ( 2013 ) woodworth , d. , et al . : mpec 2013-n02 ( 2011 )