the perturbative nonrenormalizability of the theory of general relativity ( gr ) is one of the greatest obstacles against proposing quantum gravity theory that should be correct at all scales .
the latter belongs to the yet - unsolved problems in the fundamental physics @xcite .
weakness of the gravitational force is another difficulty preventing the quantum effects of gravity to be detectable at the currently - probed scales .
the recent detection of the gravitational waves is strictly limited to pulsars in binary star systems @xcite .
accordingly , the resulting gravitational waves become very great .    at a gaussian fixed point , gr does not satisfy the perturbative renormalizability because the gravitational coupling ( @xmath0 ) is a dimensioned quantity ( @xmath1 ) and accordingly the ultraviolet ( uv ) divergences are uncontrollable . by introducing higher terms of the scalar curvature ,
the action shall be cured . on the other hand
, this introduces an additional problem , for instance , to the uniqueness @xcite . from a minimal approach that will be elaborated below , which assures perturbatively renormalizable uv completion of gr as quantum field theory ( qft ) for @xmath2-dimensional metric field , one looses some properties of gr and creates additional degrees - of - freedom .
such an approach benefits from the success of qft in describing all forces ( except gravity so - far ) and the standard model of the elementary particles .
the scalar field theory , that was proposed by lifshitz in order to explain the quantum critical phenomena in condensed matter @xcite , inspired horava to propose a theory for the quantum gravity with an anisotropic scaling in uv .
this new approach is thus known as horava - lifshitz gravity ( hlg ) @xcite , where the dynamical critical exponent ( @xmath3 ) makes the theory non - relativistic .    in analogy to critical systems , the uv fixed point has been treated @xcite .
the hlg approach assures causal dynamical triangulations @xcite , renormalization group approaches based on asymptotic safety @xcite and symmetries of gr .
this approach was discussed in an extensive amount of literature @xcite , where problems such as , internal consistency and compatibility with observations are treated , as well .
the new dimensionless coupling @xmath4 is assumed to approach unity in the infrared ( ir ) limit with the argumentation that at the hlg - charactering parameter @xmath5 , hlg theory ca nt be reduced to gr .
for the sake of completeness , we mention that this is not always the case for all hlg models @xcite .
recent astrophysical observations obviously indicate that the universe expands with an accelerating rate , probably due to dark energy @xcite . besides inflation , a series of symmetry - breaking phase transitions , where topological defects may have been formed @xcite , were conjecture to take place . on the other hand
, the hlg - type approach @xcite proposes that the modification in gravity is responsible for the accelerated expansion of the universe .
furthermore , various modified theories for gravity can be implemented @xcite .
they successfully unify various inflation models with the late - time acceleration and cosmological observations .
but when endorsing any proposal about ( non)singular big bang due to controversial quantum gravity approaches such as generalized uncertainty principle ( gup ) and modified dispersion relation ( mdr ) @xcite should be seen as unwarranted claims .
the big bang singularity is a fundamental problem , which should be tackled by a good quantum gravity theory , such as hlg or even better .
the present paper is organized as follows . the cosmological evolution governed by a simple version of hlg
shall be introduced .
the horava - lifshitz gravity is introduced in section [ sec : hlg ] .
the field equations in friedmann - lemaitre - robertson - walker ( flrw ) cosmology are elaborated in section [ sec : fe ] .
section [ sec : nondet ] and section [ sec : det ] are devoted to horava - lifshitz gravity with non - detailed and detailed balance conditions , respectively . for the different equations - of - state ( eos ) , different solutions for the resulting differential equations shall be outlined . an extensive comparison between scale factor in cosmological radiation , matter , @xmath6cdm , de sitter
, chaplygin gas @xcite and quantum chromodynamical ( qcd ) eos shall be introduced .
the results and discussion are given in section [ sec : disc ] .
to section [ sec : conc ] the final conclusions are assigned .
the basic assumptions of horava - lifshitz gravity are antisotropic uv scalings between minkowskian space and time @xmath7 where @xmath3 is the dynamical critical exponent given in uv , @xmath8 is constant by which the scaling is performed and to the critical exponent the value @xmath9 is assigned . because the anisotropic scaling implies a preferred time coordinate , the usual @xmath2-dimensional lorentzian metric ca nt be the only fundamental structure on a pseudo - riemannian ( lorentzian ) manifold .
the original theory postulated @xmath3 to take the value assuring dimensionless gravitational coupling , which in turn guarantees power - counting renormalizability . at @xmath10 ,
the symmetries of gr are recovered . in order to single out a spacial ( temporal ) coordinate in a differentiable manifold , a codimension foliation on the manifold as a basic structure of the theory
was proposed @xcite .
the foliation structure fulfils local galilean invariance but implies impossible full - diffeomorphism invariance in gr .
the latter means the existence of an undo one - to - one mapping @xmath11 and its inverse besides their differentiabilities .
thus , gr may be considered as an emerging in an infrared fixed point @xcite .
the hlg is a projectable approach minimizing the number of independent couplings in the potential and adopting an extra principle to contract the potential , i.e. detailed balance conditions . in constructing the theory , a new set of symmetries
should be introduced .
imposing invariance under foliation - preserving diffeomorphisms , @xmath12 and @xmath13 would define field content , which paradoxically not always the one in gr and would make the theory violates the lorentz invariance principle . on the other hand , for arnowitt
, deser and misner [ adm ] formalism @xcite ( gr with global time foliation , @xmath14-dimensional metric of the spatial hypersurfaces gives all the dynamics ) , the metric is given as @xmath15 where @xmath16 , and @xmath17 ( @xmath18 ) being lapse function ( shift @xmath14-vector ) , which are the gauge fields of the diffeomorphism group .
the lapse variable is taken to be just time - dependent so that the projectability condition holds and by using the foliation - preserving difeomorphisms @xmath19 can be fixed to unity .    the variables in eq .
( [ adm1 ] ) are dynamical so that under the above scaling @xmath20 in terms of the metric eq .
( [ adm1 ] ) , the ricci scalar reads @xmath21 where @xmath22 , @xmath23 is the extrinsic curvature of spatial slices , @xmath24 is the spacial scalar curvature , and @xmath25 is a unit vector perpendicular to a hypersurface of constant time .
@xmath26 is the covariant derivative on the spatial slice .
the hlg action is decomposed of kinetic , potential and matter parts .    *
the kinetic action can be built with at most two time - derivatives of the metric .
this is a universal and most general invariant term .
the new dimensionless coupling @xmath4 , whose role in the theory is still a subject of debate , differentiates the hlg kinetic term from the gr one @xmath27 when combining this kinetic term with the line element and the anisotropic scaling of space and time , the scaling dimension of the gravitational coupling can be determined .
* the potential action term @xcite ( gr contains only two derivatives ) is affected by the choice between projectable and non - projectable theories @xcite @xmath28 .
\label{eq : termp}\end{aligned}\ ] ] the coefficients are given in just three parameters , @xmath29 , @xmath30 and @xmath31 @xcite @xmath32 where @xmath33 and @xmath30 are einstein coupling and cosmological constant , respectively .
the lower - order variables in expression ( [ eq : termp ] ) ; the cosmological constant and the ricci scalar , match the ones in gr . *
the matter part of the action @xcite is given as @xmath34 where @xmath35 is the lagrangian density of matter fields @xmath36 .      in ref .
@xcite , a simple scenario for the scale - invariant quantum fluctuations was introduced .
this was based on hlg without an additional scalar degree - of - freedom .
accordingly , it is believed that the inflation is self - consistently existing in the early universe . for this simple scenario ,
the detailed balance conditions were not necessary but the inflation scenario itself may be still or no longer needed , for instance , serious horizon problems remain unsolved , such as monopole and domain walls .
they still require inflation with slow - roll conditions .
later on many authors seem to disagree with ref .
the inflation was studied in hlg without the projectability conditions @xcite .
but , opposite to the proposal of ref @xcite ,
the linear scalar perturbations equations of the flrw universe are derived for a single scalar field .
a master equation of the perturbations has been specified for a particular gauge . the power spectrum and the spectrum index of the comoving curvature perturbations
have been determined .
it was noticed that the perturbations remain scale - invariance , and hlg without the projectability conditions is consistent with all current cosmological observations .
this is another solid support for adding scalar field(s ) .    in framework of nonrelativistic hlg with the projectability conditions and an arbitrary coupling constant @xmath4 ,
the inflation has been studied @xcite .
accordingly , the flrw universe without specifying the gauge is necessarily flat . but by adding a single scalar field , it was noticed that both metric and scalar field become strongly coupled and almost identically oscillating in sub - horizon regions . in super - horizon regions ,
the comoving curvature perturbation remains constant although the flrw perturbations become adiabatic .
furthermore , the perturbations of the slow - roll parameters , for instance , both scalar and tensor are found scale - invariant .
concrete tuning the coupling coefficients makes the spectrum index of the tensor perturbation identical as that in gr .
but the ratio of scalar to tensor spectra can be similar to that from gr and seems to depend on the spatial higher - order derivative terms .
furthermore , the emergence of finite - time future singularities has been studied in ref .
it was found that such singularities can be cured by adding a higher - order spatial derivative term .
recently , the general formulas for the inflationary power spectra of scalar and tensor are driven in the presence of a scalar field @xcite .
we conclude that in simple scenarios for scale - invariant cosmological perturbations with hlg , the inflation is not necessarily guaranteed .
for a cosmological context , we recall flrw metric , which is an exact solution of the einstein s field equations of gr , @xmath37 where @xmath38 is the scale factor and @xmath39 . accordingly , the homogeneous and isotropic metric
is given as @xmath40,\end{aligned}\ ] ] with the curvature constant @xmath41 , @xmath42 and @xmath43 represents open , flat or closed universe , respectively .
in hlg with non - detailed balance conditions and by varying @xmath19 and @xmath44 , the friedman equations can be extracted @xcite , @xmath45 , \label{13 } \\ ( 3 \lambda - 1 ) \left(\dot{h } + \frac{3}{2 } h^{2}\right ) & = & - \frac{\kappa^{4}}{4 } p - 3 \left[-\frac{\sigma}{6 } - k \frac{\gamma}{3 a^{2}}+ 2 k^{3 } \frac{9\ , \sigma_{1 } + 3\ , \sigma_{2 } + \sigma_{3}}{a^{6}}\right ] .
\label{14}\end{aligned}\ ] ] by substituting with the variables given in eqs .
( [ eq : variables ] ) , ( [ 13 ] ) and ( [ 14 ] ) , the freidmann equations read @xmath46   - \dfrac{\kappa^4\ , \mu^2\ , \lambda_{\omega}^2}{16(3 \lambda -1)^2}\right\}\ ; a + \dfrac{k^2\ , \kappa^4\ , \mu^2}{32(3 \lambda - 1)^2 } a^{-3}. \label{21b}\end{aligned}\ ] ]    based on equations of state ( eos ) describing the cosmological background geometry , this differential equation can be solved . from the continuity equation @xmath47 and with the eos @xmath48
, we get @xmath49 and @xmath50 , where @xmath51 - in natural units - is a dimensionless quantity .
it is related to the speed of sound squared within the system of interest .
then , eq . ( [ 21b ] ) can be expressed as @xmath52      equation ( [ 20 ] ) can be rewritten as & = & - h^2 - p - l + , [ 20pert ] where @xmath53 $ ] and @xmath54 $ ] . assuming perturbation in eq .
( [ 20pert ] ) leads to + h & = & - h^2 - 3 h h-
p - l + - 2 .
[ 20pert2 ] for an eos in which @xmath55 and for the first - order perturbation , we get h & = & - 3 h h - 2 .
[ 20pert3 ] however from the continuity equation , eq .
( [ eq : rhodot ] ) , & & = -3 ( 1 + ) . then h & = & m h + n , + & = & p h + q , with @xmath56 , @xmath57 $ ] , @xmath58 , and @xmath59 .
thus , h - ( m + + q ) h - ( - m + p n - q m ) h & = & 0 , for @xmath60 , the term with the first - order perturbation reads @xmath61 is apparently positive .
this gives a stable solution irrespective to the third term in left - hand side .
in describing the cosmic geometry , we wanted to analyse systematically the effects of implementing various equations - of - state .    * when the cosmological geometry is dominantly filled with matter ( dust approximation ) , @xmath62 , eq .
( [ 19 ] ) can be rewritten as + @xmath63 + when assuming that the cosmological constant @xmath30 and einstein coupling @xmath64 are negligibly small , eq . ( [ eq : add0 ] ) gets the solution @xmath65 where @xmath66 $ ] and @xmath67 $ ] .
deducing scale factor does not restrict the applicability of the present analysis .
other cosmological quantities such as hubble and deceleration parameters are good examples . * for cosmological geometry filled with radiation ( perfect fluid approximation ) , @xmath68 , eq .
( [ 19 ] ) becomes @xmath69 which has the following solution @xmath70 where @xmath71 - k^2 \kappa^4 \mu^2/[32(3 \lambda - 1)^2]$ ] and @xmath72 $ ] and @xmath73 .
* for cosmological geometry filled with cold dark - matter , @xmath74 , eq . (
[ 19 ] ) is written as , @xmath75 the solution looks similar to the one in eq .
( [ solsim ] ) @xmath76 * for cosmological geometry filled with vacuum ( dark ) energy ( cosmological constant ) , @xmath77 , then eq . (
[ 19 ] ) can be written as @xcite , @xmath78 then , the solution reads @xmath79 * for cosmological geometry filled with chaplygin gas @xcite , which is a hypothetical ( non - baryonic ) substance satisfying @xmath80 , where @xmath81 is positive and @xmath82 ranges between @xmath42 and @xmath43 , i.e. negative pressure , and does not cluster , gravitationally .
it was named after sergey chaplygin ( 1869 - 1942 ) , and expresses energy of quantum vacuum , and gravity modification .
it might signal extra - dimensions and is conjectured as a candidate for dark energy .
the energy density and pressure , respectively , read @xmath83^{1/(1+\alpha ) } , \label{eq : chplp}\\ p & = & -a\ , \left[a+\frac{b}{a^{3(1+\alpha)}}\right]^{-\alpha/(1+\alpha ) } , \label{eq : chple}\end{aligned}\ ] ] where @xmath84 is an integration constant .
+ in flrw cosmology , the chaplygin equations of state , eqs .
( [ eq : chplp ] ) and ( [ eq : chple ] ) , are substituted in eqs .
( [ 19 ] ) and ( [ 21b ] ) to give @xmath85 \right.{\nonumber}\\ & & \left .
- \dfrac{\kappa^4\ , \mu^2\ , \lambda_\omega^2}{16(3\lambda-1)^2}\right\}\ , a + \dfrac{k^2\ , \kappa^4\ , \mu^2}{32(3\lambda-1)^2}\ , a^{-3}. \label{eq : daachapa0}\end{aligned}\ ] ] + when ignoring the taylor terms with orders @xmath86 , then @xmath87 + the possible solution for this differential equations depends on the choices of @xmath81 and @xmath82 . here , we elaborate different cases : * * first for a modified chaplygin gas @xcite with @xmath88 and by assuming that @xmath89 , then eq . ( [ eq : daachapa1 ] )
is reduced to @xmath90 - \dfrac{\kappa^4\ , \mu^2\ , \lambda_{\omega}^2}{16(3\ , \lambda-1)^2}\right\}\ , a + \dfrac{k^2\ , \kappa^4\ , \mu^2}{32(3\lambda-1)^2}\ , a^{-3 } , \label{eq : addchapb1 } \\
\dot{a}^2 & = & \left[\dfrac{\kappa^2\ , a^{3/4}}{6(3\lambda-1 ) } - \dfrac{\kappa^4\ , \mu^2\ , \lambda_{\omega}^2}{16(3\ , \lambda-1)^2}\right ] a^2 + \left[\dfrac{\kappa^2\ , b a^{-1/4}}{8(3\lambda-1 ) } - \dfrac{k^2\ , \kappa^4\ , \mu^2}{16(3\lambda-1)^2 } \right ] a^{-2 } { \nonumber}\\ & + & \dfrac{\kappa^4\ , k \mu^2 \lambda_{\omega}^2}{8 ( 3 \lambda - 1)^2},\end{aligned}\ ] ] which has a solution @xmath91,\end{aligned}\ ] ] where @xmath92 , { \nonumber}\\ b_3 & = & \left[\dfrac{\kappa^2\ , b a^{-1/4}}{8(3\lambda-1 ) } - \dfrac{k^2\ , \kappa^4\ , \mu^2}{16(3\lambda-1)^2 } \right ] .
{ \nonumber}\end{aligned}\ ] ] * * for matter , the pressure vanishes or @xmath93 .
( [ eq : daachapa1 ] ) is reduced to @xmath94 which has a solution @xmath95 * for cosmological geometry filled with qcd matter , the eos can be deduced from the recent lattice qcd simulations @xcite , where @xmath96 ranges between @xmath97 and @xmath98 . for sake for completeness , we highlight the importance of precise estimation of eos from heavy - ion collisions @xcite . by substituting the deduced values of @xmath51 in eq .
( [ 21a ] ) , the resulting differential equation turns to be unsolvable . for simplicity , the first term in eq .
( [ 21a ] ) can be approximated , @xmath99\ ; a^{-3 } - \dfrac{\kappa^4 \mu^2 \lambda_w^2}{16(3 \lambda -1)^2}\ ; a , \\
\dot{a}^2 & \approx & \left [ \frac{\kappa^2\ , \rho_0}{16 ( 3 \lambda -1 ) } - \dfrac{k^2 \kappa^4 \mu^2}{32(3 \lambda - 1)^2}\right]\ ; a^{-2 } - \dfrac{\kappa^4 \mu^2 \lambda_w^2}{16(3 \lambda -1)^2}\ ; a^2 .
\end{aligned}\ ] ] similar to eq .
( [ eq : awm1 ] ) , the solution reads @xmath100}{2\ , \sqrt{b_2}\ ; { \cal a}},\end{aligned}\ ] ] where @xmath101^{1/2}$ ] .
the detailed balance condition is a technical trick proposed to reduce the couplings .
furthermore , it was argued that hlg with detailed balance conditions has to be broken in order to enable the theory to be compatible with the observations @xcite .
but , according to horava s remarks , its restrictive power remains useful @xcite .
there is a possible connection between detailed balance conditions and the entropic origin of gravity @xcite .
the detailed balance condition is often believed to be abandoned with the aim to obtain an ultraviolet stable scalar field in the theory .
but because of its several attractive features , we wanted it to be implemented in constructing the potential instead of effective field theory .
many authors wanted to improve it by adding extra term to the action of the theory that softly violates the detailed balance conditions @xcite as an attempt to assure a more realistic theory in its infrared - limit .
it was shown @xcite that the properties of the extra scalar degrees - of - freedom are the reasons why detailed balance conditions do nt lead to gr in the infrared - limit .    from perturbation study of detailed balance conditions
@xcite , two different strong coupling problems have been identified .
the first one explains why the gr solutions are typically not recovered .
the second is not necessarily associated with detailed balance conditions but refers to breaking of diffeomorphism invariance , which is required for anisotropic scaling in uv .
the authors of ref .
@xcite claimed to uncover an additional mode in the detailed balance conditions satisfying an equation of motion that is of first - order in time derivative .
but , the proposed mode leads to very fast exponential instabilities at short distances and becomes strongly coupled at an extremely low cutoff scale .
for hl gravity with detailed balance conditions , the field equations read @xcite @xmath102 the continuity equations @xmath103 can be substituted in and the friedmann equation can be derived from the action @xmath104 where @xmath105 , \\ \ell_1 & = & \sqrt{g }
n \left[\dfrac{\kappa^2 \mu^2 ( 1- 4\lambda)}{32(1 - 3\lambda ) } r^2 - \dfrac{\kappa^2}{2 w^4 } \left(c_{ij } - \dfrac{\mu w^2}{2 } r_{ij } \right ) \left ( c^{ij } - \dfrac{\mu w^2}{2 } r^{ij } \right ) \right ] , \end{aligned}\ ] ] and @xmath4 , @xmath64 , @xmath29 , @xmath51 and @xmath30 are constants , and @xmath106 being ricci scalar ( tensor ) .
@xmath107 is the cotton tensor , @xmath108 .
when comparing hlg @xmath109 to that of gr in adm formalism , then the speed of light and gravitational constant , respectively , reads @xmath110 when assuming loss of generality , @xmath111 . in solving eqs .
( [ 26 ] ) and ( [ 27 ] ) four cases are likely to occur :    * if @xmath112 , and @xmath113 @xmath114 where @xmath82 is constant of integration ( arbitrary constant ) . *
if @xmath115 , and @xmath116 @xmath117 * if @xmath118 @xmath119 * if @xmath115 , and @xmath113 @xmath120 with @xmath121 .    the solutions of the differential equations ( [ 26 ] ) and ( [ 27 ] ) depend on the eos charactering the cosmological background geometry .
this shall be elaborated in the section that follows .      at finite cosmological constant
, the friedmann equations read @xcite @xmath122.\end{aligned}\ ] ] with continuity equation , this equation becomes @xmath123 due to various equations - of - state , this differential equation can be solved .
* at @xmath62 , eq . ( [ eq : bdhrgadd ] ) can be reduced to @xmath124 then , the solution reads @xmath125 where @xmath126 $ ] and @xmath127 $ ] .
* at @xmath68 , we obtain @xmath128 which can be solved as follows .
@xmath129 , \label{eq : dtl2}\end{aligned}\ ] ] where @xmath130 and @xmath131   + k^2/[\lambda_{\omega } ( 3\ , \lambda-1)]$ ] an d @xmath132 . * at @xmath74 , it results is @xmath133 , \label{eq : dtldt3}\end{aligned}\ ] ] with the solution @xmath134 , \label{eq : dtl3}\end{aligned}\ ] ] where @xmath135 $ ] .
* at @xmath77 , we get @xmath136 with the solution @xmath137 , \label{eq : dtl4}\end{aligned}\ ] ] * for the generalized chaplygin gas @xcite , it becomes @xmath138 , \label{eq : dtl5}\end{aligned}\ ] ] where @xmath139 $ ] . + in a matter - dominated universe , i.e. @xmath140 or @xmath93 , then @xmath141 and @xmath142 similar to eq .
( [ eq : dtldt1 ] ) , @xmath30 is assumed to be negligibly small @xmath143 where @xmath144 $ ] .
* for cosmological geometry filled with qcd matter , the eos can be deduced from the recent lattice qcd simulations @xcite . as given in fig .
[ fig : peqcd ] , @xmath96 ranges between @xmath97 and @xmath98 . for simplicity , the second term in eq .
( [ eq : bdhrgadd ] ) can be approximated , @xmath145\ , a^{-3 } + \frac{\lambda_{\omega}}{3\ , \lambda-1}\ , a. \end{aligned}\ ] ] similar to eq .
( [ eq : dtl2 ] ) , the solution reads @xmath146}{2\ , \sqrt{d_2}\ ; { \cal b}_7 } , \label{eq : dtl2}\end{aligned}\ ] ] where @xmath147^{1/2}$ ] and @xmath130 and @xmath148 $ ] .
based on studying a @xmath4-dependent version of the friedmann equations , it was concluded that hlg becomes an attractive gravity theory if the generalized wheeler - dewitt metric @xmath149 has an indefinite signature .
but , it becomes a repulsive gravity theory when the metric gets a positive sign .
these results are also found in our calculations . the various equations - of - state lead to different behavior of the scalar factor as an example about the cosmological quantities .    in fig .
[ fig:1 ] , the cosmic time ( @xmath150 ) is given in dependence on the scale factor ( @xmath151 ) . in these calculations
all parameters are kept fixed , @xmath152 and @xmath153 , @xmath154 and @xmath155 .
correspondingly , neither @xmath150 nor @xmath151 is given in physical units .
we compare between hlg with non - detailed and detailed balance conditions and gr gravity for various equations - of - state .
it is worthwhile to notice that flrw which is characterized by gr gravity is non - singular for four types of equations of state ; @xmath68 ( a ) , @xmath62 ( b ) , @xmath74 ( c ) and @xmath77 ( d ) .
the singularity in the flrw cosmology , which can be characterized by hlg with non - detailed and detailed balance conditions , varies due to the equations of state .
such conclusions are only valid for the given parameters .
furthermore , the dependence of @xmath150 on @xmath151 apparently varies with the equation of state , as well .
again , such conclusions are explicitly valid for the given parameters .
in other words , absence or presence of singular solutions and even the dependence of @xmath150 on @xmath151 strongly depend on the equation of state and the gravity itself .
the non - singular solutions is obvious , as @xmath30 is finite .
the top - left panel shows the results at @xmath68 ( a ) , i.e. radiation - dominated ( perfect fluid approximation ) flrw universe .
the hlg with non - detailed balance conditions , eq .
( [ eq : ndwp13 ] ) leads to a very rapid increase in @xmath151 , almost a parabolic shape . at a certain value of @xmath150 ,
the scale factor became @xmath150-independent , i.e. @xmath151 diverges .
the results from hlg with detailed balance conditions , eq .
( [ eq : dtl2 ] ) , and gr , eq .
( [ eq : w13std ] ) , are also compared with .
the bottom - left panel shows the results at @xmath62 ( b ) , i.e. matter - dominated ( dust approximation ) flrw universe .
again , the hlg with non - detailed , eq . ( [ eq : add0a ] ) , and detailed balance conditions , eq .
( [ eq : bhlg1 ] ) , lead to a rapid increase in @xmath151 ; an almost a linear dependence . also here , the gr gravity , eq .
( [ eq : egw1 ] ) , has a non - singular solution .
the bottom - right panel presents the results at @xmath77 ( d ) , i.e. vacuum or dark energy or finite cosmological constant .
we find that the validity of hlg with non - detailed balance conditions is very limited .
also , hlg with detailed balance conditions shows a very rapid change in @xmath151 against @xmath150 .
the gr gravity results in an inflationary phase followed by a slower expansion .    despite the non - singular solutions ,
we notice that the expansion depicted in the right - hand panel is faster than the one illustrated in the left - hand panel .
again , the universe seems to start up its evolution from a very large scale factor .
figure [ fig:2 ] presents the results from chaplygin ( a ) and qcd ( b ) equations of states .
the results of the chaplygin gas strongly depends on the parameter @xmath81 .
the hlg with detailed balance conditions results in an almost linear dependence of @xmath151 and @xmath150 .
the non - singular solutions appear in both hlg approaches and equations of state .
we conclude that the difference between the chaplygin and qcd results is due to non - baryonic and baryonic equations - of - state , respectively .
while the present article was under review , we have completely conducted and published a research paper in which the proposed solutions are confronted to recent planck and bicepii observations @xcite . with single scalar field potentials describing power - law and minimal - supersymmetrically extended inflation , we have derived possible modifications in the friedmann equations . for various eos
, the dependences of the tensorial and spectral density fluctuations ( and their ratio ) on the inflation field are characterized .
the tensorial - to - spectral density fluctuations are calculated with varying spectral index .
it was found that , they decrease when moving from hlg with non - detailed balance conditions , to friedmann gravity , to hlg without the projectibility conditions , and to hlg with detailed balance conditions .
such a pattern remains valid for various eos and different inflation potential models .
the calculations fit well with the recent planck observations @xcite .
even , when studying the consequences of the quantum fluctuations on our understanding of landau - raychaudhuri equations , we found that cosmic eos play an essential role @xcite .
a theory for quantum gravity with an anisotropic scaling in uv based on a scalar field theory was proposed by horava , who was inspired by lifshitz theory for the quantum critical phenomena in condensed matter physics .
horava benefited from the success of qft in describing all forces ( except gravity ) and the standard model for the particle physics .
we study the impacts of various equations - of - state on flrw cosmology based on horava - lifshitz gravity .
we also implement the so - called hlg with detailed balance conditions , which was thought as a technical trick to reduce couplings in the theory . from the possible connection between detailed balance conditions and the entropic origin of the gravity
, the constructed potential can be used instead of adopting an effective field theory . both are compared with the gr gravity .    in this work ,
we have studied the dependence of scale factor ( @xmath151 ) and cosmic time ( @xmath150 ) in the context of horava - lifshitz gravity in early universe by using different equations - of - state .
we have compared the results with the gr gravity .
remarkable differences from what is predicted by gr are found .
also , a noticeable dependence on the equations of state is observed in singular and non - singular big bang .
these observations are precautionarily corresponding to the fixed parameters @xmath152 and @xmath153 , @xmath154 and @xmath155 .
the results presented in present paper are strongly depending on these parameters , for instance , the non - singular solutions can be to a large extend understood due to finite cosmological constant , @xmath30 .
the present work presents a systematic analysis for the equations of state characterizing the cosmic geometry and proposes a link between the conventional gravity by gr and its horava - lifshitz counterparts .
the latter are power - counting renormalizable quantum theories which assure causal dynamical triangulations , renormalization group approaches based on asymptotic safety and symmetries of gr .
the horava - lifshitz gravity introduces a new set of symmetries imposing invariance under foliation - preserving diffeomorphisms .
various types of equations of state are implemented .    in future works , we intend to extend this study to cover various topics , such as comparison with other theories such as the gravity s rainbow @xcite , and determining the possible impacts of the generalized uncertainty principle @xcite .
also , we plan to study third quantization of this theory @xcite , such as in a simple cosmological model @xcite .
we conclude that hlg can explain various epochs in the early universe .
furthermore , hlg might be able to reproduce the entire cosmic history with and without singular big bang .
it intends to study the stability of the given flrw solutions and then analyse the universe evolution scenarios , especially that hlg is conjectured to be a candidate for a quantum field theory of gravity , despite its yet - unresolved problems .
for eos @xmath48 and friedmann equations with finite cosmological constant , the continuity equation reads @xmath156 accordingly , we get following solutions @xmath157}{3\ , \sqrt{b}\ , { \cal c}_1 } , \label{eq : egw1 } \\
\mathrm{at}\ ; \omega=\frac{1}{3},\;\;\ ; t & = & \frac{1}{\sqrt{b } } \ln \left[a^2 - \dfrac{2a}{b } \right ] , \label{eq : w13std } \\ \mathrm{at}\ ; \omega=-\frac{1}{3},\;\;\ ; a(t ) & = & c_1 \exp\left(\sqrt{b}\ , t\right ) + c_2 \exp\left(-\sqrt{b}\ , t\right ) , \\ \mathrm{at}\ ; \omega=-1:\;\;\ ; a(t ) & = & c_3 \exp\left(\sqrt{c}\ , t\right ) + c_4 \exp\left(-\sqrt{c}\ , t\right ) , \\
\mathrm{and}\ , \mathrm{for}\ , \mathrm{generalized}\ , \mathrm{chaplugin}\ , \mathrm{gas},\ ; a(t ) & = & c_5 \exp\left(\sqrt{d}\ , t\right ) + c_6 \exp\left(-\sqrt{d}\ , t\right ) , \ ] ] where @xmath158 , @xmath159 , @xmath160 , and @xmath161 , @xmath162 .
@xmath163 are integration _ arbitrary _ constants .
from the recent lattice qcd simulations @xcite , we illustrate in fig . [
fig : peqcd ] , the pressure in dependence on the energy density . as discussed in ref .
@xcite , the ultimate goal of the high - energy experiments is first - principle determination of the underlying dynamics in the strongly interacting matter . determining the thermodynamical quantities ,
which are well - known tools to describe nature , degrees of freedom , decomposition , size and even the overall dynamics controlling evolution of the medium from which they are originating , is therefore very essential . for the thermodynamical pressure , an approximate attempt utilizing the higher - order moments of the particle multiplicity seems to be promising @xcite . on the other hand , determining the thermodynamical energy density and even relating the bjorken energy density to the lattice energy density depends on lattice qcd at finite baryon chemical potential and first - principle estimation of the formation time of the quark - gluon plasma ( qgp ) .
the energy density can be deduced from the derivative of the free energy with respect to inverse temperature .
this would explain the need to implement another variable different than the one responsible for the derivatives of the high - order moments .
e. thomas and a. r. zhitnitsky , _  topological susceptibility and contact term in qcd . a toy model  _ , phys .
d * 85 * , 044039 ( 2012 ) .
t. zhu , f .- w .
shu , q. wu , and a. wang , _  general covariant horava - lifshitz gravity without projectability condition and its applications to cosmology  _ , phys .
d * 85 * , 044053 ( 2012 ) .
d. vernieri and t p sotiriou , _  horava - lifshitz gravity : detailed balance revisited  _ , phys .
d * 85 * , 064003 ( 2012 ) .
m. saravani , n. afshordi , and r. b. mann , _  dynamical emergence of universal horizons during the formation of black holes  _ , phys .
d * 89 * , 084029 ( 2014 ) .
f. kheyri , m. khodadi , and h. r. sepangi , _  horava - lifshitz cosmology , entropic interpretation and quark - hadron phase transition  _ , annals phys . * 332 * , 75 ( 2013 ) .
r. loll and l. pires , _  role of the extra coupling in the kinetic term in hoava - lifshitz gravity  _ , phys . rev .
d * 90 * , 124050 ( 2014 ) .
a. de felice and s. tsujikawa , _  inflationary gravitational waves in the effective field theory of modified gravity
d * 91 * , 103506 ( 2015 ) .
m. cruz , m. gonzalez - espinoza , j. saavedra , and d. vargas - arancibia , _  scalar perturbations of two - dimensional horava - lifshitz black holes  _ , eur .
j. c * 76 * , 75 ( 2016 ) .
a. maciel , _  quasilocal approach to general universal horizons  _ ,
d * 93 * , 104013 ( 2016 ) .
m. chaichian , a. ghalee , and j. kluson , _  restricted f(r ) gravity and its cosmological implications  _ , phys .
d * 93 * , 104020 ( 2016 ) .
j. ambjorn , j. jurkiewicz , and r. loll , _  quantum gravity as sum over spacetimes  _ , _ lect .
notes phys . _
* 807 * , 59 - 124 ( 2010 ) .
d. f. litim , _  fixed points of quantum gravity  _ , _ phys . rev
. lett . _ * 92 * , 201301 ( 2004 ) .
m. reuter and f. saueressig , _  renormalization group flow of quantum gravity in the einstein - hilbert truncation  _ , _ phys .
d _ * 65 * , 065016 ( 2002 ) .
s. weinfurtner , t. p. sotiriou , and m. visser , _  projectable horava - lifshitz gravity in a nutshell
_ , _ j. phys .
ser . _ * 222 * , 012054 ( 2010 ) .
t. p. sotiriou , _  horava - lifshitz gravity : a status report  _ , _ j. phys .
ser . _ * 283 * , 012034 ( 2011 ) .
m. visser , _  status of horava gravity : a personal perspective  _ , _ j. phys .
ser . _ * 314 * , 012002 ( 2011 ) .
luis piresa , master thesis on _  horava - lifshitz gravity : hamiltonian formulation and connections with cdt  _ , institute of theoretical physics at utrecht university , ( 2012 ) .
abdel nasser tawfik and abdel magied diab , _
`` generalized uncertainty principle and recent cosmic inflation observations '' _ , _ electron . j. theor .
phys . _ * 12 * , 9 - 30 ( 2015 ) .
f. kheyri , m. khodadi , and h. r. sepangi , _  horava - lifshitz early universe phase transition beyond detailed balance  _ , _ eur .
j. c _ * 73 * , 2286 ( 2013 ) .    s. nojiri and s. d. odintsov , _
`` introduction to modified gravity and gravitational alternative for dark energy '' _ , _ int .
phys . _ * 4 * , 115 ( 2007 ) .
m. chaichian , sh .
nojiri , s. d. odintsov , m. oksanen , and a. tureanu , _ class .
_ * 27 * , 185021 ( 2010 ) ; erratum - ibid .
* 29 * , 159501 ( 2012 ) .
a. tawfik and a .diab , _  generalized uncertainty principle : approaches and applications
* 23 * , 1430025 ( 2014 ) .
a. kamenshchik , u. moschella , and v. pasquier , _  an alternative to quintessence
_ , _ phys .
b _ * 511 * , 265 ( 2001 ) .
t. p. sotiriou , m. visser , and s. weinfurtner , _  phenomenologically viable lorentz - violating quantum gravity  _ , _ phys .
_ * 102 * , 251601 ( 2009 ) .
t. p. sotiriou , m. visser , and s. weinfurtner , _  quantum gravity without lorentz invariance
_ , _ jhep _ * 0910 * , 033 ( 2009 ) .
p. horava .
_ `` quantum gravity at a lifshitz point '' _ , _ phys .
d _ * 79 * , 084008 ( 2008 ) .
d. giulini and c. kiefer , _  wheeler - dewitt metric and the attractivity of gravity
_ , _ phys .
a _ * 193 * , 21 ( 1994 ) .
shinji mukohyama , _  scale - invariant cosmological perturbations from horava - lifshitz gravity without inflation
_ , _ jcap _ * 0906 * , 001 ( 2009 ) .
t. zhu , y. huang , and a. wang , _  inflation in general covariant hoava - lifshitz gravity without projectability  _ , _ jhep _ * 1301 * , 138 ( 2013 ) . y. huang , a. wang , and q. wu , _
 inflation in general covariant theory of gravity
