the physics of many 1 + 1 dimensional lattice models at criticality is captured in the continuum limit by two - dimensional conformal field theory ( 2d cft ) . as a result
, one can numerically extract the conformal data , for instance scaling dimensions of the primary operators from the scaling of the two - point correlation functions in the critical lattice model , see @xcite for the standard reference . on the other hand , measures of entanglement proved to be useful quantities for exploring the critical points of many - body systems .
for example , in the ground state of a critical chain , the entanglement entropy of a block of spins has a universal logarithmic scaling with the size of the block that is proportional to the central charge @xcite .
this provides an efficient numerical way to obtain the central charge of the cft that governs the critical point , see refs .
@xcite for review .
it is then natural to ask if entanglement measures in excited states can be used to extract more cft data , such as for instance the modular s or t matrices or quantum dimensions , numerically .    in rational cfts local elementary excitations
are catalogued into finite number of conformal families containing primary operators and their descendants .
the simplest excited states can then be obtained by inserting local cft operators at some spatial points . in such states
, one can study how the local operator changes the structure of entanglement in the ground state .
more precisely , it is possible to compute the time evolution of the change in rnyi entropies for a reduced density matrix of a single interval due to the operator insertion @xcite  see ref .
@xcite for some results in various cft setups and ref .
@xcite for entanglement in a related class of globally excited states . in 2d cft
this analysis can be preformed analytically and rnyi entropies detect an increase in entanglement equal to the logarithm of the quantum dimension of the conformal family @xcite .
having such a clear and elegant prediction from the cft , it is then natural to wonder if and how the logarithms of quantum dimensions are reproduced on the lattice . in this article
we initiate such program for the simplest case of the critical ising chain .
the advantage of the ising model is that it is exactly solvable and the action of a family of local operators can be efficiently simulated for large system sizes . the main subtle point of this analysis is the identification between the cft operators and their lattice counterpart .
in fact , there are only few models where such map is well established
see for instance the discussion in @xcite  and a general belief is that a given lattice operator corresponds in the continuum to a primary operator plus its descendants . in the case of the two - point functions these extra contributions from descendants lead to corrections that are suppressed as higher powers with the distance .
in this work , given a well established identification for the ising model operators , we will be able to check the contribution form this non - unique identification to the physics of entanglement propagation .
we will see that for truly local operators on the lattice  like the ising spin  we recover the cft answer , but for operators with non - local support  like ising energy
the subleading contributions modify the leading answer and lead to a mismatch .
the computations of rnyi entropies in cft are done using the replica method that , for excited states , boils down to calculation of correlation functions on complicated riemann surfaces . even for states locally excited by
more than a single operator such objects are notoriously difficult to compute analytically and features of entanglement measures in this class of states remain unexplored .
similarly , measures of distance between quantum states like , e.g. , relative entropy for locally excited states require the access to higher - point correlators @xcite . in this work
, we will further explore the ising model to shed a new light in these directions by numerically performing the time evolution of the relative entropy , as well as rnyi entropies in more general states excited by multiple local operators .
this paper is organized as follows . in section [ sec : cft ]
we summarize the relevant results from the two dimensional cft . in section [ sec : num ] we present our main numerical results for the evolution of entanglement in the critical ising model in states excited by single local operators . in section [ sec : gen ] we consider the evolution away from the critical point , as well as more general operator excitations . in section
[ sec : rel ] we present the evolution of the relative entropy after local operator excitation in this model .
finally , we conclude and present the details of our numerical approach in the appendix [ numer ] .
in this section we briefly review the existing results for the evolution of rnyi entropies in locally excited states in 2d cfts and then show some details of the computation for the ising cft .
finally , we discuss the minor modifications that appear for the cfts on the cylinder that , in the following sections , we will be comparing to numerics from the periodic chain .
let us start with a 2d cft on the real line and a state excited by a local operator @xmath0 at distance @xmath1 from the entangling interval @xmath2 $ ] , as presented pictorially on fig .
[ fig : schematic ] .
the density matrix is given by ( t)&=&ne^-ite^-(0,-l)|00|^(0,-l)e^-e^it & & n(w_2,|w_2)|00| ^(w_1,|w_1 ) , [ denmat ] where the insertion points of the local operators are defined as & & w_1=i(-it)-l ,  w_2 = -i(+it)-l , + & & |w_1=-i(-it)-l ,  |w_2=i(+it)-l .
[ points ] the factor of @xmath3 is the uv regulator for the local operators and we take @xmath4 at the end of the computation .
the normalization @xmath5 ensures that the trace of the density matrix is equal to 1 .
+ using the replica trick , we can compute how a family of rnyi entropies @xmath6 , indexed by integer @xmath7 , changes due to the local operator insertion . the answer to this question
is expressed in terms of the logarithm of the ratio @xcite s_a^(n ) , [ eq : delta ] where the correlator in the numerator is computed on the @xmath7-sheeted surface @xmath8 with cuts on each copy corresponding to interval @xmath9 , and the two - point function in the denominator is on a single sheet @xmath10 with an interval cut @xmath9 . for a detailed derivation and further illustrative explanations
see @xcite .    in 2d cft
, one can apply a conformal map from @xmath8 to a complex plane and evaluate the correlators explicitly .
it turns out that the answer is universal and the increase in the rnyi entanglement entropies is equal to a constant that is the same for all the members of a conformal family , i.e. primary operators and their descendants @xcite . in rational cfts
this constant is equal to the logarithm of the quantum dimension of the local operator @xcite which is defined as d_a=,[qdim ] where @xmath11 denotes the elements of the modular s - matrix of the cft , see e.g. ref .
@xcite .    in this work
we focus on the 2d ising model which is the @xmath12 minimal model with three primary operators : the identity @xmath13 , the energy @xmath14 with conformal dimensions @xmath15 @xmath16 and the spin @xmath17 with @xmath18 @xmath19 , where @xmath20 marks the total scaling dimensions .
+ the modular s - matrix of the ising model is given by s= (    ccc 1 & 1 & + 1 & 1 & - + & - & 0 +    ) , so the three quantum dimensions are d_=d_=1,d_=. this way , at criticality , only excitations with primary @xmath17 can non - trivially change the entanglement in the vacuum state , and for all the rnyi entropies we have @xcite s^(n)_a=. the standard , chiral ( anti - chiral ) descendants are obtained by either acting on the primary operators with chiral ( anti - chiral ) derivatives @xmath21 ( @xmath22 ) or taking the operator product expansion ( ope ) of the primary operators with the energy momentum tensor .
such descendants increase the entropies by the same amount as the primaries . if we however act with the linear combination of the two derivatives , there is an additional contribution to the entropy equal to @xmath23 @xcite .
we will see this in case of the spatial derivative @xmath24 acting on @xmath25 , which we consider in sec .  [
sec : num ] c.    .,width=264 ]    in order to compare the cft results with numerics we have to take into account finite size of the system , @xmath26 . in the cft computation it enters through the invariant cross - ratios .
let us , for simplicity , consider the change in the second renyi entropy @xmath27 that requires the correlator on two cylinders .
the correlators entering are computed using a composition of the conformal map @xmath28 from each cylinder to the plane with a cut and the uniformization map @xmath29 .
after some standard cft manipulations , the change in the entropy can be written as s^(2)_a=-(|z(1-z)|^4h_^_(z,|z)),[secren ] where @xmath30 is the canonical 4-point function on the complex plane with operators @xmath31 inserted at @xmath32 and the cross - ratios are defined as z=,|z= , with @xmath33 , @xmath34 and similarly for @xmath35 . from the conformal map , we can also show that @xmath36 and @xmath37 ( similarly for @xmath35 ) .
once we plug the insertion points of the operators , in the @xmath4 limit , the cross - ratios become periodic functions of time , as shown on fig .
[ fig : cr1 ] .
the difference with the cft on the infinite line is that in one cycle of time equal to @xmath26 , both @xmath38 and @xmath35 reach their maximal value of @xmath39 .
more precisely , @xmath40 in the time @xmath41=[l+\alpha n , l+l+\alpha n]$ ] and zero outside , whereas @xmath42 inside @xmath43=[n-(l+l)+\alpha n , n - l+\alpha n]$ ] and zero outside , for integer period @xmath44 . in order to extract the increase in the second rnyi entropy analytically in the @xmath4 limit , we simply take the limit of @xmath45 or @xmath46 in in the appropriate time intervals .
let us focus on the explicit example of the ising model . for the @xmath17 operator , from the fusion rule @xmath47
, the correlator can be decomposed as @xcite ^_(z,|z)=(c^1_)^2|^_(1|z)|^2+(c^_)^2|^_(|z)|^2 , with the three - point coefficients @xmath48 and @xmath49 , as well as the conformal blocks ^_(1|z)=,^_(|z)=. we can then check that that the non - zero contribution in the two time intervals comes from the identity block and is equal to @xmath50 , in accordance with the modular s - matrix elements .
this behavior is naturally explained from the quasi - particle picture where left and right moving sets of quasi - particles propagate from the insertion point of the operator and , on the circle , there are two time intervals where either left or right particles are inside the entangling interval @xmath9 .    on the other hand , for the @xmath14 excitation , using the fusion @xmath51 , we can write the correlator as ^_(z,|z)=(c^1_)^2|^_(1|z)|^2 , with @xmath52 and conformal block ^_(1|z)=. clearly , inserting the cross - ratios and taking @xmath4 yields @xmath53 for all times .
this suggests that the quasiparticles produced by @xmath14 are in a product state .
after this short review , below we perform the numerical analysis and compare how the cft predictions are reproduced on the discrete chain .
we consider quantum ising model in a transverse magnetic field on a 1d chain of @xmath26 spins-@xmath54 described by the hamiltonian h=-_n=1^n , [ eq : spin_ising ] where @xmath55 are the standard pauli matrices acting on @xmath7-th spin and we assume periodic boundary conditions @xmath56 .
the model is critical for @xmath57 and unless stated otherwise , in this article we set @xmath58 .
the model can be mapped onto the free fermion system using the jordan - wigner transformation and our numerical results are performed in such a setup .
we refer to the appendix [ numer ] for details . after the mapping
, the hamiltonian can be diagonalized as @xcite h=_k _ k ( _ k^_k - 12 ) , [ eq : h_ising ] where @xmath59 are the fermionic annihilation operators .
the ground state of @xmath60 is the vacuum state annihilated by all @xmath59 and the dispersion relation reads _
k = 2 , [ eq : epsilon_k ] where the quasi - momenta @xmath61 take discrete values @xmath62 $ ] .    even at the critical point , for @xmath58 ,
the dispersion relation is not strictly linear for all values of @xmath63 $ ] due to the discrete nature of the system , approaching linear behaviour only in the limit of long - wavelength @xmath64 . as a result
, the velocity of the quasiparticles depends on the momentum @xmath61 , especially for larger values of @xmath65  in contrast with what is the case for cft . in order to account for that , as well as for comparisons with cft
, we simultaneously consider the model with linearized dispersion relation fixing the velocity of quasiparticles , h_lin=_k ^lin_k ( _ k^_k - 12 ) , where @xmath66 and @xmath59 are the annihilation operators diagonalizing ising hamiltonian in eq .  .
in our conventions , @xmath67 is the velocity of quasiparticles at the critical point for that hamiltonian in the long - wavelength limit of @xmath68 , but we will present our plots appropriately rescaled for comparisons with cft .
next , we study the evolution of rnyi entanglement entropies in states locally excited by operators on the lattice .
more precisely , we excite the ground state of the critical ising model with local operators @xmath69 , which can have support on more then one lattice site , and perform a unitary time evolution .
we then numerically calculate the entropy of a block @xmath9 of @xmath70 consecutive spins at a distance @xmath1 from the excitation for different times and subtract from it the entropy of the block with no excitation .
this setup is schematically illustrated in fig .
[ fig : schematic ] .
the lattice operators that have @xmath17 and @xmath3 fields as their leading contributions in the continuum are @xcite ( n)&=&^x_n , + ( n)&=&^x_n ^x_n+1- ^z_n .
moreover , to confront the cft predictions for other members of a given conformal family , we also consider the simplest descendant , namely the spatial discrete derivative of the @xmath17 field d(n)=^x_n+1 - ^x_n .
that corresponds to descendant @xmath71 .
+ in the following sections , we discuss the results for classes of states locally excited by these three different operators .
the results for @xmath72 excitation are collected in fig .
[ fig : x ] . column ( a ) shows the change in entropy for fixed @xmath70 and @xmath1 .
we observe the expected plateau appearing at a time when @xmath73  and subsequent plateaus when the signal enters the block @xmath9 from the opposite side or after making some number of circle around the chain .
the plateau is however oscillating slightly and then vanishing in a long tail for @xmath74 .
the mean values in the first plateau are @xmath75 and @xmath76 , close to the expected value of @xmath77 .
both the oscillations and the tail visible for the first plateau are mostly independent of the system size .
when we use linearized ising hamiltonian @xmath78 ( dashed lines ) in place of ising hamiltonian @xmath60 ( solid lines ) both the oscillations and the tails disappear and fully periodic structure consistent with the cft prediction is recovered .
this supports the natural interpretation that the tail is related to different velocities of excited quasiparticles which naturally appear for local hamiltonian on a chain and in this case are smaller then @xmath79 .
notice also that the dashed line has smooth edges , what can be interpreted as manifestation of the nonzero @xmath3 on the lattice , compare with fig .
[ fig : cr1 ] . finally , for @xmath60 , the subsequent plateaus appearing for longer times are visibly shifted up , as the signal is slowly dissipating .
column ( b ) shows the change of entropy when both @xmath1 and @xmath70 are proportional to the system size @xmath26 and the time is properly rescaled .
this validates that the obtained value of @xmath80 is independent of the block size  in contrast to the subtracted background entropy of the block without excitation which in this case grows logarithmically with @xmath81 .
additionally , we observe that after the rescaling the tails for different @xmath81 are collapsing  at least up to corrections which are not visible in this scale and for short enough times  which suggests that the characteristic time scale at which the tails are disappearing is proportional to the block size .
finally , in column ( c ) we place the excitation symmetrically with respect to the block , so that the quasiparticle with the same absolute momentum traveling left and right should be entering the block at the same time . even in this setup
@xmath80 acquires large non - zero value when the fastest quasiparticles traveling with @xmath79 reach block @xmath9 from both sides .
this shows that the simplest interpretation valid for global ( translationally invariant ) quench that only the pairs of quasiparticles with opposite momenta @xmath82 contribute to the entanglement is no longer valid for local excitation which breaks translational invariance .
this signal almost disappears when @xmath78 is used and all the quasiparticles enter or leave the block in the same instance of time .
we refer for a consistent discussion in the context of local quenches to @xcite and evolution of the negativity to @xcite .
the results for @xmath83 excitation are presented in fig .
[ fig : e ] . for fixed @xmath1 and @xmath70 in column ( a )
we obtain non - zero signal with sharp peaks when the signal is first entering and leaving the block , which then disappears in a long tail .
when the evolution is governed by @xmath78 the tails are not present and we recover the periodic structure of plateaus , which however acquire non - zero values of @xmath84 , and @xmath85 .
interestingly , in this case , there are still additional sharp peaks when the signal is entering or leaving the block , which are similar to the structure obtained from cft with small but finite @xmath3 .
we are however not able to explain the value of @xmath80 in the plateaus as the logarithm of the quantum dimension of @xmath14 that should be zero .
we conclude that the sub - leading contributions to @xmath83 in the continuum turn out to be more relevant than for @xmath72 . at the same time
the observed results suggest that the quasiparticles with large @xmath61 significantly contribute to the observed value of entanglement , which might be related here with the fact that @xmath83 is supported on two lattice sites .    the stability of the obtained signal is validated in column ( b ) where we set @xmath70 and @xmath1 proportional to system size @xmath26 , and after properly rescaling of the time we observe the collapse of @xmath80 for different values of @xmath26 .
finally , in ( c ) we show that the non - zero value of @xmath80 obtained in the plateaus for @xmath78 can be also recovered from the evolution governed by the original ising hamiltonian @xmath60 when the block @xmath9 is placed just next to excitation .
this way difference in velocity of excited quasiparticles turn out to be unimportant for the first plateau as ( almost ) all right - moving quasiparticles are able to enter the block and the situation resembles that for @xmath78 .
this suggests that similar strategy can be used in general spin chains , which can not be mapped onto system of free fermions which we use to construct @xmath78 .
such systems can be conveniently simulated using the toolbox of matrix product states ( mps ) @xcite , where , for instance , specific algorithms to study the dynamics of localized excitations in infinite systems have been put forward @xcite .
the results for the discrete spatial derivative of @xmath72 are collected in fig .
[ fig : dx ] .
we note that , on the contrary to the two previous cases , the operator @xmath87 is not unitary , so even at @xmath88 it changes the entropy of a block supported on sites other then @xmath7 and @xmath89 .
we however see that this effect is relatively small and vanishing with increasing @xmath1 .    in column ( a ) , for fixed @xmath90 and @xmath70 , we observe that @xmath80 does not form plateaus and then disappears in a long tails , which suggests that the slower quasiparticles with large @xmath61 significantly contribute to the observed signal .
this is further corroborated by using @xmath78 which allows to recover the periodic structure of plateaus with the values at the peaks close to @xmath91 predicted by cft .
similar observation also hold in column ( b ) for @xmath1 and @xmath70 proportional to the system size @xmath26 .
finally , in column ( c ) , similarly to the situation in the previous section , we observe that for local ising hamiltonian @xmath60 we are able to recover the structure of the first plateau with the peak value @xmath92 if the block @xmath9 is placed just next to the excitation and @xmath93 .
in this section we present numerical results that , in principle , could be reproduced from cft , but in practice the analytical computations become very difficult . in
such cases numerics is a great tool for understanding the phenomenology of entanglement evolution and we explore it below . since we successfully recovered the cft predictions for the spin @xmath17 operator , we will mostly consider states excited by applying @xmath94 to the lattice sites , but we believe that the universal features of our analysis remain valid for other conformal families .
+ we begin with evolution of entanglement of a block of spins after exciting the ground state by local operator but in the non - critical model . having the diagonalized ising hamiltonian for any value of the parameters , we can study the evolution of the rnyi entropies away from the critical point . in fig .
[ non - crit ] we plot the evolution of the second rnyi entropy after acting with @xmath72 for a few values of @xmath95 around the critical one .
clearly , in each case the entropy only changes after time of order @xmath1 , however the clear plateau with the logarithm of the quantum dimension only appears at criticality .
this , in principle , might serve as a proxy for seeing the critical point with local excitations but in general the value of the plateau might be a more complicated expression in terms of the quantum dimensions of the model .
figure  [ non - crit ] also indicates that the critical hamiltonian leads to the smallest increase of entanglement of the block after we excite a vacuum by a local operator .
based on our limited analysis , it is not obvious if this lower bound is universal and satisfied by general families of hamiltonians that have a critical point in their parameter space .
it would be very interesting to provide further checks using other available methods
like e.g. mps  in order to support this observation or possibly find counter - examples .
+     excitation for different values of the magnetic field @xmath95 .
the distance of excitation from the block is ( a ) @xmath93 and ( b ) @xmath96 .
@xmath97 , @xmath98.,width=529 ]    next , by using the map to free fermions we can also study more general excited states by acting with multiple local operators on the critical chain . on the other hand ,
the computations in cft using the replica trick become very cumbersome and were only done for two excitations in @xcite .
this part will then serve as a collection of new predictions for the evolution of the rnyi entanglement entropies in ( rational ) cfts .
we begin with states where few operators were inserted on sufficiently separated sites and in some distance to the entangling block .
we chose the block to be large enough so that there is a time where all the excited quasiparticles are inside the entangling region . from numerical results
it is clear that in such excited states the entanglement entropy is a sum of the quantum dimensions of the operators , see fig .
[ steps ] for a case of 3 and 5 @xmath17-excitations .
the time evolution of entanglement in such states is very similar to periodic quench studied , e.g. , in ref .
@xcite .     on ( a ) @xmath99 sites , @xmath100 and ( b ) @xmath101 sites , @xmath102 .
results for ising hamiltonian ( solid lines ) and linearized ising hamiltonian ( dashed lines ) .
block @xmath9 supported on sites @xmath103 ( @xmath104 ) .
, width=529 ]    based on these observations , we conjecture that in rational cfts , and more generally cfts for which the quasi - particle picture remains a good effective description of the dynamics of entanglement , for @xmath106 such excitations the change in maximal contribution to the entanglement rnyi entropies is simply s^(n)_l=^m_i=1d_i with quantum dimension of the @xmath107th operator @xmath108 ( can be different or the same ) .
we verify numerically , using @xmath78 , that if we use some combination of operators @xmath72 , @xmath109 and @xmath110 placed in the setup considered in fig .
[ steps ] where the insertion points are sufficiently separated , then the total increase of the entropy of the block @xmath9 indeed is a sum of contributions from single excitations discussed in section [ sec : num ] .
a similar result was reported in free quantum field theories in @xcite but it would be very interesting to formulate at least a necessary conditions for the validity of this formula in arbitrary interacting 2d cft and we leave this as an open future problem .
+ another class of excited states that we consider is defined by acting with local operators on all the sites of the chain | _ g > = ^n_i=1 _ i| 0 > .
they could be thought of as a version of a global quench @xcite and have recently been employed in large @xmath111 holographic cfts as states dual to the matter collapsing to a black hole @xcite .    .
block size @xmath112 and the system size @xmath105 . , width=302 ]    as we can see , the evolution looks qualitatively similar to the global quench in the finite size system @xcite , but the final value in that case is given by the entropy density times the length of the interval . here , the evolution of entanglement entropies can still be interpreted in terms of the quasiparticles propagating from the lattice sites to the left and right .
this picture suggests that the maximal value of the entropies is equal to @xmath113 , where the factor of two comes from the fact that at each site we can have two ( left and right ) quasiparticles . on the other hand ,
the maximal value of the entropy of the block of @xmath70 spins is attained by the density matrix with all @xmath114 eigenvalues equal and is also equal to @xmath115 .
apparently these two numbers coincide for @xmath116 and the value of the plateaux for @xmath117 is equal to this maximal value minus the ground state entropy , i.e. @xmath118 attains its maximal possible value , which can be seen in fig .
[ allsx ] for @xmath78 .
if ising hamiltonian @xmath60 is used instead , the maximal value of the plateaux is not reached and the clean periodic structure of revivals is obscured , showing the relevance of coherent dynamics of quasiparticles for obtaining those effects . in other words , the memory effects for the quasiparticles in the critical model seem to be suppressed by entanglement between quasiparticles with different momenta and for late times entanglement measures saturate .
i would be interesting to compare how this behaviour changes for local operators with different quantum dimensions and for different densities of excitations leading possibly to a saturation , so we leave it as another open problem .
moreover , according to @xcite , in the large @xmath111 cfts , entanglement entropy saturates at the thermal value with an effective temperature and the time for returning into the initial value ( poincare recurrence ) is expected to be exponential in the central charge .
it would be also interesting to further explore what happens in between these two regimes and how our memory effects are corrected once we consider states in chaotic or non - local toy models for black holes as for instance in ref .
[ sec : rel_ent ] we finish our numerical explorations with evolution of the relative entropy at criticality . as far as we are aware , numerical analysis of the relative entropy in critical systems has been much less explored than the rnyi entropies . nevertheless , it is an important tool for understanding the notion of the distance between quantum states in field theories and plays an interesting role in uncovering the features of holographic cfts @xcite .
the relative entropy is defined for two reduced density matrices @xmath119 and @xmath120 as s(|)=()- ( ) .
[ eq : rel_ent ] in general excited state of a 2d cft , the second term makes it very hard to compute analytically . even for locally excited states the replica method requires the knowledge of the correlation function of @xmath121 operators in order to continue to @xmath122 ,
see , e.g. , @xcite . on the other hand , if we compare excited states @xmath119 with @xmath120 obtained in the vacuum , given that the reduced density matrix @xmath120 can be written as the exponent of a known modular hamiltonian , @xmath123 , we can express the entropy as s(|)=h_m - s^(1 ) .
[ rre ] the expectation value of the vacuum modular hamiltonian @xmath124 is computed in the excited state and @xmath125 denotes the difference of the von - neumann entropies of the two density matrices , see @xcite for more details .
now , in 2d cfts , the expectation value of the vacuum modular hamiltonian for an interval of length @xmath70 in a state locally excited by a primary operator is universal .
namely , it follows from the ope of the stress tensor with the primary operator , and can be computed from ( see e.g. @xcite ) h_m=_ldx ( l - x)xt_00(x ) .
moreover , as we argued in section [ sec : cft ] , the change of the entanglement entropy of the block is equal to the logarithm of the quantum dimension of the primary operator . from these two results we evaluate the relative entropy in our cft setup and we can possibly compare it with numerics .    .
colors indicate different block sizes @xmath126 ( blue ) , @xmath127 ( red ) , @xmath128 ( green ) .
results for ising hamiltonian ( solid lines ) and linearized ising hamiltonian ( dashed lines ) .
distance of the entangling block from the excitation ( a ) @xmath93 , ( b ) @xmath96 .
@xmath105 . inset ( c ) shows the maximum of relative entropy from panel ( b ) for linearized ising hamiltonian , where the linear fit @xmath129 .
, width=529 ]    in the following , we simply take @xmath120 to be the reduced density matrix of a block of @xmath70 spins in the ground state and @xmath119 as the density matrix of the block for the state locally excited by operator @xmath72 in a distance @xmath1 from the block .
we find an elegant expression for the relative entropy in terms of fermionic covariance matrices for gaussian @xmath119 and @xmath120 @xmath130 and we refer to the appendix [ numer ] for details .
we only notice here that as @xmath131 and @xmath132 , or equivalently @xmath119 and @xmath120 , can not be simultaneously diagonalized , finite numerical precision limits the calculation of the second term in the above expression only to relatively small block sizes .
we present the results in fig .
[ fig : rel_entropy ] .
we observe that the signal strongly changes depending if @xmath60 or @xmath78 is governing the time evolution . for the ising hamiltonian @xmath60 the maximum of the relative entropy for given
block size @xmath70 quickly disappears with the distance between the block and the excitation  suggesting that modes with high momenta contribute significantly to the observed value .
if linearized ising hamiltonian @xmath78 is used instead , the signal does not dissipate , with both the width and the value at the peak being proportional to @xmath70 .
one can verify that , for non - zero @xmath3 , the numerical evolution of the relative entropy is consistent with the cft computation with the expectation value of the stress tensor in our locally excited state .
moreover , it is possible to compute in cft the value of the peak of the relative entropy .
the maximum of the relative entropy is universal in the small @xmath3 limit ( s(| ) ) e_l+o(),[fl ] where the energy due to the operator insertion is @xmath133 .
numerically , this is shown on the inset in fig .
[ fig : rel_entropy ] where we see the linear growth of the peak of the relative entropy with the length of the interval .
our numerical results in the ising chain are consistent with the energy @xmath134 close to the slope @xmath135 fitted in fig .
[ fig : rel_entropy](c ) .
let us finally compare our formula with the first - law for entanglement entropy for a family of the nearby equilibrium states @xcite that reads @xmath136 with @xmath137 .
this relation is a consequence of the vanishing relative entropy .
interestingly , we observe an analogous first - law like relation for the maximal value of the distance between our two quantum states as measured by the relative entropy . moreover , the maximal value is universal and contained in @xmath138 .
the details of the cft analysis for this class of locally excited states will appear elsewhere @xcite .
we have shown that , by using local operators on the critical lattice , one is able to extract further non - trivial cft data  like quantum dimensions  numerically .
we have successfully done so for the spin operator @xmath17 and its derivative descendant .
however , we also saw that , even after linearization of the dispersion relation , the lattice energy operator gave a non - zero contribution to the entropy of a block .
this fact might be explained by the ambiguity in the identification between the lattice and the cft operators or a non - trivial contribution from the  tail " of operators in the continuum and deserves further investigation .
our analysis could be naturally extended to other cfts with operators that enjoy known lattice counterparts .
a natural setup might be the three - state potts model in which @xcite recently analyzed the lattice realization of the local operators in the parafermionic cft .
clearly , it would be much more interesting to explore analytically what is the contribution from a general lattice operator and what kind of information can be extracted from the increase in rnyi entropies .
in fact there has been a lot of work on extracting local cft operators from the multiscale entanglement renormalization ansatz ( mera ) @xcite .
it would be interesting to apply those developments in the studies of entanglement evolution after local operator excitations .
the evolution of entanglement measures after local excitations away from the critical point appears to be relatively unexplored .
as we saw , the critical behavior appears to be very special  characterized by formation of a clear plateaux  and might be used as a smoking gun of a critical point .
moreover , the contribution to the rnyi entropies appears to be the smallest for the ising hamiltonian with critical parameters what might be a sign of a general bound .
we hope that our analysis will serve as a starting point for more complicated systems and will help to uncover other unknown universal phenomena in the propagation of entanglement .
last but not least , we performed the evolution of the relative entropy between the vacuum and a locally excited state .
interestingly , this distance measure shows universal features analogous to the first law of entanglement and its maximum  maximal quantum distance  is proportional to the change in the energy with an effective temperature . exploring this relation in cfts or free field theories where explicit computations are under control opens a new interesting path for investigation . +
* acknowledgements * we would like to thank john cardy , paul fendley , masahiro nozaki , tokiro numasawa , tadashi takayanagi , kento watanabe , luca tagliacozzo and alvaro veliz - osorio for correspondence and discussions .
we would also like to thank xueda wen for comments sharing some of his unpublished results .
we acknowledge support by the swedish research council ( vr ) grant 2013 - 4329 ( p.c . ) and narodowe centrum nauki ( ncn , national science center ) under project no .
2013/09/b / st3/01603 ( m.m.r ) .
p.c . would like to thank yukawa institute for theoretical physics for hospitality and support during the `` quantum information in string theory and many - body systems '' workshop where some of this work was performed .
we would like to thank the organizers of the tensor network summer school in ghent where this work was initiated .
_ hamiltonian. _ the ising hamiltonian is diagonalized in a standard way @xcite by mapping it onto a system of free fermions using the jordan - wigner transformation @xmath139 where @xmath140 are fermionic annihilation operators . for convenience ,
we introduce majorana fermions @xmath141 , @xmath142 which are hermitian , unitary and satisfy canonical anticommutation relations
@xmath143_+ =   2 \delta_{m , n}$ ] . the ising hamiltonian with periodic boundary conditions can be rewritten as @xmath144 above , @xmath145 are projectors on the subspaces with respectively even and odd number of fermions , where the parity operator @xmath146 commutes with @xmath147 .
the hamiltonians in both subspaces can be expressed in terms of fermionic operators as @xmath148 differing only at the boundary term .
this is accounted for by enforcing the boundary conditions : antiperiodic @xmath149 for @xmath150 , and periodic @xmath151 for @xmath152 . for @xmath153 , which we use in this work
, the ground state of @xmath147 belongs to the subspace with even parity for any value of @xmath26 , see e.g. ref .
@xcite .
it is convenient to introduce @xmath154 , where @xmath155 is a column vector composed of operators @xmath156 and @xmath157 are @xmath158 hermitian matrices . in the following ,
we employ the matrix notation whenever is serves to simplify the notation . in each parity subspace the system is solved independently by a canonical transformation to a new base of majorana fermions , @xmath159 , where @xmath160 are real and orthogonal matrices .
the hamiltonian then reads @xmath161 expressed above in terms of standard annihilation operators @xmath162 .
this can be done analytically by a subsequent fourier and bogoliubov transformations @xcite , resulting in @xmath163 given by eq .  .
the fermionic modes are indexed by quasi - momenta consistent with the respective boundary conditions , i.e.   @xmath164 $ ] for @xmath165 , and @xmath166 $ ] for @xmath167 .
this procedure is equivalent to bringing @xmath168 into the canonical form @xmath169 .
we can now formally introduce the linearized ising hamiltonian by @xmath170 defined as @xmath171 , where the linearized dispersion relation @xmath172 .
the ground state of the system , @xmath173 , which is the starting point for our further considerations , is the even parity state annihilated by all annihilation operators diagonalizing @xmath174 , namely @xmath175 .
all the information about this state is encoded in the matrix @xmath176 .
_ local operators and time evolution. _ we consider local operators which in the fermionic language read @xmath177 up to the irrelevant phase factors .
we simulate the action of those operators on the ground state by employing the heisenberg picture . as all operators in eq .
are expressed as a product of linear combinations of operators @xmath156 , this can be done iteratively .
starting with the simplest case of operator @xmath178 , which is unitary by construction , we define the operator in the heisenberg picture as ^o_1_n = _ 1^_n _ 1
= - ( -1)^_nm_1 _ n = - _ j ( -1)^_nm_1 [ u^+]_nj ^+_j , where a new matrix @xmath179_{nj } =    - ( -1)^{\delta_{nm_1 } } [ { \rm u}^+]_{nj}$ ] describes the transformation to the fermionic base in which the initial state is the vacuum  note that operator @xmath180 is applied to the ground state of the ising hamiltonian @xmath181 in eq .  .
subsequent application of unitary operators @xmath182 for @xmath183 is obtained similarly as @xmath184 , and @xmath185_{nj } =    - ( -1)^{\delta_{nm_l } } [ { \rm u}^{o_{l-1}}]_{nj}$ ] for @xmath186 , corresponding to the action of the operator @xmath187 on the initial ground state .    in a general case
when @xmath188 , with @xmath189 being a row vector of coefficients , we have to generalize the above procedure .
notice that @xmath180 , in principle , does not have to be unitary .
as we are only interested in the action of @xmath180 on the ground state , @xmath173 , we introduce unitary operator @xmath190 , such that @xmath191 up to normalization . for @xmath190 to be unitary ,
it is enough to require that @xmath192 , where @xmath193 is a row vector of real coefficients .
we can now define the operators in the heisenberg picture ^o_1 = _ 1^ _ 1 = u^+ ( _ 1^^+ _ 1 ) = u^+ ( 2 - ) ^+ , defining new transformation matrix @xmath194 .
the vector @xmath193 is found as @xmath195 _ { 2n-1 } - i [ \vec v    \mathrm{u}^+ ] _ { 2n})$ ] and @xmath196 _
{ 2n-1 } - i [ \vec v    \mathrm{u}^+ ] _ { 2n})$ ] .
this follows from the condition @xmath197 , where the two operators may differ by operators @xmath198 annihilating the vacuum state @xmath173 .
finally the time evolution is simulated in the heisenberg picture as @xmath199 $ ] , which leads to @xmath200 , where @xmath201 , and the matrix @xmath202 for the even ( odd ) parity of the excited state @xmath203 .
notice that operator @xmath109 is preserving the parity , and @xmath72 , @xmath110 are changing it .
finally , the time evolution with the linearized hamiltonian is obtained by using matrices @xmath204 instead of @xmath205 above .    _ calculating entropies. _ the entropy of a block of consecutive spins is calculated in a standard way @xcite , as all the information about the reduced density matrix of a block is encoded in the @xmath206 covariance matrix @xmath207_{m , n = 2l_0 + 1,2l_0 + 2 , \dotsc,2l_0 + 2l}\ ] ] supported on @xmath70 consecutive lattice sites @xmath208 .
it is found as a submatrix of @xmath209 , with @xmath210 being the correlation matrix in the canonical ( vacuum ) base .
rnyi entropy of the block is then found as , s^()_l = _ j=1^l , where @xmath211 are the eigenvalues of @xmath212 . for @xmath213
, the von - neumann entropy reads s^(1)_l = - _ j=1^l = - ( c_l c_l ) , where the second equation provides the direct expression in term of covariance matrix @xmath212 .
relative entropy can be similarly found similarly as @xmath214 as the matrices @xmath215 and @xmath132 ( calculated as the correlation matrix in the ground state ) can not be diagonalized simultaneously , it becomes difficult to numerically calculate the second term in the above equation for large blocks .
in that case @xmath132 has many eigenvalues approaching @xmath216 , falling below numerical precision , which prevents the calculation of the logarithm .
for this reason , in sec .
[ sec : rel_ent ] , we calculate residual entropy only for small enough block where this problem does not yet occur .    finally , in order to derive the non - standard second term on the r.h.s . of eq .
we employ the fact that @xmath120 is gaussian and can be expressed as @xmath217 , where the diagonal form in terms of fermionic annihilation ( creation ) operators @xmath218 ( @xmath219 ) can be directly obtained from the correlation matrix @xmath220 @xcite .
we introduce here @xmath221 as a column vector consisting of majorana fermions @xmath222 in the entangling block .
then @xmath223 , where @xmath224 is a column vector consisting of @xmath225 , and @xmath226 is a unitary matrix diagonalising @xmath227 \mathrm{u}^{\vartheta \dagger}$ ] .
this allows to compute @xmath228 .
using matrix notation this equals @xmath229 \left [ \bigoplus_{j = 1}^l \begin{pmatrix } \log q_j & 0 \\   0   & \log ( 1-q_j ) \\
\end{pmatrix } \right ]   \right )   =   \tr \left ( \mathrm{u}^{\vartheta \dagger } c_l^\rho \mathrm{u}^{\vartheta }   \left [ \bigoplus_{j = 1}^l   \begin{pmatrix } \log q_j & 0 \\ 0   & \log ( 1-q_j ) \\
\end{pmatrix } \right ] \right )   =    \tr \left ( c^\rho_l \log c^\vartheta_l \right)$ ] , which gives eq .  .
calabrese p and cardy j 2016 j. stat .
mech . 064003 nozaki m , numasawa t and takayanagi t 2014 phys .
112 * 111602 nozaki m 2014 jhep * 1410 * 147 caputa p , nozaki m and takayanagi t 2014 ptep * 2014 * 093b06 + caputa p , simon j , stikonas a and takayanagi t 2015 jhep * 1501 * 102 + asplund c t , bernamonti a , galli f and hartman t 2015 jhep * 1502 * 171 + caputa p , simon j , stikonas a , takayanagi t and watanabe k 2015 jhep * 1508 * 011 + shiba n 2014 jhep * 1412 * 152 + guo w z and he s 2015 jhep * 1504 * 099 + guo w z 2015 arxiv:1510.07142 + nozaki m , numasawa t and matsuura s 2016 jhep * 1602 * 150 + caputa p , numasawa t and veliz - osorio a 2016 arxiv:1602.06542 + sivaramakrishnan a 2016 arxiv:1604.00965 + david j r , khetrapal s and kumar s p 2016 arxiv:1605.05987 + caputa p , nozaki m and numasawa t 2016 phys . rev .
d * 93 * 105032 + nozaki m , and watamura n , arxiv:1606.07076 + zhou t , arxiv:1607.08631 alcaraz f c , berganza m i and sierra g 2011 phys .
lett . * 106 * , 201601 + berganza m i , alcaraz f c and sierra g 2012 j. stat .
* 1201 * p01016 + v.  alba , m.  fagotti and p.  calabrese , 2009 j. stat . mech . * 0910 * , p10020 + palmai t 2014 phys .
b * 90 * 161404 + sheikh - jabbari m m and yavartanoo h 2016 arxiv:1605.00341 + chen b and wu j q 2016 arxiv:1605.06753 + taddia
l , ortolani f and palmai t 2016 arxiv:1606.02667 + giusto s , and russo r , 2014 phys .
d * 90 * , no .
6 , 066004    he s , numasawa t , takayanagi t and watanabe k 2014 phys .
d * 90 * 041701 caputa p and veliz - osorio a 2015 phys .
d * 92 * 065010 chen b , guo w z , he s and wu j q 2015 jhep * 1510 * 173 mong r s k , clarke d j , alicea j , lindner n h and fendley p 2014 j. phys .
a * 47 * 452001 lashkari n 2016 phys .
* 117 * 041601          zauner v , ganahl m , evertz h g , nishino t 2015 j. phys . : condens .
matter * 27 * 425602 + phien h
n , vidal g and mcculloch i p 2012 phys .
* 245107 + milsted a , haegeman j , osborne t j and verstraete f 2013 phys .
b * 88 * 155116        anous t , hartman t , rovai a and sonner j 2016 jhep * 07 * 123 mezei m and stanford d 2016 arxiv:1608.05101 magan j m 2016 arxiv:1601.04663 + jansen a and magan j m 2016 arxiv:1604.03772 n.  lashkari , j.  lin , h.  ooguri , b.  stoica and m.  van raamsdonk , arxiv:1605.01075 sarosi g and ugajin t 2016 jhep * 07 * 114 blanco d d , casini h , hung l y and myers r c 2013 jhep * 1308 * 060 cardy j and tonni e arxiv:1608.01283 caputa p _ in preperation _    bhattacharya j , nozaki m , takayanagi t and ugajin t 2013 phys .
lett . * 110 * 091602 giovannetti v , montangero s and fazio r 2008 phys .
* 101 * 180503 + pfeifer r
n c , evenbly g and vidal g 2009 phys .
a * 79 * 040301 + evenbly g , corboz p and vidal g 2010 phys .
b * 82 * 132411 + evenbly g and vidal , g 2013 quantum criticality with the multi - scale entanglement renormalization ansatz . in : avella ,
a. , mancini , f. ( eds . ) strongly correlated systems : numerical methods .
springer , heidelberg , arxiv:1109.5334        peschel i 2003 j. phys . a : math.gen . * 36 * l205 + vidal g , latorre j i , rico e and kitaev a 2003 phys .
rev . lett . *
90 * 227902 + jin b - q and korepin v e 2004 j. stat . phys .
* 116 * 79 + peschel i and eisler v 2009 j. phys .
a : math . theor .
* 42 * 504003


