progress in single molecule techniques has enabled observations of single steps in many motor proteins , and accurate measurement of the distribution of dwell times , i.e. , the periods of little or no motion between steps .
examples are the forward and backward steps of processive molecular motors like kinesin @xcite , myosin v @xcite , cytoplasmic dynein @xcite or rna polymerase @xcite , and stepwise rotations in atp synthase @xcite and the flagellar motor @xcite . observations of steps and dwell time distributions offer a route to gain insight into the microscopic mechanism and detailed motion of such systems , beyond what is available through knowledge of average turnover rates alone .
dwell times are examples of first passage times , which have been extensively studied in the theory of random walks @xcite .
an important modeling step is therefore to formulate a first passage problem that describes the experimental situation .
close examination of step trajectories with high time resolution from several motor proteins @xcite reveals that the steps are very rapid events compared to typical dwell times . in a discrete - state description , it is therefore reasonable to identify a step with a single transition , and the direction of the step with the direction of that transition @xcite .
the identification of steps with single transitions follows naturally from the assumption that each state has a well - defined average position , and is a basic assumption in this paper .
motor proteins are cyclic enzymes , but contrary to earlier assumptions @xcite , the steps and waiting times obtained in many single - molecule experiments are not well described in terms of the cycle completions in enzyme kinetics .
this was recently demonstrated for models where a single cycle accounts for both forward and backward steps @xcite . in this case , the average number of forward and backward steps differ from the average number of completed forward and backward cycles , and the average dwell times between steps differ from the average cycle completion times
. the differences can be significant even in conditions where almost no backward steps occur .
the basic reason for this difference is that a motor with one or more intermediate states per cycle is in a different state immediately after forward and backward steps .
therefore , an experimental trajectory of forward and backward steps does not give explicit information about completed forward and backward cycles .
the different states obtained just after forward and backward steps also implies that consecutive step directions are correlated @xcite .
this prediction will be confirmed in a later section , where we analyze stepping data from @xcite for a flagellar motor , and find the clear step - step correlations shown in ( b ) .
in contrast , consecutive cycle completions are statistically independent @xcite .
the observation that steps are correlated has important implications for the interpretation of stepping experiments , and also motivates further theoretical study of stepping statistics that goes beyond the assumption @xcite of independent steps and dwell times . in this paper
, we extend the theory of @xcite in two ways .
first , we consider the distributions of conditional dwell times instead of mean values .
second , we consider a larger class of models , including motor detachments , substeps , and multiple pathways .
our main results are a distribution symmetry for conditional dwell times , and a simple relation between the dissipated free energy per cycle , @xmath0 , and certain conditional stepping probabilities . for a large number of models relevant to describe reversible motor proteins
, we get    [ mainresult ] @xmath1    here , @xmath2 and @xmath3 are the probability density functions for the conditional dwell times between two consecutive forward and backward steps respectively , @xmath4 are the corresponding integrated probability functions , @xmath5 is the probability that a forward step is followed by another forward step , and @xmath6 that a backward step is followed by another backward step .
the different types of dwell times are illustrated in , and simply means that the conditional dwell times @xmath7 and @xmath8 are random variables with equal distributions .
however , does not say anything about the distributions of dwell times between steps of different directions , which in general have different distributions .     could be position of a linear motor , net rotation of a rotary motor , or the net number of a substrate or product in an enzymatic reaction , and changes stepwise in steps of size @xmath9 .
the staircase line is an idealized running average .
]    equation holds also at finite average velocity . for example
, a forward - moving motor will take mostly forward steps , but might eventually produce two consecutive backward steps .
the dwell times between such @xmath10 and @xmath11 step pairs have equal distributions , although the probability of observing two consecutive backward steps might be very small for a motor with strong forward bias .    in practice , the need to observe a significant number of @xmath10 and @xmath11 events to test or apply can be an experimental challenge .
for example , atp driven motors like kinesin or myosin v typically have @xmath12 _ in vivo _ @xcite .
this means that @xmath11 events are very rare ( see ) , unless an external load is applied @xcite .    to extract qualitative information out of an experimental test of dwell time symmetry
, one needs a characterization of the class of models which satisfy .
we try to formulate a general characterization of this class .
this means that some of the models within this characterization will not be realistic descriptions of biological systems . on the other hand ,
a large number of candidate models can be excluded if dwell time symmetry is not observed in some system .
we find two sufficient assumptions for to hold , which we call _ strong coupling _ and the _ bottleneck property_. within the discrete state modeling framework that we use , these assumptions are easy to formulate , but can also be disposed of @xcite .
the bottleneck property is an assumption about the model topology .
it means that both forward and backward steps correspond to transitions to or from a single state in the mechanochemical cycle , which we call the bottleneck state .
this ensures that the state of the system is uniquely determined after each observed step , independent of previous step directions .
the bottle - neck property alone results in a particularly simple form for the step - step correlations , given in , which can be tested experimentally .
the assumption of strong coupling , defined mathematically in , is related to microscopic reversibility .
physically , it means that there is no futile free energy dissipation in the mechanochemical cycle .
in addition to tight coupling in the usual sense @xcite , i.e. , a one - to - one correspondence between fuel consumption and forward steps , it includes the assumption that backward steps is tightly coupled to synthesis of fuel . as a result ,
motors with strong coupling are effectively one - dimensional systems , in the sense that only one reaction coordinate , e.g. , position , is needed to describe their operation .
strong coupling might not necessarily hold for all motors .
for example , the forward steps of myosin v are tightly coupled to atp hydrolysis , while backward steps independent of atp binding was recently reported @xcite .
another example is kinesin , for which atp binding during steps in both directions has been reported @xcite . on the other hand ,
strong coupling might be possible in some rotary motors , at least under certain conditions @xcite .
several recent works @xcite use hidden markov models to estimate kinetic parameters directly from experimental trajectories , and thus in principle utilize all information in the data and bypass the difficulties associated with step detection .
however , this approach does not replace the need for a theoretical understanding of the capabilities and limitations of the underlying stochastic models .
our results should provide useful guidance when applying these techniques .
equation corresponds to a similar result in cycle kinetics @xcite . using the symbol `` @xmath13 '' to denote cycle time properties
, one has @xmath14 where @xmath15 are the probability density functions for waiting times before forward ( @xmath16 ) or backward ( @xmath17 ) cycle completions , @xmath18 are the relative frequencies of completed forward and backward cycles respectively , and @xmath19 are one - way cycle fluxes @xcite .
equation is based on similar assumptions as @xcite .
however , the waiting times in do not describe the experimental dwell times observed in stepping experiments on motor proteins @xcite . distinguishing between cycle completion times and dwell times results in significant differences when modeling and interpreting stepping data .
one example of such a difference is the possibility to calculate @xmath0 from observations of forward and backward steps . in this case , using the cycle completion result in on stepping data can give large systematic errors in the estimated @xmath0 , which can be avoided if is used instead .
a related example concerns how to interpret the ratio of observed forward and backward steps , which has been measured for kinesin over a range of forces and atp concentrations @xcite . the force dependence of this ratio can be described as roughly proportional , where @xmath20 nm is the step length , @xmath21 is the applied load , and @xmath22 is a numerical factor significantly smaller than unity .
the dissipated free energy per cycle depends on applied load as @xmath23 this means that predicts @xmath24 if steps are assumed to correspond to completed cycles .
in contrast , gives no general reason to even expect an exponential behavior of the ratio of forward and backward steps @xcite , and suggests that @xmath25 is a more relevant quantity to study . as mentioned above , kinesin might not satisfy strong coupling , in which case one must look at more complicated models . in any case , making the distinction between steps and completed cycles is clearly important in order to interpret the experiments correctly .
the rest of the paper is organized as follows . in the next section
, we sketch a derivation of our results for the models studied in refs .  @xcite .
we also discuss step - step correlations .
after that , we generalize our results to a broader class of models , and discuss detachments and multiple pathways in the mechanochemical cycle .
we also generalize to the case of several detectable substeps . to illustrate our results ,
we then analyze stepping data for the flagellar motor , before the concluding discussion .
detailed derivations are given in the appendices .
states per cycle , @xmath26_l,[2]_l,\ldots,[n-1]_l$ ] . the model can be viewed as a biased random walk on the periodic lattice of states .
numbers in brackets represent different states , and indices different cycles .
state @xmath27_l$ ] denotes the same structural state as @xmath28_l\equiv [ k]_{l+1}$ ] , but with position or angular orientation @xmath29 ( for linear and rotary motors respectively ) shifted by the step length @xmath9 .
arrows indicate the allowed transitions between neighbor states with transition rates @xmath30 , @xmath31 .
fat solid arrows denote the major conformational changes that produce observable steps like those in .
the other transitions are indicated by dashed arrows . in a motor protein , these transitions could for example be binding or release of ligands , or small conformational changes that are hidden in the experimental noise . ]    in this section , we derive and some related results in the simplest case , i.e. , a sequential model .
consider a model with a single cycle consisting of a sequence of @xmath32 states ( see ) , where we defined @xmath27_l$ ] to denote state @xmath33 in cycle @xmath34 , associated with position @xmath35 .
this is a basic model in enzyme kinetics @xcite and has been used to describe motor proteins like kinesin @xcite , @xcite and f@xmath36-atpase @xcite . by construction , it satisfies both strong coupling and the bottleneck property .
we denote the forward ( backward ) transition rates from state @xmath37 to adjacent states by @xmath30 ( @xmath31 ) as indicated in .
the transition rates are positive and periodic , i.e. , @xmath38 , @xmath39 , and possibly functions of external loads and concentrations of various species in the surrounding solution . for the purpose of this discussion ,
the rates are assumed to be arbitrary positive constants , some of which can be tuned experimentally .
we assume that transitions between different cycles , @xmath40_l\rightleftharpoons[0]_{l+1}$ ] , produces observable forward ( @xmath16 ) or backward ( @xmath17 ) steps .
if the model describes an enzyme , the observable step could be the release or uptake of a product or substrate molecule .
counting steps is different from counting cycle completions , which can be illustrated with the following thought experiment .
consider the sequential model in with @xmath41 states per cycle , and a trajectory where the motor goes through the states @xmath26_0\to[1]_{-1}\to[0]_0\to[1]_0\to[0]_1 $ ] .
this completes one forward cycle from @xmath26_0 $ ] to @xmath26_1 $ ] , but produces three steps : one backward ( @xmath26\to[1]_{-1}$ ] ) , followed by two forward ( @xmath26_{-1}\to[0]_0 $ ] , and @xmath42_0\to[0]_1 $ ] ) .
due to such events , neither the number of steps nor the dwell times between steps are accurately described in terms of cycle completions .
a different treatment is needed .
the time evolution of the system is a random walk on the periodic one - dimensional lattice of states , where the average velocity may be nonzero . to simplify the notation
, we now use @xmath37 to denote a state in any cycle , with the convention that @xmath37 and @xmath43 are equivalent states in different positions .
the probability @xmath44 to be in state @xmath37 at time @xmath45 evolves according to a master equation @xcite , in this case @xmath46 state @xmath37 has free energy @xmath47 , and according to detailed balance , the free energy difference between two adjacent states are related to the transition rates through @xmath48 . from the periodicity of the rates we find @xmath49 the direction of the average drift is positive ( to the right in fig.[seqpath ] ) if the dissipated free energy per cycle , @xmath0 , is negative , and zero if @xmath50 @xcite .    tracking
the position produces a series of forward ( @xmath16 ) and backward ( @xmath17 ) steps separated by random dwell times , as sketched in .
our aim is to describe the statistics of such a trajectory , i.e. , the fluctuations of the dwell times and the number of forward and backward steps , which can then be compared with experimental data .    following ref .
@xcite , we introduce the pairwise splitting probabilities @xmath5 , @xmath51 , @xmath52 , and @xmath6 , where @xmath5 and @xmath51 is the probability that a forward step is followed by a forward step or a backward step , respectively , and similar for @xmath53 .
the splitting probabilities satisfy @xmath54 similarly , we introduce random variables @xmath7 , @xmath55 , @xmath56 , and @xmath8 for the conditional dwell times , where @xmath7 is the dwell time between two consecutive forward steps , and so on , as illustrated in .
once the proper first passage problem to describe the dwell times has been formulated , explicit expressions for the splitting probabilities and mean conditional dwell times can be computed for arbitrary @xmath32 @xcite , using standard methods @xcite .
we now describe the first passage problem for dwell times in sequential models , which we then use to derive .
the dwell time symmetry is due to a one - to - one mapping between each @xmath11 event , i.e. , two consecutive backward steps and the dwell time between them , and a corresponding @xmath10 event .
this is the key observation for the generalization of to more complex models in later sections .
we model the splitting probabilities and probability distribution functions for the dwell times as a first escape problem from the interval of states @xmath57 .
one way to approach this problem is the approach with absorbing boundaries @xcite , which means solving a reduced master equation for the states @xmath57 , with absorbing boundaries at both ends @xcite . from the assumption that steps are produced by the transition @xmath40_l\rightleftharpoons[0]_{l+1}$ ]
, it follows that immediately after a @xmath16 step , the system is in the state @xmath58 ( @xmath59 ) , and just after a @xmath17 step , it is in the state @xmath60 ( @xmath59 ) .
we write the reduced master equation in matrix form , , where the matrix @xmath61 has elements @xmath62 and @xmath63 .
the element @xmath64 is the transition rate from state @xmath37 to state @xmath65 .
the probability functions are given by the outgoing probability current , @xmath66 where @xmath67 indicates an initial condition @xmath68 , given by the direction the of the previous step ( see ) .
normalization of the distributions requires that @xmath69      to derive , we show that the taylor series of @xmath70 and @xmath71 are identical up to a factor @xmath72 , from which follows .
the actual calculation is given in appendix a , and we now discuss some of its consequences .
first , it is interesting to note that the periodicity of the model is not necessary for the dwell time symmetry , only to get the simple relation between @xmath5 , @xmath6 and @xmath0 in .
second , the @xmath73 probability distribution has a simple closed form for ( periodic ) sequential models . as we show in appendix b , @xmath74 where the @xmath75 are the eigenvalues of the matrix @xmath61
this is the distribution of a sum @xmath32 of independent exponential random variables with mean values @xmath76 .
explicit expressions for all four conditional dwell time distributions for the case @xmath77 are given in ref .
as expected from the calculations in appendix b , there is no simple relation between @xmath78 and @xmath79 .
equation is only valid if the eigenvalues @xmath75 are distinct .
this is the generic situation , since there are no symmetries or other reasons to expect degeneracies .
however , is valid also in the degenerate case .
finally , it is worth noting that the sequential models in this section can also be obtained as a discretization of overdamped one - dimensional diffusion @xcite in an arbitrary potential @xmath80 between two points a and b ( see ) . in this case ,
our results mean that the waiting time between last touch at a and first touch at b , and the reverse waiting time between last touch at b and first touch at a , have equal distributions .
this pair of first passage problems was previously studied by @xcite .
the fact that the motor is in a different state immediately after steps in different directions means that the step directions might be correlated . in this section
, we formulate this observation mathematically , and point out some experimentally relevant consequences .
let @xmath81 and @xmath82 be the probabilities that step @xmath33 in a trajectory is a forward or backward step respectively . in the simplest case of no substeps or detachments , the definitions of the pairwise splitting probabilities give @xmath83 equation describes the sequence of forward and backward steps as a markov chain @xcite , and the normalization constraints in leave two independent parameters in the @xmath84 transition matrix .
the average frequencies @xmath85 of forward and backward steps have been measured for several motor proteins @xcite . in our model ,
those frequencies are given by the stationary distribution of @xcite , @xmath86 in contrast to the ratio @xmath87 of the number of forward and backward cycles in , the ratio @xmath88 of forward and backward steps is in general not equal to @xmath89 .
looking back at , we see no general reason to expect the ratio of forward and backward steps to depend exponentially on the applied load .
several experiments have divided the dwell times according to the direction of the following step @xcite .
our model gives the forward and backward dwell time probability density functions as @xmath90 since @xmath78 and @xmath79 are in general not related in a simple way , one should not expect equal dwell times before forward and backward steps , in contrast to the case of forward and backward cycle completion times @xcite .
equation implies that step directions are correlated , and gives the step - step correlation function @xmath91 where @xmath92 indicates the direction of step @xmath93 , and @xmath94 is an eigenvalue of the transition matrix in @xcite .
the other eigenvalue is one .
note that @xmath95 since @xmath96 .
in addition , it is reasonable to expect @xmath97 , since a forward step following a backward step can be accomplished by a single transition , while the whole cycle must be passed before a forward step is followed by another forward step .
inserting @xmath97 in , and using the normalization in , we get @xmath98 apparently , we should expect negatively correlated steps , @xmath99 .
uncorrelated steps ( @xmath100 ) occur in systems with @xmath101 , i.e. , simple random walks , as well as some special cases of the extended models in the next section .
an upper bound for the number @xmath102 of @xmath11 events is useful to estimate the regime where our results can be applied in practice .
a rough estimate is given by @xmath103 hence , the number of steps @xmath104 required to get adequate statistics grows at least exponentially with decreasing @xmath0 .
beyond the simplest descriptions in terms of a sequential model for the dominating pathway ( if there is one ) , more complicated situations with multiple pathways are possible @xcite .
a general ( coarse - grained ) description of a motor protein would keep track of both the position and the consumption of fuel molecules , which requires the model to include an effective chemical coordinate in addition to the position of the motor @xcite .
are the dwell time symmetry properties of sequential models valid for more general cases as well ?
as we will see next , the answer is yes for a large class of models , in which the motion is still effectively one - dimensional .    with the extended models , we look for general properties that imply dwell time symmetry , even if not all models with these properties are biologically relevant .
we find that the strong coupling and bottleneck assumptions mentioned in the introduction are sufficient .
this means that if dwell time symmetry is observed not to hold for a particular system , then one can conclude that the strong coupling or bottleneck properties are absent . to model such a system , one must go beyond the extended models presented below , for example along the lines discussed in refs .
the special situation without strong coupling but with the bottleneck property is illustrated below for an example model ( fig .
[ n2example2d ] ) .
the extension of our results to detachments and observable substeps in later subsections is directly motivated by experimental observations of such events @xcite .      in this section
, we describe a large class of extended models which display dwell time symmetry in a little more detail , and indicate how the dwell time symmetry is derived .
we use chemical kinetics , but write a master equation with arbitrary transition rates @xmath105 from state @xmath37 to @xmath65 , @xmath106 non - zero `` diagonal '' rates @xmath107 are allowed , for example to describe irreversible detachments of motors from their tracks @xcite , or other events that can be filtered out experimentally .
we also need to specify the condition of strong coupling .
one way is to demand that the transition rates around any closed loop @xmath108 of transitions satisfy @xmath109 this property is automatically fulfilled for sequential models , and makes it possible to define a free energy @xmath47 for each state @xmath37 , up to an additive constant . from detailed balance , we have @xmath110 for the free energy difference along the transition @xmath111 . hence , says that the sum of free energy differences along any closed loop is zero , which means that free energy dissipation that produces no net motion , e.g. , futile atp hydrolysis , is ruled out .
( note that the mechanochemical cycle , say , from a state @xmath112 to @xmath113 , is not a closed loop in our description . )
we also assume that no transition is irreversible , i.e. , , although the transition rates may be arbitrarily small .    to reflect the cyclic operation of motor proteins , the transition rates are periodic with period @xmath32 , and
there is a well defined free energy per period , @xmath114    an important step in the derivation of is to conclude that the state of the motor immediately after a step is independent of earlier steps . in a sequential model
, this occurs because there is only one step - producing transition per cycle .
this condition can be relaxed somewhat .
it is sufficient to assume the bottleneck property mentioned in the introduction , i.e. , that all transitions corresponding to a step are either to or from a single state , the bottleneck state , as illustrated in .
as the term indicates , the system must visit the bottleneck state each time it goes through the cycle .
whether the bottleneck assumption holds must be determined for each system separately .
for example , deviations from the step - step correlations predicted in means that the assumption is not valid .
although the models considered in this section are more general than the sequential ones , several properties of the sequential models are retained .
these include strong coupling , periodicity , and that steps ( roughly ) correspond to one transition in the enzymatic cycle .
they therefore describe an effectively one - dimensional motion , where the motion along the spatial and chemical reaction coordinates @xcite are tightly coupled to each other for motion in both directions .
hence , backward and forward motion proceeds in opposite directions along the same reaction paths , just as for the sequential models .
states per cycle .
all states inside the gray area belong to the same cycle .
the bottleneck states are @xmath115 .
black arrows denote transitions coupled to a mechanical step of length @xmath9 , and dashed arrows indicate transitions within a cycle , which are undetectable .
it is not possible to go from one cycle to the next without either leaving or arriving at a bottleneck state . in general ,
transitions between any states within a cycle are allowed , and bottleneck states might have a transition to any state in the neighbor cycle . ]
dwell time symmetry for the extended models with detachments can be derived using the same methods as for sequential models , and further details on that derivation are given in appendix c. just as in the sequential case , periodicity is only necessary to establish , while the dwell time symmetry , , follows also without assuming periodicity .
examples of the non - periodic case are the subcycles in , which shows a model with two observable substeps per cycle .
each step in the example satisfies the bottleneck property , but both bottleneck states are in the same subcycle . deriving dwell time symmetry for the different subcycles proceeds as in appendix c , with slightly different initial conditions .
as mentioned above , detachments of a motor from its track are observed in several systems @xcite .
we model detachments by introducing death rates @xcite as extra terms in the diagonal elements of the master equation and @xmath61 matrix of eqs .   and .
the presence of detachments affect the long time behavior of the model @xcite , as well as dwell time distributions and pairwise splitting probabilities .
conditional detachment probabilities also have to be added to eqs .   and .
however , in the derivations of dwell time symmetry in appendix c , the transition rates always enter as ratios of forward and reverse rates , e.g. , @xmath116 in , so that diagonal elements ( @xmath117 ) always cancel . hence , remains valid even in the presence of detachments .       and
subcycle 1 has @xmath119 states , and subcycle 2 has @xmath120 states .
both steps have the bottleneck property , although the bottleneck states @xmath58 and @xmath121 are both in subcycle 1 . if strong coupling is satisfied , both subcycles display dwell time symmetry . ]    for some motors , the full forward and backward steps are divided into observable substeps @xcite . in this section ,
we show how the relation between splitting probabilities in can be generalized , if strong coupling holds and all substeps satisfy the bottleneck property .    to start with , consider the system in , which produces substeps @xmath122 and @xmath123 during each cycle , and assume that each sub - cycle displays dwell time symmetry . from now on ,
we denote a forward step of length @xmath122 with @xmath124 , and similarly for the other substeps .
keeping track of only the @xmath118-steps , we can analyze the system as in previous sections , using four splitting probabilities @xmath125 and dwell times @xmath126 , which satisfy .
on the other hand , if we keep track of and discriminate between all substeps , we could instead measure eight splitting probabilities : @xmath127 , @xmath128 , @xmath129 , @xmath130 , and so on .
summing over all paths between two consecutive @xmath131-steps , we get @xmath132 and similarly @xmath133 this means that the relation between free energy and splitting probabilities for a system with only one step per cycle generalizes to @xmath134 note that we summed over all paths between two @xmath131 steps without explicit reference to all paths that start with a @xmath131 step and end with something else .
this means that is valid also for a system with detachments , or with several parallel pathways .    for systems which go through substeps @xmath122 , @xmath118  @xmath135 in each cycle
, we can use the same argument to relate the splitting probabilities of the full analysis ( all substeps included ) to one where the first substep ( @xmath122 ) is ignored : @xmath136 iterating this transformation to ignore substeps @xmath137 as well , we find the following relation for a cycle with @xmath138 visible substeps : @xmath139 in a complicated system with parallel pathways where a cycle can be completed using different sequences of substeps , holds for every such sequence separately .
states per cycle , with two parallel pathways that produce a step in the @xmath29 direction .
solid arrows are step - producing transitions , while dashed arrows indicate hidden transitions that do not produce an observable step .
shaded boxes are bottleneck states . ]    in this section , we consider a small but non - trivial model with 2 states per cycle , to present an explicit example that illustrates the results of the previous sections .
the model has two parallel pathways as sketched in , and the steady state velocity and effective dispersion are known exactly @xcite .
the bottleneck property is satisfied , and strong coupling is equivalent to @xmath140    computing the dwell time distributions and splitting probabilities ( see appendix d for details ) , we get @xmath141 the eigenvalues @xmath142 are given in appendix d , and @xmath143 . for this model , is obviously equivalent to the dwell time symmetry of . for the other two pairs of steps , we get @xmath144 we see that @xmath78 and @xmath79 differ both from each other and from the @xmath10 and @xmath11 distributions , independent of whether holds or not .
the step directions are anti - correlated , since @xmath145     and @xmath93 respectively .
shaded boxes are bottleneck states .
projecting out the chemical coordinate as explained in the text , this model can be brought to the same form as in . in general ,
the resulting effective one - dimensional model does not satisfy the strong coupling assumption of . ]
the example illustrates the logic of our results .
strong coupling and the bottleneck property together are sufficient conditions for the dwell time symmetry .
the bottleneck property alone is sufficient for the markov - chain description of step directions in , which remains valid even if is not satisfied .    a possible reason for why the model in might not satisfy is illustrated in
this is again an @xmath41 state model , but with an independent `` chemical '' reaction coordinate ( the superscript @xmath93 on the states ) , e.g , the number of hydrolyzed atp molecules @xcite .
for example , the path @xmath26_l^m\leftrightharpoons [ 1]_l^m\leftrightharpoons [ 0]_{l+1}^{m+1}$ ] could be a step driven atp hydrolysis , while the reaction @xmath26_l^m\leftrightharpoons [ 0]_{l+1}^m$ ] could describe an atp independent step .
this model can be transformed to the 1d model in through the projection @xmath146_l=\sum_m [ j]_l^m$ ] , which leaves the rates unaffected .
however , the open - ended reaction path @xmath26_l^m\leftrightharpoons [ 1]_l^m\leftrightharpoons [ 0]_{l+1}^{m+1}\leftrightharpoons [ 0]_l^{m+1}$ ] is transformed into the closed loop @xmath26_l\leftrightharpoons[1]_l\leftrightharpoons [ 0]_{l+1}\leftrightharpoons [ 0]_l$ ] . since strong coupling places no constraints on the rates in the original model in
, there is no reason to assume that holds for the projected model . in the effective 1d description ,
this loop then becomes a slip loop , which dissipates free energy without producing net motion .
hence , the simple step - step correlations of is retained , but not the dwell time symmetry .    in systems where it is not possible to measure the `` chemical position '' with single molecule precision ,
the resulting violation of dwell time symmetry is a useful test for such projected slip loops , if the bottle - neck property is satisfied .
this property can in turn be ruled out , for example if the step - step correlations are more complicated than predicted in .
the most interesting candidates for such tests seem to be the rotary motors driven by ion fluxes , e.g. , the f@xmath147 part of atp synthase , or the bacterial flagellar motor , since it is difficult to measure the ion flow with single molecule precision .
we are not aware of any experiments on atp synthase under conditions that produces both forward and backward steps .
data from a flagellar motor will be analyzed in the next section .
@xmath148    in this section , we apply our theoretical results to analyze stepping data from a recent experiment with a chimaeric flagellar motor @xcite .
as we will see , the data is consistent with the predicted dwell time symmetry and step - step correlations within the experimental uncertainty .
we also estimate the free energy per cycle , and show that the estimate based on cycle completions can have a significant systematic error , compared to .
the flagellar motor propels many swimming bacteria , by driving the rotation of flagellar filaments .
each filament is driven at its base by a transmembrane rotary motor , powered by ion flux ( na@xmath149 in this case @xcite ) down an electrochemical gradient across the cell membrane .
the motor is about 45 nm in diameter , and is believed to contain 13 torque - generating units @xcite . in this experiment , only one unit was active , and the flagellum is expected to switch between 26 distinct orientations per turn , corresponding to a step length of @xmath150 @xcite .
we analyzed one trajectory with close to zero net velocity , and one with finite velocity .
aided by chung - kennedy filtering @xcite to enhance the steps , and the step - finding algorithm described in ref .
@xcite , the raw data was converted to stair - case stepping data , as shown in ( a ) . upon close inspection
, the trajectory seems to be divided into intervals where the step lengths are consistent with each other , but out of phase with steps outside .
this apparent drift might be due to dynamical exchange of the stator units @xcite that anchor the motor to the cell .
we identified intervals of consistent stepping by inspection ( see ( a ) ) , and excluded the drift from our analysis by treating inconsistent steps as detachments , i.e. , each interval of consistent stepping was treated as an independent run .
the result is summarized in tab .
[ stepdir ] , and we now proceed to compare step - step correlations and dwell time distributions from trajectory a with the theoretical predictions .
trajectory b contained too few steps to make such comparisons meaningful .    .
stepping is consistent with the theoretical step length in intervals interrupted by drift events ( arrows ) , and adjacent intervals are out of phase with each other .
( b ) step - step autocorrelation function as defined in .
the measured correlation function ( squares ) is consistent with the theoretical prediction ( circles ) using conditional splitting probabilities from tab .
[ stepdir ] .
symbols are the same size as the estimated standard deviations .
( c ) distribution functions @xmath151 for the four different dwell times . ]
the step - step autocorrelation function for trajectory a is shown in ( b ) .
the prediction @xmath152 of is valid for trajectories without detachments . in order to compare theory and experiments , the splitting probabilities in tab .
[ stepdir ] must therefore be renormalized by a factor @xmath153 to account for the fact that intervals of consistent stepping ( by definition ) contain no drift events .
this gives the theoretical prediction @xmath154 a least squares fit to the experimental correlations gives @xmath155 , in agreement with , but inconsistent with uncorrelated step directions .
having confirmed that the step - step correlations are consistent with , we go on to test the dwell time symmetry . the empirical distribution functions @xmath151 are plotted in ( c ) for trajectory a. error bars are standard deviations from bootstrap estimates @xcite . using the kolmogorov - smirnov test @xcite , we conclude , with 95% confidence , that @xmath56 have different distribution than @xmath7 and @xmath8 , and that @xmath55 have different distribution than @xmath8
. however , the test could not detect significant differences between the other three pairs ( @xmath73 , @xmath156 , and @xmath157 ) .
this is again consistent with the theoretical predictions , and also shows that the statistics is good enough to detect differences between distributions that are not equal .
it is also interesting to note from ( c ) that @xmath7 and @xmath8 are longer than @xmath55 and @xmath56 on average .
this is reasonable , since @xmath157 events in principle only requires a single transition , while the system must go through a complete cycle to complete a @xmath10 or @xmath11 event .
the ionic motive force was not measured independently @xcite , so it is not possible to compare the free - energy estimate of with an experimental value .
however , we can compare our result with that based on , which was used in ref .
the results are summarized in tab.[dgestimates ] .
as expected , the free energy per cycle is close to zero in trajectory a , which has almost no net velocity , while trajectory b clearly has a finite free energy per cycle to drive the rotation .
also note the differences in the estimates based on and , which show that the cycle - completion estimate can give rise to significant systematic errors in estimated free energies if applied to stepping data .
the statistical properties of steps in reversible molecular motors are not the same as that of cycle completions in enzyme kinetics .
to interpret stepping trajectories correctly , this difference must be taken into account .    in this paper , we have extended the theory for dwell times in sequential models @xcite , and derived a symmetry relation for the conditional dwell time distributions for a large class of experimentally relevant models , including parallel pathways , visible substeps , and detachments .
in contrast to the statistics of cycle completions @xcite , the dwell time distributions and splitting probabilities of steps depend on the step directions of both steps before and after the dwell period .    the dwell time symmetry is a consequence of strong coupling , i.e. , tight coupling for both forward and backward steps , and a bottleneck property of the underlying kinetic scheme .
this means that the motion of the system is essentially a one - dimensional random walk in the ( two - dimensional ) space of spatial and chemical reaction coordinates .
this makes it possible to use our results to infer information about the coupling and efficiency of a system from kinetic data , without detailed assumptions about the underlying kinetic scheme . in this respect , it is similar to the relation @xmath159 between the randomness parameter @xmath160 and the number of steps @xmath32 in the mechanochemical cycle of molecular motors @xcite .
an alternative to the discrete formalism used here , is to include continuous spatial degrees of freedom , as is common for example when modeling ratchet - type motor mechanisms @xcite . since such models can be discretized @xcite , it is in principle possible to inquire about strong coupling and the bottleneck property in continuous models as well .
an example of this is diffusion in one dimension ( see fig .  [ contlim ] ) , which strictly satisfies strong coupling and the bottleneck property .
an interesting problem for further research is to quantify how well the dwell time symmetry is preserved in system with small deviations from these assumptions .
we analyzed stepping data from single motor subunits in a flagellar motor @xcite .
the data seems to be consistent with the predicted dwell time symmetry and step - step correlations , although only one trajectory with almost zero velocity had enough steps to make such a comparison meaningful .
the form of the step - step correlations is consistent with an underlying kinetic scheme that satisfies the bottleneck property .
so , does the dwell time symmetry in fig .
[ flagellapc ] indicate that the flagellar motor subunits are strongly coupled to the driving ion flow ?
another possibility is that the steps are equilibrium fluctuations , as indicated by the low @xmath0 estimated from the stepping probabilities . in equilibrium
, there is no free energy change associated with ion transport , which means that dwell time symmetry could be obtained also for a system with loose coupling . to say anything about the coupling in the flagellar motor
, the dwell time symmetry must be tested in a regime with finite ( and constant ) velocity and ionic driving force .
this is in principle a question of observing more steps in such conditions , but a systematic way to identify and separate steps from drift events would also be useful .
the step - step correlations predicted in might be present also in systems that violates the dwell time symmetry .
deviations from this form would say something about the topology of the kinetic pathways in such systems , and it would probably be useful to apply our analysis also to kinesin and myosin v.    an accurate theoretical model is often crucial in order to correctly interpret experiments on such complex systems as motor proteins . many previous theoretical works on steps and dwell times in molecular motors @xcite derive or assume descriptions where consecutive steps or cycles are statistically independent of each other .
as we have demonstrated , this does not apply to the flagellar motor of ref .
we have also presented a large class of simple models where easily accessible quantities like the dwell times and step directions are correlated .
we expect our results to be of practical use in both data analysis and design of experiments .
in particular , further experimental efforts are motivated in order to detect correlations and collect significant statistics for both forward and backward steps .
we now give a detailed derivation of for a sequential model .
the initial condition is the state immediately after a step , and can be written @xmath161 to derive , we compute @xmath162 and @xmath163 using and the initial conditions in .
this gives @xmath164 where the summation over @xmath165 goes from @xmath58 to @xmath60 .
the same calculation for @xmath166 yields @xmath167 the products of matrix elements in eqs.([seqffdiff]-[seqbbdiff ] ) correspond to @xmath168-step paths between states @xmath58 and @xmath60 , plus the extra escape step . for sequential models ,
the shortest such path is @xmath32 steps , so the first @xmath60 derivatives are zero . for @xmath169 , we note that there is a one - to - one correspondence between the nonzero terms in and . for each term including a path from @xmath58 to @xmath60 in , there is a corresponding term for the reverse path from @xmath60 to @xmath58 in .
the ratio of two corresponding non - zero terms is @xmath170 going back to for the elements of @xmath61 , and using @xmath171 by periodicity , we see that @xmath172 for all pairs of corresponding terms .
moreover , @xmath173 , as there is no transition which produces two steps at once .
therefore , @xmath174 for all @xmath175 .
the underlying reason is that the sum of free energy changes is the same along all possible paths going forward one cycle from @xmath58 .    since @xmath176 and @xmath166 are integrals of the solutions of the reduced master equation , which is a finite system of ordinary differential equations with constant coefficients , they are smooth functions which have taylor series . hence , together with the normalization in implies .
a similar correspondence between paths was used in ref .
@xcite to derive .
the periodicity is necessary to get the simple relationship between @xmath5 , @xmath6 and @xmath0 in , but not to establish the dwell time symmetry in . without the periodicity ,
is replaced by @xmath177
to derive explicit expressions for @xmath178 and @xmath179 for sequential models , we first note that the occupation probabilities @xmath44 in are the solution of a system of linear ordinary differential equations with constant coefficients , namely the elements of @xmath61 in . if the eigenvalues @xmath180 of @xmath61 are non - degenerate , the solutions have the form @xmath181 . note that the eigenvalues have negative real parts , to guarantee that the system eventually leaves the interval .
we make the ansatz @xmath182 where the term 1 ensures proper normalization .
we need @xmath32 equations to determine the coefficients @xmath183 .
these are @xmath184 for @xmath185 .
the case @xmath186 was argued just before , and @xmath187 follows from the argument between eqs .   and .
the resulting system of equations for @xmath183 is of the vandermonde type , @xmath188 solving with cramer s rule and setting @xmath189 we arrive at .
the distributions of @xmath55 and @xmath55 can be computed with the same ansatz , but the derivatives @xmath190 and @xmath190 are non - zero for @xmath191 , so the results are more complicated and there is in general no simple relation between @xmath78 and @xmath79 .
in this appendix , we derive for a periodic extended model like the one in , generalizing the first passage problem as we go along .
the derivation proceeds much as for the sequential model , by solving a reduced master equation for the states @xmath57 , @xmath192 with absorbing boundaries at @xmath193 and @xmath32 , and allowing @xmath194 to model detachments .
note that the sum in have at most @xmath195 terms , due to the bottleneck property .
after a forward step , the system is in state @xmath58 , while it could in principle be anywhere in the cycle after a backward step .
the initial conditions , describing the distribution of states just after a @xmath67 step , are therefore @xmath196 since there are several escape transitions for each step , for the dwell time probability functions is generalized to @xmath197 we now proceed to compute derivatives of @xmath198 and @xmath199 , and note that by eqs .
( [ micrev]-[gendg ] ) , the free energy difference @xmath200 is independent of path : @xmath201 @xmath202 from the periodicity of @xmath105 expressed in , we get @xmath203 , which proves and thereby .
the derivation for the case when backward steps ( instead of forward steps ) end in the bottleneck state is analogous .
note that our definition of the bottleneck property means that there are no transitions that produce two steps at the same time .
this rules out discontinuities at @xmath204 in the distribution functions @xmath151 , and justifies our use of taylor expansions .
in this appendix , we solve the example model in fig .
[ n2example1d ] in some detail .
in addition to deriving eqs .
( [ ffbbexample]-[n2xgamma ] ) , this also illustrates a method that can easily be generalized to larger systems and implemented on computer algebra systems . for @xmath205 states per cycle ,
the eigenvalues @xmath75 must be calculated numerically . with this exception
, the method gives explicit expressions for the dwell time distributions in the real time domain , that can be directly compared to experimental dwell time histograms .
the first escape problem for the model in is governed by a reduced master equation as described in , with the matrix @xmath206 the eigenvalues of @xmath61 are @xmath207 and satisfy the relations @xmath208    the initial conditions after forward and backward steps are given by , and the conditional dwell time distributions are given by : @xmath209 using the ansatz in together with eqs.([aic]-[apdfdef ] ) to compute @xmath210 and @xmath211 , we get the following systems of linear equations : @xmath212 solving this and using @xmath213 , we eventually arrive at the probability density functions given in eqs.([ffbbexample]-[fbbfexample ] ) .    finally , we need expressions for the splitting probabilities , which we derive using the adjoint equations @xcite .
let @xmath214 , with @xmath215 , be the probability that a system starting in state @xmath37 will next produce a @xmath67 step . starting with @xmath216 , we note that from state @xmath58 , a backward step can be accomplished either by the next transition being a backward step , which has probability @xmath217 , or by jumping to state @xmath218 and from there eventually get a backward step .
the last possibility has total probability @xmath219 , so that @xmath220 next , we multiply with @xmath221 and rearrange the terms . applying the same procedure for the other @xmath214 , we finally get @xmath222 to get the pairwise splitting probabilities in eqs.([ffbbexample]-[n2xgamma ] ) , we first solve for the @xmath214 , and then weight them according to the initial conditions of : @xmath223    we are grateful to richard berry and yoshiyuki sowa for helpful comments about the stepping data and flagellar motors in general , to michael e. fisher , denis tsygankov and r. dean astumian for stimulating discussions , and to hong qian for sharing ref.@xcite with us prior to publication . financial support from the royal institute of technology , the wallenberg foundation ( m.l . ) and the swedish research council , grant 2003 - 5001 , ( m.w . )
is gratefully acknowledged .