_ , _ jcap _ * 1210 * , 010 ( 2012 ) .
e. elizalde , s. nojiri , s. d. odintsov , and d. saez - gomez , _  unifying inflation with dark energy in modified f(r ) horava - lifshitz gravity  _ , _ eur .
j. c _ * 70 * , 351 ( 2010 ) .
r. kase and s. tsujikawa , _  effective field theory approach to modified gravity including horndeski theory and hoava - lifshitz gravity  _ , _ int .
* 23 * , 1443008 ( 2015 ) .
a. d. polyanin and v. f. zaitsev , _  exact solutions for ordinary differential equations  _ , ( chapman & hall / crc , london , 2003 ) .
s. borsanyi , z. fodor , ch .
hoelbling , s. d. katz , s. krieg , and k. k. szabo , _  full result for the qcd equation of state with 2 + 1 flavors  _ , _ phys .
b _ * 730 * , 99 - 104 ( 2014 ) . abdel nasser tawfik , _  an estimate of the thermodynamic pressure in high - energy collisions  _ , _ int
* 30 * , 1550027 ( 2015 ) .
e. p. verlinde , _  on the origin of gravity and the laws of newton
_ , _ jhep _ * 1104 * , 029 ( 2011 ) .
chien - i chiang , je - an gu , and pisin chen , _  constraining the detailed balance condition in horava gravity with cosmic accelerating expansion
_ , _ jcap _ * 1010 * , 015 ( 2010 ) .
h. lu , j. mei , and c. n. pope , _  solutions to horava gravity
_ , _ phys .
* 103 * , 091301 ( 2009 ) .
a. tawfik and a. diab , _
`` emergence of cosmic space and minimal length in quantum gravity : a large class of spacetimes , equations of state , and minimal length approaches '' _ , _ indian j. phys . _ , * 90 * , 1095 ( 2016 ) .    y .- w .
kim , s. k. kim , and y .- j .
park , _  thermodynamic stability of modified schwarzschild - ads black hole in rainbow gravity  _ , eur .
j. c * 76 * , 557 ( 2016 ) . s. gangopadhyay and a. dutta , _  constraints on rainbow gravity functions from black hole thermodynamics  _ , europhys .
115 * , 50005 ( 2016 ) .
r. garattini and e. n. saridakis , _  gravity s rainbow : a bridge towards hoava - lifshitz gravity
j. c * 75 * , 343 ( 2015 ) .
mir faizal and barun majumder , _
 incorporation of generalized uncertainty principle into lifshitz field theories  _ , annals phys .
