the essential point of stochastic resonance @xcite is the following : if a stochastic resonator is subjected to an external influence , commonly referred to as the signal , its response exhibits signal features , which are most pronounced at a certain level of noise present in the system .
the fascinating aspect of stochastic resonance is therefore that , counter - intuitively , an increasing noise level does not steadily deteriorate the transmission of the signal through the resonator . on the contrary
, noise may rather be used to optimise the transmission process @xcite .
the implications of these effects are intensively investigated . in biological systems , for example , stochastic resonance apparently plays a role in the neural transmission of information @xcite . from a technical point of view
one might possibly see the emergence of a novel type of detector which incorporates an optimal amount of noise to perform best . in this
regard superconducting quantum interference devices ( squids ) have been studied @xcite .
they may be made to detect weak magnetic fields reasonably well without the usual costly shielding from environmental noise .    in recent years
it was shown that the performance of a single stochastic resonator can be enhanced , if it is embedded into an ensemble of other stochastic resonators which are properly coupled @xcite@xcite@xcite@xcite@xcite : compared to being operated isolated the response of the resonator to the signal increases within the coupled ensemble .
however , if the coupling becomes too strong this response was found to deteriorate again .
thus in an ensemble of stochastic resonators the coupling strength turns out to be a second design parameter : apart from the noise level it can be tuned to achieve an optimal performance of the embedded resonator .
this effect was called _
array - enhanced stochastic resonance _ @xcite .
it will possibly find technical exploitation and might also be relevant to biological systems , for example , coupled neurons .
the aim of the present paper is a further investigation of the phenomena related to array - enhanced stochastic resonance . to this end
we study a simple prototype system of coupled two - state stochastic resonators under periodic modulation .
the model we propose allows to analytically calculate the weak - signal limits of two prominent stochastic - resonance quantities , spectral power amplification ( spa ) and signal - to - noise ratio ( snr ) , respectively . following the concept of array - enhanced stochastic resonance
both quantities refer to the response of a resonator as part of the coupled ensemble .
the impact of the coupling on both response measures will be studied .
there is a close link of the present analysis to the work presented in @xcite and @xcite , which to our knowledge are the most comprehensive in the field .
however , compared to @xcite , where the spa of a system of globally coupled bistable elements was studied analytically , too , the present approach offers a more refined description , since it allows to avoid a mean - field approximation .
compared to @xcite , on the other hand , where the snr of a chain of locally coupled elements was investigated in a ( rather extensive ) simulation , the snr can now be obtained analytically . for a quick review of the results
the reader is referred to the 3rd and 4th passage of the summary section .    from a conceptual point of view
the model proposed here can be seen as an extension of the two - state model developed by mcnamara and wiesenfeld @xcite to study stochastic resonance in noisy bistable systems . instead of considering individual two - state elements as in @xcite these elements
are now arranged in a chain .
a simple next - neighbour interaction is introduced which brings the model close to the system studied in @xcite .
the interaction is chosen in a way that the resulting evolution of the elements is given by glauber s stochastic ising model @xcite .
a detailed description of this approach is given in the next section which also provides the necessary background .
we note that the connection to the glauber model is made for mathematical convenience .
we are not concerned with the observation of stochastic resonance in ising systems which was the central theme in @xcite and @xcite .
our purpose will also lead to the assumption of arrhenius - type transition rates of uncoupled unmodulated elements which is unusual in the context of ising systems .
first results of the present analysis had been published in @xcite .    a further interesting approach to arrays of stochastic resonators is c entered around a response of a more collective nature .
here stochastic resonance is studied in the summed output of @xmath0 resonators which need not necessarily be coupled @xcite@xcite@xcite .
recently , it was shown analytically that the snr of this output approaches the input snr , if a sufficiently large number of uncoupled resonators is used @xcite .
whether in some way coupling may still be beneficial for the collective response of the present model will be the subject of a future study .
only the summed output of a very large number of resonators shall be briefly discussed here .
it will be shown that in this case the snr always deteriorates under coupling .
this indicates that the coupling - induced improvement of the performance of stochastic resonators associated with array - enhanced stochastic resonance is a local rather than a global effect .
an early theory of stochastic resonance in noisy bistable systems with one ( generally continuous ) variable was worked out by mcnamara and wiesenfeld @xcite .
they studied the effect of a periodic modulation of these systems in a way which is independent of the precise dynamics involved .
to this end the behaviour of the bistable system was approximated by a random telegraph process . for simplicity
this process was taken to be symmetric , randomly switching between two states @xmath1 .
for convenience we shall consider here @xmath2 .
the probabilities @xmath3 to find the process in state @xmath4 satisfy @xmath5 and their time - evolution is governed by @xmath6    where @xmath7 denotes the rate of the transition @xmath8 .
these rates must be extracted from the precise dynamics at hand . to build a theory of stochastic resonance
they were assumed to be approximately given by @xcite    @xmath9\ ] ]    with @xmath10 being the small parameter of the theory .
obviously , at this level the periodic modulation of the bistable system with ( angular ) frequency @xmath11 and phase @xmath12 is taken into account by the cosine term , only .
this simple model allows to calculate the long - time limit of the power spectrum of the random telegraph process averaged over a uniformly distributed phase @xmath12 .
the spectrum consists of a continuous part and a delta function at modulation frequency @xmath11 .
the continuous part of the spectrum is the lorentzian of the unperturbed process times a frequency - dependent prefactor .
the latter is close to 1 and governs the modulation - induced transfer of broad - band power to the delta peak in the spectrum .
from this spectrum mcnamara and wiesenfeld obtained their central theoretical result with respect to stochastic resonance : an exact analytical expression for the signal - to - noise ratio ( snr ) of the two - state process .
this snr is defined as the ratio of the weight of the delta function to the continuous part of the spectrum at modulation frequency . in view of the present work we shall neglect the signal - induced suppression of the continuous part of the spectrum and consider the linear response approximation of the mcnamara - wiesenfeld result with respect to the small parameter @xmath10 .
doing so one finds for the snr    @xmath13    if the approximation ( [ mwrate ] ) is performed on a particular noisy bistable system subject to a periodic modulation , the parameters @xmath14 and @xmath10 become functions of noise strength and modulation amplitude , respectively .
the dependence of the snr on the noise intensity can then be studied and the occurrence of stochastic resonance may be established for that particular system .    as an example
the overdamped double - well system was given in @xcite .
the model equation reads    @xmath15    where @xmath16 is gaussian white noise with @xmath17 and @xmath18 . using a modified kramers formula for the transition rates valid for sufficiently low modulation frequencies mcnamara and wiesenfeld found    @xmath19 { } \delta = \frac{a}{d}.\ ] ]    the central idea of
the present work is to pass from the single bistable element to a set of coupled bistable elements by extending the two - state model of mcnamara and wiesenfeld . to this end
we consider a chain of two - state elements which for convenience is taken to be of infinite length .
if the elements interact the simple gain - loss balance ( [ mwmaster ] ) has to be modified to describe the evolution of the set of probabilities @xmath20 to find the chain in a particular configuration @xmath21 at time @xmath22 . introducing a formal operator @xmath23 defined for any function @xmath24 by @xmath25
the new gain - loss balance reads    @xmath26    with the initial condition @xmath27 the @xmath20 are interpreted as the transition probabilities @xmath28 of a markovian process with infinitely many discrete components .
the statistics of the process is fully determined by ( [ master ] ) .    to introduce interactions
we assume that the transition rates ( [ mwrate ] ) of an element depend on the states of its next neighbours .
a simple choice for this coupling controlled by a parameter @xmath29 is    @xmath30,\ ] ]    where @xmath31 is the mcnamara - wiesenfeld rate ( [ mwrate ] ) . with positive @xmath29 neighbouring elements
prefer to be in the same state , whereas they tend to be in opposite states , if @xmath29 is negative .
both tendencies grow with growing coupling strength @xmath32 . to avoid negative transition rates @xmath33 has to be required .
we note that the assumed type of coupling is thus able to model ferromagnetic - type interactions for which a coupling - induced improvement of the performance of stochastic resonators was found in @xcite and @xcite . for positive @xmath29
the present model is particularly close to the model simulated in @xcite , since both are linear arrays of next - neighbour - coupled elements . for convenience we shall often use the terms ferromagnetic or antiferromagnetic coupling instead of coupling with positive or negative @xmath29 .    of course
, other types of transition rates ( [ transrate ] ) could be considered as well .
the advantage of the present choice is that most of the relevant stochastic properties of the resulting model are already known .
they were studied by glauber who introduced this model as a stochastic form of the ising model @xcite .
( we note that in the glauber model the term @xmath34 in ( [ transrate ] ) is replaced by a general time - dependent parameter @xmath35 . ) as in the mcnamara - wiesenfeld model it is again possible to find an analytical expression for the snr in leading order of the modulation parameter @xmath10 which is the subject of the next section .
the respective analysis is essentially an exploitation of glauber s work .    it may be possible to reduce the dynamics of some coupled systems to the two - state model given by ( [ master ] ) and ( [ transrate ] ) and hence to directly employ the snr formula derived in the next section
however , this approach is certainly more involved than the rate expansion ( [ mwrate ] ) needed to make use of the mcnamara - wiesenfeld snr .
we shall not pursue this more general aspect of the model in further detail .
instead we are interested in devising a simple prototype system which allows to study the impact of the coupling on the stochastic - resonance effect . to this end
we are looking for simple assumptions on the dependence of the model parameter @xmath14 , @xmath29 , and @xmath10 on some noise intensity and signal amplitude , respectively .
we proceed in two steps :    first , we retain the glauber s original relation of the present model to the ising model given by the hamiltonian    @xmath36    this relation is defined by an assumption on the probability @xmath37 to find the @xmath38-th element in state @xmath39 , if all other elements are fixed : in an adiabatic limit the @xmath37 of both models are to be identical .
this implies that    @xmath40    holds which allows to establish a relation between the parameters of both models .
one finds @xmath41 and @xmath42 , respectively @xcite .
the role of the noise intensity is thus played by the temperature @xmath43 . to skip unnecessary parameters we set @xmath44 and @xmath45 . in the present paper
we assume that @xmath46 , which shall be referred to as the signal , is given by    @xmath47    for weak amplitudes @xmath48 this results in @xmath49 .
( in section [ srsnr ] we shall also consider periodic signals ( [ forcing ] ) with an additional noisy background . )
second , a plausible assumption has to be made on the temperature dependence of the time - scale parameter @xmath50 which is not effected by requiring ( [ detailedbal ] ) .
it is assumed that it retains the qualitative dependence ( [ mwparam ] ) found by mcnamara and wiesenfeld for the double - well system , i.e. , @xmath51 with properly scaled temperature @xmath43 .
this dependence is typical for many rate processes .
an appropriate temperature dependence of @xmath14 is also necessary to turn the elements into stochastic resonators for vanishing coupling . for simplicity
@xmath52 is considered .
the prototype system we are looking for is thus specified by ( [ master ] ) and ( [ transrate ] ) together with    @xmath53 { }    \gamma = \tanh \left(\frac{2j}{t } \right ) , \makebox[0.3 cm ] { }    \delta = \frac{h_0}{t}.\ ] ]    we emphasise that this system is not an approximation of the ising system ( [ hamiltonian ] ) in the same way as the double - well system ( [ doublewell ] ) is approximated by ( [ mwmaster]),([mwrate ] ) , and ( [ mwparam ] ) , respectively : the ising system and the stochastic glauber model are based on totally different dynamical concepts .
nevertheless the devised prototype system appears to capture in a simple way the essence of the effect of coupling and forcing on coupled continuous - variable systems like the double well .    to motivate this ,
let @xmath54 denote the height of the barrier to be surmounted in order to escape from state @xmath4 .
in a very simple approximation the effect of next - neighbour coupling and forcing on the barrier height could be modelled as    @xmath55    @xmath56    comparing this to ( [ hamiltonian ] ) one finds that within this simple picture the effect of coupling and forcing on the barrier heights is the same as their effect on the energy levels of the ising model .
( we note that for coupled double - well systems the chosen approximation would only hold for weak coupling . for strong coupling the systems
may no longer be bistable and the entire concept of a simple - minded two - state approximation breaks down . )    within an adiabatic limit and with some noise strength @xmath43 the transition rate @xmath57 of the barrier system would roughly read    @xmath58    from here rates of the form ( [ transrate ] ) together with parameters ( [ params ] ) can be obtained in two steps .
first , the factor related to @xmath59 is kept as the parameter @xmath14 with a respective choice of @xmath59 .
second , the factor related to @xmath60 , which describes the effect of forcing and coupling , is replaced by the respective term of ( [ transrate ] ) , which is easier to handle . a connection between both types of rates is again made by requiring that they lead to identical stationary distributions in the sense of ( [ detailedbal ] ) .
a particular close relation between the glauber - dynamics model and the barrier system should thus arise within the regime of quasistatic response .
nevertheless we shall not restrict the following investigation of the simple prototype system to the quasistatic regime . whether or not the suggested model provides insights which are also relevant to more realistic settings and possibly helpful in the design of coupled stochastic resonators
remains to be seen .
the aim of this section is to calculate the snr of the response of a two - state element embedded into the chain of coupled elements given by ( [ master ] ) and ( [ transrate ] ) .
this calculation will be done in terms of the parameters @xmath14 , @xmath29 , and @xmath10 and is essentially based on glauber s work .
we shall only take account of terms in leading order of the modulation amplitude @xmath10 .
the resulting snr will then be an extension of the simplified mcnamara - wiesenfeld result ( [ mwsnr ] ) . at the end of the section the general result will be applied to the prototype system specified by ( [ params ] ) .
in addition to the snr we shall also give an analytical expression for the so - called spectral power amplification ( spa ) .
following glauber we start by calculating the averages @xmath61 . from ( [ master ] ) one derives @xmath62 .
it results in    @xmath63 \nonumber \\ & & +   \left[1 - \frac{\gamma}{2 } ( r_{k-1,k } + r_{k , k+1})\right ] \delta \cos ( \omega t + \phi)\end{aligned}\ ] ]    with @xmath64 .
a closed set of equations for the @xmath65 is found by linearising ( [ mean ] ) in @xmath10 . in this case
the @xmath66 can be taken from the unperturbed model ( @xmath67 ) .
after long times these @xmath66 read @xmath68 with @xmath69 @xcite .
the resulting long - time - limit set of equations can be further simplified .
since all elements are forced uniformly , a particular element can not be distinguished from any other , once the initial distribution has been forgotten .
hence , all elements will have identical statistics and the indices in ( [ mean ] ) can be skipped . together with the previous assumption one
obtains    @xmath70    obviously , the long - time dynamics of the average @xmath71 is identical to the long - time dynamics of the average of an uncoupled element with rescaled relaxation rate @xmath72 and rescaled modulation term @xmath73 .
the resulting long - time limit of the averaged state simply reads @xmath74 @xmath75 { } \tan \psi = - \frac{\omega}{\alpha(1-\gamma)},\ ] ]    where @xmath76 is the response to static signals @xmath77 :    @xmath78    formally this static response might even exceed @xmath79 .
we conclude that in this case @xmath10 is not sufficiently small and the linearisation of ( [ mean ] ) is no longer justified . later on it
will be shown that this does not amount to a considerable restriction .
the power spectrum of an element is determined from the average @xmath80 . like in @xcite
its @xmath22-dependence is removed by averaging @xmath81 over an uniformly distributed initial phase @xmath12 of the modulation term in ( [ transrate ] ) . introducing the correlation function @xmath82 [ \sigma_k(t+\tau)-\langle \sigma_k(t+\tau ) \rangle ] \rangle $ ]
one finds    @xmath83    the second term on the r.h.s
. contributes a delta function with weight @xmath84 at signal frequency @xmath11 to the ( one - sided ) power spectrum as defined below .
the first term , on the other hand , forms the continuous part of the spectrum .    for the purpose of calculating the snr in leading order of the modulation amplitude
@xmath10 it is sufficient to approximate this continuous part by the power spectrum of the unperturbed model .
this implies that we neglect any possible signal - induced transformation of the continuous part of the spectrum as we did to obtain the simplified mcnamara - wiesenfeld result ( [ mwsnr ] ) . for the unperturbed model the long - time correlation function @xmath85 is not effected by the average @xmath81 .
it also no longer depends on the index of the element .
according to glauber @xmath85 reads    @xmath86    with modified bessel functions @xmath87 . using the relation @xcite    @xmath88    which holds for @xmath89 and @xmath90
the power spectrum @xmath91 is calculated .
for the one - sided spectrum defined by @xmath92 at @xmath93 we eventually obtain    @xmath94    with @xmath95 and @xmath96 .    now the snr of the response of a two - state element embedded into the chain can be calculated , which is the central result of this section .
it reads    @xmath97    for vanishing coupling it reduces to the simplified mcnamara - wiesenfeld result ( [ mwsnr ] ) .
in addition to the effect of the coupling on the snr we also wish to study its impact on the delta peak in the power spectrum . in this regard the spectral power amplification ( spa ) is a convenient measure ( cf .
e.g. jung in @xcite ) .
it is defined as the ratio of power contained in the signal peaks of output to input spectrum , i.e. , the ratio of the weights of the respective delta functions .
because of its dependence on the input spectrum the spa depends on the precise dynamics to be modelled by the two - state chain .    for the prototype system specified by ( [ params ] ) the spa reads @xmath98 , or explicitly    @xmath99    @xmath100    where @xmath101 is the spa of static signals . obviously , the spa does no longer depend on the signal amplitude @xmath102 . for convenience
we also remove the @xmath102-dependence of the snr by considering    @xmath103    for the prototype system the parameter @xmath104 involved in the unperturbed spectrum ( [ spectrum ] ) simplifies to @xmath105 @xcite .    both quantities ,
spa ( [ spa ] ) and rescaled snr ( [ snr ] ) , only depend on three parameters : temperature @xmath43 , coupling strength @xmath106 , and signal frequency @xmath11 .
restrictions arise from @xmath48 , which lead to a simple expression for @xmath10 , as well as from the linearisation of ( [ mean ] ) , which results in an upper bound on @xmath10 as discussed in connection with ( [ meanstatic ] ) .
together one finds    @xmath107    it implies that finite @xmath102 place a lower bound on @xmath43 and an upper bound on @xmath106 .
however , these restrictions are rather weak : @xmath102 can be made arbitrarily small because its size is immaterial within the present weak - signal approximation .
we are now in a position to study the impact of the coupling on the response measures spa and snr of a single resonator embedded into the chain
. this will be the subject of the following two sections .
the spa ( [ spa ] ) has an unique maximum over temperature @xmath43 and coupling parameter @xmath106 for any time - dependent signal ( [ forcing ] ) .
the maximum spa is obtained for a frequency - dependent value @xmath108 which is always positive . in other words : a properly tuned ferromagnetic - type coupling yields the best spa performance of the two - state resonator element which is embedded into the chain .    for our simple model
this maximum can be studied analytically .
the partial derivatives of the spa with respect to @xmath43 and @xmath106 are found to vanish at pairs @xmath109 given by    @xmath110    and    @xmath111    respectively .
for @xmath108 to be real and finite , @xmath112 has to hold . with this restriction ( [ simple2 ] )
has a unique solution @xmath113 for any given @xmath114 .
the corresponding @xmath108 is found via ( [ simple ] ) . within the given range of @xmath115 one
can easily check on ( [ simple ] ) that @xmath108 is indeed always positive .
an inspection of the second derivatives of the spa finally reveals that there is a maximum at @xmath109 .
its peak value is implicitly given by    @xmath116    from ( [ simple ] ) , ( [ simple2 ] ) , and ( [ spamax ] ) the following features of the maximum can be derived which are shown in fig.1 : tuning the signal frequency @xmath11 from very large to vanishing small values the temperature @xmath115 falls from @xmath117 to @xmath118 . at the same time , the coupling strength @xmath108 as well as the peak height @xmath119 increase from vanishing small to very large values .
in addition , an analysis of the curvatures @xmath120 and @xmath121 , both expressed in terms of @xmath115 , shows that the sharpness of the maximum grows as its height increases .
this indicates that to achieve optimal performance the system parameters @xmath43 and @xmath106 have to be tuned with increasing accuracy as the signal frequency decreases .
fig.1 also includes a comparison to the spa of an uncoupled element . at given @xmath114
this spa has a maximum located at a temperature @xmath122 , implicitly given by    @xmath123    on the interval @xmath124 , whereby its peak value is found to be @xmath125 .
it can be shown from ( [ simple2 ] ) and ( [ simple3 ] ) that @xmath115 always exceeds @xmath122 .
the fact that @xmath108 was found to be never zero also implies that @xmath119 always exceeds @xmath126 .
however , looking at fig.1 it is obvious that at frequencies @xmath127 both compared quantities of the coupled element approach those of the uncoupled one and @xmath108 approaches zero .
the coupling - induced increase in the spa is thus vanishingly small at sufficiently high frequencies of the signal .
this comparison allows to distinguish between two different types of spa behaviour : at low frequencies the spa is enhanced under ferromagnetic coupling whereas at high frequencies it is basicly not , although a tiny increase still occurs .
both situations are illustrated at selected frequencies in fig.2 and fig.3 , respectively .
the low frequency spa exactly reproduces the qualitative effect found analytically by jung et al .
@xcite in a system of globally interacting elements .
the high frequency behaviour , on the other hand , where the spa does practically not increase under coupling , has to our knowledge not been reported before .    in both cases fig.2 and fig.3
clearly demonstrate that at fixed coupling parameter @xmath106 stochastic resonance occurs : the spa has a maximum over temperature @xmath43 which , unfortunately , can not be established analytically .
fig.2(bottom ) shows that while ferromagnetic coupling improves the spa at any fixed temperature @xmath43 this improvement is lost , if the coupling becomes too strong .
there is thus an optimal coupling strength @xmath128 for every given temperature , which can even be exactly calculated at any @xmath114 .
one finds    @xmath129    which is always positive .
since the spa does not have a further extremum over @xmath106 , it implies that antiferromagnetic coupling always decreases the spa .
moreover , @xmath128 increases with growing temperature t and decreases as the signal frequency grows . in fig.2 the spa at the optimal coupling strength is included ( dash - dotted curves ) . in fig.3(bottom ) , where the high - frequency spa is illustrated , the coupling - induced improvement of the spa is hardly detectable and @xmath128 is almost zero .    considering the spa of time - independent signals ( @xmath130 ) given by ( [ spa_static ] )
one finds that it has a maximum over temperature for antiferromagnetic coupling @xmath131 , only , whereby it decreases as the coupling strength @xmath132 grows . for @xmath133
the spa increases with increasing @xmath106 as well as with decreasing @xmath43 .
it formally diverges for @xmath134 and @xmath135 , respectively . in both cases
the weak - signal limit ( [ range ] ) breaks down .    in general ,
the spa ( [ spa ] ) is given by its static value ( [ spa_static ] ) times a dynamical factor @xmath136 ( cf.([meanresult ] ) ) . here
@xmath137 is the ratio @xmath138 of signal frequency to long - time relaxation rate , which also governs the phase shift @xmath139 in ( [ meanresult ] ) .
with growing @xmath137 the elements gradually loose their ability to follow the signal : the spa weakens and the phase shift grows .
this effect occurs , for example , if the signal frequency @xmath11 increases .
subsequently , the spa decreases with growing @xmath11 .
the impact of @xmath137 is also responsible for the occurrence of stochastic resonance and optimal coupling in the spa . here
a decrease of the long - time relaxation rate @xmath72 plays the crucial role : the elements dynamics slows down as @xmath106 increases or as @xmath43 decreases .
thus this slow - down occurs whenever the static spa grows .
hence , the dynamical factor always counteracts the static spa as @xmath43 or @xmath106 is changed .
eventually , the increase of the monotonous static spa is outperformed by the decrease of the dynamical factor which results into a maximum of the spa over @xmath43 and @xmath106 , respectively .
in other words , stochastic resonance and optimal coupling occur .    as illustrated in fig.2
the spa can be seen as a transition between two limits , the static spa given by @xmath140 ( @xmath141 , dashed curves ) and @xmath142 ( @xmath143 , dotted curves ) , respectively . since both limits intersect at @xmath144 , i.e. at @xmath145 , the plots nicely show the well - known approximate matching of time scales at the spa peak .
this matching does not only occur over @xmath43 ( fig.2 , top ) , but over @xmath106 ( fig.2 , bottom ) , too .
if the signal frequency is changed , the peaks of the spa shift whereby all curves share the static spa as a limit .
this results into the qualitative frequency dependence found by mcnamara and wiesenfeld for the delta function part of the spectrum of the double - well system ( [ doublewell ] ) .    roughly speaking
, an improvement of the spa under ferromagnetic coupling only occurs at @xmath146 , where the impact of the static spa is not yet outperformed by the dynamical factor .
this explains why there is no improvement at signal frequencies @xmath147 : since under ferromagnetic coupling @xmath148 holds , @xmath146 can not be fulfilled in this case . at @xmath149 ( fig.3 )
the spa is already within line width of the limit @xmath150 , where the impact of the dynamical factor prevails .
this limit is basicly given by the modulation term in ( [ simplemean ] ) .
its non - monotonous temperature dependence is known to be no longer associated with a matching of time scales @xcite .
the static response , which is thus the origin of the coupling - induced enhancement of the spa , can be interpreted as follows . without coupling and for a given fixed signal the two - state elements prefer to be in the state of lowest energy ( [ hamiltonian ] ) or highest barrier ( [ deltau ] ) , respectively . since the signal is homogeneous , this state is the same for all elements . without a signal but under ferromagnetic coupling neighbouring elements
tend to be in the same state , too , although they do not favour a particular state @xmath151 or @xmath152 .
taking now signal and coupling together , both effects add up . the tendency to find the elements in the state favoured by the signal grows and hence @xmath76 and the static spa grow compared to the uncoupled case .    with antiferromagnetic coupling , on the other hand ,
neighbouring elements prefer to be in opposite states .
this counteracts the effect of the signal and leads to a decrease in the spa .
finally , the decrease of the static spa with growing temperature directly follows from ( [ detailedbal ] ) . there
the imbalance in the distribution of probability between the two states , which dependent on the view taken either differ in energy or barrier height , decreases as the system heats up : the mean response @xmath76 weakens due to increasing fluctuations .
the investigation of the snr presented here basicly relies on numerical evaluations of ( [ snr ] ) .
an explicit calculation of the power spectrum ( [ spectrum ] ) already yields a rather complicated expression which does not lend itself to a detailed analytical study .    at finite signal frequencies
@xmath114 the snr displays the same qualitative dependence on temperature @xmath43 and coupling parameter @xmath106 as shown for the spa in fig.2 and fig.3 , respectively . at low frequencies the snr is enhanced under ferromagnetic coupling whereas at high frequencies it is not . due to this close similarity we omitted the respective plots for the snr .
we note , however , that this similarity is by no means a trivial result , since the spectrum ( [ spectrum ] ) is itself a non - monotonous function of @xmath43 and @xmath106 .
the low - frequency snr qualitatively reproduces the snr behaviour found in a chain of next - neighbour - coupled overdamped double - well systems simulated in @xcite .
it is the first analytical confirmation of this behaviour .
we note that this correspondence occurs although the simulation in @xcite was performed with strong forcing while the present model is studied within a weak - signal limit . as for the spa
, the high - frequency behaviour of the snr has to our knowledge not been reported before .
if the elements are not coupled , their snr ( [ mwsnr ] ) does not depend on the signal frequency .
this was already found in @xcite , if the signal - induced reduction of the continuous part of the spectrum is neglected , as ( [ mwsnr ] ) clearly shows . under coupling this frequency independence of the snr
for ferromagnetic coupling numerical evaluations of ( [ snr ] ) predict a decrease of the snr with growing signal frequency . only in the limit of low and high frequencies
the snr approaches constant values . in both cases
the snr formula ( [ snr ] ) simplifies significantly .
one finds    @xmath153    @xmath154    where @xmath155 is the simplified mcnamara - wiesenfeld snr ( [ mwsnr ] ) of uncoupled elements .
@xmath156 is the static snr while @xmath157 represents the leading order term of its high - frequency expansion .
fig.4 illustrates the typical qualitative behaviour of the snr as a function of the signal frequency . at various frequencies
it shows a set of snr curves over temperature @xmath43 and coupling parameter @xmath106 , respectively .
the dashed curves represent @xmath156 and @xmath157 , respectively .
the apparent fact that the snr decreases with growing signal frequency implies that its static value can not be exceeded at any other signal frequency .
then it follows immediately from ( [ snrlim1 ] ) that the coupling - induced enhancement of the snr possesses an upper limit : it can not be better than a factor four .
this is a rigorous result for the general chain model and not limited to the prototype system , where a special choice of the parameters was made . due to the simplicity of the chain model
this result may well reflect a limit of more general nature for the improvement of a weak - signal snr .
fig.5 shows the static snr over temperature @xmath43 .
it differs qualitatively from the static spa which according to ( [ spa_static ] ) has a monotonous behaviour .    turning briefly to antiferromagnetic coupling
one finds that the continuous part of the spectrum is insensitive to the sign of the coupling parameter @xmath106 .
now even the static snr decreases under coupling . as shown in fig.4(bottom ) it is found below the high - frequency snr curve , i.e. , here the snr increases with growing signal frequency . since at medium frequencies
again a transition occurs between both curves similar to the transition in fig.4(top ) , a second maximum in the snr over @xmath43 may emerge in this situation .
the investigation of the impact of the coupling on the snr can also be extended to the more natural situation , where the signal itself is embedded into noise .
the question to address is , whether the improvement of the snr , which was established for independent internal noise sources , can still be found with additional coherent external noise .
( this problem does not occur for the spa , where only the height of the signal peak is of interest . )    to study this case it is assumed that the input spectrum consists of the previous signal peak described by @xmath158 and a noise part @xmath159 . within the weak - signal limit it was shown that the signal peak of the input spectrum leads to a respective peak in the output spectrum at signal frequency , only : no additional peaks at multiples of that frequency occur . therefore , any additional contribution to the continuous part of the output spectrum at signal frequency can only arise from @xmath160 .
in analogy to ( [ meanresult ] ) one finds that this additional contribution is given by @xmath161
.    together with ( [ snr ] ) the resulting snr @xmath162 reads , again in units of @xmath163 , @xmath164 it can be expressed in terms of the snr @xmath165 ( cf.([snr ] ) ) and the input snr @xmath166 written in units of @xmath163 , respectively .
one finds @xmath167    @xmath162 is thus a steadily growing function of the previously studied snr @xmath165 , whereby it can not exceed the input snr @xmath168 .
( the latter is a meanwhile well - known result of linear response theory @xcite which is in fact the limit we are taking here . )
since @xmath168 is constant at fixed signal frequency @xmath11 , a coupling - induced increase of @xmath165 will lead to an increase in @xmath162 , too .
the improvement of the snr under coupling is thus preserved with external coherent noise . from ( [ noisysnr2 ] )
one can easily show that the maximum enhancement of @xmath162 compared to the uncoupled chain reaches the previously found factor four for vanishing input noise , only .
in general , this factor is smaller approaching 1 , if the input noise is so strong that the input snr and hence @xmath162 go to zero .    at the end
we shall make a brief excursion to the collective response of the chain , which was mentioned in the introduction .
we consider the following sum which involves the states of @xmath0 elements :    @xmath169    it is not difficult to show that the spa of this new quantity is simply @xmath170 times the previously studied spa ( [ spa ] ) . for the snr , however , this new situation is completely different .
previously the continuous part of the spectrum was determined from the auto - correlation function @xmath171 of an element alone ( cf .
( [ correl1 ] ) ) .
now this spectrum will involve cross - correlation contributions @xmath172 of different elements @xmath173 , too .
the new contributions to the continuous part of the spectrum change the qualitative behaviour of the snr under coupling .
we shall demonstrate this in the limit @xmath174 .
for this case glauber calculated the unperturbed spectrum @xcite which reads in its one - sided version    @xmath175    inserting this spectrum into ( [ snr ] ) , multiplied by @xmath170 due to the mentioned increase of the spa , one finds for the new snr per element @xmath176    which , of course , again reduces to the simplified mcnamara - wiesenfeld snr ( [ mwsnr ] ) for vanishing coupling .
clearly , the snr is now a decreasing function of the coupling parameter and the sign of the latter is no longer important .
hence , there is a drastic difference in the snr of single and collective response with respect to coupling .
of course , we have so far merely investigated the limits of ( [ global0 ] ) , @xmath177 and @xmath174 , respectively .
we expect that the snr of the collective response of only a few elements still increases under coupling @xcite .
a detailed analysis of the present model with respect to these collective effects will be the subject of a further investigation .
in this paper we extended the two - state model of stochastic resonance introduced by mcnamara and wiesenfeld to a chain of infinitely many coupled two - state elements which are periodically modulated .
the interaction of the elements was chosen in a way that the chain evolves according to glauber s stochastic ising model .
in analogy to the work of mcnamara and wiesenfeld and based glauber s results analytical expressions for the signal - to - noise ratio ( snr ) and the spectral power amplification ( spa ) have been obtained in the limit of weak modulations . here
both quantities refer to the response of a single element as part of the chain .    instead of approximating the dynamics of a particular coupled system by the chain model
, we used the latter to build a prototype system which hopefully captures the essential features of an entire class of coupled stochastic resonators .
to this end additional assumptions were made on the dependence of the model parameters on some noise intensity and signal amplitude , respectively .
the prototype system was used to study the effect of the coupling on the response of a single resonator embedded into the chain .
the results show that array - enhanced stochastic resonance occurs for ferromagnetic - type coupling in spa and snr .
the qualitative features of the effects reproduce those previously found in coupled stochastic resonators . the simplicity of the chain model allowed for a detailed analytical investigation of the spa .
for the snr the observed effects have been confirmed analytically for the first time .
in addition , it was found that the improvement of the single - resonator response compared to the response of the uncoupled resonator still occurs , if the signal is embedded into noise .
for the snr this improvement was shown to be limited by a factor four which is reached for vanishing input noise , only .
a brief excursion into the collective response of @xmath0 resonators , on the other hand , disclosed that coupling can not improve the snr , if @xmath0 is very large
.    a closer look at the mechanisms behind the effects revealed that in the present model an improvement of the stochastic resonators under coupling is associated with the regime of quasistatic response .
since the model studied here possesses an upper bound to the time scale of its dynamics , the desired improvement is essentially restricted to sufficiently slow signals . for the spa
this improvement is based on a stronger tendency of the two - state elements to align in parallel , if signal and ferromagnetic coupling act together , compared to this tendency caused by the signal alone .
the reason why the improvement of the spa is lost , if the coupling is too strong , was found to be the slow - down of the system dynamics under coupling : it simply prevents the resonators from responding quasistatically
.          for review - type articles on stochastic resonance see , e.g. : + p. jung , phys .
* 234 * , 175 ( 1993 )
. + f. moss , d. pierson , d. ogorman , int .
4*(6 ) , 1383 ( 1994 ) .
+ k. wiesenfeld , f. moss , nature * 373 * , 33 ( 1995 ) .
+ a. bulsara and l. gammaitoni , physics today , march 1996 , 39 ( 1996 ) .
+ l. gammaitoni , p. hnggi , p. jung , and f. marchesoni , rev .
* 70 * , 223 ( 1998 ) + a bibliography on stochastic resonance can be found at http://www.pg.infn.it / sr/.    j. k. douglas , l. wilkens , e. pantazelou , and f. moss , nature * 365 * , 337 ( 1993 ) .
+ s. m. besrukov and i. vodyanoy , nature * 378 * , 362 ( 1995 ) .
+ j. e. levin and j. p. miller , nature * 380 * , 165 ( 1996 ) .
+ p. cordo , j. t. ingils , s. verschueren , j. j. collins , d. m. merfeld , s. rosenblum , s. buckley , and f. moss , nature * 383 * , 769 ( 1996 )
. + j. j. collins , t. t. imhoff , and p. grigg nature * 383 * , 770 ( 1996 ) .
r. rouse , siyuan han , and j. e. lukens , appl .
phys . lett . * 66 * , 108 ( 1995 ) .
+ a. d. hibbs , a. l. singsaas , e. w. jacobs , a. r. bulsara , j. j. bekkedahl , and f. moss , j. appl .
* 77 * , 2582 ( 1995 ) .
+ a. d. hibbs and b. r. whitecotton , in _ applied nonlinear dynamics and stochastic systems near the millennium _
, ed . by j. b. kadtke and a. bulsara , ajp conference proceedings , vol .
411 , p.227 ( aip , woodbury , new york , 1997 ) .
a. bulsara , g. schmera , phys . rev .
* e47 * , 3734 ( 1993 ) .
+ a. bulsara , a. maren , and g. schmera , biol . cybern . * 70 * , 145 ( 1993 ) .
+ m. e. inchiosa , a. r. bulsara : phys . lett .
* a200 * , 283 ( 1995 ) and phys . rev .
* e52 * , 327 ( 1995 ) .