Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: quantum walks have been a particularly fruitful field of research in quantum information going back to ideas from feynman @xcite , meyer @xcite , and aharonov @xcite . one focus of recent work on quantum walks is their application to spatial search algorithms . it is well - known that grover s algorithm @xcite for searching on unstructured databases offers a quadratic improvement in search time over classical models . grover s algorithm is not designed to search physical systems where only local operations are possible and quantum walk algorithms have been employed for this purpose . in developing such algorithms , the major consideration is the search time attempting to match the quadratic improvement over the classical case offered by grover s algorithm . the first spatial search algorithm to do this , using the standard model of quantum walks , was developed by shenvi , kempe and whaley @xcite for searches on a hypercube . while there exist discrete - time searches on @xmath0-dimensional cubic lattices which are faster than classical searches for @xmath1 @xcite , effective continuous - time quantum searches only exist for @xmath2 @xcite or else they require additional memory ( in the form of spin degrees of freedom ) in order to improve their search time in lower dimensions @xcite . in @xcite , we demonstrated that effective searches over two - dimensional lattices may be achieved in an arguably simpler way which does not require extra degrees of freedom , and could , therefore , be viewed as more efficient . this is achieved through the choice of a different lattice , specifically , a honeycomb lattice which is the underlying lattice structure of carbon atoms in the material graphene . the association with graphene is important as , although we first study a purely theoretical problem in quantum information , the use of a graphene lattice also offers a potential physical realisation . we thus envisage using the quantum walk and quantum search algorithm framework to investigate the effect of perturbations on the dynamics on graphene and other carbon structures . this offers not only the possibility for demonstrating two - dimensional continuous - time quantum searching , but also paves the way for looking for novel effects in the material graphene .    in this paper , we offer a detailed account of the theory and numerical results thus expanding on the findings as presented in @xcite . generalisations of the results in @xcite have been given in @xcite , where a approach to solving the problem has been taken through encoding extra degrees of freedom into crystal lattices . recent experiments on artificial microwave graphene @xcite have shown that additional site perturbations , as discussed in sec . [ sec : alternativemarking ] , can be used to create a search protocol . searching on honeycomb lattices in discrete - time setting has been considered in @xcite ; however , no improvement over discrete time searches on cubic lattices was found . the paper is structured as follows : in section  [ sec : quantumwalksandsearching ] we give an introduction to the formalism of continuous - time quantum walks and the construction of quantum walk search algorithms . we also explain why previous search algorithms struggled on lower - dimensional lattices and why graphene offers a solution . section  [ sec : graphene ] will detail the relevant properties of graphene and set - up the notation we shall use throughout . sections  [ sec : threebondperturbation ]  &  [ sec : threebondperturbation_analysis ] contain our main analytical results , where we detail the specifics of our search algorithm and offer an analysis of the search running time and success probability . we shall then show in section  [ sec : communication ] how this search protocol can be adapted to demonstrate novel communication setups . sections  [ sec : alternativemarking ]  &  [ sec : alternativenanostructures ] contain numerical work demonstrating the possibility of using alternative methods of marking to create search behaviour and other carbon nanostructures . we conclude with a review and discussion of our results in section  [ sec : discussion ] . continuous - time quantum walks ( ctqw ) , first defined in @xcite , are the quantum analogue of continuous - time markov chains . they are defined purely on the state space , that is , the hilbert space @xmath3 , spanned by the states @xmath4 which represent the @xmath5 site of the lattice . thus , the time - evolution of such systems is defined by the schdinger equation @xmath6 where @xmath7 is the probability amplitude at the @xmath8 vertex of a system described by the state vector @xmath9 , and @xmath10 is the hamiltonian describing the connectivity of the lattice . note that we are using a dimensionless description setting e.g. @xmath11 . the dynamics of a quantum walk over a network is defined by the nature of the interaction between connected sites . therefore , the hamiltonian is generally constructed from the adjacency matrix of the underlying lattice . the adjacency matrix @xmath12 is defined as @xmath13 typically , the hamiltonian is chosen as @xmath14 . the parameter @xmath15 determines the coupling strength between connected sites and the parameter @xmath16 is an on - site energy that only enters the dynamics in a trivial way and thus can be set to a desired value . if @xmath17 and @xmath16 is equal to the valency of the lattice @xmath18 is the discrete laplacian . this form of hamiltonian is closely related to the tight - binding model for condensed matter systems @xcite . as first explained in @xcite , a quantum walk is transformed into a search protocol by introducing a localised perturber state , forming an avoided crossing in the spectrum of the search hamiltonian between an unperturbed eigenstate and the localised perturber . thus , initialising the system in the unperturbed eigenstate involved in the crossing and allowing the system to evolve in time , one finds that the system rotates into the localised perturber state .    the first ctqw search over @xmath0-dimensional cubic lattices @xcite introduced the localised perturber state as a projector onto a single site @xmath19 , resulting in the search hamiltonian @xmath20 here , @xmath21 is a parameter governing the strength of interactions between sites in the lattice . this parameter is chosen carefully @xcite such that the perturber state @xmath22 is brought into resonance with the ground state of the unperturbed hamiltonian , the uniform superposition @xmath23 . let us assume for the moment that all other unperturbed eigenstates are energetically sufficiently separated from the ground state  we will come back to this assumption later . then there will be an avoided crossing of two eigenvalues corresponding to the perturber state and the ground state . perturbation theory estimates that the energy splitting of these resonant states is proportional to the overlap of the localised perturber state @xmath22 and the uniform ground state @xmath23 , which scales as @xmath24 and that the corresponding eigenstates are of the form @xmath25 . by preparing the system initially in @xmath23 and allowing it to evolve for a set period of time @xmath26 , a measurement of the system will result in the state @xmath22 being measured with high probability . this explains the speed - up of the quantum walk search . the grover algorithm and its speed - up can be understood in analogous terms and one can show that all other states are indeed energetically very well separated .    in the present case this separation only holds above a critical dimension @xmath27 . a simple argument for the critical dimension can be obtained by comparing the scaling behaviour of the energy splitting @xmath28 with the energy separation of the first excited state from the ground state in the unperturbed lattice . for a cubic lattice we have a quadratic dispersion relation which allows us to estimate the energy separation @xmath29 for the first excited state in a @xmath0-dimensional lattice ) . for @xmath30 one then has @xmath31 as @xmath32 and a detailed analysis indeed proves that the quantum walk search works with optimal speed - up @xcite . at the critical dimension @xmath33 the two energy scales scale in the same way  the detailed analysis shows that the search still works but the dynamics is more complicated due to the interference of excited states . while there is still a speed - up it is only almost optimal ; the optimal search time @xmath34 gets multiplied with a logarithmically increasing factor . for the experimentally relevant regimes of 2- and 3-dimensional cubic lattices , all states participate in the dynamics as any avoided crossing gets dissolved completely as the number of sites grows . as a result no speed - up over classical searches is found . the above estimate offers a simple way how to reduce the critical dimension . if one can construct a search around a uniformly distributed state at an energy @xmath35 where the dispersion relation is conic ( linear in @xmath36 ) , then in @xmath0 dimensions @xmath37 . in @xcite this was implemented using a modified dirac hamiltonian . however , this requires the addition of a spin degree of freedom , essentially a doubling of memory , with the added complication that it is not immediately clear how such a system would be physically realised .    instead , our solution here is to change the lattice from a square ( cubic ) to a graphene lattice . this change of topology automatically implies the first step in our solution , as one of the important electronic features of graphene is the conic dispersion relation around the dirac energy . this arises naturally from the tight - binding descriptions of graphene . note that graphene has the critical dimension @xmath38 . indeed our construction as presented in @xcite has an almost optimal speed - up with logarithmic corrections that need to be evaluated in a detailed analysis that goes beyond the simple perturbative description given above . graphene is a single layer of carbon atoms arranged in a honeycomb lattice . the lattice is bipartite with two sublattices , labelled @xmath39 and @xmath40 , and a unit cell containing two carbon atoms . the spatial and reciprocal lattices are shown in figure  [ fig : both_lattices ] . the primitive vectors describing the lattice are @xmath41 , such that the position of a unit cell in the lattice is given by @xmath42 . we use dimensionless units in space where the distance between nearest neighbor sites is @xmath43 . the reciprocal lattice shows two important points , the dirac points @xmath44 and @xmath45 at the two inequivalent corners of the brillouin zone . , translation vectors @xmath46 and unit cell ( dashed lines ) . right : reciprocal lattice with basis vectors @xmath47 , symmetry points @xmath48 , @xmath44 , @xmath45 , @xmath49 and first brillouin zone ( hexagon).,title="fig : " ] , translation vectors @xmath46 and unit cell ( dashed lines ) . right : reciprocal lattice with basis vectors @xmath47 , symmetry points @xmath48 , @xmath44 , @xmath45 , @xmath49 and first brillouin zone ( hexagon).,title="fig : " ] [ fig : reciprocal_lattice ]    the energy spectrum of electrons in graphene was first derived by wallace @xcite when considering the band structure of graphite using a tight - binding hamiltonian @xmath50 where @xmath51 is the on - site energy ( which we will identify as the energy of the dirac points ) and @xmath15 the hopping strength ( both dimensionless in our setting ) . the tight - binding model for graphene and the derivation of the solution are well - known @xcite and give rise to the dispersion relation @xmath52 shown in figure  [ fig : infinite_graphene_lattice ] for an infinite graphene lattice .    as there are two atoms per unit cell the spectrum has two branches , the upper branch being the conduction band and the lower the valence band , which meet at the corners of the brillouin zone , the @xmath44-points . the energy at the @xmath44-points is @xmath16 which we name the dirac energy . it is around these points , that the behaviour of the spectrum is conical , that is , @xmath53 with a reduced density of states , a necessary feature for the creation of the search dynamics .    as the lattice possesses a translational symmetry the hamiltonian can be solved using linear superpositions of bloch functions over both sublattices . as a basis we use the orthonormal states @xmath54 to denote states on either the @xmath39 or @xmath40 sublattice in the cell at position @xmath55 . for the majority of what follows ( except in section  [ sec : alternativenanostructures ] ) , we will focus on finite - sized lattices with assumed periodic boundary conditions along the axes of both primitive vectors so that the topology of our lattice is a torus ; that is , our wavefunction is of the general form @xmath56 . our boundary conditions imply that the state vector must satisfy @xmath57 where @xmath58 denote the period of the lattice . thus , the wavefunctions on the torus take the bloch function form @xmath59    \ , , \label{eq : wavefn}\end{aligned}\ ] ] where @xmath60 is the number of sites in the lattice , @xmath61 is the momentum , and @xmath62 is a relative phase contribution dependent upon whether the state belongs to the conduction or valence bands ; it can be calculated by explicitly working through the tight - binding model . application of the periodic boundary conditions results in the following quantised momenta @xmath63 where @xmath64 and @xmath65 . in the following and whenever we consider quantum walk dynamics on a torus , the number of cells in each direction is generally chosen to be the same , that is , @xmath66 . this choice is purely for simplifying the notation ; alternative torus dimensions are possible as are other choices of boundary conditions corresponding to alternative carbon structures ( e.g.  nanotubes or a graphene sheet ) as will be demonstrated in section  [ sec : alternativenanostructures ] . for our choice of torus dimensions , we find there are momenta equal to the @xmath44-points and , consequently , eigenstates with energies equal to the dirac energy when both @xmath67 and @xmath68 are some multiple of 3 .    in fact , using the quantised momenta in eq . obtained for periodic boundary conditions , we find that there are four degenerate eigenstates with an energy that coincides exactly with the dirac energy when @xmath67 and @xmath68 are both multiples of 3 . these four states , known as the dirac states , can be constructed to live only on one of the sublattices , and are given by @xmath69 where @xmath70 for states on the @xmath39-sublattice or @xmath71 for states on the @xmath40-sublattice . in this section we will describe how a site is marked and will derive our optimal search starting state . our approach will be to analyse the system s spectrum and its dynamics in a reduced hamiltonian model involving only the relevant states from the avoided crossing . in the next section we will then validate this approach with a more detailed analysis , using the results obtained here as an initial guide . as already established , we introduce the localised perturber state in a region of the spectrum with a conic dispersion relation and thus a low density of states . simply altering the on - site energy in eq.([hgamma ] ) as done in @xcite does , however , not work here . as the on - site energy @xmath16 and the dirac energy are equal , one finds that the perturbation only interacts with the dirac states in the limit of zero perturbation strength , returning the unperturbed lattice . therefore , we choose an alternative perturbation method : namely , we modify the coupling strength between the site we wish to mark and its neighboring vertices . we focus in this section on changing the coupling to all three neighboring vertices of a particular site equally . our choice of perturbation matrix @xmath72 to mark the @xmath39-type vertex @xmath73 is then @xmath74 where the state @xmath75 is the symmetric superposition over the three neighbors of the perturbed site @xmath76 this leaves us with the search hamiltonian @xmath77 where @xmath21 is a free parameter . in what follows , we always set the on - site energy @xmath78 . considering our search hamiltonian , we can see that setting @xmath79 corresponds to a coupling strength of @xmath80 from the perturbed vertex and its nearest - neighbors ; our perturbation essentially removes the marked vertex @xmath73 from the lattice . note that vacancies are a common , naturally occurring , defects in graphene lattices @xcite .    in order to establish the critical value of @xmath21 , we numerically calculate the spectrum of @xmath81 as a function of @xmath21 , plotted in figure  [ fig:3bond_spectrum ] for a torus of dimensions @xmath82 . as @xmath72 is a rank-2 perturbation , we see in figure  [ fig:3bond_spectrum ] two perturber states working their way through the spectrum to an avoided crossing around @xmath83 when @xmath84 , that is , when the perturbed vertex is removed from the underlying lattice . 12x12_3p_spectrum_colour ( 54.25,26.75 )    as well as establishing a critical value for @xmath21 , we also find from this figure the states involved in the avoided crossing : the four dirac states at the dirac energy and the two perturber states which form our perturbation matrix @xmath72 . however , with further consideration , we may reduce the number of states involved further . as we have removed the site @xmath85 from the lattice it can no longer interact with the rest of the lattice and the corresponding state drops out . also , by direct calculation one finds that the @xmath40-type dirac states do not interact with the perturbation , that is , @xmath86 . thus , they remain an eigenstate of the search hamiltonian and do not interact with the perturbation . we are left with three states taking part in the avoided crossing : @xmath87 . we reduce our search hamiltonian in this three state basis at the critical point @xmath84 to obtain the following reduced hamiltonian describing the local dynamics at the avoided crossing @xmath88 with @xmath89 and @xmath90 . this reduced hamiltonian has eigenvalues @xmath91 , and eigenvectors @xmath92    using the eigenvectors of the reduced hamiltonian , we can construct a search starting state which is a superposition of dirac states    @xmath93    allowing our search starting state @xmath23 to evolve under the reduced hamiltonian we find    @xmath94    so that our system rotates from our dirac superposition @xmath23 to a state which is localised on the neighbors of the perturbed vertex @xmath75 in a time @xmath95 . thus , we find a @xmath96 search time , a polynomial improvement over the classical search time .    we plot in figure  [ fig : search_prob ] a numerically calculated search for a lattice with @xmath97 sites . the system has been prepared in @xmath23 and allowed to evolve under the full search hamiltonian from eq . ( [ eq : full_search_hamiltonian ] ) at the critical value @xmath84 . one finds that the system localises on the neighbors of the marked vertex and , for this particular system size , the probability for a measurement of the system to return one of the neighbors is around @xmath98 at its peak . this is orders of magnitude higher than the average probability to measure other vertices in the lattice , @xmath99 , which for this system size is around @xmath100 . the localised state interacting with the spectrum under the full hamiltonian has a tail into the rest of the lattice ; this leads to a loss of probability to be found at the neighboring vertices below @xmath101 . it is worth emphasising that the marked vertex plays no role in the search - it is removed from the lattice . rather , we force the system to localise on the nearest - neighbors , essentially making the marked vertex conspicuous by its absence . this feature differentiates our algorithm from previous searches , which have generally concentrated on localising directly on the marked vertex . an approach considering the neighboring vertices was detailed in @xcite , however , the focus was here still on localising on a single site but with classical post - processing steps considering the tail of the localised state extending into the neighborhood around the marked vertex . at this point , there are a number of issues which need to be mentioned . the first is that our starting state in eq . ( [ eq : starting_state ] ) contains information about the marked vertex in the form of the relative phase between the dirac states . however , this phase can only take three different values . therefore , we have three possible optimal starting states for an @xmath39-type perturbation and the same number for @xmath40-type perturbations ; there are thus in total six possible optimal states . as not all of these states are orthogonal , we only find an increase of necessary runs for a successful search by a factor of 4 ; this additional overhead is independent of @xmath60 and , therefore , does not affect the overall time complexity of the search . the particular representation of the dirac states in eq . ( [ eq : k_state ] ) was chosen in such a way as to make the calculations and conceptualisation of the system dynamics easier . in reality , constructing starting states which exist on one sublattice only , experimentally is extremely difficulty ; rather , in an experiment , one is likely to excite a superposition of dirac states , reducing the success probability , on average , by a factor of @xmath102 . our search is based on the conic dispersion relation in the spectrum and the @xmath103 scaling between successive energy levels in the linear regime , giving rise to the @xmath96 search time found in our reduced model . however , previous searches in continuous - time using a modified dirac hamiltonian @xcite or in discrete - time @xcite have found logarithmic corrections to the search time . as we have seen in eq . ( [ eq : search_time ] ) the search time is related to the energy gap at the avoided crossing . in the inset of figure  [ fig:3bond_spectrum ] we have numerically calculated the scaling of the gap in the spectrum of the full search hamiltonian and we indeed find evidence for a logarithmic correction . the lack of a logarithmic term in the search time for our reduced model is due to the neglecting contributions from the rest of the spectrum .    in what follows we establish a more accurate estimate of the running time and success probability by working with the full hamiltonian . we also justify the findings of our reduced model : while this model is insufficient to estimate the finer details , such as logarithmic correction terms , the accurate calculations show that our search algorithm indeed takes place mainly in a two - dimensional subspace spanned by the dirac states and a state localised on the neighbors of the marked vertex . we follow here closely our derivation as presented in @xcite which builds on ideas given in @xcite for a similar a similar calculation for regular rectangular lattices . we focus first on the search success amplitude , that is , @xmath104 where @xmath105 is the search time , that is , the time our search probability reaches a maximum , and our search hamiltonian @xmath106 is given by eq . ( [ eq : full_search_hamiltonian ] ) with @xmath84 , that is,@xmath107 here , @xmath108 , @xmath109 are the eigenstates and eigenenergies of @xmath106 . we assume , without loss of generality , a specific starting state @xmath110 where the marked vertex is chosen such that @xmath111 leading to @xmath112 in the following , we suppress the sublattice superscript as the analysis is the same regardless on which sublattice the perturbation lives . we also denote @xmath113 the positive eigenenergies of @xmath114 from eq . ( [ eq : graphene_unperturbed_eigenenergies ] ) and set @xmath83 .    for eigenstates @xmath115 with eigenenergies @xmath116 that are not in the unperturbed graphene spectrum ( @xmath117 for all points @xmath61 in the dual lattice ) , we may rewrite eq . ( [ eq : specific_search_hamiltonian ] ) in the form @xmath118 where @xmath119 and the phase of @xmath115 is chosen such that @xmath120 .    at this point , we can remove several states from the summation in eq . ( [ eq : amplitude ] ) . note that the basis state associated with the marked vertex , @xmath121 , is itself an eigenstate of the search hamiltonian with @xmath122 , and also , @xmath123 so that @xmath121 does not contribute to the time - evolution . we may also remove eigenstates of @xmath124 which are at the same time eigenstates of @xmath114 with the same energy . this can be seen in the following way . we first consider an unperturbed eigenstate @xmath125 such that @xmath126 . let us assume that there is an eigenvector @xmath115 of the search hamiltonian with the same eigenenergy @xmath109 , that is , @xmath127 . considering the matrix element @xmath128 , we find @xmath129 . as @xmath121 is an eigenvector of the search hamiltonian we know that @xmath130 . this leaves us with @xmath131 . as the unperturbed eigenstate @xmath125 is simply a bloch state we know @xmath132 . thus , we obtain @xmath133 . it is then clear that eigenstates of the search hamiltonian whose eigenenergies remain in the spectrum of @xmath134 do not play a role in the time - evolution of the search . we now derive an eigenvalue condition for those perturbed eigenvalues which are in the spectrum of @xmath124 . using the orthogonality of eigenstates @xmath135 , eq . ( [ eq : eigenstate_representation ] ) leads to @xmath136 by expressing @xmath137 in terms of the eigenstates of the _ unperturbed _ walk hamiltonian @xmath138 , we may write this as a quantization condition @xmath139 with @xmath140\ , , \ ] ] were @xmath60 is the total number of sites in the lattice . ( [ fe ] ) is written as two summations to incorporate the fact that the spectrum of @xmath114 as well as @xmath124 is symmetric around @xmath141 . choosing @xmath115 to be normalised @xmath142 , eq . ( [ eq : eigenstate_representation ] ) also implies @xmath143 which allows @xmath144 to be rewritten as @xmath145 we may now rewrite the amplitude in eq .   in the form @xmath146 where we have used the adjoint of eq . ( [ eq : eigenstate_representation ] ) . note again that eigenstates of @xmath124 , which are also in the spectrum of @xmath114 , do not contribute to the time - evolution of the search ; it may be seen from the definition of @xmath147 in eq . ( [ fe ] ) that @xmath148 , where @xmath109 is in the spectrum of @xmath114 . as we saw in the previous section , the avoided crossing is formed by two perturbers approaching symmetrically either side of the dirac states with energy @xmath83 . thus , we concentrate on evaluating the contribution to the time - evolution from these perturbed eigenstates of @xmath124 either side of the dirac point and we label these states @xmath149 . in order to evaluate these contributions we calculate the perturbed eigenenergy @xmath150 and also derive a leading - order expression for @xmath151 ( since the spectrum is symmetric , we have @xmath152 ) .    using the definition of @xmath147 in eq . ( [ fe ] ) , we estimate @xmath153 by separating out the dirac points , where @xmath154 , from the summations and then taylor expanding the remaining terms at @xmath141 , to find @xmath155 where the sums @xmath156 are given by @xmath157^{n } }      +       \frac{1}{\left[-\epsilon\left(\underline{k}\right)\right]^{n } }    \right ]    \,.\ ] ] as the unperturbed spectrum is symmetric only those @xmath156 with even @xmath68 are non - zero ; thus from now on we focus only on @xmath158 where @xmath159 . we stated earlier that the spectrum of graphene is well - approximated by a conic dispersion relation around the dirac points . thus , it is from around these points that the major contributions to the @xmath158 sums arise . by applying the linear approximation from eq . ( [ eq : linear_dispersion ] ) and the momenta quantum numbers from eq . ( [ eq : quantised_momenta ] ) we may approximate the @xmath158 summations as @xmath160 \nonumber \\     & + & o(1 ) \ , . \label{eq : i_epstein}\end{aligned}\ ] ] here @xmath161 is the epstein zeta - function @xcite @xmath162 for a real positive definite real symmetric @xmath163 matrix @xmath164 . for our purposes we use @xmath165 which describes the spectrum close to the dirac points . the linear approximation around both dirac points @xmath166 and @xmath167 is the same and , thus , the matrices @xmath168 and @xmath169 are equal .    our summation in eq . ( [ eq : i_epstein ] ) is over the rectangular regions @xmath170 and @xmath171 of the lattice @xmath172 . both are centered on @xmath173 and have side lengths proportional to @xmath174 , however the center , @xmath173 , corresponding to the relevant dirac point , is omitted from the summation . convergence of the epstein zeta function is well - known for @xmath175 @xcite , and leads to the bounds @xmath176 for @xmath177 . a sharp estimate for @xmath178 is given in appendix  [ appendix_a ] . in order to calculate an estimate for @xmath150 , we truncate the taylor expansion of @xmath147 in eq . ( [ eq : exp ] ) at the @xmath179 term ( the first term in the summation ) , and apply the eigenvalue condition @xmath180 . that is , we solve @xmath181 and obtain the approximation @xmath182 . this solution also leads to the estimate @xmath183 . we consider next whether our solution for @xmath150 lies inside the radius of convergence of the taylor expansion . we note that each term in @xmath184 is smaller than the corresponding term in @xmath185 for @xmath186 , and so it follows that @xmath187 for @xmath186 . this property of the epstein zeta function and the estimate from eq . ( [ eq : ikest ] ) imply @xmath188 i.e.  the infinite sum in eq . ( [ eq : exp ] ) converges for @xmath189 . thus , for large @xmath60 , our solution lies within the radius of convergence of our taylor expansion in eq . ( [ eq : exp ] ) .    in order to show that the dominant contributions to the search come from @xmath149 , we need to establish the leading order error term . all the @xmath158 sums are positive so that , given the sign of all the terms in the taylor expansion in eq . ( [ eq : exp ] ) , the true value of @xmath190 has to be smaller than the estimate we have obtained . thus , we write the true value of @xmath150 as @xmath191 with @xmath192 . one may rewrite @xmath193 , using the true value of @xmath150 , to give @xmath194 we follow the same arguments as used for the calculation of the radius of convergence to obtain an upper bound for the summation . this together with the already established fact that @xmath195 is inside this radius for sufficiently large @xmath60 , we get the following inequality : @xmath196 so @xmath197 . thus , we obtain @xmath198 it also follows that @xmath199    this shows that our perturbed eigenstates @xmath149 have an @xmath200-overlap with the starting state and are , therefore , the relevant states to be considered in the time - evolution of the algorithm . using the definitions of @xmath115 and @xmath144 in eqs . ( [ eq : eigenstate_representation ] ) , ( [ eq : ra ] ) , the inner product of the starting state and the perturbed eigenvectors can be expressed as @xmath201 where @xmath202 is the overlap of the starting state with the marked vertex state . applying our previous approximations for @xmath150 and @xmath151 , that is , for our perturbed eigenstates closest to the dirac point , we find @xmath203 thus , our starting state is indeed a superposition of the perturbed eigenstates , @xmath149 as assumed in the previous section . this makes it possible to investigate the running time and success amplitude of the algorithm in more detail , that is , @xmath204 it is clear from our earlier results for @xmath205 and @xmath179 that our algorithm localises on the neighbor state @xmath75 in time @xmath206 with probability amplitude @xmath207 . this confirms the logarithmic correction observed numerically and displayed in the inset of fig . [ fig:3bond_spectrum ] . it has been demonstrated in @xcite that discrete - time search algorithms can be modified to create a communication protocol . we show here that a communication setup can be established also in the continuous - time search algorithm with minor changes due to the subtleties of our search . we use the same unperturbed walk hamiltonian as before , that is , @xmath208 and the same type of perturbation matrix as in eq . ( [ eq:3bond_perturbation_matrix ] ) , but now at two different sites in the lattice . explicitly , our communication hamiltonian is @xmath209 where the perturbation matrices @xmath210 mark the source and target sites , respectively . as prescribed in @xcite , the communication protocol operates by preparing a state localised on the source perturbation , in our case on the neighbors of the source site , and allowing the system to evolve under our communication hamiltonian . the system then localises on the neighboring vertices of the target site . the effect of the two perturbation matrices @xmath210 amounts to disconnecting the two vertices from the underlying lattice .    as noted before , our search algorithm has several different optimal starting states depending on the position of the marked vertices . therefore , we divide our communication analysis into three different cases : the two sites i ) are on the same sublattice and have the same optimal search starting state ; ii ) are on the same sublattice but _ do not _ share the same optimal search starting state ; iii ) are on different sublattices . both of the first two cases can be treated by applying the reduced hamiltonian method demonstrated earlier , but communication between different sublattices is not tractable by this method and so we focus on numerics in this case .    for signal transfer between two sites on the same sublattice , we will assume , without loss of generality , that our communication takes place between two sites on the @xmath39-sublattice , @xmath211 ; here , the subscript @xmath212 or @xmath213 denotes the source or target vertices respectively . using arguments as employed in sec . [ sec : threebondperturbation ] , we can reduce the number of relevant states . ultimately , the search dynamics takes place in the subspace spanned by the basis @xmath214 , the @xmath39-type dirac states , and @xmath215 , the uniform superpositions over the neighbors of the source and target vertices . this basis leads to the reduced hamiltonian    @xmath216    where @xmath217 and @xmath218 .    in the first case considered we assume that there are two perturbations on the graphene lattice located at the points @xmath219 and @xmath220 , chosen such that @xmath221 . this implies that a search for either vertex , using our search algorithm from secs . [ sec : threebondperturbation ]  &  [ sec : threebondperturbation_analysis ] , would use the same optimal starting state . as the phases are equal , we drop the subscript on the perturbation coordinates in what follows . fixing our phases and diagonalising the reduced hamiltonian , we find it has eigenvalues @xmath222 and @xmath223 with eigenvectors    @xmath224    using these eigenstates , we may rewrite the source neighbor state as @xmath225 placing the system in the source state @xmath226 and allowing the system to evolve under the reduced hamiltonian , one finds @xmath227 we note that , in the last line above , the term in the brackets involving only the dirac states is actually the optimal search starting state for both perturbed vertices , defined in eq . ( [ eq : starting_state ] ) . the system thus oscillates between the states localised on the neighbors of the perturbed vertices , @xmath226 and @xmath228 , in a time @xmath229 , via their optimal search starting state . figure  [ fig : equivalent_communication_torus ] shows the system evolving under the full communication hamiltonian , eq .  ( [ eq : communication_hamiltonian ] ) . the initial state used for the time evolution shown in figure  [ fig : equivalent_communication_torus ] is the true localised state , that is , we run the quantum search with a single perturbation located at vertex @xmath219 until it reaches maximum success probability , and then apply the second perturbation to the vertex @xmath220 . the figure confirms the behaviour expected from the reduced model calculation . the communication mechanism essentially works in the same way as the quantum search algorithm , as it can be viewed as one marked vertex ` finding ' another . the initial localised source state decays back towards the search starting state , and the system then searches for the target state .    in our second case of signal transfer we analyse the behaviour of a communication system where we have two perturbed vertices on the same sublattice , but with the restriction that @xmath230 . as such , the two marked sites can not interact via the same search starting state . however , the optimal search starting states are not orthogonal so that signal transfer is still possible , but the resulting interference effects make the analysis slightly more complicated . we rewrite the coordinates of the target in terms of the source , @xmath231 and @xmath232 . again , fixing our phases and diagonalising the reduced hamiltonian in eq . ( [ eq : reduced_communication_hamiltonian ] ) , we find it has eigenvalues @xmath233 and @xmath234 with eigenvectors @xmath235 using these eigenstates , one may write the source perturbation as @xmath236 the full expression for the time - evolution is rather cumbersome ; therefore , we only show the terms and prefactors we are interested in , namely @xmath237 and @xmath238 @xmath239{\left| \ell_{s } \right\rangle } \label{eq : simulated_nonequivalent_communication } \\    & -    \frac{1}{2}\left[\cos\left(\lambda^{+}_{\sqrt{3}}t\right)-\cos\left(\lambda^{+}_{1}t\right)\right]{\left| \ell_{t } \right\rangle }    \nonumber \ , .\end{aligned}\ ] ] we can see here that the prefactors do not depend upon the coordinates of either the source or the target sites ; the transport signal between sites on the same sublattice but with different optimal search starting state is thus independent of the position of the source or target site .    in figure [ fig : both_nonequivalent_communication ] we show the system evolved under the full communication hamiltonian . again , the initial state is the true localised state on the nearest - neighbors of the source vertex , obtained by running the search algorithm using one marked vertex until it reaches its peak success probability . the time - evolution is radically different to the previous communication case ; the behaviour here is erratic with uneven peaks of probability at the two perturbations involved in the protocol . however , there are still significant probability revivals . the transport behaviour from eq . ( [ eq : simulated_nonequivalent_communication ] ) , calculated using the reduced hamiltonian , is also shown in figure  [ fig : both_nonequivalent_communication ] . the probability at time @xmath240 has been scaled to match that shown in figure  [ fig : both_nonequivalent_communication ] . our calculated behaviour has the same signal pattern as the numerically calculated behaviour from the full hamiltonian , although over a shorter period of time . as our reduced model only makes use of the dirac states and the perturber states , we lose the contribution to the time - evolution from the rest of the spectrum giving rise to logarithmic corrections as discussed in sec . [ sec : threebondperturbation_analysis ] . this leads to differences in the overall time scales for the reduced model and the full hamiltonian . this also supports our findings that signal transfer between all non - equivalent vertices on the same sublattice is the same . the communication protocol is again set up by a search mechanism in reverse , where one vertex finds another . the slightly erratic behaviour that emerges here is due to interference between the two separate search mechanisms which interact due to the non - zero inner product of the three possible optimal search starting states for vertices on the @xmath39-sublattice .    in our final case of signal transfer , we consider communication between sites on different sublattices which we can not treat in the reduced hamiltonian model . it has been demonstrated previously that perturbations to one sublattice do not interact with the dirac states which live on the other sublattice . therefore , when attempting to reduce the communication hamiltonian as done before , we merely find that it decouples into two non - interacting reduced hamiltonians , each describing a search protocol on one sublattice . the interaction between the two sublattices is here facilitated due to interactions via the bulk of the spectrum . in what follows , we focus on numerics and inspect the behaviour for systems involving the same source perturber but different targets . the numerics show a finite number of different signal patterns , two examples are shown in figure  [ fig : both_differentsublattice_communication ] . we can see from these examples that the communication mechanism takes place over a much longer timescale with a superimposed oscillatory dynamics . the increase in timescale can be attributed to the weak nature of the interaction between the two perturbations due the fact that the localised states live mainly on one sublattice . the signal pattern in figure  [ fig : both_differentsublattice_communication ] thus decouples into a fast oscillation between one site , say the source site , and the delocalised lattice state and a slow time scale on which a small amount of probability amplitude escapes into the other sublattice due to the weak interaction . this process continues until the recurrence probability at the target perturbation reaches the same peak as the initial localised source state , and then the behaviour reverses . in the previous sections , we have discussed a method of marking a vertex through modifying the hopping potential from the perturbed site to all three of its nearest - neighbors . here , we will discuss alternative methods of marking a vertex which still keep the perturber interaction in the conic dispersion region of the spectrum . this can be done by altering the hopping potential in different ways as discussed above ; this approach is necessarily a rank-2 perturbation and these two perturber states must meet at the dirac energy . we will discuss several ways of applying such a perturbation in the next paragraph . another approach is based on coupling extra sites to the lattice ; this set - up will be discussed at the end of this section . focussing on hopping potential perturbations , we may consider perturbing the bonds to any number of nearest - neighbors to any strength . we have seen previously that a symmetric three - bond perturbation successfully creates a search . here , we demonstrate a search based on perturbing the hopping potential from a given site to only one of its nearest - neighbors . that is , we use the same search hamiltonian as in eq . ( [ eq : full_search_hamiltonian ] ) with the perturbation matrix @xmath241 and eigenstates @xmath242 note that it no longer makes sense to speak of a single marked vertex . rather , our perturbation marks both vertices @xmath243 simultaneously ; thus , it may be more accurate to view our perturbation matrix as marking a single cell of the lattice .    for this type of perturbation , the avoided crossing used to generate search behaviour is not necessarily at @xmath84 , see figure  [ fig : singlebond_spectrum ] , where the spectrum of our search hamiltonian as a function of @xmath21 for a @xmath244 cell torus is shown . as stated before , our single bond perturbation is also a rank-2 matrix , and so again there are two perturber states approaching the spectrum from the negative and positive regions of the spectrum . inspecting the region around the dirac energy , see inset in figure  [ fig : singlebond_spectrum ] , we see that avoided crossings are formed by four states in total , the two blue ( solid ) curves and the two red ( dashed ) curves . at @xmath245 there is an exact crossing between the red and blue curves indicating that these states are orthogonal . singlebond_spectrum ( 54.26,26.75 )    we can make this avoided crossing picture clearer by breaking down the eigenstates of the search hamiltonian in terms of the symmetries of the lattice . in particular , we use the rotation operator @xmath246 , which is a rotation of @xmath247 about the mid - point of the perturbed bond . considering the action of @xmath248 on the marked vertices and the eigenstates of the perturbation matrix , one finds @xmath249 , @xmath250 and @xmath251 . the inset of figure  [ fig : singlebond_spectrum ] shows the states nearest to the dirac energy , where the same colour indicates the same @xmath248 parity . it now becomes clear that the states near the dirac point actually form two avoided crossings , one for each parity with a minimum energy gap at @xmath245 .    without going into details , we can again find optimal starting states using a reduced hamiltonian approach ( one for each of the avoided crossings ) ; the evolution of the system using these optimal states is shown in figure  [ fig : singlebond_search ] for a @xmath244 cell lattice . we find there is a significant localisation on the marked vertices and their nearest - neighbors , with the probability of being found on either of the marked vertices peaking at around @xmath252 . we can also see that the probability of being found on each of the nearest - neighbors peaks at around @xmath253 , resulting in a total probability of being found on the marked vertices and their nearest - neighbors of approximately @xmath254 . the success probability fluctuates , as the probability amplitude oscillates between the marked vertices and their nearest - neighbors . this is due to the probability amplitude being constrained in the local area by the increased hopping potential between the two marked vertices . although our demonstration is for a @xmath244 cell lattice , numerical investigations show that the search behaviour remains the same as the lattice size increases with critical value fixed at @xmath245 . we also find that the gap at the two avoided crossings scales as @xmath255 , in the same way as for the three - bond perturbation shown in figure  [ fig:3bond_spectrum ] . as the search time is inversely proportional to the energy splitting at the avoided crossing , it also gives an estimate of the running time of the search @xmath256 . we now turn to coupling additional sites to the lattice as a way of introducing a perturbation as also considered in @xcite ; such a treatment is in many ways closer to experimental realisations as defects due to additional add - on atoms are quite common in graphene @xcite . we focus here on the idealised case where an additional site is coupled to a single lattice vertex . we use the perturbation matrix @xmath257 for coupling the additional site @xmath258 to the @xmath39-type vertex @xmath73 and @xmath21 is a free parameter related to the on - site energy of the additional vertex . the choice of the coupling terms is such that the binding energy between the additional site and the lattice vertex is the same as the internal couplings in the lattice . thus , our search hamiltonian is of the form @xmath259 in figure  [ fig : singlesite_spectrum ] , the spectrum of the search hamiltonian is given as function of @xmath21 . as our free parameter @xmath21 only changes a single term , our spectrum only has a single perturber state . one finds a clear avoided crossing around the dirac energy ( @xmath141 ) when @xmath260 , that is , when the on - site energy of the additional site matches the on - site energy of the lattice vertices . singlesite_spectrum ( 11.75,28.25 )    in the reduced hamiltonian picture , we write eq . [ eq : additional_site ] in terms of a basis consisting of the dirac states and our perturber state , @xmath258 . only three states are involved in the search , that is , @xmath261 , and @xmath10 reduces , up to a different prefactor , to the same @xmath262 matrix found in eq . ( [ redh ] ) for the three - bound perturbation . thus , our previous reduced hamiltonian analysis holds for this case and the optimal search starting states are the same ; the time - evolution for the additional site search is shown in the inset in figure  [ fig : singlesite_spectrum ] . the only differences between this case and the three - bond marking are that the system localises on the additional site , and also the change in prefactor leads to a search time of @xmath263 .    we note that other , more realistic types of perturbations involving additional sites coupling to a lattice site and its neighbors , can be shown to result in effective search protocols . also , the single additional site perturbation described here can , like the three - bond perturbation , be used to setup a communication protocol , with the same signal patterns and behaviours as found previously . the initial state is in this case completely localised on the additional site . so far we have described the dynamics of searches on graphene lattices with periodic boundary conditions , that is , graphene on a torus . here we consider more realistic boundary conditions such as nanotubes and graphene sheets . for the dynamics on nanotubes , we move to periodic boundary condition along one axis only and impose dirichlet boundary conditions along the other directions . the properties of nanotubes are well - known @xcite , and it has been shown that the band structure of armchair nanotubes , that is , nanotubes with armchair boundaries , always allows for an energy at the dirac energy regardless of the nanotube diameter ; we will focus on these types of nanotubes . we are interested in searching on finite length nanotubes , where the band structure becomes discrete ; the length of the nanotube are in addition chosen such that there exists an eigenenergy at the dirac energy . an example of the cell we use to construct finite armchair nanotubes is shown in figure  [ fig : armchairnanotube_example ] . we choose the finite length of the nanotube to be along the horizontal axis and we close the underlying graphene lattice into a nanotube along the vertical axis . we denote the basis states of sites in our nanotube as @xmath264 , where @xmath67 indicates the @xmath265 @xmath266-type vertex in the horizontal direction in the @xmath267 cell . it is simple to see that , due to the periodic boundary conditions along the circumference of the tube , the vertical component of the eigenstates must be bloch states and the horizontal component of the amplitudes must be sinusoidal in nature . that is , the eigenstates of the finite nanotube are standing waves along its length .    by working through the tight - binding model for this system , one can show that @xmath268 is the only point where there exists an eigenenergy equal to the dirac energy . it is also possible to show that the spectrum includes dirac points when the number of sites along the length of the nanotube @xmath269 , where @xmath270 is an integer . as there is only one potential dirac point for finite armchair nanotubes , it follows that there are only two dirac states ( one from the bonding and the anti - bonding regions of the spectrum ) . as we have @xmath271 at the dirac point , the bloch wave around the circumference of the nanotube is simply a uniform superposition . another important feature of the dirac states on the nanotube , which we have not encountered previously , is the existence of nodal points where the amplitude of the eigenstate is 0 . we find these nodal points occur at every third site along the horizontal axis . we focus on applying our three - bond perturbation from sections  [ sec : threebondperturbation ]  &  [ sec : threebondperturbation_analysis ] to the armchair nanotube . in our new labelling our perturbation takes the form @xmath272 we have assumed that we are perturbing an @xmath39-type vertex on an even horizontal coordinate , so that the perturbed site and its nearest - neighbors remain within one nanotube cell . while it is easy to re - write this perturbation matrix for other sites we restrict ourselves to this form for simplicity . note even in this form it can still be expressed in terms of a marked site , @xmath273 , and a state which lives on the neighbors , @xmath274 . our search hamiltonian for the nanotube is the same as in eq . ( [ eq : full_search_hamiltonian ] ) . inspecting the spectrum of the search hamiltonian as a function of the free parameter @xmath21 where the perturbation is not located on a nodal point of a dirac state , one finds an avoided crossing around the dirac energy when @xmath84 . this spectrum is very similar to the one shown in figure  [ fig:3bond_spectrum ] for the search on the torus and will be omitted . if the perturbation is located on a nodal site the picture is different , however ; one finds that there is an exact crossing at the dirac energy when @xmath84 and searching is not possible . recall that the effect of the perturbation at the critical point is to completely remove the site from the lattice , if the site already has zero amplitude then the perturbation will not interact with the dirac state . similar to the previous section , we numerically reduce the full search hamiltonian in a basis consisting of the two dirac states and the state living purely on the neighbors , @xmath75 . the perturbed site , @xmath273 , is decoupled from the lattice . through this process , we find two possible starting states , one for each sublattice . the starting states are weighted superpositions of the dirac states over the sublattice containing the marked vertex . the nodal points are excluded from this treatment . figure  [ fig : armchair_search ] shows the system evolution for two searches using the numerically found optimal starting states . one search has the marked site located in the centre of the nanotube , the other has the marked site positioned near the edge of the nanotube . the dimensions of the nanotube have been chosen so that there are 320 sites , comparable to the searches shown in earlier sections . comparing these searches to the three - bond perturbation search on the torus in figure  [ fig : search_prob ] , we see that there is a marked difference in behaviour and success probability induced by relaxing the periodic boundary conditions along one axis .    using a perturbation located near the edge , the lower image in figure  [ fig : armchair_search ] , we see a reduction in success probability by a factor of 2 - 3 and strong fluctuations in peak height when compared to the search on the torus . searching in the bulk , shown in the lower image of figure  [ fig : armchair_search ] , displays behaviour closer to searches on a graphene torus , but again with significant fluctuations at each peak . we propose that this effect is due to the reflection of probability amplitude from the edges of the nanotube . this is supported by the changes in the interference pattern in the signal as the perturbation is moved across the lattice . numerics demonstrate that additional site perturbations and communication protocols can also be used on nanotube lattices . we now move to working on graphene sheets , that is , removing the periodic boundary conditions along both axis . a detailed account of a specific version of this set - up together with experimental results has been presented in @xcite for graphene sheets consisting of armchair boundaries only . note that such a configuration can not be achieved on rectangular sheets ; the simplest configuration has the form of a parallelogram , see @xcite for further details . the advantage of such a configuration is that the form of the boundary does not admit so - called ` edge - states ' . in the following , we will look at the influence of these edge - states in more detail by considering rectangular graphene sheets , such as in figure [ fig : finite_lattices_example ] . in addition to armchair edges , here along the vertical axis , these sheets have boundaries formed by bearded edges ( left ) and zigzag edges ( right ) along the horizontal axis . these boundaries support edge states , i.e.  states localised along these edges , with an eigenenergy close to the dirac energy @xcite . in the following , we will investigate , how the existence of these edge states influences the search . as in the case of the finite armchair nanotube , the imposition of dirichlet boundary conditions at an edge generates sinusoidal eigenstates . as a result , there are no extended ( bulk ) eigenstates at the dirac energy due to the inability to equate the quantised momenta with the necessary points in @xmath61-space . we thus need to find other extended eigenstates in the sea of edge states near the dirac point to undertake a search in this set - up . we mark sites using the triple - bond perturbation such as in eq . ( [ eq : full_search_hamiltonian ] ) . throughout this section we choose the dimensions in terms of primitive cells of the graphene sheets , @xmath275 ; a bearded edge sheet thus consists of @xmath276 sites and a zigzag sheet is formed of @xmath277 sites . the spectra of the search hamiltonians for the bearded lattice is shown in figure  [ fig : spectrum_bearded_3p ] constrained to the energy region of interest . the spectrum for zigzag edges is very similar and is not show here . one can clearly see an avoided crossing around the dirac energy at @xmath84 , the critical value for this type of perturbation . there are several states very close to the dirac energy ; these are the edge states which are all non - degenerate . therefore , it is not possible for us to construct a superposition of degenerate eigenstates which is optimal for searching . rather , our initial starting state must be a single eigenstates of the unperturbed hamiltonian @xmath278 .    using the search hamiltonian from eq . ( [ eq : full_search_hamiltonian ] ) and fixing @xmath84 , we proceed by allowing the system to evolve after being prepared in one of the unperturbed eigenstates . one finds for both types of lattice , that the edge states near the dirac energy fail to produce localisation behaviour ; the probability at the neighboring vertices of the marked site do not rise above noise levels . rather , search behaviour only begins to emerge when we use the first delocalised , non - edge state as our initial state . we also note that we _ can not _ search for sites near the zigzag or bearded edges , where the edge states exist . only as we move further into the bulk of the lattice or along an armchair - type edge , the localisation behaviour returns . note that the edge states can be viewed as a kind - of one - dimensional system and , as we saw from the scaling argument towards the end of sec . [ sec : quantumwalksandsearching ] , one - dimensional systems imply an energy spacing between successive energy levels of @xmath279 . thus , the perturber state interacts with many states in a dense part of the spectrum and the search fails .    in figure [ fig : bearded_searches ] we show the search behaviour arising from marked sites placed at different positions on a bearded graphene lattice ( zigzag lattice types display similar behaviour ) . the results are similar to those found for the nanotube searches . there is a slight increase in variation of signal pattern with position , and an increase in success probability maxima as we move towards the centre of the lattice from either of the bearded edges . we have shown that continuous - time quantum search can be done effectively on a two - dimensional graphene lattice without the use of internal degrees of freedom . this is achieved by making use of the conical ( linear ) dispersion relation in the graphene spectrum . the search succeeds in time @xmath280 with probability @xmath281 . this is the same time complexity found in @xcite for discrete - time searches and in @xcite for continuous - time searches . to boost the probability to @xmath282 , @xmath283 repetitions are required giving a total time @xmath284 . amplification methods @xcite may be used to reduce the total search time further . our main result focusses on perturbations which involve altering the hopping potential from a marked site to all three of its nearest - neighbors equally . we have also demonstrated other types searches based on perturbing the hopping potential in a cell and the adding extra sites . we have shown that search mechanisms can be utilised for the purposes of signal transfer . our findings point towards applications in directed signal transfer , state reconstruction , or sensitive switching . this opens up the possibility of a completely new type of electronic engineering using single atoms as building blocks of electronic devices . our results demonstrate that a range of nanostructures constructed from graphene could be used to this end . we thank julian bhm , ulrich kuhl and fabrice mortessagne for discussions on how the ideas of @xcite could be adapted to more physical systems , which lead to some of the work on additional perturbation types and graphene flakes . we give here some details regarding the derivation of eq . ( [ eq : i2est ] ) and the logarithmic divergence of @xmath179 . it is clear that the dominant contributions to the @xmath158 summations come from the vicinity of the dirac points . approximating the spectrum close to the dirac points , one has @xmath286 + o(1 ) \ .    \end{split}\ ] ] here the sums over integers @xmath287 and @xmath288 is over a rectangular region @xmath170 of the lattice @xmath172 which is centered at @xmath173 and has side lengths proportional to @xmath174  the center @xmath173 , corresponding to the relevant dirac point , is omitted from the sum . as stated in the main text , for @xmath289 the corresponding sums converge which proves eq .  ( [ eq : ikest ] ) . + for @xmath290 we will establish constant @xmath291 and @xmath246 such that @xmath292 which then directly leads to eq . ( [ eq : i2est ] ) . to establish @xmath291 note that because each term in the sum ( [ logineq ] ) is positive its value decreases by restricting it to a square region @xmath293 , @xmath294 which is completely contained in @xmath170 . up to an error of order one the sum over a square region can in turn be written as a sum over eight terms of the form @xmath295 for fixed @xmath287 we can find @xmath296 such that @xmath297 we may choose @xmath298 for some constant @xmath299 , so @xmath300 which diverges as @xmath301 .    establishing @xmath246 and the corresponding logarithmic bound from above , and thus proving eq . ( [ eq : i2est ] ) , follows the same line by first extending the sum to a square of side length @xmath302 that completely contains @xmath170 and then establishing @xmath303 with @xmath304 . d. aharonov , a. ambainis , j. kempe , and u. vazirani . quantum walks on graphs , in _ proceedings of the thirty - third annual acm symposium on theory of computing _ , ( stoc 01 , acm , new york , ny , usa , 2001 ) , pages 50 - 59 . l. k. grover . a fast quantum mechanical algorithm for database search , in _ proceedings of the twenty - eighth annual acm symposium on theory of computing _ , ( stoc 96 , acm , new york , ny , usa , 1996 ) , pages 212 - 219 . a. ambainis , j. kempe , and a. rivosh . coins make quantum walks faster . in _ proceedings of the sixteenth annual acm - siam symposium on discrete algorithms _ , ( society for industrial and applied mathematics , philadelphia , pa , usa , 2005 ) , pages 1099 - 1108 . j.  boehm , m.  bellec , f.  mortessagne , u.  kuhl , s. barkhofen , s.  gehler , h.j . stckmann , i.  foulger , s. gnutzmann , and g.  tanner . experiments on ` quantum ' search and directed transport in microwave artificial graphene , _ phys . _ , * 114*:0 110501 , 2015 . a. ambainis , a. backurs , n. nahimovs , r. ozols and a. rivosh . search by quantum walks on two - dimensional grid without amplitude amplification . in _ theory of quantum computation , communication , and cryptography _ , ( springer berlin heidelberg , berlin , germany , 2013 ) , pages 87 - 97 .
Abstract:

Answer: quantum walks have been very useful in developing search algorithms in quantum information , in particular for devising of spatial search algorithms . however , the construction of continuous - time quantum search algorithms in two - dimensional lattices has proved difficult , requiring additional degrees of freedom . here , we demonstrate that continuous - time quantum walk search is possible in two - dimensions by changing the search topology to a graphene lattice , utilising the dirac point in the energy spectrum . this is made possible by making a change to standard methods of marking a particular site in the lattice . various ways of marking a site are shown to result in successful search protocols . we further establish that the search can be adapted to transfer probability amplitude across the lattice between specific lattice sites thus establishing a line of communication between these sites .