* 357 * , 49 ( 2015 ) .
abdel nasser tawfik and eiman abou el dahab .
_  corrections to entropy and thermodynamics of charged black hole using generalized uncertainty principle  _ , int .
a * 30 * , 1550030 ( 2015 ) .
abdel nasser tawfik and abdel magied diab , _  black hole corrections due to minimal length and modified dispersion relation
a * 30 * , 1550059 ( 2015 ) .
s. robles - perez , y. hassouni , and p. f. gonzalez - diaz , _  coherent states in the quantum multiverse
b * 683 * , 1 ( 2010 ) .
s. gielen , d. oriti and l. sindoni , _  homogeneous cosmologies as group field theory condensates  _ , jhep * 1406 * , 013 ( 2014 ) .
mir faizal , _  multiverse in the third quantized horava - lifshits theory of gravity  _ , mod . phys . lett .
a * 27 * , 1250007 ( 2012 ) .


quantum information protocols have been demonstrated experimentally in both the discrete - variable ( dv ) and so - called continuous - variable ( cv ) settings .
dv quantum information protocols employ qubits  @xcite and qudits  @xcite , and cv quantum information protocols regard continuously parameterized canonical position states as the logical elements analogous to qubits for the dv case  @xcite .
cv quantum information is experimentally appealing because sophisticated squeezed light experiments have led to claims of successful quantum information protocols such as teleportation  @xcite , key distribution  @xcite , and memory  @xcite , but the theoretical status of cv quantum information is challenged by unresolved issues concerning quantum error correction  @xcite , non - distillability  @xcite , no - go theorems for quantum computation  @xcite , and the absence of full security proofs for key distribution .
cv information processing has also been studied for classical models , including the now named blum - shub - smale machine @xcite and continuous turing machines @xcite .
these models are of background relevance to the research into cv quantum information and are referenced here for contextual purposes .    in this paper , we establish a sound theoretical framework for studying quantum algorithms and apply this framework to study the cv analogue of the early dv quantum algorithm , known as the deutsch - jozsa ( dj ) algorithm  @xcite . the problem solved by dj algorithm is the following .
[ problem : balancedorconstant ] given a function @xmath0 that is promised to be either _ constant _ ( @xmath1 takes the same value everywhere ) or _ balanced _
( @xmath1 takes the value 0 on exactly half the inputs ) , determine whether @xmath1 is constant or balanced .
the best classical algorithm requires @xmath2 evaluations in the worst case .
if error is tolerated , for any integer @xmath3 , to achieve an error of at most @xmath4 , any probabilistic algorithm requires a number of evaluations that is at least of order  @xmath5  @xcite .
if the function is accessible on a quantum computer as a quantum oracle , then the dj algorithm is exact and requires just one evaluation to solve the problem .
our focus here is on the cv analogue of the dj algorithm , and we are inspired by the braunstein and pati formulation  @xcite of the cv dj algorithm ; however , our work differs from theirs in that ours relies only on logical states that are elements in the hilbert space and thus provides a strict cv version of the dj problem .
we introduce a particular model for the computation of the dj problem in a cv setting . within the constraints of this model
, our analysis shows that the cv dj algorithm is necessarily probabilistic and its performance must therefore be compared to the classical case where bounded error is tolerated and not to the classical deterministic case .
we choose the dj algorithm for the following reasons .
two types of quantum algorithms dominate the field , those that implement a version of the hidden subgroup problem and those that use a version of grover s search algorithm  @xcite .
an early example of the former is the deutsch - jozsa ( dj ) algorithm  @xcite , which is amongst the oracle class of problems  @xcite that have been important in demonstrations of quantum speed - ups . finally the cv dj algorithm has a head start in the work of braunstein and pati so our analysis can build on their concepts  @xcite .
our paper is presented as follows . in sec .
[ sec : background ] , we review the dj algorithm . although this algorithm is well known , our review serves as a foundation for careful construction of the cv version of this algorithm .
furthermore we compare the dj algorithm s performance against both deterministic and probabilistic strategies , especially because the cv case can only be properly compared against probabilistic strategies because the cv dj algorithm can never be deterministic .
our description of the dv dj algorithm comprises three steps so that these steps can be discussed separately during the construction of the cv analogue .
our approach emphasizes a recasting of the dv dj algorithm in that we do not need the target qubit .
this approach leads to an easier adaptation to the cv case . in sec .
[ sec : background ] , we review the formalism of rigged hilbert spaces ( rhs )  @xcite since our cv algorithm , as well as any other cv quantum algorithms , must work in a rhs .
this will have implications when we discuss the limitations of error inherent in our cv dj algorithm in sec .
[ sec : cvrep ] and in sec .
[ sec : bounding ] .    in sec .
[ sec : cvrep ] , we adapt the dj problem to the cv case and develop the cv dj algorithm through the same three fundamental steps of the algorithm .
we pay particular attention to the challenge of encoding a finite @xmath6-bit string into functions over the real numbers .
overcoming this challenge enables us to recognize that perfect encoding results in the inability to determine if the encoding is of a constant string or balanced in a single execution of the algorithm .
we show that this probabilistic nature of the algorithm is the result of an uncertainty relation between the continuous representation and its fourier - transform dual representation .    in sec .
[ sec : bounding ] , we determine an upper bound on the query complexity of the cv dj algorithm .
we note that because the cv dj algorithm is shown to be probabilistic , its performance can only logically be compared to the classical probabilistic algorithm and not to the classical deterministic algorithm .
we conclude that the formalism presented herein is applicable to a wide range of oracle identification problems in a cv setting .
we cast the dj problem into the class of ` oracle identification problems ' in subsec .
[ subsec : oracle ] .
we then review deterministic algorithms in subsec .
[ subsec : deterministic ] and probabilistic algorithms in subsec .
[ subsec : randomized ] . in subsec .
[ subsec : quantumdj ] , we analyze an alternative representation of the quantum dj algorithm that uses @xmath7 qubits instead of the traditional @xmath8 qubits . in subsec .
[ subsec : cv ] , we present a primer on the rigged hilbert space and close with a discussion of the concepts required to transition from discrete variables to continuous variables .      the dj problem is an identification problem in which we are given a function from some candidate set @xmath9 of functions .
the candidate set @xmath10 is the disjoint union of two collections of functions , and our task is to determine which of the two collections the function @xmath1 is drawn from .    [
problem : oip ] consider the set of all functions from @xmath7 bits to one bit , @xmath11 .
let @xmath12 and @xmath13 be disjoint subsets of @xmath14 .
given some oracle @xmath15 with the promise that either @xmath16 or @xmath17 , determine the index @xmath18 such that @xmath19 .
for @xmath20 , we impose lexicographic order on the @xmath6-bit strings of @xmath21 .
we can then specify any function @xmath22 by writing all its @xmath6 function values in a list @xmath23 of length  @xmath6 .
the @xmath24 bit @xmath25 in the list is 1 if @xmath1 takes the value  1 on the @xmath24 bit - string of @xmath21 .
there are @xmath26 functions from @xmath7 bits to one bit , and thus our candidate set has cardinality upper bounded by  @xmath27 . in the following , we often write @xmath22 to denote the function that corresponds to the @xmath6-bit string  @xmath28 .
we are interested in finding an efficient strategy to identify the property of whether @xmath1 belongs to set @xmath12 or to set @xmath13 without necessarily determining @xmath1 itself . in the dj case ,
the property we are interested in is whether @xmath1 is balanced or constant  @xcite .
the cost of the algorithm is the number of queries made to the oracle .    with the promise of balanced or constant functions ,
there are far fewer than @xmath26 functions .
the number of balanced and the number of constant functions is readily ascertained from the binomial theorem applied to power sets .
the strings @xmath28 of length @xmath6 that correspond to the constant functions are the string consisting only of  1s and the string consisting only of  0s .
there are thus just two constant functions .
the strings @xmath28 of length @xmath6 that correspond to balanced functions are the strings in which exactly half of the bits are  0 and half are  1 .
there are thus precisely @xmath29 balanced functions .      on a classical turing machine , problem  [ problem : balancedorconstant ]
can be solved deterministically .
a deterministic algorithm corresponds to submitting queries in the form of @xmath7-bit inputs and obtain the one - bit output for each query .
there are  @xmath6 unique input strings , but the promise of balanced versus constant functions implies that only @xmath30 are required to determine whether the given function is balanced or constant , with certainty .
the reason that fewer than @xmath30 queries is insufficient is that only @xmath31 queries may reveal all output bits being the same , suggesting a constant function , whereas the remaining @xmath31 outputs could all be the opposite of the first @xmath31 queries .      in subsec .
[ subsec : deterministic ] , we saw that fewer than @xmath30 queries is insufficient for a deterministic algorithm , but that case seems highly unlikely .
more formally , fewer than @xmath30 queries will identify most of the balanced functions as non - constant in much fewer than @xmath30 queries .
here we ask the question about how many queries are required if we are prepared to tolerate a small number of errors .
in fact a probabilistic algorithm achieves an exponentially small error of @xmath4 with a number of queries that is only linear in  @xmath5  @xcite . to understand how a probabilistic algorithm can help ,
consider that , although a single query with a random input provides no information , two queries with two random inputs can be highly informative .
if the output from the second query differs from the first output , then the function is proved not to be constant and therefore must be balanced .
if , on the other hand , the second output is the same as the first , then the outcome is not certain , but the more times the outputs are the same , the more confident one can be about the function being constant .
-qubit `` control '' state , and the lower line represents the 1-qubit `` target '' state.,width=321 ]    we calculate the probability of successfully determining whether the given function  @xmath22 is balanced or constant . a lower bound on the success probability  @xmath32 for  @xmath5 queries can be achieved by examining a sampling - without - replacement strategy , which is expressed as @xmath33 here the equality is calculated assuming sampling without replacement and shows dependency on  @xmath6 , whereas the inequality in eq .
( [ eq : classdjran ] ) is based on sampling with replacement and is independent of @xmath6 .
the failure probability @xmath34 declines exponentially in  @xmath5 , the number of queries .    in subsec .
[ subsec : quantumdj ] , we study the quantum dj algorithm next where we show that the problem can be solved with a single query independent of @xmath6 .
although this exponential speed - up is impressive when compared to the classical deterministic approach , it is much less so when compared to the classical probabilistic approach .
the quantum dj algorithm has been shown to solve problem  [ problem : balancedorconstant ] in a single query  @xcite .
the quantum dj algorithm is usually studied via its corresponding quantum circuit .
we present a standard circuit version  @xcite in fig .
[ fig : tradqdjcircuit ] .    the state represented by the lower line in fig .
[ fig : tradqdjcircuit ] is referred to as the _ target _ qubit . in order for easier adaption of this circuit to the cv setting , we choose an alternative , and equivalent , circuit formulation  one without the target state .
we take this approach to avoid some of the difficulties the target state introduces in @xcite . the unitary operator associated with the oracle function
changes slightly in this alternative circuit .
we discuss these differences before proceeding with analysis of the circuit .
this simpler algorithm without the target qubit is given in fig .
[ fig : targetlessqdjcircuit ] .     defined in eq.[eq : qdjtargetlessuf],width=8 ]    .
[ fig : targetlessqdjcircuit ]    oracle application is the critical part of the algorithm .
the oracle construct originally proposed by dj is expressed , for @xmath35 and @xmath36 , as @xmath37 this construction yields a matrix representation for the  @xmath38 as a permutation matrix , hence always unitary  @xcite .
with respect to the ordered basis @xmath39 the unitary matrix @xmath38 can be expressed in the following insightful form @xmath40 with @xmath41 the @xmath42 @xmath43 operator in this case . here
@xmath38 is a @xmath44 matrix , which results from there being @xmath45 strings ( the arguments of @xmath1 ) and an additional target qubit .
the operator @xmath38 can also be expressed in the alternative ordered basis @xmath46 as @xmath47 with @xmath48 the @xmath49 identity operator .
furthermore the operator @xmath50 is expressed as the @xmath51 matrix @xmath52 and thus provides a reduced representation for  @xmath38 .
it is apparent that the operator  @xmath50 acts on a @xmath49 subspace of  @xmath38 since @xmath53 we make the assumption that if we have the oracle @xmath38 , we we also have the oracle @xmath50 .
we thus conclude that the construction employing both control and target qubits is not strictly necessary .
that is , one could construct this algorithm employing the @xmath7-qubit control state only .
apparently the choice of representation simply depends on the nature of the actual physical implementation .
we now present a step - by - step analysis of the alternative circuit presented in fig .
[ fig : targetlessqdjcircuit ] .
we shall analyze the cv circuit in the same steps for cross reference and comparison .
we use the hat notation @xmath54 in order to emphasize that this analysis is of the algorithm presented in fig .
[ fig : targetlessqdjcircuit ] , which employs @xmath7-qubit states and not of that presented in fig .
[ fig : tradqdjcircuit ] , which employs @xmath55-qubit states . the @xmath7-qubit input state of the circuit in fig .
[ fig : targetlessqdjcircuit ] is a string of qubits prepared in @xmath56 .
the next step in state preparation is to place the state @xmath57 into an equal superposition of all computational basis states @xmath58 for  @xmath59 the single qubit hadamard operator .      given the definition of the reduced operator @xmath50 defined in eq .
[ eq : qdjtargetlessuf ] , its effect on the equal superposition of basis states expressed in the state @xmath60 is to effectively encode the @xmath6-bit string @xmath28 unitarily into the state @xmath61 .
we express this as @xmath62 which is a convenient representation .
we shall show that this representation naturally extends to the cv setting .
measurement proceeds by first undoing the superposition created during the state preparation step .
this is achieved through the application of the operator @xmath63 , which modifies the state after oracle application @xmath64 the resultant state is @xmath65 we rewrite eq .
( [ eq : qdjstate3 ] ) with the operator @xmath66 expressed in terms of a recursive definition as follows @xmath67 given @xmath68 the combination of eq .
( [ eq : qdjtargetlessstate3recursive ] ) and eq .
( [ eq : hadamardsinglequbit ] ) allows us to see that all of the rows ( and columns ) of the operator @xmath66 have an equal number of positive and negative ones except for the first row , which consists entirely of plus ones .
it is this feature that permits the constant and balanced functions to be distinguished in a single measurement .
for the two constant cases , eq .
( [ eq : qdjtargetlessstate3recursive ] ) may be expressed as @xmath69 as only the first row does not result in amplitude cancellation of the @xmath45 constant amplitude components of the state @xmath61 .
each of the balanced functions result in the amplitudes of the state @xmath61 having an equal number of positive and negative ones .
this feature coupled with action of the operator @xmath66 results in the first component of the state @xmath70 having zero amplitude for all the balanced functions .
we express this result as @xmath71 where we use the symbol  @xmath72 to represent that the non - zero value(s ) will land on the other @xmath73 components depending on which of the @xmath29 balanced functions the oracle is set to .
it is interesting to note that the number of rows in the state @xmath74 potentially having a non - zero value is @xmath73 whereas the number of balanced functions is exponential in @xmath6 .
this means that many of the balanced states can be expressed as real - valued mixtures of the computational basis states with the condition that the amplitude of the first component is always zero .    for the final measurement step ,
we employ the projection operator  @xcite defined for @xmath75 as follows @xmath76 we are only concerned with the first component as discussed above , so for the constant cases we have @xmath77=\left\langle\hat{\psi}_{3\text{c}}\right|m_{(0\cdots 0 ) }          \left|\hat{\psi}_{3\text{c}}\right\rangle=1,\ ] ] and for all balanced cases we have @xmath78=\left\langle\hat{\psi}_{3\text{b}}\right|          m_{(0\cdots 0)}\left|\hat{\psi}_{3\text{b}}\right\rangle=0\ ] ] as required .
we have completed the study of the quantum dj algorithm in a form that allows us to adapt readily to the cv setting .
our strategy will be to construct a cv algorithm analogous to that shown in fig .
[ fig : targetlessqdjcircuit ] and whose operator representation is given by @xmath79 this approach is simpler , and we can worry about whether or not an implementation will require target states when a particular implementation is considered . before delving into the cv algorithm , we present some background cv information .      the transition from dv to cv quantum information requires an extension of hilbert spaces to rigged hilbert spaces  @xcite , which allows the use of position states  @xmath80 with @xmath81 but restricts dual states to so - called ` test functions ' .
an inner product between position states and test functions is meaningful , but the inner product between two position states leads to the dirac relation @xmath82 , which must be treated carefully . for  @xmath7 the size of problem  [ problem : balancedorconstant ] , the target - less quantum dj algorithm requires @xmath7 qubits , which requires a hilbert space of size @xmath20  @xcite .
the hilbert space for cv problems seems quite generous in this respect as it is infinite - dimensional .
in fact the cv hilbert space is congruent to the space of square - integrable complex functions over the real field @xmath83  @xcite .
a function @xmath84\rightarrow \mathbb{c}$ ] , for @xmath85\subset\mathbb{r}$ ] , is in @xmath83 if @xmath86 the inner product of two functions  @xmath87 is @xmath88 with positive definite norm and distance metric defined by @xmath89 respectively .    typically in cv quantum information discourse ,
the position states  @xmath80 are introduced as a basis set of the hilbert space with each  @xmath80 an eigenstate of a position operator  @xmath90 , with  @xmath81 .
unfortunately the state  @xmath80 does not exist in the hilbert space ; this problem is evident in the standard inner product @xmath91 as @xmath92 is not a proper function , position states are not proper states .
fortunately the position states are correct as a representation ; for example @xmath93 is the position representation of test function  @xmath1 within the context of the rigged hilbert space . also eq .
( [ eq : cv_deltafunction ] ) is meaningful in the context of distribution theory .
a rigged hilbert space is a pair @xmath94 such that @xmath95 is a hilbert space and @xmath96 is a vector space that is included by a continuous mapping into @xmath95 : @xmath97 .
elements of  @xmath96 are referred to as ` test functions ' , and the dual to @xmath96 is @xmath98 , for  @xmath99 dual to @xmath100 and  @xmath101 comprising generalized functions , or ` distributions ' .
the inner product @xmath102 is in @xmath103 $ ] for any @xmath104 and for any @xmath105 @xcite .
note that the adaptation of the dv dj algorithm to the cv regime needs to be done in the context of a computational problem .
here the relevant problem is still problem  [ problem : balancedorconstant ] , and the notion of the oracle remains unchanged .
thus , in the cv case , our task is still to determine whether the function  @xmath22 belongs to the set of constant functions or to the set of balanced functions .
we begin by giving a strategic overview in order to convey the key concepts of our approach to developing a cv computation model .
we follow this with a subsection giving some preliminary definitions allowing us to set the stage for detailed analysis .
we then proceed with a step - by - step analysis of our cv dj algorithm .
although we are now working with cv , instead of dv , quantum information , the computational problem to be solved remains problem  [ problem : balancedorconstant ] .
in other words , we want to learn whether the function  @xmath22 is constant or balanced with as few oracle queries as possible .
another way to think of this is that we wish to determine the index @xmath106 such that @xmath107 .
we now give a conceptual overview of our model for cv quantum computation of the dj problem , which we follow later with a rigourous treatment .    in our model of cv quantum computation , we will use the continuous position and momentum variables of a particle . for @xmath108
, we use the particle s position wave function ,  @xmath109 , to describe where the particle is concentrated and the particle s momentum wave function ,  @xmath110 , to describe its momentum distribution .
the position and momentum wave functions are fourier transform pairs , and the relationship between the particle s position and its momentum is governed by heisenberg s uncertainty principle .
-bit string in a region of momentum extending from  @xmath111 to  @xmath112 using the  @xmath113 ,  @xmath114 example .
note that each of the bits @xmath115 are uniquely represented.,width=8 ]    there are many position and momentum wave function pairs on which we could base our computational model .
we select our particular pair as follows .
first , we wish to encode the unknown @xmath6-bit string , @xmath28 , in the momentum domain .
we do so because encoding in the momentum domain is the continuous analogue of the discrete case , where encoding is performed on an equal superposition of computational basis states .
second in order to fix one of the degrees of freedom of the problem , we want each of the bits comprising the string @xmath28 to be unambiguously represented in the momentum space . by unambiguous we mean that each of the bits are represented by equal - sized , non - overlapping , contiguous regions in the momentum space .
since we want each of the @xmath6-bits comprising the string to be represented unambiguously , we naturally think of each bit as being manifested by a finite - width square pulse whose position in momentum space represents the bit position in the string @xmath28 and whose magnitude represents the bit value .
continuing along this line of reasoning to the representation of the entire string @xmath28 , we can imagine we have a region of momentum extending from @xmath111 to @xmath112 .
all the contiguous momentum pulses within this region thus have  width " @xmath116 , and for @xmath117 , the @xmath118 momentum pulse is centred at position @xmath119 and takes on value @xmath120 .
we illustrate this concept in fig .
[ fig : cvdjalgencodeconcept ] for a particular  @xmath113 case .
is that of a sinc function .
( b ) the fourier transform of the position state wave function is a  pulse " function in the momentum domain , which acts as the encoding  substrate " .
( c ) the @xmath6-bit string  @xmath121 modulates this momentum  substrate " .
( d ) the inverse fourier transform of the encoded  square wave " produces a  generalized " sinc function whose infinite position extent necessitates an optimal measurement  window " parameterized by  @xmath122.,width=8 ]    the picture that thus emerges is that each of the @xmath26 possible strings may be represented by a uniquely shaped  square wave " having extent @xmath123 comprising @xmath6 pulses each of width @xmath124 and having magnitude @xmath125 . with the encoding concept clear
, we conceptually illustrate the four key stages of the algorithm in fig .
[ fig : cvdjalgoverview ] .
we begin with a position wave function centered at  @xmath126 and illustrated in fig .
[ fig : cvdjalgoverview ]  ( a ) .
note that this position wave function is a sinc function since sinc / pulse functions are fourier transform pairs . in fig .
[ fig : cvdjalgoverview ]  ( b ) , we present the momentum wave function , a pulse function , which acts as the  substrate " into which the @xmath6-bit strings are encoded . in fig .  [
fig : cvdjalgoverview ]  ( c ) , the pulse function is encoded with the particular @xmath6-bit string  @xmath121 . finally ,
the inverse fourier transform of this  square wave " is presented in fig .
[ fig : cvdjalgoverview ]  ( d ) .
since the inverse fourier transforms of finite pulses in the momentum domain have infinite extent in the position domain , we need to limit the extent of our measurement to @xmath122 .    in summary
, we see that our algorithm will need the parameters @xmath6 , @xmath127 and @xmath92 .
we note that as @xmath6 gets large , the individual pulse width associated with a single bit becomes small appearing to pose a limit on the maximum value of @xmath6 .
we will return to this issue once we have determined the relationship between @xmath127 and @xmath92 .
there are many potential models for quantum computation in a cv setting .
we have chosen to study one where we unambiguously encode an n - bit string into the continuous momentum variable of a particle . within the constraints of this model
, we will show that the cv dj is necessarily probabilistic and prove an upper bound on the query complexity of the cv dj problem .
we speculate that we ca nt do better than this .
for example if the momentum / position pair are described by gaussian / gaussian functions , as would be the case for the physically meaningful states of quantum optics , imperfect encoding of the @xmath6-bit string in the momentum domain will result in increased position error .
whether or not this will in turn impact the
big oh " representation of the query complexity requires further research as does a general proof of a lower bound .
the challenge will be to show that another strategy can do better than the model described herein .
we now proceed to formalize some of the concepts presented in the previous subsection . here
we describe a ` natural ' way of encoding a finite - dimension , @xmath6-bit string in a continuous domain .
we define the following function , along with its fourier dual , to help us achieve this end .    for @xmath128 , the ` top hat ' function @xmath129\\0,&\text{if}\quad p\notin[p_0-p , p_0+p]\end{array }              \right.\end{aligned}\ ] ]
will be especially useful in bridging the gap between dv and  cv quantum information because @xmath130 so the state @xmath131 is , in some sense , a momentum eigenstate @xmath132 in the limit @xmath133 .
the inverse fourier transform of the function @xmath134 is @xmath135 where @xmath136 defines the position of the sinc function .
the limit of @xmath109 as @xmath127 goes to @xmath137 yields @xmath138 .
the position eigenstate @xmath139 is likewise formed in the limit @xmath140 .
now imagine we want to sum a contiguous string of  pulses " described by the top hat function ( [ eq : tophat ] ) with all pulses having width @xmath141 and the @xmath142 pulse having complex amplitude @xmath143 .
this results in the composite function @xmath144 this function can also be used as a basis @xcite of cv kets in dirac notation as @xmath145 thus allowing us to encode quantum information in the cv domain .
note that @xmath146 is the complex amplitude for real - valued @xmath147 .
this affords a consistent way of encoding a discrete wave function over a continuous domain .    before proceeding with a formal analysis of the algorithm
, we give an overview of our proof strategy .
the oracle is either set to one of two constant strings or to one of @xmath148 balanced strings . a string and its complement
have indistinguishable probability distributions , so there are a total of one constant probability distribution plus @xmath149 balanced probability distributions representing the possible oracle settings . in order to simplify the analysis
, we wish to replace this exponential number of balanced probability distributions with a single  worst - case " balanced probability distribution .
thus we seek a particular balanced string ( and its complement ) whose probability distribution is most likely to  fool " us into concluding it is a constant string .    intuitively , the balanced strings that have the fewest number of changes between adjacent bits in the interval @xmath150 $ ] will be the most  constant like " of the balanced strings .
there are no balanced strings with zero changes - this is the key feature that separates the constant strings from the balanced strings .
there is however , a single pair of balanced strings having only one change .
these strings exhibit the feature that the first @xmath31 bits are constant and the second @xmath31 bits are the complement of the first .
we call these strings the anti - symmetric balanced ( asb ) strings
. one of these two strings is illustrated in fig .
[ fig : cvdjalgoverview ]  ( c ) .
note that all other balanced strings have more than one change .
our proof strategy begins by making the assumption that the asb case is the  worst case " of all balanced cases .
we use this assumption to determine the optimum value of @xmath92 , which is the extent of our measurement in the position domain and is illustrated conceptually in fig .
[ fig : cvdjalgoverview ]  ( d ) . given this optimum value of @xmath92 ,
we then prove by induction that the worse balanced case is indeed the asb case .
our strategy is to create a cv analogue of the alternative formulation of the discrete dj algorithm presented in fig .
[ fig : targetlessqdjcircuit ] .
the cv extension of this is presented in fig .
[ fig : targetlesscvdjcircuit ] .
our construction of a cv dj algorithm employs some of braunstein and pati s techniques  @xcite and avoids the pitfalls .
in particular , we employ position states as a logical representation ( states in  @xmath101 ) analogous to the discrete computational basis states .
encoding is not , however , into the position states but rather into test functions  @xmath151 with @xmath152 .
furthermore we employ a fourier transform to operate as a cv version of the dv hadamard transform ( extending the hadamard transformation to the cv case is not unique  @xcite ) .    for  @xmath72 the canonical position and @xmath147 the canonical momentum , the fourier transform maps a function @xmath109 to its dual @xmath110 according to  @xcite @xmath153 such that @xmath154 and @xmath155 note that we make use of the momentum variable @xmath147 as the fourier dual of the position variable @xmath72 .
the function @xmath156 can be a test function in @xmath96 and @xmath109 is the inner product of @xmath156 with the position in @xmath101 : @xmath157 .
the momentum state @xmath158 is the fourier transform of  @xmath80 and @xmath159 .    with these concepts in order
, we now proceed through the cv dj algorithm analogous to the three steps in the dv dj algorithm .
the function notation  @xmath109 and @xmath110 is more convenient here rather than the dirac notation in the previous section .
we have argued previously that we need the fourier transform of the input state to be the top hat function defined in eq .
( [ eq : tophat ] ) .
we add several conditions that do not take away from the generality of the solution .
first , we want the top hat to have zero phase , which gives @xmath160 and to be centred at @xmath161 .
second , we want the pulse to have extent @xmath123 .
this gives us the simplest form of the sinc function for the initial state @xmath162 we note that the limit of @xmath163 as @xmath164 gives a @xmath165 .
thus we can think of the quantity @xmath127 as playing the role of the standard deviation in a gaussian distribution .
the final step in state preparation is to perform the fourier transform , which yields the top hat function with extent @xmath123 @xmath166\\0,&\text{if}\,\,p\notin[-p , p].\end{array }              \right.\ ] ] this function forms the raw substrate , which will be ` modulated ' by the individual @xmath6-bit strings @xmath28 .
we perform encoding by partitioning the real numbers representing momentum into non - overlapping , contiguous and equal - sized bins . in this digital - to - analogue strategy ,
the width of each @xmath147-bin is  @xmath124 , and @xmath167\\0,&\text{otherwise . }
\end{array }          \right.\ ] ] the oracle encodes the index  @xmath28 into the function @xmath22 as follows : @xmath168 where the factor  @xmath169 serves to modulate the phase of the top hat function according to the bit value .
consider the case @xmath170 ; hence @xmath171 .
as one case , the function corresponding to the four - bit string @xmath172 is @xmath173 the only two four - bit strings yielding constant functions would be @xmath174 , for which the function is identically unity over the whole domain  @xmath175 $ ] , and @xmath176 , for which the function is identically @xmath177 over  @xmath175 $ ] .
four cases are presented in fig .
[ fig : cvdjsqwaveneq2 ] .
we refer to the function @xmath178 as the  lowest - order " antisymmetric balanced wave as it has just one zero crossing in @xmath175 $ ] .
( a )  @xmath179 , ( b )  @xmath180 , ( c )  @xmath114 , and ( d )  @xmath181.,width=8 ]    in the limit that @xmath182 with  @xmath127 fixed , @xmath183 for @xmath184 the midpoint of the @xmath185 bin .
the limit @xmath182 thus gives a prescription for approaching a continuous variable representation where the @xmath28 index seems to approach a continuum ; however this limit yields a countable , rather than uncountable , set  @xmath186 , and the finite domain @xmath175 $ ] has important ramifications on the nature of the functions corresponding to fourier transforms of  @xmath187 .
we express the state after encoding as @xmath188 where we observe the `` modulating '' effect of the encoded string @xmath189 on the momentum `` substrate '' @xmath110 .    in the context of the digital - to - analogue strategy , the constant functions are analogous to direct current ( dc ) signals and the balanced functions to alternating current ( ac ) signals .
the number of zero - crossings corresponds to frequency information , and the question of whether the output is balanced or constant is essentially a problem of querying whether there is a non - zero frequency component of the output signal . as noted previously ,
the asb function has the lowest frequency component .
we now proceed to analyze the measurement stage .
we have the strings @xmath152 encoded into the momentum state  ( [ eq : encodedmomentum ] ) .
the next step prior to the final measurement is to take the inverse fourier transform of this pulse train . for @xmath115 the @xmath142 bit of @xmath28 , this is expressed as @xmath190 the expression given in eq .
( [ eq : cvdjphaser1 ] ) can be simplified to yield @xmath191where we have defined @xmath192 we see that the magnitude of an individual generalized sinc function , @xmath193 , is determined by a vector sum of @xmath6 phasors , which is modulated by a particular @xmath6-bit string @xmath28 .
note that the phasors , @xmath194 , are equiangular divisions of the angular interval @xmath195\nonumber,\ ] ] and they exhibit the pairwise complex conjugate property @xmath196 .     defined in eq .
( [ eq : cvdjphaser2 ] ) for @xmath197 ( a ) @xmath198 phasors range between @xmath199 in steps of @xmath200 , and(b ) @xmath201 phasors range between @xmath202 in steps of @xmath203.,width=7 ]    in fig .
[ fig : cvdjsincphasers ] , we present the phasors for @xmath197 with @xmath198 and @xmath201 to illustrate these features . note that the phasors are added constructively or de - constructively depending on the phase of the angles , which results from the term @xmath120 .
this effect defines the magnitude of the resulting sinc function .
we note that the only functions with @xmath204 are the two constant sinc functions .
this is clear given that @xmath205 for the two constant cases , and @xmath206 since for all balanced cases , the latter sum always resolves to @xmath207 .
this feature of the set of @xmath208 sinc functions represented by @xmath209 defined by eq .
( [ eq : cvdjphaser2 ] ) implies the strategy for measurement that will distinguish between the constant and balanced cases .    in order to refine this strategy , we focus on two cases .
the first of these cases is for the two constant functions for which eq .
( [ eq : cvdjphaser2 ] ) gives the probability distribution @xmath210 where we have @xmath211 the second case deals with the two balanced functions having the lowest ` frequency ' content , which occurs when the first @xmath31 bits and the last @xmath31 bits have opposite values .
we think of these two balanced functions as having the lowest frequency content since eq .
( [ eq : fz - encoding ] ) has a single zero crossing in the interval @xmath175 $ ] for these two balanced strings only .
all other balanced strings have more than one zero crossing and thus higher frequency content . for this pair of balanced functions , which we call the antisymmetric balanced ( asb ) functions , eq .
( [ eq : cvdjphaser2 ] ) gives the probability distribution @xmath212 where we have @xmath213 note that of the @xmath214 balanced functions , there are many that are also antisymmetric about the midpoint .
however , we reserve the term asb for these two _ lowest - order _ antisymmetric balanced functions .
we will use these cases to bound the success probability of distinguishing between the constant and all balanced cases .
we first illustrate the concept of frequency in the following example .
again consider the case @xmath215 , @xmath170 ; hence @xmath171 .
as one case , the function corresponding to the four - bit string @xmath172 is @xmath216 this function corresponds to the @xmath113 asb function .
the probability distributions for the four distinct @xmath113 cases are presented in fig .
[ fig : cvdjsincftnsneq2 ] .
we clearly see that of the three balanced cases , the @xmath113 asb function has probability peaks closest to @xmath160 .     for ( a )  @xmath179 , ( b )  @xmath180 , ( c )  @xmath114 , and ( d )  @xmath181 .
we clearly see that of the three balanced cases ( b ) through ( d ) , the asb function ( b ) has probability peaks closest to @xmath160.,width=8 ]    our measurement strategy is to measure the probability distribution in a small band around the position @xmath160 parameterized by @xmath122 .
the cv analog of the projection operator given in eq.([eq : qdjmeasure])is defined as @xmath217 where @xmath218[eq : cvdjdoorway]due to the symmetry of the sinc functions about @xmath136 , we set @xmath219 and @xmath220 .
we now need to determine the optimal value of @xmath92 that will maximize our ability to distinguish between the constant and balanced cases .
we will determine the optimum value of @xmath92 by first assuming that the probability distribution @xmath221 given by eq .
( [ eq : cvdjstate3pdfabal ] ) dominates all other balanced probability distributions in the region @xmath222 $ ] . after using this assumption to determine a value for the optimal delta , we will state and prove a theorem justifying our assumption . as an illustration that our assumption is true for the @xmath113 case , we plot the four distinct cases in fig .  [ fig : cvdjoptimalmeasurelimits ] .
the ability to effectively distinguish between two random events is proportional to the separation of the individual probabilities of occurrence .
thus we need to select @xmath92 such that we get as much separation between the constant distribution and the asb distribution as possible . given this concept we can think that when we make a measurement we are distinguishing between two events , the probabilities for which we define as follows @xmath223=   e^{\delta}_{-\delta}\left(\mathcal{p}_{\text{\tiny{c}}}(x)\right )   \label{eq : cvdjprobconst},\ ] ] and @xmath224=   e^{\delta}_{-\delta}\left(\mathcal{p}_{\text{\tiny{asb}}}(x)\right )   \label{eq : cvdjprobasb}.\ ] ]    , the optimal value of @xmath225 .
this graph shows that only the constant and the antisymmetric balanced functions significantly contribute to probability between @xmath226 .
, width=6 ]    we can determine the optimum value of @xmath92 by maximizing the expression @xmath227 .
it suffices to find the value of @xmath92 for which @xmath228 , which may be expressed as @xmath229this occurs where @xmath230 for @xmath231 , which gives a global maximum at @xmath232 .
it is interesting to think of this result as an uncertainty relationship @xmath233we shall return to this concept in our discussion in the conclusion .
we have determined the optimum value for @xmath92 based on our assumption that for @xmath234 the balanced probability distribution @xmath235 dominates all other balanced probability distributions .
we now proceed to prove this assumption .    in order to proceed with the proof
, we define a set @xmath96 of @xmath5 pairwise conjugate angles with @xmath236 .
note that @xmath6 is not restricted to being equal to @xmath45 for the purpose of this proof .
also for the purpose of this proof , we set @xmath215 and incorporate @xmath72 into the definition of @xmath237 for @xmath238 .
we let @xmath239 with @xmath240 and note the pairwise conjugate property @xmath241 . now consider @xmath242 where @xmath243\mapsto\pm 1 $ ] subject to the balanced condition @xmath244 , then    max @xmath245 occurs under the specific balanced conditions @xmath246 and @xmath247 which we refer to as the asymmetric balanced functions ( asb ) .
proof is done by induction on @xmath5 .
we begin with the base case @xmath248 , @xmath249 .
this case is trivial since the only balanced cases are the two asb cases represented by the strings @xmath250 .
we proceed with the base case for @xmath251 , @xmath113 .
this case is a little more involved .
we begin by labelling the angles and phasors as shown in fig .
[ fig : cvdjoangledefintionsforproof ] .     base base .
note that the effect of varying @xmath72 over @xmath252 $ ] simply focuses or expands the double angles @xmath253 and @xmath254 proportionally .
, width=6 ]    there are @xmath255 balanced cases , and we have to consider the strings @xmath256 .
since the latter three are complements of the first three , we have to consider only three vector sums .    with reference to fig .
[ fig : cvdjoangledefintionsforproof ] , we have @xmath257 .
we simplify and express the resultant along with the three other cases as    1 .
@xmath258 2 .
@xmath259 3 .
@xmath260 .
clearly @xmath261 .
we use the trigonometric identities,@xmath262 and @xmath263 .
we note that @xmath264 and @xmath265 for all @xmath6 and the specified range of @xmath72 .
since @xmath266 for @xmath267 , we conclude that @xmath268 and thus @xmath269 for @xmath270 .
this proves that the theorem is true for the @xmath251 , @xmath113 base case .
we are now ready to prove the inductive step .
we consider two cases .
case ( i ) assumes every pair is balanced . by this
we mean that @xmath271 . by inspection
, this gives the same result as @xmath272 and @xmath273 for the @xmath251 , @xmath113 case .
case(ii ) assumes that case(i ) is not true and is proved by induction . since case(i ) is not true , there must exist two non - balanced pairs for which @xmath274 and @xmath275 . as an illustration in the @xmath276 , @xmath197 case
, the balanced string @xmath277 has this property .
the inductive step is @xmath278 where @xmath279 is maximized for the @xmath251 , @xmath113 base case .
only when @xmath280 itself is maximized is equality achieved and the total sum maximized .
this occurs for the asb strings .
we have established that we can bound the probabilities of determining whether an unknown function is balanced or constant in a single query in a cv setting . in the next section
will determine an upper bound for the query complexity of a cv algorithm in terms of success probability in terms of the number of queries .
before we bound the query complexity , we make some important observations regarding the comparison between the discrete dj algorithm and the cv dj algorithm .
first , we note that probability distributions @xmath281 and @xmath235 defined by eqs .
( [ eq : cvdjstate3pdfconst ] ) and  ( [ eq : cvdjstate3pdfabal ] ) respectively , are in @xmath282 , the hilbert space of @xmath83 functions over the interval @xmath283 $ ] .
this implies that since we are measuring over a finite interval , the cv dj algorithm is necessarily probabilistic .
furthermore , we noted that @xmath127 and @xmath92 are related by the uncertainty relation given in eq .
( [ eq : cvdjuncertainty ] ) .
this leads to the conclusion that even in the limit of the improper delta function @xmath138 , the cv dj algorithm remains probabilistic .
this conclusion is contrary to that made in  @xcite .
second , we compare the operator descriptions of the dv dj and the cv dj , which we express as @xmath284 the first equation represents the quantum dj algorithm operator expression given in eq .
( [ eq : qdjtargetlessoperator ] ) .
the second equation is the analogous cv operator expression determined by concatenating the steps of the previous section .
there is a high degree of similarity between these two expressions , but there are mathematical subtleties .    1 .
the cv position state @xmath163 is not a perfect analog to the computational basis state @xmath56 except in the limit .
however , this limit creates a state that is not in the rhs we argued is necessary for consistency  @xcite .
the continuous fourier transform is not equal to the cv extension of the hadamard operator in a cv parameterized system with a finite hilbert space .
it is however , a convenient extension when the hilbert space is infinite .
3 .   finally ,
the diagonal operator @xmath50 given by eq .
( [ eq : qdjtargetlessuf ] ) has each entry taking on the value @xmath285 dependent on the value of @xmath189 .
the cv analogue to this operator is the function @xmath286 , where each of the @xmath6 partitions of the real interval @xmath175 $ ] similarly take on the value @xmath125 .
we now determine numerical values of the probabilities determined in eqs .
( [ eq : cvdjprobconst ] ) and  ( [ eq : cvdjprobasb ] ) .
we can readily calculate the probability of detecting if the function is constant @xmath287 where the sine integral is given by @xmath288 note this probability depends only on the product @xmath289 . if the function is the lowest - order antisymmetric balanced ( asb ) @xmath290[eq : cvdjmeasureabalprb]for @xmath291 , the numerical values of these two probabilities are @xmath292 and @xmath293 given this probabilistic nature of the cv dj algorithm , we need to develop a strategy to bound the error probability
. we will employ the technique sometimes called _ probability amplification _
@xcite .
our strategy will be to make @xmath5 repetitions of the cv dj algorithm where we assume that the oracle is set to the same function for each of the repetitions .
each repetition ends with a measurement . from this sequence of measurements we want to determine whether the unknown function is balanced or constant with high probability .
an error of @xmath294 can be achieved by making @xmath295 repetitions of the cv dj algorithm .
we will adopt the convention that when we make a query to the cv dj algorithm we either detect something ( algorithm returns a 1 ) , or we do not ( algorithm returns a 0 ) .
we can thus treat multiple queries as a sequence of bernoulli trials @xcite .
we assume that we have set our measurement limits to the optimal @xmath226 .
the two events we are trying to uncover are the constant cases where , for ease of calculation we set the probability of detecting something is @xmath296 , and the balanced cases where the probability detecting something is @xmath297 . note that we have set the probabilities to these rational numbers for illustrative purposes and to simplify the calculation .
we can make this arbitrary setting , and we will get the same result as long as the probabilities are bound from @xmath298 by a constant .    if each measurement is based on an independent preparation of the state @xmath163 , then each of the queries are independent . after a series of @xmath5 queries
, we can use the chernoff bounds of the binomial distribution to amplify the success probability  @xcite .
the simplest ( but somewhat weak ) chernoff bound on the lower tail is given by  @xcite as @xmath299<\text{e}^{-\frac{\mu \epsilon^2}{2}},\label{eq : cvdjchernoffweakl}\end{aligned}\ ] ] and on the upper tail as @xmath300<\text{e}^{-\frac{\mu \epsilon^2}{4}},\label{eq : cvdjchernoffweaku}\end{aligned}\ ] ] where @xmath301 is the expected mean of the resulting binomial distributions after @xmath5 queries , and @xmath302 is the relative distance from the respective means .
first , we bound the lower tail corresponding to the distribution of the constant case for which we have @xmath303 . here
we set @xmath304 , which expresses the probability for the value being less than half way between the two means as @xmath305<\text{e}^{-\frac{m}{24}}$ ] .
similarly , we bound the upper tail for the balanced case for which we have @xmath306 .
here we set @xmath307 , which expresses the probability for the value being greater than half way between the two means as @xmath308<\text{e}^{-\frac{m}{16}}$ ] . clearly the success is worse for the lower tail allowing us to bound the success probability of the cv dj algorithm after @xmath5 queries as @xmath309\geq1-\text{e}^{-\frac{m}{24}}.\label{eq : cvdjsuccessprob}\end{aligned}\ ] ] this gives an error probability that is @xmath310 as required .
we note that this is of the same order as the exponentially good success probability we have for the classical probabilistic approach given by eq .
( [ eq : classdjran ] ) . also note that this query complexity is independent of the value of @xmath6 .
we have made no attempt to obtain a tighter bound preferring to show only that the success probability of the cv dj algorithm is of the same order of that of the classical probabilistic approach to solving the dj problem .
in this paper we have presented a rigorous framework for the analysis of the dj oracle identification problem in a cv setting . the rigged hilbert space ( rhs )
affords a consistent transition from the traditional discrete hilbert space to the cv setting .
our framework allows us to define a consistent way of encoding @xmath6-bit strings into functions over the real numbers .
we have used this framework , and the selection of the sinc / pulse fourier transform pair , to prove that a cv implementation of the dj algorithm can not provide the exponential speed - up of its discrete quantum counterpart .
additionally , we have presented a bounded - error , upper bound on the query complexity of the dj problem within the constraints of our model .
the lack of speed - up results from an uncertainty principle between the ability to encode perfectly in a continuous representation and the subsequent inability to measure perfectly in the fourier - dual representation .
this uncertainty relationship is manifest in eq .
( [ eq : cvdjuncertainty ] ) , which relates @xmath127 , the encoding extent , to @xmath92 , the measurement extent . a natural extension of this work would be prove a lower bound perhaps exploring the techniques along the lines of @xcite from the perspective of different fourier transform pairs .
this uncertainty relationship appears to be a natural feature of the cv setting , but it could also be used to advantage .
there is likely to be oracle function symmetries that are particularly suited to different cv settings .
for example in sec .
[ sec : cvrep ] , we showed that balanced functions with a higher number of zero crossings create sinc functions with frequency components further away from @xmath136 .
it appears that an oracle identification problem designed to separate balanced functions according to frequency separation could be implemented in a cv setting and possibly provide advantage over classical or discrete quantum settings .
furthermore , it would be interesting to classify the balanced functions from the perspective of different coherent states  @xcite in cv parameterized settings of both finite and infinite dimensions .
the former would naturally involve the study of the coherent spin systems  @xcite .
furthermore the use of squeezed spin states should be studied  @xcite .
infinite dimension systems would naturally involve the study of implementations involving the coherent states of quantum optics  @xcite .
additionally , we have set up this framework in a manner that should allow any oracle identification problem to be analyzed in a similar manner in the cv setting .
an implementation of a discrete quantum oracle , for example  @xcite , requires a unitary operator representing the oracle .
provided we can create a diagonal representation of this oracle along the lines of @xmath50 given in eq .
( [ eq : qdjtargetlessoperator ] ) , our framework will naturally extend to it .
of course we need to be able to create an implementation of these oracles and that remains an important open question .
other avenues of the extension of this framework include cv implementations of other hidden subgroup problems .
the solution of simon s problem  @xcite in this setting would be an obvious starting point as would the exploration of a cv implementation of shor s algorithm  @xcite . additionally , the cv framework could be extended to include analysis of noisy oracles along the lines of  @xcite .
in closing we note that the transition from a discrete quantum information setting to a cv setting has many subtleties .
in particular the improper delta functions must not be used . limiting behaviour can be explored , but only if the limits are taken from the perspective of functions defined in the rigged hilbert space .
we appreciate financial support from the alberta ingenuity fund ( aif ) , alberta s informatics circle of research excellence ( _ _ i__core ) , canada s natural sciences and engineering research council ( nserc ) , the canadian network centres of excellence for mathematics of information technology and complex systems ( mitacs ) , and general dynamics canada .
ph is a scholar and bcs is an associate of the canadian institute for advanced research ( cifar ) .
d. akamatsu , y. yokoi , m. arikawa , s. nagatsuka , t. tanimura , a. furusawa , and m. kozuma , `` ultraslow propagation of squeezed vacuum pulses with electromagnetically induced transparency '' , * 99 * , 153602 ( 2007 ) .                                  c.  m.  caves , `` physical resources , entanglement , and the power of quantum computation '' , lecture given in squint summer retreat university of southern california , 2005 july 7 .
reference downloaded 03 dec .
2008 .
m.  adcock `` the classical and quantum complexity of the goldreich - levin problem with applications to bit commitment '' , _ master of science thesis , department of computer science , university of calgary ( 2004)_.                    r.  cleve , w.  van  dam , m.  nielsen , and a.  tapp , `` quantum entanglement and the communication complexity of the inner product function '' , lecture notes in computer science * 1509 * ( springer - verlag ) , pp .  61 - 74 ( 1999 )
.        m.  adcock and r.  cleve , `` a quantum goldreich - levin theorem with cryptographic applications '' _ proc .
19@xmath311 international symp
aspects comp .
( stacs 2002 ) _
, h. alt and a. ferreira , eds . , lecture notes in computer science * 2285 * ( springer - verlag ) , pp .  323 - 334 ( 2002 ) .