coupled chaotic systems are recently calling much attention due to the verification that they may be useful to the understanding of natural systems in a variety of fields as in ecology @xcite , in neuroscience @xcite , in economy @xcite , and in lasers @xcite .
it has been verified that despite of the higher dimensionality of a coupled chaotic system , the coupling among the elements might make them to synchronize @xcite , reducing the dynamics of the system to a few degrees of freedom .    in this work ,
we focus our attention in the phenomenon of phase synchronization ( ps ) , which describes the appearance of a phase synchronous behavior between two interacting chaotic systems @xcite , i.e. , given two chaotic systems , their phase difference remains bounded , despite of the fact that their amplitudes may be uncorrelated .
this phenomenon is particularly interesting since it can arise from a very small coupling strength .
its presence was reported in a variety of experimental systems .
it was experimentally demonstrated in electronic circuits @xcite , and latter in electrochemical oscillators @xcite .
it was found in plasma @xcite , in the chua s circuit @xcite , and there were also found evidences of phase synchronization in communication processes in the human brain @xcite .    to detect ps in a real - time experiment , one has to measure the phase of the chaotic trajectory @xcite .
however , the phase is not always an easily accessible information . to overcome this difficulty , it is important to understand fundamental properties of phase synchronous systems , that could be experimentally easily verified . for chaotic systems that are phase synchronized with a periodic forcing @xcite , it was reported that a stroboscopic map of the trajectory was a subset of this attractor and occupies only partially the region delimitated by a projection of the attractor .
this property was used to detect in a real - time experiment phase synchronization between the chua s circuit and a sinusoidal forcing @xcite .
this approach of detecting phase synchronization through the stroboscopic map can be extended for coupled chaotic oscillators , in a formal way .
the stroboscopic map is generalized to the * conditional poincar map*. given two oscillators , at least one being chaotic , the conditional poincar map is constructed by collecting points in one oscillator at the moment at which some event happens in the other one . if the set of discrete points generated by this conditional map does not visit any arbitrary region of a especial projection of the chaotic attractor , we call this set a * p - set*. this property of the conditional poincar map is called non - transitivity @xcite , i.e. , an initial condition under the conditional poincar map does not visit everywhere in a subspace of the attractor . alike
the stroboscopic maps of oscillators that are in phase synchrony with a forcing , the conditional poincar maps of coupled chaotic oscillators , in ps , also only partially occupy a projection of the attractor .    in this work ,
we show how the conditional poincar map can be used to detect ps , without actually having to measure the phase . for phase - coherent oscillators , a special type of p - set , that we call ps - set ( phase synchronization set ) ,
conversely , its existence also implies ps .
we illustrate our findings and ideas with numerical and experimental analyzes in the forced chua s circuit , and the coupled rsller oscillator @xcite .
further , we extend these results to non - phase coherent attractors .
finally , we also introduce a phase of a chaotic trajectory to be a quantity related to the amount of rotation of the tangent vector .
this definition can be used to chaotic attractors , independently whether they have phase - coherent or non phase - coherent dynamics .
this work is organized as follows : in sec .
[ sec : phase ] , we define a way to measure the phase of a chaotic flow , and discuss the relation between the average returning time and the average angular frequency . in sec .
[ sec : ps ] , we discuss the conditions for ps states and , in sec . [
sec : chua ] , we describe the phenomenon of ps in the forced chua s circuit .
we introduce the notion of a conditional poincar map in sec .
[ cpm ] and the p - sets ( as well as the ps - sets ) in sec .
[ p - sets ] . in sec .
[ sec : bsc ] , we show how ps can be found by the detection of these sets in the forced chua s circuit and in sec .
[ sec : bsr ] , for the coupled rssler oscillator .
further , in sec . [ extension ] , we discuss the extension of these ideas to a class of non - coherent oscillators . in sec .
[ sec : cr ] , we make some remarks and the conclusions of this work . in appendix
a , we formally introduce the conditional poincar map and the p - set , and in appendix b , we show that for coherent dynamics the ps - sets exist if , and only if , there is ps .
in other words , ps implies ps - sets and vice - versa .
the phase of a chaotic attractor in a projection @xmath0 ( a subspace ) is defined to be the amount of rotation of the tangent vector in this projection , and can be represented by an integral function of the type    @xmath1    with @xmath2 being an infinitesimal displacement of the tangent vector of the flow , from time @xmath3 to time @xmath4 , and @xmath5 .
note that in eq .
( [ fase ] ) , we are measuring the amount of rotation , per unit time , of a projection of the tangent vector of the flow , on the same subspace @xmath0 where the phase is defined .
we call this subspace @xmath6 .
the attractor @xmath7 , projected on the subspace @xmath6 is regarded as @xmath8 .
the instantaneous angular frequency of the trajectory in @xmath8 , named @xmath9 is given by @xmath10 .
so , from eq .
( [ fase ] ) , @xmath11 and , the average angular frequency @xmath12 is @xmath13 @xmath14 represents the average . equation ( [ natural_frequence_2 ] ) can be put into the form @xmath15 = @xmath16 .
whenever a poincar section can be defined , the average period of the chaotic attractor on the subspace @xmath6 is calculated by @xmath17 where @xmath18 , and @xmath19 represents the time at which the trajectory in the subspace @xmath6 makes the @xmath20-th crossing with this poincar section .
we introduce @xmath21 to be the average displacement of the phase for a typical period as @xmath22 with @xmath23 being the phase associated to the subspace @xmath6 at the moment that the @xmath24-th crossing between the trajectory and the poincar section happens .
thus , we can put eq .
( [ natural_frequence_2 ] ) as @xmath25 for the forced chua s circuit , the subspace @xmath26 is defined by a suitable projection of the circuit variable .
we have that @xmath27 .
so , eq . ( [ natural_frequence_1 ] ) can be written as @xmath28 .
for the coupled rssler oscillator , the quantities in eq .
( [ natural_frequence_1 ] ) can be calculated in two subspaces , one associated to the variables of one rssler , the subspace @xmath26 , and the subspace @xmath29 , associated to the variables of the other rssler system . as shown in @xcite , @xmath30 might slightly differs from @xmath31 , and thus , @xmath32 .
( [ natural_frequence_2 ] ) and ( [ natural_frequence_1 ] ) relate the average period , the average angular frequency , and the phase of a chaotic trajectory .
this shows that the average period ( recurrence ) and the average angular frequency are intimately connected in phase coherent chaotic systems , and both these quantities can be calculated from the phase .
having defined phase , ps exists whenever the following condition is satisfied    @xmath33    the minimal bound for the phase difference , @xmath30 , in terms of the phase as defined by eq .
( [ ps ] ) , was theoretically estimated in ref .
@xcite . equation ( [ ps ] ) means that the phase difference between the two coupled systems is always bounded , and @xmath34 is a rational constant @xcite .
also , @xmath35 in this work , we will consider cases for @xmath34=1 .
otherwise , a simple change of variables could eliminate this constant from eq .
( [ nc ] ) . for the forced chua s circuit @xmath36 , with @xmath37 , representing the angular frequency of the forcing .
there is ps , if @xmath38=@xmath39 .
therefore , @xmath40=@xmath41 . for the coupled rssler ,
if ps exists , @xmath42 = @xmath43 , @xmath44=@xmath45 , and @xmath46=@xmath47 , and therefore , we could have in the right term of eq .
( [ ps ] ) @xmath48 , instead of @xmath49 .
the circuit is represented by :    @xmath50    where @xmath51 , @xmath52 , and @xmath53 represent , respectively , the tension across two capacitors and the current through the inductor ( see @xcite for more details ) , @xmath39 and @xmath54 are the angular frequency and the amplitude of the forcing , respectively
. the piecewise - linear function , @xmath55 , is given by :    @xmath56\ ] ]    where we have chosen the parameters @xmath57=0.1 , @xmath58=0.574 , @xmath59=1 , @xmath60=1/6 , @xmath61=-0.5 , and @xmath62=-0.8 , such that we obtain a rssler - type attractor , for @xmath54=0 .    to calculate the phase of the chaotic trajectory
, we first define the subspace @xmath26 to be given by the pair of variables @xmath63 , and then , we use eq .
( [ fase ] ) . in fig .
[ ps_fig6 ] , we show the difference between the phase of the chaotic circuit [ as calculated by eq .
( [ fase ] ) ] and the phase of the forcing , @xmath64 . in ( a ) , the phase difference is bounded and the average period of the chaotic attractor , defined as the average recurrence time of trajectories that cross the section @xmath65 , is equal to @xmath66=3.57015 , which is equal to 1/f , since @xmath67=0.2801 .
the average angular frequency can be calculated using eq .
( [ natural_frequence_1 ] ) , which gives us @xmath68=1.75992 . or , from eq .
( [ natural_frequence_2 ] ) , which gives us @xmath38=1.75992 .
note that the average growing of the phase [ calculated by eq .
( [ fase ] ) ] for a typical average period is @xmath69 which is @xmath31 . in fig .
[ ps_fig6](b ) , we have that @xmath70=3.57006 which is different from @xmath41 , since @xmath67=0.279 .
so , in ( b ) there is no ps , and consequently inequality ( [ ps ] ) is not satisfied .    in ref .
@xcite , we detected @xmath71 experimentally in the forced chua s circuit , whenever stroboscopic maps could be constructed for a time interval equal to @xmath72 , such that this map , projected into the same subspace considered to calculate the phase , does not occupy the region occupied by the attractor projected on the same subspace .    to understand this technique ,
we assume that @xmath73 , the phase of the chaotic trajectory , is the angle ( on the lift ) described by the vector position of this trajectory ( rssler - like attractor ) , and @xmath74=@xmath75 the phase of the forcing . if eq .
( [ ps ] ) is satisfied at any time , then it is satisfied at multiples of the period of the forcing @xmath76 .
so , we get @xmath77 , which means that a stroboscopic map has to be concentrated in an angular section smaller than @xmath31 . the stroboscopic map , which is already a subset of the chaotic flow , projected into the same subspace considered to calculate the phase , does not occupy the whole region occupied by the attractor projected in this same subspace .
another property of the stroboscopic map is that points in it are mapped into it by looking at the trajectory after a time interval given by @xmath78 , so , it is a subset that is recurrent to itself .    using this technique , and for the same parameters as @xcite , we show in fig .
[ ps_fig1](a ) the experimental synchronization region for the forced chua s circuit in the parameter space @xmath79 .
the triangular shaped region represents parameters for which the stroboscopic map has the points concentrated in an angular section smaller than 2@xmath80 .
the bump at the bottom right side of the ps region is due to non - synchronous states that present a long bounded phase difference . by a typical time interval within which the experiment is realized , which was of the order of 40,000 cycles , the system seemed to be phase synchronized , i.e.
, localized stroboscopic maps were found .
in fact , we have detected these maps even for observation times corresponding to 150,000 cycles , in the region of the bump .    in a short ,
the bump region is an extended structure in the circuit parameter space , that presents intermittent behavior in the phase difference @xcite , but with a long laminar regime , even for parameters far away from the border between the ps and the non ps region .
this intermittency differs from the usual one , observed in the transition to ps , in the fact that this latter happens very close to the border between the ps and the non ps region .
the reason for this intermittency is due to the presence of a periodic window , close to the region of the bump .
simulation is shown in fig .
[ ps_fig1](b ) , where black points represent perturbing parameters for which eq .
( [ ps ] ) is satisfied , with @xmath81 defined in eq .
( [ fase ] ) .
one sees that the @xmath71 region resembles a triangle .
the triangular shaped region , denoted by the light gray dashed line , represents the region where the system is not phase synchronized , but the phase difference remains bounded for a long time interval , that might be longer than 100,000 cycles of the systems .
so , we reproduce numerically the same atypical intermittency observed experimentally , i.e. , long laminar regime in the phase difference , for parameter regions away from the border between ps and non ps states .
this happens associated with a periodic window , as the one shown in fig .
[ ps_fig1](b ) .
the shape of the synchronization region in fig .
[ ps_fig1](b ) is equivalent to the region in the experiment , constructed by detecting the stroboscopic maps contained within a small angular section .
this proposes an equivalence between the existence of a recurrent subset and the verification of eq .
( [ ps ] ) . inside the synchronization region
, a stroboscopic map appears like in fig .
[ ps_fig7](a ) , where the light gray points represent the attractor , and the dark filled circles , the map . outside of the @xmath71 region , there are parameter sets for which the stroboscopic maps do not occupy the region occupied by the attractor , at the projection in which the phase is calculated .
as one sees in fig .
[ ps_fig7](b ) [ for the parameters represented by the plus symbol in fig .
[ ps_fig1](b ) ] there is a region of the projected attractor ( pointed by the arrow ) for which the stroboscopic map never visits .
the difference between the stroboscopic map that appears while there is ps in ( a ) , and the stroboscopic map that appears while there is not ps in ( b ) , is that points of the stroboscopic map in ( a ) are all concentrated in an angular section of the attractor .
as it will be further classified , the stroboscopic map in ( a ) is a ps - set and in ( b ) is a p - set .    as a way to better characterize the ps phenomena in the chua s circuit
, we calculate the lyapunov exponents .
as usually expected , the transition to ps is associated to one of the exponents becoming smaller than zero .
since there is already one exponent that is smaller than zero , ps induces the creation of a second stable direction in the chua s circuit .
previous to the transition to phase synchronization , this exponent was zero . in fig .
[ ps_fig8 ] , we show the exponent ( and the error bars ) associated to ps in black and the quantity @xmath82 ( in gray , for a fixed amplitude of @xmath54=0.0015 ) with respect to the frequency . in the region that the exponent becomes smaller than zero , @xmath38=@xmath39 , satisfying eq .
( [ nc ] ) .
the finding of maps of the attractor that appear as localized structures implies ps .
the conditional poincar map introduced in this chapter as a generalization of the stroboscopic map , is an efficient way of revealing the existence of such special mappings .
the stroboscopic map technique defined for periodically driven chaotic systems was explained in the previous sections .
to extent the idea of stroboscopic map in coupled chaotic oscillators we came up with the conditional poincar map , which is a map of the flow , constructed by observing it for specific times at which events occur in the subsystems @xmath8 .
an event is considered to happen when the trajectory of one subsystem crosses a poincar section .
so , given two oscillators @xmath83 and @xmath84 ( at least one being chaotic ) , with trajectories in the subspaces where the phase is defined , the conditional poincar map of @xmath83 is the trajectory position at the moment that a series of equal events happens , in @xmath84 .
analogously , the conditional poincar map of @xmath84 is the trajectory position at the moment that a series of equal events happens , in @xmath83 . in the case of
periodically forced chaotic systems , an event may be defined to happen whenever the forcing reaches a specific value , and the conditional poincar map is the usual stroboscopic map , because the time interval between two events is constant . in coupled chaotic systems ,
the time interval between two successive events is no longer constant .
we define a time series of events @xmath85 , by the following rule :    * @xmath86 represents the time at which the @xmath20-th crossing of the trajectory of @xmath84 occurs in a poincar plane .
* @xmath87 represents the time at which the @xmath20-th crossing of the trajectory of @xmath83 occurs in a poincar plane .
the discrete set of points observed at the times @xmath85 is called set @xmath88 .
this set , projected at the subspaces @xmath6 ( where the phase is calculated ) , is named @xmath89 .
the conditional poincar map is represented by @xmath90 .
so , we say that @xmath89 is the set of points generated by @xmath90 .
the next step is to define when @xmath89 can be regarded as either a p - set or a ps - set , this last set implying phase synchronization .
the @xmath89 set is a p - set , if it does not completely fulfill the projection @xmath8 of the attractor . in other words ,
a discrete set @xmath89 is considered to be a p - set , if for balls of radius @xmath91 , centered in all points of the attractor projection @xmath8 , one does not find points of @xmath89 inside of all these balls .
if @xmath89 completely fulfill @xmath8 , we say that these two sets are equivalent ( and we represent this by the symbol @xmath92 ) .
more details , see appendix a.    so , a p - set exists , if the conditional poincar map is not transitive on @xmath8 @xcite .
that is , the flow , observed by the times for which the conditional poincar map is defined , does not visit arbitrary regions of @xmath8 .
note that however , the attractor is chaotic , and therefore , the chaotic set is always transitive through the flow .
so , given a set of initial conditions , its evolution through the flow eventually reaches arbitrary open subsets of the original chaotic attractor .
we can classify three relevant types of sets generated by the conditional poincar map    type - a : :    @xmath89 is equivalent to    @xmath8    ( @xmath93 ) .
the conditional    poincar map @xmath90 is transitive in    @xmath8 .
p - set : :    @xmath89 is not equivalent to    @xmath8    ( @xmath94 ) . the    conditional poincar map is not transitive in    @xmath8 .
ps - set : :    @xmath89 is p - set , with the additional condition    that it is localized in the vicinity of the poincar section chosen to    define the events .    in the following , we comment each case :      _
these sets appear whenever there is no ps_. if two non - identical coupled chaotic systems ( topologically similar ) are not phase synchronized , the chaotic trajectories do not make correlated events in both subspaces @xmath26 and @xmath29 . as a consequence while the trajectory is positioned at the specified poincar plane , at the subspace @xmath26 , the trajectory in the subspace @xmath29 is everywhere in this subspace , making the set @xmath89 to be equivalent to @xmath8 .
an interesting illustration is the case of two uncoupled equal chaotic systems , but with different initial conditions .
as we construct the conditional poincar maps , they will be a type - a set , since also the distance between the trajectories in the two oscillators are sensitive to the initial conditions , and will diverge exponentially .      _
these sets constitutes a necessary , but not sufficient , condition to describe ps . _    in periodically forced chaotic systems
, they might appear when there is not ps .
as we already mentioned , the points in these sets are not localized in special spots of the attractor projection . as a consequence
, the domain of the absolute difference between the time at which the same number of events happen in both oscillators has a broad character .
_ ps - sets imply ps and vice - versa_. they exist , if and only if , there is phase synchronization , as shown in appendix b , and illustrated with the examples in sec .
[ sec : chua ] , and throughout this work .
another important point is that the ps - set provides a real - time detection that can be easily experimentally implemented ( sec .
[ sec : chua ] ) , and easily constructed from a data set .
ps - set implies ps because the difference between the time at which the @xmath24-th event happens in both oscillators is small , which means that the time difference @xmath95 is smaller than a finite constant value . as a consequence
, the points in the conditional poincar map of one oscillator are confined around the poincar section chosen to define the events .
therefore , the detection of a ps - set can be done by observing this characteristic of the conditional poincar map .    for rssler - like oscillators ,
in which the trajectory spirals around an equilibrium point , the ps - set is confined within an angular section .
having found the ps - set , we can study properties of these sets that give us the level of organization and coherence of the oscillators    a ps - set , @xmath89 , is said to have length 1 , if the set is constructed by the time series of the same events . defining the event to be given by the crossing of the trajectory to a poincar plane , the corresponding time series of events @xmath96 is given by @xmath97 .
for this ps - set , points in @xmath89 are mapped in @xmath89 , after one application of the conditional poincar map .
a ps - set is said to have length 2 , if it is obtained by a time series of two different events .
so , a ps - set can be constructed from more complex series of events .
we construct a length-2 basic set using a time series of events @xmath98 given by @xmath99 @xcite . as an example , for the perturbed chua s circuit , the length-2 basic set is constructed by a stroboscopic map that collects points every half period of the forcing .
a length-2 basic set is assumed to be composed by two other subsets , named minimal sets , the subsets @xmath100 and @xmath101 .
they have the property that if a point @xmath102 is such that @xmath103 , this point , iterated by the conditional poincar map , goes to the minimal set @xmath101 , and if @xmath102 is such that @xmath104 , this point , iterated by the conditional poincar map , goes to the minimal set @xmath101 .
thus , points in @xmath101 are mapped to itself after 2 applications of the conditional poincar map .
the minimal set @xmath101 is said to be disjoint to the set @xmath105 , if they do not intersect , i.e. , @xmath106 .
for some systems that present a strong phase - coherent state , as the ones here studied , in which the instantaneous trajectory velocity does not differ too much from the average velocity on typical orbits , it is possible to find length-2 basic set , with disjoint minimal sets , when ps is present .    for a general case , we do not expect to find a length-2 ps - set with disjoint minimal sets .
as an example , one can think of a spiking - firing oscillator , phase synchronized with a periodic forcing . due to the fact that the spiking - firing dynamics has a fast and a slow mode , the conditional poincar maps might overlap .      here , we explain through a diagram the possible emerging sets from the conditional poincar map .
@xmath107^{\mathcal{p}_j } \ar[ld]_{\mathcal{p}_j \circ t } &   \\
\mathcal{d}_j \ar[rd ] & & \mathcal{x}_j \ar[ld]\\   &   \mathcal{d}_j \stackrel{?}{\equiv } \mathcal{x}_j \ar[ld]_{\txt{yes } } \ar[rd]^{\txt{no } } & \\ type - a &   & p - set \ar[d ] \\ &    ps \ar@{<=>}[r ]   &   ps - set } \ ] ]    starting from the chaotic attractor @xmath7 , the set @xmath88 is constructed from the conditional poincar map , represented by @xmath108 , we project @xmath88 and @xmath7 into the subspace @xmath6 , obtaining the sets @xmath89 and @xmath8 , respectively .
we classify the set @xmath89 into type - a set or p - set , by checking whether the conditional poincar map is transitive in @xmath8 , i.e. , by verifying the equivalence between the sets @xmath8 and @xmath89 .
then , if the p - set is localized in the vicinity of the poincar section where the events occur , the p - set is a ps - set , which means that ps is present .
for applying our formalism to the periodically forced chua s circuit , the event times are @xmath109 , with @xmath110 representing the forcing period .
the time series for the length-2 basic set is given by @xmath111 .
everywhere inside the @xmath71 region , we find length-1 [ as an example , see fig . [ ps_fig7](a ) ] and length-2 [ as an example , see fig .
[ ps_fig3](a ) ] * ps*-sets , this latter with disjoint minimal sets . in fig .
[ ps_fig3](a ) , the application of the conditional poincar map in the minimal set @xmath112 leads to the minimal set @xmath113 .
both sets form a length-2 ps - set .
note that both @xmath112 and @xmath113 can be regarded as a length-1 ps - set . to our numerical precision
, we have checked that there are not basic sets beneath the synchronization region tip . which means that * type - a * set is present for very small but finite amplitude forcing . outside the @xmath71 region
, there is a * p - set * , i.e. , a non - transitive conditional poincar map on the chaotic attractor projection .
so , @xmath114 . in this case
, there is not ps .
more examples of length - q ps - sets in the phase synchronous forced chua s circuit can be seen in @xcite .
we can use the formalism of the conditional poincar map to study the appearance of ps in coupled chaotic systems , as the two coupled rssler oscillators given by :    @xmath115    with @xmath116 , and @xmath117 .
the index denotes systems 1 and 2 .
the subspaces @xmath6 are defined by @xmath6=@xmath118 . in a coupled chaotic system , @xmath119 ( @xmath120 ) does not increase uniformly , but it is given by the time the trajectory crosses the poincar plane @xmath121=0 ( @xmath122=0 ) . for these times , and using the time series @xmath98 , we construct the minimal sets @xmath112 and @xmath113 ( @xmath123 and @xmath124 ) .
the parameters are @xmath125=0.01 and @xmath126 = 0.001 .    in fig .
[ ps_fig3](b ) , we show a length-2 basic set with disjoint minimal sets .
the application of the conditional poincar map in the minimal set @xmath123 leads to the minimal set @xmath124 , and vice - versa .
these two sets form together a length-2 ps - set , but each one separated can be regarded as a length-1 ps - set .
the characteristic of these ps - sets is that they appear as localized structures around the poincar section chosen to define the events .
as one can see , the set @xmath123 is localized in the neighborhood of the line @xmath127 , where the poincar section is chosen .    in fact ,
this length-2 ps - set with disjoint minimal sets ( as well as a length-1 ps - set ) is found everywhere in the ps region , as shown in fig .
[ ps_fig4 ] . in it ,
filled squares represent parameters for which these special ps - sets are found , and empty circles parameters for which ps exists .
a ps - set of length - q is detected using in eq .
( [ practice_basic_sets ] ) @xmath128 , from which we can check whether @xmath129 occupies ( * type - a * discrete set @xmath88 ) the whole space occupied by @xmath130 .
the set @xmath131 in eq .
( [ practice_basic_sets ] ) is constructed assuming squares of size @xmath132 , in points of the set @xmath128 . in fig .
[ ps_fig4 ] , we show a case for @xmath133=2 .    in fig .
[ ps_fig40](a ) , we show the attractor in the subspace @xmath26 , the subset @xmath130 in gray , and the discrete set @xmath134 in dark empty circles . in ( b ) , we show a magnification of the box in ( a ) . note that neighborhoods of arbitrary points in the trajectory of @xmath130 ( gray ) always contain a point of the discrete set @xmath134 ( dark empty circles ) .
so , the set @xmath130 is equivalent to the set @xmath134 and , therefore , @xmath134 is not a ps - set , and therefore , it does not exist ps . differently from the chua s circuit ,
the rssler coupled system presents no p - sets for parameters outside of the ps region .
the approach presented in this paper can be extended for non - coherent attractors . as noticed in ref .
@xcite , attractors that present non - coherent phase motion in the phase space may present a coherent motion in the space of the velocities , i.e. @xmath135 , which is the case of the funnel attractor @xcite . in this case , the extension of our approach to the non - coherent phase motion is straightforward . instead of defining a conditional poincar map in the phase space , we analyze the dynamics in the velocity space , in which the phase is coherent , and therefore , all the ideas introduced herein can be applied in this new space .    to some chaotic attractors
, it might not be possible to define a poincar section to construct the conditional poincar map .
however , one can still construct these conditional maps by defining different types of events .
for example , in coupled neurons , this event can be chosen as being a beginning or the ending of the bursts , which can be well defined by the crossing of the trajectory with some given threshold .
a chaotic set is always transitive through the flow .
so , given a set of initial conditions , its evolution through the flow eventually reaches arbitrary open subsets of the original chaotic attractor .
however , a stroboscopic map of the flow , whose generalization is here called conditional poincar map , might not possess the transitive property .
that is , given a set of initial conditions , its evolution through the conditional poincar map * _ might _ * not reach arbitrary open subsets of the chaotic attractor .
the introduction of the term `` conditional '' in the map nomenclature comes from the unconventional and non rigidly way we adapt the established definition of a stroboscopic poincar map . for coupled chaotic oscillators , this conditional map is constructed based on events , conveniently chosen at the same subspaces where the phase of a chaotic system is defined . in addition , the application of this map through the flow , which results in a discrete set , is inspected not in the whole phase space , but also in the same subspaces where the phase is defined .
if phase synchronization exists the conditional map generates special discrete sets , named ps - sets .
the contrary is also true .
i.e. , a ps - set implies phase synchronization .
this was illustrated in the periodically forced chua s circuit , in the coupled rssler oscillator , and in other more general ( topologically equivalent ) coupled chaotic systems .
the ideas introduced here provide an efficient way of detecting phase synchronization , without having actually to measure the phase of the chaotic trajectory .
indeed , this detection can be done in experiments in a real - time , as it was done here in the perturbed chua s circuit .
it is worthy to say that the ps - sets are robust under small additive noise that could corrupt the data in an experiment .
this is so , because a small additive noise does not interfere much with the time that the trajectory crosses the poincar section , but just deviates the timing of the crossing .
then , the ps - sets still remains .
this robustness against the noise is an important property in order to apply these ideas in experiments , as done in this work .
we have also introduced the phase as a quantity that measures the velocity of rotation of a projection of the tangent vector along the trajectory .
this definition can be used to arbitrary flows , independently whether or not they present a coherent or non - coherent phase dynamics .
finally , our formalism of the conditional poincar map can be used in coupled maps ( or perturbed ) to detect synchronous behavior ( not * phase * synchronization ) between the systems ( or between the forcing ) , for the case where one does not find full synchronization between the maps .
one particular example where that happens is in the periodically forced logistic equation @xcite or for a system of coupled logistic maps @xcite , where one can find a finite number of synchronous chaotic subsets , the basic sets .
* acknowledgment * research partially financed by fapesp , alexander von humboldt foundation , and cnpq .
next , we present some formalism using elements of ergodic theory in order to introduce the conditional poincar map and the @xmath89 sets .    given a flow @xmath136 , we call @xmath7 the chaotic attractor
. a subset of the attractor is its projection on the subspace @xmath6 , that we call @xmath8 .
we assume @xmath137 .
the notion of stroboscopic map of a flow can be generalized to any temporal translation on the trajectory .
we define the temporal translation to be a transformation represented by @xmath138 , called conditional poincar map .
the initial condition @xmath139 is iterated under the temporal transformation @xmath90 , to the point @xmath140 . applying the transformation @xmath90 in a typical trajectory , for the time sequence @xmath141 , @xmath142 , @xmath143 , give us the points @xmath144 , @xmath145 , and @xmath146 .
so , @xmath145 is the point @xmath144 integrated by the flow for a time interval given by @xmath147 .
applying the transformation @xmath90 , for an infinite series of @xmath141 , that is @xmath148 , gives us a discrete set that we call @xmath88 .
the projection of @xmath88 on the subspace @xmath6 is named @xmath89 .
now , we introduce the notion of transitivity .
let us assume we have a chaotic set @xmath149 .
let us choose two disjoint subsets @xmath150 and @xmath151 in @xmath149 .
so , @xmath152 . there is a transformation @xmath153 that generates the set @xmath149 , with the property that @xmath154 and @xmath155 .
so , clearly the transformation @xmath153 and its @xmath156-fold application ( @xmath157 ) can always place an arbitrary initial condition , belonging either in @xmath150 or @xmath151 , anywhere in the set @xmath149 .
this property makes the transformation @xmath153 to be transitive in the set @xmath149 .
however , the conditional poincar map might no longer be transitive .    for example , in the case of the @xmath158-fold of the transformation @xmath108 , i.e. , @xmath159 .
given an initial condition in the set @xmath150 , its iteration by @xmath159 will never reach the set @xmath151 .
therefore , the conditional poincar map @xmath159 is not a transitive transformation in the set @xmath149 . in the case of flows , we have seen that the temporal transformation represented by @xmath90 , whenever ps is present , does not place arbitrary points of the attractor everywhere in the subsets @xmath8 of the chaotic attractor .
therefore , if ps is present the transformation @xmath90 is not transitive to this subset of the attractor
.    the conditional poincar map @xmath90 is topologically transitive in a set @xmath149 @xcite if for any two open sets @xmath160 , @xmath161    to detect whether the conditional poincar map is transitive in some subset , we introduce the notion of equivalence .
two sets @xmath149 and @xmath150 are ( not ) equivalent [ equivalence is represented by the symbol * @xmath92 * , and non equivalence by * @xmath162 * ] if they ( do not ) occupy the same neighboring space . in a more general way :    two sets @xmath149 and @xmath150 are equivalent , @xmath163 , if @xmath164 , a set @xmath151 can be constructed by the union of open sets @xmath165 , open @xmath166 volumes centered at @xmath167 with length @xmath168 , such that @xmath169 @xmath170 @xmath171 , and @xmath172 if @xmath173 .
the set @xmath88 is a p - set if @xmath174 .
if @xmath89 can be decomposed into a collection @xmath175 of subsets of @xmath176 , with @xmath177 , such that a point in @xmath178 iterated by the conditional poincar map goes to @xmath179 , we refer to each minimal set @xmath180 as a _
recurrent decomposition_. in the particular case we have @xmath181 ( constant ) this minimal set is a
_ periodic decomposition_. the number of sets @xmath133 is called _ length _ of the decomposition @xcite .    if @xmath182 is transitive on @xmath8 , then @xmath89 @xmath183 .
if @xmath182 is non transitive in @xmath8 , then a subset @xmath149 can be constructed , such that @xmath184 and @xmath185 .
proofs of propositions 1 and 2 can be done by using eq .
( [ transitive ] ) and the definitions .    definition 3 can be understood by the ps - set of length-1 and length-2 considered in this work .
if the length-1 ps - set is constructed by the time series @xmath97 , the length-2 is obtained by the time series given by @xmath99 @xcite . for the particular case of a periodically forced system , @xmath186=@xmath41 , and so , the length-1 ps - set is the standard stroboscopic map , that collects point every period of the forcing .
the length-2 ps - set is constructed by a stroboscopic map that collects points every half period of the forcing . for the length-2 ps - set , we have two minimal sets , the sets @xmath100 and @xmath101 , with the property that if @xmath187 , this point iterated by the conditional poincar map goes to @xmath101 and @xmath101 goes to @xmath100 , under the conditional poincar map .
to check whether @xmath94 , we do the following , for @xmath188 it exists @xmath189 such that @xmath190 where @xmath131 is a open ball of radius @xmath91 centered at the point @xmath167 .
@xmath91 is a small positive value .
in this appendices , we show that a ps - set exists if , and only if , phase synchronization exists , which in other words ps - set implies ps and ps implies ps - sets .    given a dynamical system @xmath191 ,
let @xmath192 be the flow and @xmath7 the attractor generated by it , we suppose that we have a chaotic dynamics now on .
let @xmath193 be the poincar section in the subspace @xmath6 , and let @xmath194 be the poincar map associated to the section @xmath193 , such that given a point @xmath195 , thus @xmath196 = @xmath197 .
from now on we use a rescaled time @xmath198 . for a slight abuse of notation
we omit the symbol @xmath199 .      to show the _ if _ in the preposition let us start by considering the time interval associated to the return of the point @xmath200 to the point @xmath201 is @xmath202 = @xmath203 , with @xmath142 being the times at which the subsystem @xmath8 crosses the poincar section @xmath193 , in the subspace @xmath6 .    as already introduced , the average return time is given by @xmath204 , and the time is rescaled , such that @xmath205
our hypothesis is that the subsystem @xmath8 has a phase - coherent oscillation , so there is a number @xmath206 for which holds @xcite : @xmath207 the number @xmath208 measures the coherence in the phase oscillation , and is linked to the phase diffusion @xcite .
this equation holds for all @xmath24 , so it implies that for a single oscillation is also true that @xmath209 . in the case of two systems that present ps , it holds @xcite :        let us introduce a new variable that measures the difference between the time interval of two events in @xmath130 and @xmath213 .
this new variable is @xmath214 , from eq .
( [ diferenca_temporal ] ) , it is true that @xmath215 .
now , we analyze one typical oscillation .
given the following initial conditions @xmath216 and @xmath217 , we evolve both until @xmath218 returns to @xmath219 . in other words , we evolve both initial conditions for a time @xmath220 .
so , @xmath221 .
analogously , @xmath222 .
so , given an initial condition in @xmath219 ( subsystem @xmath130 ) evaluated by the time initial conditions return in the section @xmath225 ( of subsystem @xmath213 ) , it returns near the section @xmath219 , and vice - versa .    for a general case
, we have to show that an initial condition , on the section @xmath219 , evolved by the flow for the time @xmath226 still remains close to this section .
in other words , we have to show that the approximation in eq .
( [ tiagoiii ] ) is valid for an arbitrary number of events @xmath24 in the subspace @xmath213 .
now , noting that @xmath227 , from our hypotheses of phase coherent dynamics @xmath228 .
the same arguments used to derive eq .
( [ tiagoiii ] ) for one oscillation , can be extended to an arbitrary number @xmath24 , so we proof the _ if _ of the proposition ( ps implies the existence of ps - sets ) .
now , to show the _ only if _ ( ps - sets imply ps ) , let us say that there is a ps - set . as a consequence , eq . ( [ diferenca_temporal ] ) is valid . then using the definition of phase coherence @xcite , and noting that in ps the average return times in a given poincar section are the same
, we see that eq .
( [ diferenca_temporal ] ) comes from the boundness of the phase , concluding the proposition .
furthermore , let us suppose that the trajectory of the oscillators are perturbed by a small perturbation , which does not destroy the phase synchronous dynamics .
the effect of the small perturbation in the time return of the trajectory to a poincar section is to deviate this time according @xmath229 , where @xmath230 represents the perturbation in system @xmath231 at the moment of the i - th event , with @xmath232 . under these hypotheses about the perturbation ,
we conclude the following result .
this result shows that a ps - set can be constructed in ps states .
moreover , for the coupled rssler - like systems , this result states that this set is confined in an angular region , which is a consequence of eq .
( [ diferenca_temporal ] ) .    to see the relation between the constant @xmath233 and the size of the ps - set
, we do the following . by the time normalization , average time interval between points in phase space
are proportional to their distances .
so , from eq .
( [ tiagoiii ] ) we write @xmath234 , with @xmath235 being the average half length of the ps - set .
a rough calculation shows that in our experiment with the chua s circuit , we have that @xmath236 , and for the coupled rssler oscillators , we have that @xmath237 .
this results completely agree with the theoretical approach done in @xcite .
a. pikovsky , m. rosenblum , and j. kurths , _ synchronization : a universal concept in nonlinear sciences _ , cambridge university press , 2001 ; s. boccaletti , j. kurths , g. osipov , d. valladares , and c. zhou , phys . rep . * 366 * ( 2002 ) 1 .
phase can be defined as a hilbert transformation of a trajectory component , as the angle of a trajectory point into a special projection of the attractor , as a function that grows @xmath31 every time the chaotic trajectory crosses some specific surface , or as the angle of the projected vector field into some subspace with respect to some rotation point .
a chaotic set is always transitive through the flow .
so , given a set of initial conditions , its evolution through the flow eventually reaches arbitrary open subsets of the original chaotic attractor .
however , the conditional poincar map , might not possesses the transitive property .
that is , given a set of initial conditions , its evolution through this map * _ might _ * not reaches arbitrary open subsets of the especial projections of the chaotic attractor , its dynamics stay confined to a subset of the attractor .
@xmath34 is considered to be a rational . however , as shown in ref .
@xcite , ps , as defined by the boundness of the phase difference was found in two chaotic systems for a finite but very large time interval , as @xmath34 approaches an irrational .
therefore , although in this work we consider @xmath34 to be rational , we should make the remark that for the especial situation as the one presented in ref .
@xcite , eq .
( [ ps ] ) can only be satisfied for a finite but large time if @xmath34 is considered to be irrational .
the intermittency observed in the phase difference is characterized as an usual intermittency , by the alternance between a laminar regime and a burst regime .
if the phase difference remains bounded in the interval @xmath238 $ ] , we say that we have a laminar regime .
if the phase difference leaves this interval , we have a burst , also known as phase slip .
as one approaches the border between the ps and the non ps region , the laminar regime in the phase difference becomes longer .
the choice of these time intervals is not unique .
it depends on what type of event one wants to identify in the system .
a length-1 basic set is constructed by measuring time intervals between the occurrence of two events of the same type .
@xmath239 length-2 basic set is constructed by measuring the time interval between the occurrence of an event @xmath239 and an event @xmath240 , and then , between the event @xmath240 and finally the event @xmath239 , and so , on .