surface plasmon polaritons ( spp ) are electromagnetic bound modes responsible for the transport of light at the interface separating a metal from a dielectric .
their ability to confine light at an air - dielectric interface offers the prospect of developing a new technology consisting of photonic nano - devices@xcite .
active research is currently focusing on the possibility of achieving control over the propagation of spps by means of optical elements that would couple or decouple light to them@xcite . in order to conceive optical elements ( lenses , mirrors , beam - splitters ) able to manipulate spps propagation
, we need to learn more about the interaction of surface plasmons with a sub - wavelength modification of the underling dielectric metal interface .
indeed the interaction of spps with surface sub - wavelength defects on a metal surface is of great interest from a theoretical standpoint.@xcite + in this article we shall study scattering of spps by a _
shallow _ surface defect .
we will consider both indentations of the metal surface ( grooves ) and protrusions on it ( ridges ) .
we shall only deal with bi - dimensional defects , which are deemed infinite in one dimension parallel to the interface ( the @xmath0-direction ) .
different aspects of this problem have been studied before with a variety of numerical techniques @xcite . here
we present a systematic comparison between the different scattering coefficients and provide both analytical expressions and qualitative explanations .
+ it must be noted that in a previous work we presented such a comparison@xcite but within an approximate numerical scheme . within that framework
it was found that ridges and grooves exhibited the same scattering , whenever they are shallow enough . here we will revise that result , which turns out to be valid only for long ( elongated ) defects .
the mistaken outcome of ref .
[ ] for short defects may be traced back to the breakdown of the assumption of small curvature in the defect geometry that was made there . in this paper
we solve the maxwell equations through a discretization method , which does not assume the previous approximation and whose accuracy depends only on the discretization mesh .
we found that , as in the previous work , long asymmetric ridges or grooves with the width much larger than the depth , do scatter very similarly .
however square shallow defects manifest a different scattering efficiency , differing in the relative radiative loss and radiation pattern .
the lack of distinction between these two cases did not emerge in the previous approximate treatment . on the whole
the problem needs to be revisited so as to : @xmath1 substantiate why the approximate result does work in the case of elongated defects , @xmath2 point out what is the correct result in the case of shallow and short symmetric defects , and @xmath3 explain qualitatively how the scattering properties of short and shallow symmetric defects are gradually transformed into the scattering properties of elongated defects , as the aspect ratio of the defect increases .
+ this paper is organized as follows . in sec .  [ s1 ]
we state the basic assumptions on the scattering system as well as the solution method . in sec .
[ coeff ] we rearrange the asymptotic expansions of the far - field to produce the scattering coefficients . namely we express the far - field and the related poynting vector in terms of the field inside the defect . still in this section
we look at an approximation for the scattering coefficients of shallow ridges . in sec .
[ ray ] we explain that , in general , we can not quantitatively represent a scatterer ( however small ) by one mesh .
we explain how we associate a small symmetric ridge or groove to a point dipole . in sec .
[ num ] we look at exact numerical results for the scattering of shallow defects of various horizontal lengths .
we analyze these results and , in the case of square defects , we associate a ridge to a vertical dipole and a groove to a horizontal dipole . in sec .
[ hv ] we produce an analytical model that explains the radiation pattern of the surface plasmons scattered by small square ridges and grooves . in sec .
[ needlesection ] we look at the solutions for the case of shallow and long defects and we present a clear - cut interpretation to support the results of the previous treatment@xcite . finally in sec .
[ transition ] we explain qualitatively that the aspect ratio of the defect determines the orientation of the field induced in a shallow defect .
the considered defects are infinite in the @xmath0-dimension and shallow with depth @xmath4 , where @xmath5 is the free space wavelength .
the defects are going to be illuminated by a monochromatic surface plasmon at normal incidence @xmath6 , associated to an impinging energy flux @xmath7 , defined and derived in appendix  [ sp ] .
therefore , only radiation into p - polarized(tm ) waves needs to be considered .
after we drop , out of symmetry , the @xmath0-dependence on the whole problem the field is expressed as : @xmath8 .
the wavevector in vacuum is : @xmath9 , where @xmath10 . the material making the slab shall be lossless silver@xcite , that is : @xmath11 .
absorption is neglected as we consider non - resonant defects with widths much smaller than the spp propagation length . + as represented in fig .
[ scheme ] , we shall be expressing the source orientation in a cartesian basis @xmath12 , and the scattered fields in a right - handed orthogonal polar basis : @xmath13 finally a question of notation : throughout we shall refer to a bi - dimensional point - source simply as a dipole , but it is meant that their emission as all of the fields are cyclical in the @xmath0-direction . as represented in fig .
[ scheme ] , each object lying in the vacuum semi - space shall be labeled by the superscript @xmath14 while any object lying in the metal shall be labeled by the superscript @xmath15 . in particular , scattering quantities related to ridges have the superscript @xmath14 while the ones related to grooves have the superscript @xmath15 .
the field within the cross - sectional area of the ridge is labeled @xmath16 and the one within that of the groove is labeled @xmath17 .
a groove is a one dimensional defect located in the metal and labeled by the index @xmath18 ]    the green tensor approach is a standard method to solve electromagnetic scattering problems @xcite .
our first task in this section , is to arrive at an explicit expression for the scattered electric far - field .
this is attained by propagating the field induced by a dipole density @xmath19 ( where @xmath20 ) inside the area of a ridge , to a point @xmath21 very far from the source . for a groove
we have the same relation between polarization and field ( except for a change of sign ) @xmath22 . to propagate the field from any of the two , we use the standard formula @xcite : @xmath23 + where @xmath24 is the green tensor for the air - metal background .
the green tensor propagates the emission of a point source at @xmath25 to the distant point of the @xmath21 .
one of the advantages of the green tensor technique is that once the fields inside the defects @xmath26 ( and thus @xmath27 ) are computed numerically the asymptotic expansions of scattered fields become analytic .
this takes us to our second task , which is making a direct connection between the orientation of the induced polarization inside the defects and the far - field radiation pattern , and in so doing define the scattering coefficients .
+ first of all , finding the scattered electric far - field @xmath28 requires the asymptotic expansions of the green tensor .
the derivation is sketched in the appendix  [ gt ] . in what follows
we give some simplifying rearrangements that will let us focus directly on the angular radiation pattern of surface defects .
the asymptotic green s tensor in the radiative zone for either a ridge or a groove can be written in a compact form as : @xmath29    & \times&e^{-i k_z^\nu z'}\;\g_{\infty}^{(\nu)}(\alpha,\r').\end{aligned}\ ] ] in such form we can factor the asymptotic scalar green function out of the dyadic part of the green tensor . from eq.([es ] ) , the direction of @xmath30 results from superposition of @xmath31 , the emission from all induced point polarization elements , or dipole density elements . yet
the direction of each contribution @xmath32 must be independent of @xmath33 .
in other words since electromagnetic waves are transverse waves in vacuum , far from their source , the field emitted by a dipole @xmath31 must be proportional to @xmath34 .
in fact using the standard asymptotic expansions ( see the appendix [ gt ] ) we can write : @xmath35\;\u_\alpha .\end{aligned}\ ] ] where for a ridge : @xmath36 and for a groove : @xmath37 the vectors @xmath38 are p - waves defined in vacuum , while @xmath39 are defined in the metal . a reminder of their expressions at normal incidence , in terms of the angle @xmath40 of fig .
[ scheme ] , is reported in the appendix [ pmodes ] , along with the expression for the fresnel reflection and transmission coefficients : @xmath41 , @xmath42 .
+ we are now in a position to write the expressions for the radiative fields . plugging eq.([asyp ] ) and eq.([inf ] ) into eq.([es ] ) we can separate the electric far field dependence into its radial and angular parts as : @xmath43 here the angular amplitude can be written as : @xmath44 where @xmath45 is the scattering coefficient into radiative - modes : @xmath46 in the last expression the scattered field in the far zone consists of a cylindrical wave , transverse to the direction of propagation @xmath47 , and with a net angular amplitude determined by the integral over the source region @xmath48 .
the latter is actually the important bit in the formula as its squared module determines the radiation pattern . as seen from eq.([angularfield ] ) this angular amplitude results from the superposition of each scattering element taken with its own amplitude , phase and optical path in analogy to how an antenna array determines its effective radiation pattern .
the radiation is given by the intensity or poynting vector in the far field .
accordingly the differential angular scattering cross - section is : @xmath49 finally , the net radiative loss @xmath50 is defined as the integrated angular radiation : @xmath51      whenever the height of the defect is small enough , typically much smaller than the wavelength of the incident light , we can make the approximation @xmath52 .
that allows for some simplification for the angular amplitude of a scattering element _ above the surface_. consider : @xmath53 hence , for shallow defects the green tensor dependence of eq.([asyp ] ) is entirely given by the exponential factors @xmath54 , for both a source in the vacuum semi - space and a source in the metal semi - space .
indeed this turns out to be a major simplification for the relative amplitude of the scattering elements in the air semi - space , which we shall perform in detail section[hv ] . + before that we need to highlight the relation between the green tensor of a defect on the metal slab and in the metal slab , under this approximation .
such relation emerges from the boundary conditions for the green s tensor at the interface , which are : @xmath55 \cdot\u_x=0&\nonumber   \\\end{aligned}\ ] ] @xmath56\cdot\u_z=0.&\nonumber\\\end{aligned}\ ] ] notice that , in the unperturbed system , space is translationally invariant in the horizontal direction @xmath57 and this is reflected in is the @xmath57-component of the vector in eq.([vector ] ) .
because of eq.([asyp ] ) and eq.([inf ] ) , we can turn eq.([gsurfacex ] ) into : @xmath58 the presence of surface charges at the interface implies , from eq.([gsurfacez ] ) , that the @xmath59-components of the vector @xmath60 on either sides of the interface have the relation : @xmath61      let us derive the scattering coefficient into surface plasmon modes . note
that , in this one dimensional problem , scattering will be into both the forward surface plasmon @xmath62 , propagating in the positive @xmath57 direction and the backwards plasmon @xmath63 propagating in the negative @xmath64 direction , as defined in the appendix [ sp ] .
the emission by a point dipole or a point polarization element must result into a plasmon final state : @xmath65 , as shown in the derivation sketched in the appendix [ gt ] .
the asymptotic greens tensor for a source _ upon _ ( @xmath14 ) or _ in _ ( @xmath15 ) the metal is : @xmath66&{\bf e } _ { spp\pm}({\bf r}).\end{aligned}\ ] ]
notice that @xmath67 complies with eq.([gsurfacex ] ) and eq.([gsurfacez ] ) .
consequently the field of the scattered plasmons are : @xmath68 \scap_{p\pm}^{(\nu)}&= &   \int_{a^{(\nu ) } } d\r'\;{\bf e } ^{*}_{spp\pm}(\r')\cdot { \bf p } ( \r')\end{aligned}\ ] ] furthermore the magnetic field related to the field scattered into spps is : @xmath69 where @xmath70 is the magnetic field of a spp , as proved in the appendix [ sp ] .
+ now , if the the source @xmath71 is produced by an incident surface plasmon field ( as is our case ) , we can define the scattering cross - section of into spps as : @xmath72 finally , we can define the total scattering cross - section , which in the lossless case is equivalent to the extinction cross - section : @xmath73
next we are going to develop solutions to point sources in a metal plane background .
however one question may be raised : how do we associate the field induced by a surface plasmon inside a ridge or a groove to a point dipole ?
the answer is the argument of this section .
+ when the field inside a defect is obtained by mesh discretization we assume that the field inside a single mesh is uniform , and deviations from the field at its center are deemed negligible . yet , in general , the field in a defect , can not be represented by the field at its center alone .
let us explain a little bit further this point .
+ for simplicity let us consider a defect in a homogenous medium with dielectric constant @xmath74 , but the argument is the same in other backgrounds . as usual@xcite ,
the field at every mesh is found by solving self - consistently a system of @xmath75 coupled equations : @xmath76where @xmath77 and @xmath78 and @xmath79 is the field at the mesh center .
@xmath80 is a term related to the depolarization of light and comes about from the quasi - static contribution of the green tensor .
@xmath81 is a correction term to the green tensor in the region of the scatterer useful to improve the accuracy of the calculation , when the inhomogeneity is discretized @xcite .
+ in practice , the number of mesh points @xmath75 is increased until the calculation converges to the required precision .
then scale variations @xmath82 of @xmath83 are properly represented in the solution . in the rayleigh limit , for a defect of area
@xmath84 so small that @xmath85 , the scatterer behaves like a point source or a point dipole and the background field ( in this case the illumination ) can be considered uniform over @xmath84 : @xmath86 .
exceptionally , for a circular defect in a homogenous medium with dielectric constant @xmath74 , the net field at any point @xmath87 converges to : @xmath88     in vacuum , and a circular defect of the same area , illuminated by a plane wave .
the dashed line represents extinction coefficient calculated with the mie theory for the circle .
the dielectric constant in the defect is @xmath89 at the wavelength of @xmath90 .
the inset represents the geometry of a discretized circle when inscribed in a square represented by 30x30 mesh points . ]
this is because for the field inside an infinitesimal ( very sub - wavelength ) circular shape is actually uniform and thus scattering by such circular defects can be described by one mesh .
in fact the extinction coefficient@xcite can be derived from the field at the center alone : @xmath91=\\ \label{xtn2}&= & a \im\left [ \delta\varepsilon \e^{*}_b\cdot\e \right]\end{aligned}\ ] ] to prove this numerically we have calculated @xmath92 for a cylinder represented by a single mesh , as in eq.([mie ] ) , and illuminated by a plane wave .
first of all we have checked that the one - mesh cross - section of eq.([xtn2 ] ) , coincides with the mie theory result .
secondly , we have subdiscretized the cylinder into square meshes as rendered in the inset of fig .  [ cylinder3 ] .
as also rendered in the figure , applying eq.([ls ] ) we found that , as the number of meshes grows , the scattering cross - section calculated by the collection of meshes eq.([xtn1 ] ) converges to the initial value of one single mesh of eq.([xtn2 ] ) .
however the field inside of a square scatterer can never be uniform if it is to satisfy real boundary conditions even in a homogenous medium or vacuum .
thus , it can not be faithfully described by one mesh .
this is illustrated in fig .
[ cylinder3 ] , which renders the extinction coefficient for a square defect of the same area as the circle .
as it turns out , the converged value is @xmath93 larger than that obtained by the one - mesh approximation .
remarkably this error is not reduced with the defect size : we obtained the same error for squares with side @xmath94 or @xmath95 .
this is just for reference in the optical range , since we found that the error actually depends on type of defect and on the dielectric constant .
+ however , even if the field is not uniform , a small defect in the rayleigh limit can be represented by a point source at the center of the mesh , with its field equal to the average field over the mesh @xmath96 .
+ indeed if the variation of @xmath97 is negligible over the area of the defect we have : @xmath98\end{aligned}\ ] ] so the object behaves as a point - dipole @xmath99 .
+ the previous results were for a homogeneous background , but they also hold for the inhomogeneous one considered in this paper .
we find that , for a defect above the surface in the optical range , the relative error is about @xmath100 , while it can reach @xmath101 for a defects below the surface . + with very small non - elongated ridges and grooves , such that @xmath102 , the equivalent point dipoles are attained by averaging the fields over the area of the defects as follows : @xmath103 accordingly if we set @xmath104 eq.([angularfield ] ) and eq.([ep1 ] ) for small non - elongated defects become : @xmath105\scap_{p\pm}^{(\nu)}(\alpha)&= & \left[{\bf e } _ { p\pm}^{(\nu)}(0)\right]^ *   \cdot\p^{(\nu)}\end{aligned}\ ] ]
as an illustration consider a square ridge and a groove of side @xmath106 .
we have calculated the scattering into radiative modes and spps without associating the defect to a point dipole but rather using eq.([sigmaout ] ) and eq.([spmas ] ) . in this case
the major task is computing the green s tensor for the plane metal surface required to attain the exact field within the surface defect .
this can be achieved following the prescriptions of ref .
+ similar numerical results for the case of shallow grooves were found in ref .
[ ] using a different computational technique .
+ the out of plane radiation pattern of a surface plasmon scattered by such defects is given in fig .
[ radpati ] .
+ calculations show that , for symmetric defects , the net radiative loss is greater for a groove than for a ridge .
this is so because , while both the scattering into spps and the radiation close to the surface ( at @xmath107 ) are similar , their radiation patterns greatly differ normal to the surface ( @xmath108 ) , where the groove radiation is maximum while the ridge radiation goes to zero .
angular radiative cross section ( eq.([sigmaout ] ) ) and surface plasmon cross section ( eq.([spmas ] ) represented by the almost horizontal lines at @xmath109 and @xmath110 ) , for square defects with 10 nm side , illuminated by a spp on silver at @xmath111 .
the scale is linear but the units are arbitrary .
each concentric line indicates an equal increment of the cross - sections , from the minimum at the the center to the maximum at the outermost . ]    ) ) and surface plasmon cross section ( eq.([spmas ] ) represented by the almost horizontal lines at @xmath109 and @xmath112 ) , for rectangular defects with @xmath113 height and 50 nm width , illuminated by a spp on silver at @xmath114 .
the scale is linear but the units are arbitrary .
each concentric line indicates an equal increment of the cross - sections , from the minimum at the the center to the maximum at the outermost . ]
the ridge radiation pattern is distributed into two lobes on either sides of @xmath108 but the groove radiation pattern forms a single lobe .
this is one of our main result and shall be analyzed in detail in the next section .
the result is not in agreement with those obtained in the approximate treatment ref .
we associate the discrepancy to the breakdown of the condition that the curvature of a short and shallow defect does not vary rapidly , used in that work .
+ notice the fraction of energy scattered into spps ,
i.e @xmath115 of eq.([spmas ] ) , is large .
the values of @xmath115 are represented by the horizontal lines of fig .
[ radpati ] ( the concentric lines indicate their amplitude in a linear scale and in arbitrary units , say , for instance from @xmath116 at the center to @xmath117 at the outermost ) . for both ridges and grooves @xmath118 and @xmath119 are roughly equal . however in the case of ridges , @xmath120 is greater than the maximum value of the scattering cross - section into radiative modes @xmath121 of eq.([sigmaout ] ) , by a factor slightly greater than @xmath122 . for grooves ,
@xmath115 is greater than the maximum value of @xmath123 by a factor slightly smaller than @xmath122 .
+ let us now keep the defects height at @xmath124 and enlarge the width @xmath125 . fig .
[ rect1 ] renders the radiation pattern for a rectangular defect of width @xmath126 nm ( @xmath124 ) .
the emergence of directivity in the out of plane radiation , is part of a transitional behavior , in which the radiation patterns tend to align and , simultaneously , one of the lobes is shrunk while the other is blown up in the ridge radiation .
notice that the scattered energy into spps exhibits the same directivity , going mainly in reflection .
eventually , if we keep enlarging the defects until they are considerably asymmetric the radiation patterns for both ridges and grooves tend to be single overlapping lobes ( see fig .
[ rect2 ] ) .
noticeably , the scattering into spps is greatly reduced .
such similarity is explainable in the approximate framework presented in ref .
[ ] which turns out to be quite acceptable in this limit of large enough defects , as we shall substantiate in sec .  [ needlesection ] .
in sec .  [ transition ] we shall account qualitatively for the transition observed in fig .
[ rect1 ] , explaining why the radiation pattern changes when the defects are enlarged .    ) ) and surface plasmon cross section ( eq.([spmas ] ) represented by the almost horizontal lines at @xmath109 and @xmath112 ) , for rectangular defects with @xmath113 height and @xmath127 width , illuminated by a spp on silver at @xmath90 .
the scale is linear but the units are arbitrary .
each concentric line indicates an equal increment of the cross - sections , from the minimum at the the center to the maximum at the outermost .
[ rect2 ] ]      the equivalence between non - elongated subwavelength defects and point dipoles gives us a chance to investigate in depth the individual radiation pattern of a single scattering element .
[ mode ] shows the averaged the field inside the 10 nm ridges and grooves , as prescribed in eq.([sqraverage ] ) and eq.([sqraverage2 ] ) .
the field induced in a groove is mainly longitudinal while the field inside the ridge is mainly transversal .
this is due to both the illumination and the polarizabilty of the scatterers .
when defects are almost symmetric their polarizibilities @xmath128 are nearly isotropic and so the induced field and the incident field are virtually parallel .
hence the field induced in a ridge and a groove are nearly parallel to the incident surface plasmon @xmath129 , which is mainly perpendicular to the plane in the vacuum semi - space and is mainly parallel to the plane in the metal semi - space .
therefore , in the rayleigh limit , a ridge scatters spps into radiative modes like a vertical dipole on the plane , while the groove scatters them into radiative modes like a horizontal dipole on a plane .
the results for grooves is in agreement with ref .
+ interestingly , we also have found numerically in fig .
[ mode ] that : @xmath130 especially at short wavelengths .
we have devised a virtual source , that can condense the orientation of the equivalent dipole representing a non - elongated symmetric ridges and grooves .
this virtual dipole is defined as : @xmath131 .
the fields inside a groove and a ridge , are respectively , represented as:@xmath132 | \overline{\e}^{r}|&\simeq & |\delta\varepsilon\,\overline{e}^r_{z}{\bf q}(90 ^ 0)|\end{aligned}\ ] ] at least as long as eq.([exz ] ) holds .
+ in reality we can see what happens by means of eq.([mie ] ) . despite the fact that this equation is only exact for a circle in a homogenous background ( as explained ) we can use it to show _ qualitatively _ the relation between the field inside the groove and the ridge , when their shapes are symmetric . if we approximate the polarizability of a ridge for that of a circle in vacuum ( whose polarizability is calculated through eq.([mie ] ) ) , so @xmath133 . if we also approximate the groove polarizability by that of a hole in a homogenous metal medium , we have : @xmath134 .
hence the field induced inside each object is : @xmath135 since these polarizabilities also have the property : @xmath136 ( the polarizability of a hole in a material is @xmath137 times larger than the polarizability of a particle of the same material and the same shape ) then @xmath138 .    ) , eq.([sqraverage2 ] ) ) for a square groove and a square ridge of 10 nm side in silver , as a function of the wavelength .
the scale is logarithmic with arbitrary units . ]
+ the symmetry of the polarizations @xmath128 and the property @xmath139 are strictly true for circular defects in homogeneous media . our numerical calculations of fig .  [ mode ] shows that , even though the field inside a ridge and a groove are quantitatively different from those of circular defects in homogenous media , the assumption that their mutual relation is preserved is in very good agreement with the exact result . because of the symmetry of the square shape , the _ averaged field _ inside
the square is very nearly parallel to the incident field .      as a corollary of the properties of the fields in a ridge and a groove @xmath140 notice that these define @xmath141 through eq.([spmas ] ) .
once @xmath50 from eq.([net ] ) and @xmath142 are determined the value of the transmission of the surface plasmon is a constrained variable : @xmath143 , at least for the lossless case@xcite .
since @xmath50 is greater for grooves than for ridges , the groove transmission is smaller .
the first part of the expression eq.([sigmaout ] ) is a pre - factor @xmath144 whereas the second part is the _ the radiation pattern _ of a point dipole : @xmath145 a groove emits like a horizontal dipole .
the angular amplitude of the field radiated by a horizontal unit dipole @xmath146 , placed close to the interface @xmath147 , is @xmath148 , and it does not matter on which side of the interface it is placed .
@xmath148 can be derived using the relations in the appendix [ key ] and the explicit result is : @xmath149 and the radiation pattern is @xmath150 .
notice @xmath148 presents a mirror symmetry about the angle @xmath151 , the normal to to the plane .
furthermore since @xmath148 never changes sign between @xmath116 and 180@xmath152 ( nor goes to zero ) , the field of a horizontal dipole has one single symmetric lobe , where the field always has the same sign .     of a horizontal point dipole at an air - metal interface .
the radiation patterns varies as the metal dielectric constant is varied .
the scale is logarithmic with arbitrary units . ]
the field intensity @xmath150 of such lobe is rendered in fig .  [ ux ] for different dielectric constants .
this radiation pattern of a groove shown in fig .
[ ux ] , is in agreement with the one represented by ref .
[ ] , obtained with a different numerical method .
notice that for @xmath153 : @xmath154that is , when @xmath137 increases this radiation pattern tends to become simultaneously isotropic and vanishing .
in fact a horizontal dipole does not radiate on a perfect conductor@xcite . on a small digression
it is interesting to notice an apparent contradiction between treatments such as ref .
[ ] , which considered that a defect in a perfect metal were equivalent to a magnetic dipole , while another work@xcite explains a defect in a real metal corresponds to an electric dipole . actually we have just reconciled the two results .
we know that a horizontal dipole on a plane tends to emit isotropically for large @xmath137 .
this means that on a first order expansion in @xmath155 , the radiation pattern of a horizontal dipole on a plane and that of a magnetic dipole in vacuum , are identical . + for finite @xmath137 the field @xmath148 of a horizontal dipole within a real metal would not be thoroughly screened , and while the pattern remains symmetric , its isotropy is disrupted parallel to the surface ( i.e. @xmath156 ) to accommodate the emergence of the surface plasmons density of states .
+ for an individual vertical dipole @xmath157 , which represents a ridge , the angular amplitude of the field is ( see appendix [ key ] ) : @xmath158 the field from a vertical dipole also goes to zero at @xmath156 for a finite @xmath159 , but since dipoles only radiate transversally , the field has a third zero at @xmath160 .
the field is antisymmetric with respect to the normal of the plane , while the intensity @xmath161 is symmetric , and is made up of the two lobes separated by a zero at 90@xmath152 , see fig .
[ symmetric ] .
yet it is important to keep in mind that the field of one lobe is in anti - phase with the field of the other .     from the virtual dipole @xmath162 at @xmath163 .
the scale is linear with arbitrary units . ]
+ unlike a horizontal dipole , the vertical dipole radiative field does not vanish for @xmath153 in fact : @xmath164 the total radiation from a vertical dipole has a larger weight than the radiation by a horizontal one , by a factor of @xmath165 .
this can be seen , in fact , from eq.([gammax ] ) if we assume @xmath159 is large , we get the following relation : @xmath166 in fig .  [ symmetric ] we represent radiation pattern of @xmath167 for the horizontal and vertical orientations respectively , @xmath168 , @xmath169 , which corresponds to our analytic analog of the emission pattern of square ridges and grooves respectively .
while we will consider an intermediate orientation in the next section , we want to remark here that , due to eq.([uguale ] ) , the radiation by both the horizontal moment @xmath170 and a vertical moment @xmath171 vanish parallel to the plane at @xmath156 in a similar manner , as illustrated in fig .
[ radpati ] + at the same time the far - field emissions of ridges and grooves become increasingly different as we approach the direction normal to the plane .
for shallow and long defects @xmath172 and @xmath173 we define the following height - averaged polarization densities and fields : @xmath174 where the last equation defines @xmath175 .
likewise for a groove we can define @xmath176 and @xmath177 through the following equation : @xmath178 notice for @xmath153 we can make the approximation @xmath179 . + the benefit of using @xmath180 is that the scattered - field coefficients for these defects in the far zone , @xmath45 and @xmath181 , are those emitted by a chain of point - dipoles on the surface over the segment @xmath125 , and set at @xmath182 and @xmath183 for ridges and grooves , respectively .
+ the scattered field angular amplitude @xmath184 from eq.([angularfield ] ) and eq.([short ] ) is obtained as : @xmath185     the averaged fields component inside of a ridge and a groove , @xmath186 , the quantity @xmath187 and the rest of the components @xmath186 , @xmath188 , @xmath189 , for rectangular defects of @xmath190 and @xmath191 .
the system is illuminated by a spp in lossless silver at @xmath192 nm .
the scale is linear with arbitrary units . ]
this holds for the scattering into surface plasmon modes as well since we have : @xmath193^*\cdot   \int_0^w dx'\;e^{\mp i k_{px}x'}\ ; { \bf \widetilde{p}}^{(\nu)}(x').\end{aligned}\ ] ] when we illuminate a shallow and long defect , with a spp , an equivalent linear density of dipole sources @xmath194 stems from how the induced fields are distorted inside the scatterer , namely by its polarizability .
when the defect is larger in the horizontal direction than in the vertical one , ridges and grooves were found to give the same scattering by an approximated rayleigh expansion@xcite .
we have an alternative first principles argument to justify the rayleigh expansion result , which is based entirely on the assumption that these defects are needle shaped .
the field induced in these defects tends to be that induced in a needle - shaped protrusion placed horizontally on the surface @xmath182 in the case of a ridge . for a groove
we have a horizontal needle - shaped cavity at @xmath183 . in such idealistic simplification
it is clear - cut to deduce the fields inside the defects from the boundary conditions .
namely the parallel component of the incident field is always continuous and equal , as in eq.([ispp1])and eq.([ispp2 ] ) : @xmath195 which preserves the continuity of eq.([gsurfacex ] ) .
however , we are generating fields which , normal to the surface , make up for the discontinuity perpendicular to the metal surface of eq.([gsurfacez ] ) .
in fact , for a horizontal needle - like ridge , the boundary conditions imposed by the continuity of the displacement vector are : @xmath196 while for a needle - like slit : @xmath197 ultimately : @xmath198 which , matched with eq.([gsurfacex ] ) and eq.([gsurfacez ] ) , yields : @xmath199 and thus the property of producing the same scattering coefficients , previously found in ref .
of course this is just an approximation , but it explains why elongated defects have similar scattering properties . in real life
the plasmon scattering by protrusions and indentations is similar because , far from the edges , a shallow but elongated defect behaves as an infinitely elongated one , as confirmed by numerical calculations . as an example we report in fig .
[ fieldlong ] a numerical calculation of the fields averaged over the height for defects of @xmath190 and @xmath124 .
this shows that eq.([needle ] ) and eq.([needle2 ] ) are quite accurate at the center of the defect , and deviate from the needle model prediction due to fringe effects at the edges .
+ it is worth mentioning that this equivalence is valid in the rayleigh limit when the defect size is much smaller than the wavelength , and may be altered at resonant wavelengths .
everything we just said for symmetric surface defects was based on the fact that their aspect ratio equals one . as the defect width
is increased , the aspect ratio becomes larger and this leads , progressively , to an asymmetric polarizability tensor .
the first effect is that the field induced is gradually less and less parallel to the incident field .
therefore a ridge would develop a non - negligible horizontal electric field component , thus ceasing to be equivalent to a vertical dipole .
likewise the groove , which in the symmetric case behaves as a horizontal dipole , gradually starts having a non - negligible vertical component as its shape is elongated .
the process goes on until we recover the case of a needle shaped defect of section [ needlesection ] .
the fields inside a defect having intermediate width , as in fig .
[ rect1 ] , are intermediate between those for the needle case and the square symmetric case . therefore in these cases defects
emit qualitatively like _ oblique dipoles _ , with the orthogonal components out of phase .
+ in order to understand better the radiation pattern by ridges and grooves we decompose the oblique dipole in its horizontal and vertical components .
+ first of all , we focus on the mechanisms involved radiation pattern for a ridge @xmath14 . from eq.([p ] ) a dipole with arbitrary orientation emits close to the surface , with a field angular amplitude : @xmath200 where @xmath201 and equals : @xmath202@xmath203 shows that the contribution to the radiative field coming from the vertical and horizontal dipole on a metal plane have a phase difference of @xmath160 .
this was already evident from eq.([uguale ] ) , when @xmath204 . such phase difference arises from the impedance of a metal plane@xcite @xmath205 .
+ the radiation pattern for a dipole with arbitrary orientation and lying above the metal , is written in our formalism as : @xmath206 .     for a point dipole : @xmath207 , lying on top of a metal surface .
the scale is linear with arbitrary units . ]
the net angular amplitude for an oblique dipole is resolved into the superposition of the angular envelope of the horizontal dipole ( shown in fig .  [ ux ] ) , with the other radiation factor @xmath208 .
this last factor contains both the orientation and phase of the field . to envisage how these combine we may develop @xmath208 into three terms . these consist in the individual emission from the horizontal and vertical dipole plus an interference term : @xmath209\cos\alpha &   \end{aligned}\ ] ] in the presence of the plane metal background
, we have that horizontal and vertical dipoles behave as individual sources but their interaction presents an intrinsic added phase difference of @xmath160 , which is due to the different interaction of a horizontal and a vertical dipole with the plane . as a result , when in phase they do not interfere , and their radiation pattern is always symmetric regardless of the orientation of the dipole .
this is the case for @xmath210 where , as in fig .
[ symmetric ] , the radiation pattern is the sum of the angular intensity of a vertical and a horizontal dipole , so that at @xmath160 there is a minimum due to the vanishing of the vertical dipole contribution , and yet never goes to zero because of the horizontal dipole contribution . nevertheless , when the dipole components are not in phase , we can get asymmetric radiation patterns and additional zeros ( to those at @xmath211 and @xmath212 ) , because the interaction term can be negative .
in such case the interaction of the horizontal radiative field ( with only one lobe ) with the vertical radiative ( with two lobes of different sign ) is responsible for an asymmetric radiation pattern and exhibits directionality .     for a point dipole : @xmath213 , lying on top of a metal surface .
the scale is linear with arbitrary units . ]
this is illustrated in fig .
[ dipoles2 ] for a dipole emission whose main contribution comes from the vertical dipole . in fig .
[ dipoles3 ] we show the radiation pattern for a dipole whose main contribution comes from the horizontal dipole radiation . + for the case of a grooves ( @xmath15 ) , the radiative angular field amplitude is , from eq.([piz ] ) : @xmath214 \delta^{(2)}(\alpha)&=&p_{x}^{(2)}+ i\ ; \frac { \cos\alpha}{\sqrt{\cos^2\alpha+|\varepsilon|}}\;p_{z}^{(2)}\end{aligned}\ ] ] where remember we have also added the approximation : @xmath215 for @xmath153 .
+ remarkably , as opposed to the the dipole emission over the surface , in the net emission from a dipole under the surface the horizontal dipole contribution has a greater weight than the vertical dipole contribution .
apart from this , all the arguments used for a dipole over the surface apply .
+ the interaction between the vertical and horizontal components of the field induced in the field generates the directional patterns of fig .
[ rect1 ] . for a ridge with length slightly larger than its height
the directional radiation is dominated by its vertical component .
[ dipoles2 ] exemplifies the effect of the interference of a dominant vertical component with a smaller but non - negligible horizontal component . for even larger aspect ratios
the contribution from the other component may be comparable .
+ likewise when a groove has a small aspect ratio it is predominantly a horizontal source interfering with a smaller vertical source .
the result is in an interference pattern that looks like the one rendered in fig .
[ dipoles3 ] . yet again this can be modified by increasing the aspect ratio .
this transition is in good agreement with fig .
11 of ref .
[ ] where , using a different numerical method , the radiation pattern of a groove was computed for different aspect ratios .
our analysis of the surface plasmon scattering by square shallow defects into radiative modes and plasmon modes , reveals that a groove scatters more of the incident energy than a ridge does .
the reflection by a symmetric ridge and a groove is similar and so is the radiative emission close to the horizontal direction .
indeed their scattering essentially differs in the vertical direction , where a groove scatterers while a ridge does not . when defects start to become longer in width we saw the polarizability gets more asymmetric .
correspondingly , since both components of the incident plasmon are out of phase , defects are equivalent to interfering horizontal and vertical dipoles on a plane , which interfere constructively in some direction , thus producing directionality in the radiation pattern .
finally when ridges and grooves are shallow and long they tend to produce the same scattering as , apart for fringe effects , their polarizability exactly counterbalances the discontinuity of the incident surface plasmon field at the air - metal interface .
the authors acknowledge financial support from the spanish ministry of science and innovation under grants no.ap2005-5185 , mat2008 - 06609-c02 and csd2007-046-nanolight.es .
the incident illumination is the field of a surface plasmon wave mode propagating in the positive x direction @xmath216 or negative x direction @xmath217 is : @xmath218 @xmath219 where @xmath220 , @xmath221 and @xmath222 .
this can , alternatively , be written as : @xmath223 for @xmath224 ; and @xmath225 for @xmath226 .
+ the magnetic field associated is continuous at the interface and equal to : @xmath227 now consider a lossless metal , characterized by a real and negative dielectric constant @xmath137 and consider a plasmon moving in the forward direction , ( the subscript + will be omitted ) .
the incident poyinting vector of the plasmon in the air side is : @xmath228 while in the metal is @xmath229 where @xmath230 .
the total poynting vector energy flux associated to a plasmon mode in a lossless metal is : @xmath231
we shall repeat , out of completeness , the explicit expression for p - waves , particularly in the far field when @xmath232 . in this case
these modes are expressed in terms of the direct space polar angle @xmath40 by noticing that @xmath233 and @xmath234 in the air semi - space and @xmath235 in the metal .
hence @xmath236      for reference , we give here the fresnel coefficients for an air metal interface . in the present treatment we only deal with the reflection coefficient for a p - wave propagating from air to metal , and this is : @xmath237 where notice that , for the sake of tidiness , we omit the superscript throughout .
+ as to the transmission coefficients the one for a wave ( 2,1 ) propagating from the metal to air is @xmath238 , while the one for a p - wave transmitted from the air medium to the metal is @xmath239 .
@xmath240 notice that the transmission coefficients are related as follows : @xmath241      the following expressions for the reflection and transmission coefficients are essential to derive eq.([gammax ] ) and eq.([gammaz ] ) : @xmath242    @xmath243 1-r_p(\alpha)=\frac{2 \varepsilon\sin\alpha } { \sqrt{\varepsilon-\cos^2\alpha}+\varepsilon \sin\alpha}\end{aligned}\ ] ]
the asymptotic expressions for the green tensor for @xmath244d scatterers are found in references@xcite .
we have already presented the derivation scheme for bi - dimensional defects in appendix b of ref .
[ ] , for a groove . as explained therein the surface plasmon green tensor and the far - field green tensor are obtained from its angular spectrum . from the relevant sommerfeld integral the surface plasmon contribution
is obtained by applying the residue theorem and the far - field green tensor instead is obtained by applying the method of the steepest descent . + for the case of the ridge we use the total green tensor of the background in the vacuum semi - space .
this can be written as the sum of the direct green tensor ( the free space green tensor ) and the indirect green tensor ( which gives the contribution due to the reflections at the metal plane interface ) .
hence @xmath245 where the spectral representation for the direct green tensor is : @xmath246 while for the indirect green tensor : @xmath247 applying the residue theorem and the steepest descent method to @xmath248 we end up with eq.([asyp ] ) and eq.([gp1 ] ) for ( @xmath14 ) . + for the groove case we need to expand the green tensor connecting a point in the metal to a point in air .
this is just : @xmath249 applying the residue theorem and the steepest descent method to @xmath250 we end up with eq.([asyp ] ) and eq.([gp1 ] ) for ( @xmath15 ) .
notice that the form of @xmath251 given in sec .
[ coeff ] , is obtained by recognizing @xmath252 .
+ one more subtlety , that might be confusing , is how we pass from the transmission coefficient @xmath253 in the integral to the transmission coefficient @xmath254 in the asymptotic form @xmath255 .
this comes about because when we apply the method of the steepest descent to the integral we get eq.([asyp ] ) with : @xmath256 where , in the last equation , we have used the identity eq.([a6 ] ) .