let a memoryless additive white gaussian noise ( awgn ) channel be described by @xmath0 , where @xmath1 is independent of @xmath2 .
if the input is imposed an average power constraint given by @xmath3 , the input distribution maximizing the mutual information is gaussian .
this is due to the fact that under second moment constraint , the gaussian distribution maximizes the entropy , hence @xmath4 on the other hand , if we use a gaussian input distribution , i.e. , @xmath5 , the worst noise that can occur , i.e. , the noise minimizing the mutual information , among noises with bounded second moment , is again gaussian distributed .
this can be shown by using the _ entropy power inequality _
( epi ) , cf .
@xcite , which reduces in this setting to @xmath6 and implies @xmath7 hence , in the single - user setting , when optimizing the mutual information as above , a gaussian input is the best input for a gaussian noise and a gaussian noise is the worst noise for a gaussian input .
this provides a game equilibrium between user and nature , as defined in @xcite , p. 263 . with these results ,
many problems in information theory dealing with gaussian noise can be solved . however , in gaussian networks , that is , in multi - user information theory problems where the external noise is gaussian distributed , several new phenomena make the search for the optimal input ensemble more complex . besides for some specific cases of gaussian networks , we still do not know how interference should be treated in general .
let us consider two users interfering on each other in addition to suffering from gaussian external noise and say that the receivers treat interference as noise .
then , if the first user has drawn its code from a gaussian ensemble , the second user faces a frustration phenomenon : using a gaussian ensemble maximizes its mutual information but minimizes the mutual information of the first user .
it is an open problem to find the optimal input distributions for this problem .
this is one illustration of the complications appearing in the network setting .
another example is regarding the treatment of the fading . over a single - user awgn channel ,
whether the fading is deterministic or random , but known at the receiver , does not affect the optimal input distribution . from
, it is clear that maximizing @xmath8 or @xmath9 under an average power constraint is achieved by a gaussian input .
however , the situation is different if we consider a gaussian broadcast channel ( bc ) . when there is a deterministic fading , using and , the optimal input distribution can be shown to be gaussian .
however , it has been an open problem to show whether gaussian inputs are optimal or not for a gaussian bc with a random fading known at the receiver , even if the fading is such that it is a degraded bc .
a reason for these open questions in the network information theoretic framework , is that gaussian ensembles are roughly the only ensembles that can be analyzed over gaussian networks , as non - gaussian ensembles have left most problems in an intractable form . in this paper ,
a novel technique is developed to analyze a class of non - gaussian input distributions over gaussian noise channels .
this technique is efficient to analyze the _ competitive situations _ occurring in the network problems described below .
it allows in particular to find certain non - gaussian ensembles that outperform gaussian ones on a gaussian bc with coherent fading channel , a two user interference channel , and it allows to disprove the strong shamai - laroia conjecture on the gaussian intersymbol interference channel .
this tool provides a new insight on gaussian networks and confirms that non - gaussian ensembles do have a role to play in these networks .
we now introduce with more details the notion of competitive situations .        consider a degraded gaussian bc with coherent memoryless fading , where the fading is indeed the same for both receivers , i.e. @xmath10 but @xmath11 and @xmath12 , with @xmath13 .
the input @xmath2 is imposed a power constraint denoted by @xmath14 . because the fading is coherent , each receiver also knows the realization of @xmath15 , at each channel use .
the fading and the noises are memoryless ( iid ) processes .
since this is a degraded broadcast channel , the capacity region is given by all rate pairs @xmath16 with @xmath17 .
the optimal input distributions , i.e. , the distributions of @xmath18 achieving the capacity region boundary , are given by the following optimization , where @xmath19 , @xmath20 note that the objective function in the above maximization is given by @xmath21 now , each term in this expression is individually maximized by a gaussian distribution for @xmath22 and @xmath2 , but these terms are combined with different signs , so there is a competitive situation and the maximizer is not obvious .
when @xmath23 , one can show that gaussian distributions are optimal .
also , if @xmath15 is compactly supported , and if @xmath24 is small enough as to make the support of @xmath15 and @xmath25 non overlapping , the optimal distribution of @xmath18 is jointly gaussian ( cf .
however , in general the optimal distribution is unknown .
we do not know if it because we need more theorems , or if it is really that with fading , non - gaussian codes can actually perform better than the gaussian ones .
we consider the symmetric memoryless interference channel ( ic ) with two users and white gaussian noise .
the average power is denoted by @xmath14 , the interference coefficients by @xmath26 , and the respective noise by @xmath27 and @xmath28 ( independent standard gaussian ) .
we define the following expression @xmath29 where @xmath30 and @xmath31 are independent random vectors of dimension @xmath32 with a covariance having a trace bounded by @xmath33 and @xmath34 , @xmath35 , are iid standard gaussian . for any dimension @xmath32 and
any distributions of @xmath30 and @xmath31 , @xmath36 is a lower bound to the sum - capacity .
moreover , it is tight by taking @xmath32 arbitrarily large and @xmath30 and @xmath31 maximizing .
now , a similar competitive situation as for the fading broadcast problem takes place : gaussian distributions maximize each entropy term , but these terms are combined with different signs .
would we then prefer to take @xmath37 and @xmath38 gaussian or not ?
this should depend on the value of @xmath26 . if @xmath39 , we have two parallel awgn channels with no interference , and gaussian inputs are optimal
. we can then expect that this might still hold for small values of @xmath26 .
it has been proved recently in @xcite , that the sum - capacity is achieved by treating interference as noise and with iid gaussian inputs , as long as @xmath40 .
hence , in this regime , the iid gaussian distribution maximizes for any @xmath32 .
but if @xmath26 is above that threshold and below @xmath41 , the problem is open .
let us now review the notion of `` treating interference as noise '' . for each user
, we say that the decoder is _ treating interference as noise _ , if it does not require the knowledge of the other user s code book . however , we allow such decoders to have the knowledge of the distribution , under which the other user s code book may be drawn .
this is for example necessary to construct a sum - capacity achieving code in @xcite , where the decoder of each user treats interference as noise but uses the fact that the other user s code book is drawn from an iid gaussian distribution .
but , if we allow this distribution to be of arbitrarily large dimension @xmath32 in our definition of treating interference as noise , we can get a misleading definition .
indeed , no matter what @xmath26 is , if we take @xmath32 large enough and a distribution of @xmath30 , @xmath31 maximizing , we can achieve rates arbitrarily close to the sum - capacity , yet , formally treating interference as noise .
the problem is that the maximizing distributions in may not be iid for an arbitrary @xmath26 , and knowing it at the receiver can be as much information as knowing the other user s code book ( for example , if the distribution is the uniform distribution over a code book of small error probability ) .
hence , one has to be careful when taking @xmath32 large . in this paper
, we will only work with situations that are not ambiguous with respect to our definition of treating interference as noise .
it is indeed an interesting problem to discuss what kind of @xmath32-dimensional distributions would capture the meaning of treating interference as noise that we want .
this also points out that studying the maximizers of relates to studying the concept of treating interference as noise or information .
since for any chosen distributions of the inputs we can achieve , the maximizers of must have a different structure when @xmath26 grows .
for @xmath26 small enough , iid gaussian are maximizing distributions , but for @xmath42 , since we do not want to treat interference as noise , the maximizing distributions must have a `` heavy structure '' , whose characterization requires as much information as giving the entire code book .
this underlines that an encoder can be drawn from a distribution which does not maximize for any value of @xmath32 , but yet , a decoder may exist in order to have a capacity achieving code .
this happens if @xmath43 , iid gaussian inputs will achieve the sum - capacity if the receiver decodes the message of both users ( one can show that the problem is equivalent to having two mac s ) .
however , if @xmath42 , the iid gaussian distribution does not maximize @xmath44 ( for the dimension 1 , hence for arbitrary dimensions ) .    in any cases
, if the gaussian distribution does not maximize for the dimension 1 , it means that iid gaussian inputs and treating interference as noise is not capacity achieving , since a code which treats interference as noise and whose encoder is drawn from a distribution can be capacity achieving only if the encoder is drawn from a distribution maximizing .
hence , understanding better how to resolve the competitive situation of optimizing is a consequent problem for the interference channel .
[ sl ] let @xmath45 , @xmath46 and @xmath1 ( independent of @xmath47 ) . for all @xmath48 i.i.d . with mean 0 and variance @xmath14
, we have @xmath49    this conjecture has been brought to our attention by shlomo shamai ( shitz ) , who referred to the strong conjecture for a slightly more general statement , where an arbitrary memory for the interference term is allowed , i.e. , where @xmath50 stands for @xmath51 .
the strong conjecture then claims that picking all @xmath52 s gaussian gives a minimizer .
however , we will show that even for the memory one case , the conjecture does not hold .
the weak conjecture , also referred to as the shamai - laroia conjecture , corresponds to a specific choice of the @xmath53 s , which arises when using an mmse decision feedback equalizer on a gaussian noise isi channel , cf .
this conjecture is investigated in a work in progress .
+ there are many other examples in network information theory where such competitive situations occur . our goal in this paper is to explore the degree of freedom provided by non - gaussian input distributions .
we show that the neighborhood of gaussian distributions can be parametrized in a specific way , as to simplify greatly the computations arising in competitive situations
. we will be able to precisely quantify how much a certain amount of non - gaussianness , which we will characterize by means of the hermite polynomials , affects or helps us in maximizing the competitive entropic functional of previously mentioned problems .
for the fading bc problem problem described in [ bc ] , we want to determine if / when the distribution of @xmath18 maximizing is gaussian or not .
[ [ ic ] ] ic ~~    for the interference channel problem described in [ ic ] , we know from @xcite that treating interference as noise and using iid gaussian inputs is optimal when @xmath54 .
we question when this coding scheme is no longer optimal .
more generally , we want to analyze the maximizers of .
we distinguish the implication of such a threshold in both the synchronized and asynchronized users setting , as there will be an interesting distinction between these two cases .
we recall how the synch and asynch settings are defined here . in the synch setting
, each user of the ic sends their code words of a common block length @xmath55 simultaneously , i.e. , at time 1 , they both send the first component of their code word , at time 2 the second component , etc . in the asynch
setting , each user is still using code words of the same block length @xmath55 , however , there might be a shift between the time at which the first and second users start sending their code words .
we denote this shift by @xmath56 , and assume w.l.o.g .
that @xmath57 . in the totally asynch
setting , we assume that @xmath56 is drawn uniformly at random within @xmath58
. we may also distinguish the cases where @xmath56 is not known at the transmitter but at the receiver , and when @xmath56 is not known at both .
note that if iid input distributions are used to draw the code books , and interference is treated as noise , whether the users are synch or asynch is not affecting the rate achievability and @xmath38 can still be defined for the totally asynch ic ] . however ,
if the users want to time - share over the channel uses , such as to fully avoid their interference , they will need synchronization .
[ ts ] time sharing over a block length @xmath55 ( assumed to be even ) with gaussian inputs refers to using @xmath37 gaussian with covariance @xmath59 and @xmath38 gaussian with covariance @xmath60 , where @xmath61 is a diagonal matrix with @xmath62 1 s and 0 s , and @xmath63 flips the 1 s and 0 s on the diagonal .
we want to determine whether conjecture [ sl ] holds or not .
our more general goal is to understand better the problem posed by any competitive situations . for this purpose ,
we formulate the following mathematical problem
.    we start by changing the notation and rewrite and as @xmath64 where @xmath65 denotes the gaussian density with zero mean and variance @xmath14 , and the functions @xmath66 are density functions on @xmath67 , i.e. , positive functions integrating to 1 , and having a well - defined entropy and second moment @xmath68 .
+ we consider the local geometry by looking at densities of the form @xmath69 where @xmath70 satisfies @xmath71 with these two constraints on @xmath72 , @xmath73 is a valid density for @xmath74 sufficiently small .
it is a perturbed gaussian density , in a `` direction '' @xmath72 .
observe that , @xmath75    we are now interested in analyzing how these perturbations affect the output distributions through an awgn channel . note that , if the input distribution is a gaussian @xmath65 perturbed in the direction @xmath72 , the output is a gaussian @xmath76 perturbed in the direction @xmath77 since @xmath78 _ convention : _ @xmath79 refers to @xmath80 , i.e. , the multiplicative operator precedes the convolution one .
+ for simplicity , let us assume in the following that the function @xmath72 is a polynomial satisfying , .
[ appro ] we have @xmath81 where @xmath82    moreover , note that for any density @xmath66 , if @xmath83 and @xmath84 , we have @xmath85 hence , the extremal entropic results of and are locally expressed as @xmath86 where 0 denotes here the zero function . if is obvious , requires a proof which will be done in section [ proofs ] .
let us define the following mapping , @xmath87 where @xmath88 denotes the space of real functions having a finite @xmath89 norm .
this linear mapping gives , for a given perturbed direction @xmath72 of a gaussian input @xmath65 , the resulting perturbed direction of the output through additive gaussian noise @xmath90 .
the norm of each direction in their respective spaces , i.e. , in @xmath88 and @xmath91 , gives how far from the gaussian distribution these perturbations are ( up to a scaling factor ) .
note that if @xmath72 satisfies - , so does @xmath92 for the measure @xmath76 .
the result in ( worst noise case ) tells us that this mapping is a contraction , but for our goal , what would be helpful is a spectral analysis of this operator , to allow more quantitative results than the extreme - case results of and .    in order to do so
, one can express @xmath93 as an operator defined and valued in the same space , namely @xmath94 with the lebesgue measure @xmath95 , which is done by inserting the gaussian measure in the operator argument .
we then proceed to a singular function / value analysis .
formally , let @xmath96 which gives @xmath97 , and let @xmath98 which gives @xmath99 denoting by @xmath100 the adjoint operator of @xmath101 , we want to find the singular functions of @xmath101 , i.e. , the eigenfunctions @xmath102 of @xmath103 : @xmath104
the following theorem gives the singular functions and values of the operator @xmath101 defined in previous section .
[ hermite1 ] @xmath105 holds for each pair @xmath106 } , \left(\frac{p}{p+v}\right)^k)\}_{k \geq 0},\ ] ] where @xmath107}(x)=\frac{1}{\sqrt{k ! } } h_k(x/\sqrt{p})\ ] ] @xmath108    the polynomials @xmath109}$ ] are the normalized hermite polynomials ( for a gaussian distribution having variance @xmath14 ) and @xmath110}$ ] are called the hermite functions . for any @xmath111 , @xmath112}\}_{k\geq 0}$ ] is an orthonormal basis of @xmath88 , this can be found in @xcite .
one can check that @xmath113 , respectively @xmath114 perturb a gaussian distribution into another gaussian distribution , with a different first moment , respectively second moment . for @xmath115 ,
the @xmath116 perturbations are not modifying the first two moments and are moving away from gaussian distributions . since @xmath117}=1 $ ] , the orthogonality property implies that @xmath109}$ ] satisfies for any @xmath118 .
however , it is formally only for even values of @xmath119 that is verified ( although we will see in section [ proofs ] that essentially any @xmath119 can be considered in our problems ) .
the following result contains the property of hermite polynomials mostly used in our problems , and expresses theorem [ hermite1 ] with the gaussian measures .
the following result contains the property of hermite polynomials mostly used in our problems , and expresses proposition [ hermite1 ] with the gaussian measures .
[ hermite2 ] @xmath120}= \frac{g_p h_k^{[p ] } \star g_v}{g_{p+v}}=\left(\frac{p}{p+v}\right)^{k/2 } h_k^{[p+v ] } , \label{it } \\ & { \gamma^{(-)}}h_k^{[p+v]}= h_k^{[p+v ] } \star g_v = \left(\frac{p}{p+v}\right)^{k/2 } h_k^{[p]}. \end{aligned}\ ] ]    last theorem implies theorem [ hermite1 ] , since @xmath121 for @xmath122    _ comment : _ the results that we have just derived are related to properties of the ornstein - uhlenheck process . _
summary : _ in words , we just saw that @xmath123 is an eigenfunction of the input / output perturbation operator @xmath93 , in the sense that @xmath124 } = \left(\frac{p}{p+v } \right)^{k/2 } h_k^{[p+v]}$ ] .
hence , over an additive gaussian noise channel @xmath90 , if we perturb the input @xmath65 in the direction @xmath109}$ ] by an amount @xmath74 , we will perturb the output @xmath76 in the direction @xmath125}$ ] by an amount @xmath126 .
such a perturbation in @xmath116 implies that the output entropy is reduced ( compared to not perturbing ) by @xmath127 ( if @xmath128 ) .
the following result states that the capacity region of a degraded fading bc with gaussian noise is not achieved by a gaussian superposition code in general .    [ ce ] let @xmath129 with @xmath2 such that @xmath3 , @xmath130 , @xmath131 , @xmath132 and @xmath133 mutually independent .
there exists a fading distribution and a value of @xmath24 for which the capacity achieving input distribution is non - gaussian .
more precisely , let @xmath22 be any auxiliary random variable , with @xmath134 .
then , there exists @xmath135 , a distribution of @xmath15 and @xmath136 such that @xmath137 is maximized by a non jointly gaussian distribution .    in the proof
, we present a counter - example to gaussian being optimal for @xmath15 binary . in order to defeat gaussian distributions ,
we construct input distributions using the hermite coordinates .
the proof also gives a condition on the fading distribution and the noise variance @xmath24 for which a non - gaussian distribution strictly improves on the gaussian one .
let @xmath138,\ ] ] where @xmath139 , @xmath140 and @xmath141 , with @xmath142 defined in below ( as explain in section [ form ] , @xmath142 is a formal modification of @xmath116 to ensure the positivity of the perturbed densities ) .
+ in other words , @xmath143 represents the gain ( positive or negative ) of using @xmath37 perturbed along @xmath116 and @xmath38 perturbed along @xmath144 with respect to using gaussian distributions . note that the distributions we chose for @xmath37 and @xmath38 are not the most general ones , as we could have chosen arbitrary directions spanned by the hermite basis to perturb the gaussian densities
however , as explained in the proof of the theorem [ calc ] , this choice is sufficient for our purpose .
[ calc ] we have for @xmath145 @xmath146^k - \frac{(a^k-1)^2}{(pa^2+p+1)^k}.\ ] ]    for any fixed @xmath14 , the function @xmath147 has a unique positive root , below which it is negative and above which it is positive .
[ inter ] treating interference as noise with iid gaussian inputs does not achieve the sum - capacity of the symmetric ic ( synch or asynch ) and is outperformed by @xmath148 and @xmath149 , if @xmath150    this theorem is a direct consequence of theorem [ calc ] .
[ f2 ] for the symmetric synch ic , time sharing improves on treating interference as noise with iid gaussian distribution if @xmath151    we now introduce the following definition .
blind time sharing over a block length @xmath55 ( assumed to be even ) between two users , refers to sending non - zero power symbols only at the instances marked with a 1 in @xmath152 for the first user , and zero power symbols only at the instances marked with a 1 in @xmath153 for the second user .
[ b2 ] for the symmetric totally asynch ic , if the receivers ( but not transmitters ) know the asynchronization delay , blind time sharing improves on treating interference as noise with iid gaussian distributions if @xmath154 where @xmath155 if the receivers do not know the asynchronization delay , blind time sharing can not improve on treating interference as noise with iid gaussian distributions if @xmath156 .    _
how to read these results : _ we have four thresholds to keep track of :    * @xmath157 is when @xmath158 . if @xmath159 , we know from @xcite that iid gaussian inputs and treating interference as noise is sum - capacity achieving . *
@xmath160 is when @xmath161 . if @xmath162 , we know from prop .
[ f2 ] that , if synchronization is permitted , time sharing improves on treating interference as noise with iid gaussian inputs .
this regime matches with the so - called moderate regime defined in @xcite .
* @xmath163 is when @xmath164 . if @xmath165 , we know from prop .
[ inter ] that treating interference as noise with iid non - gaussian distributions ( opposites in @xmath166 ) improves on the iid gaussian ones .
* @xmath167 is when @xmath168 . if @xmath169 , we know from prop . [ b2 ] that , even if the users are totally asynchronized , but if the receivers know the asynchronization delay , blind time sharing improves on treating interference as noise with iid gaussian inputs .
if the receivers do not know the delay , the threshold can only appear for larger values of @xmath26 .
the question is now , how are these thresholds ranked .
it turns out that @xmath170 and if @xmath171 , the above inequality reads as @xmath172 .
this implies the following for a decoder that treats interference as noise . since @xmath173
, it is first better to time share than using non - gaussian distributions along @xmath166 .
but this is useful only if time - sharing is permitted , i.e. , for the synch ic .
however , for the asynchronized ic , since @xmath174 , we are better off using the non - gaussian distributions along @xmath166 before a gaussian input scheme , even with blind time - sharing , and even if the receiver could know the delay .
we notice that there is still a gap between @xmath157 and @xmath160 , and we can not say if , in this range , iid gaussian inputs are still optimal , or if another class of non - gaussian inputs ( far away from gaussians ) can outperform them . in @xcite , another technique ( which is related to ours but not equivalent ) is used to find regimes where non - gaussian inputs can improve on gaussian ones on the same problem that we consider here .
the threshold found in @xcite is equal to 0.925 for @xmath171 , which is looser than the value of 0.680 found here .    finally , the following interesting and curious fact has also been noticed . in theorem [ calc ] , we require @xmath145 .
nevertheless , if we plug @xmath175 in the right hand side of theorem [ calc ] and ask for this expressions to be positive , we precisely get @xmath176 , i.e. , the complement range delimited by @xmath157 . however , the right hand side of theorem [ calc ] for @xmath175 is not equal to @xmath177 ( this is explained in more details in the proof of theorem [ calc ] ) .
indeed , it would not make sense that moving along @xmath113 , which changes the mean with a fixed second moment within gaussians , would allow us to improve on the iid gaussian scheme .
yet , getting to the exact same condition , when working on the problem of improving on the iid gaussian scheme , seems to be a strange coincidence .
we show in section [ proofs ] that conjecture [ sl ] does not hold . we provide counter - examples to the conjecture , pointing out that the range of @xmath178 for which the conjecture does not hold increases with snr .
the hermite polynomial corresponding to @xmath179 is @xmath117}=1 $ ] and is clearly not a valid direction as it violates . using the orthogonality property of the hermite basis and since @xmath117}=1 $ ]
, we conclude that @xmath109}$ ] satisfies for any @xmath118 .
however , it is only for @xmath119 even that @xmath109}$ ] satisfies . on the other hand , for any @xmath180
, we have that @xmath109 } + \delta h_{4k}^{[p]}$ ] satisfies , whether @xmath119 is even or not ( we chose @xmath181 instead of @xmath182 for reasons that will become clear later ) .
now , if we consider the direction @xmath183}$ ] , is not satisfied for both @xmath119 even and odd . but
again , for any @xmath180 , we have that @xmath184 } + \delta h_{4k}^{[p]}$ ] satisfies .
hence , in order to ensure , we will often work in the proofs with @xmath185 } + \delta h_{4k}^{[p]}$ ] , although it will essentially allow us to reach the performance achieved by any @xmath186}$ ] ( odd or even ) , since we will then take @xmath187 arbitrarily small and use continuity arguments .
_ convention : _ we drop the variance upper script in the hermite terms whenever a gaussian density with specified variance is perturbed , i.e. , the density @xmath188 always denotes @xmath189 } ) $ ] , and @xmath190 always denotes @xmath191 } $ ] , no matter what @xmath14 is .
same treatment is done for @xmath89 and @xmath192 .
now , in order to evaluate the entropy of a perturbation , i.e. , @xmath193 , we can express it as the entropy of @xmath194 minus the divergence gap , as in , and then use lemma [ appro ] for the approximation .
but this is correct _ if _ @xmath195 has the same first two moments as @xmath65 .
hence , if @xmath72 contains only @xmath116 s with @xmath128 , the previous argument can be used .
but if @xmath72 contains @xmath113 and/or @xmath114 terms , the situation can be different .
next lemma describes this .
[ corr ] let @xmath196 and @xmath197 } =     \begin{cases }        b(h_k^{[p ] } + \delta h_{4k}^{[p ] } ) , & \text{if } b \geq 0 ,   \\
b(h_k^{[p ] } - \delta h_{4k}^{[p ] } ) , & \text{if } b < 0 .
\label{tild }     \end{cases}\end{aligned}\ ] ] we have for any @xmath198 , @xmath199 , @xmath200 @xmath201    finally , when we convolve two perturbed gaussian distributions , we get @xmath202 + { \varepsilon}^2 g_a h_j \star g_b h_k.$ ] we already know from theorem [ hermite2 ] how to describe the terms in @xmath74 , what we still need is to describe the term in @xmath203 .
we have the following .    [ new ] we have @xmath204 } \star g_b h_l^{[b ] } = c g_{a+b } h_{k+l}^{[a+b]},\ ] ] where @xmath205 is a constant depending only on @xmath206 and @xmath207 .
in particular if @xmath208 , we have @xmath209 .
we start by reviewing the proof of , as it brings interesting facts .
we then prove the main result .
_ proof of : _
+ we first assume that @xmath73 has zero mean and variance @xmath14 . using the hermite basis , we express @xmath72 as @xmath210}$ ] ( @xmath72 must have such an expansion , since it must have a finite @xmath88 norm , to make sense of the original expressions ) . using , we can then express as @xmath211 which is clearly negative .
hence , we have proved that @xmath212 and is maximized by taking @xmath213 .
note that we can get tighter bounds than the one in previous inequality , indeed the tightest , holding for @xmath166 , is given by @xmath214 ( this clearly holds if written as a series like in ) .
hence , locally the contraction property can be tightened , and locally , we have stronger epi s , or worst noise case .
namely , if @xmath215 , we have @xmath216 and if @xmath217 , @xmath65 is outperformed by non - gaussian distributions .
now , if we consider the constraint @xmath218 , which , in particular , allows to consider @xmath219 and @xmath220 , we get that if @xmath221 , @xmath222 and if @xmath223 , @xmath65 is outperformed by @xmath224 for some @xmath225 . it would then be interesting to study if these tighter results hold in a greater generality than for the local setting .
_ proof of theorem [ hermite2 ] : _ + we want to show @xmath226 } \star g_v = ( \frac{p}{p+v})^{k/2 } h_k^{[p ] } , \\ & g_p h_k^{[p ] } \star g_v=(\frac{p}{p+v})^{k/2 } g_{p+v } h_k^{[p+v ] } , \ ] ] which is proved by an induction on @xmath119 , using the following properties ( appell sequence and recurrence relation ) of hermite polynomials : @xmath227 } ( x)&= & \sqrt{\frac{k+1}{p } } h_{k}^{[p ] } ( x ) \\ \frac{\partial}{\partial x } \left ( g_p ( x)h_k^{[p]}(x ) \right ) & = & -\sqrt{\frac{k+1}{p } } g_p(x ) h_{k+1}^{[p]}(x ) . \ ] ]    _ proof of theorem [ ce ] : _ + we refer to as the mu - rate .
let us first consider gaussian codes , i.e. , when @xmath18 is jointly gaussian , and see what mu - rate they can achieve . without loss of generality , we can assume that @xmath228 , with @xmath22 and @xmath229 independent and gaussian , with respective variance @xmath230 and @xmath231 satisfying @xmath232 . then , becomes @xmath233 now , we pick a @xmath136 and look for the optimal power @xmath231 that must be allocated to @xmath229 in order to maximize the above expression .
we are interested in cases for which the optimal @xmath231 is not at the boundary but at an extremum of , and if the maxima is unique , the optimal @xmath231 is found by the first derivative check , which gives @xmath234 since we will look for @xmath136 , @xmath24 , with @xmath235 , previous condition can be written as @xmath236 we now check if we can improve on by moving away from the optimal jointly gaussian @xmath18 .
there are several ways to perturb @xmath18 , we consider first the following case . we keep @xmath22 and @xmath229 independent , but perturb them away from gaussian s in the following way : @xmath237 } ( u ) + \delta h_{4 } ) ) \label{udef}\\ & p_{v_{\varepsilon } } ( v ) = g_r(v ) ( 1-{\varepsilon}(h_3^{[r]}(v ) - \delta h_{4 } ) ) \label{vdef}\end{aligned}\ ] ] with @xmath238 small enough .
note that these are valid density functions and that they preserve the first two moments of @xmath22 and @xmath229 .
the reason why we add @xmath239 , is to ensure that is satisfied , but we will see that for our purpose , this can essentially be neglected . then , using lemma [ hermite2 ] , the new distribution of @xmath2 is given by @xmath240 } - { \varepsilon}\left(\frac{r}{p}\right)^{\frac{3}{2 } } h_3^{[p ] } ) + f(\delta)\end{aligned}\ ] ] where @xmath241 } +    \left(\frac{r}{p}\right)^{\frac{4}{2 } } h_4^{[p]})$ ] , which tends to zero when @xmath187 tends to zero .
now , by picking @xmath242 , we have @xmath243 hence , by taking @xmath187 arbitrarily small , the distribution of @xmath2 is arbitrarily close to the gaussian distribution with variance @xmath244 .
we now want to evaluate how these hermite perturbations perform , given that we want to maximize , i.e. , @xmath245 we wonder if , by moving away from gaussian distributions , the gain achieved for the term @xmath246 is higher than the loss suffered from the other terms .
using theorem [ hermite2 ] , lemma [ appro ] and lemma [ corr ] , we are able to precisely measure this and we get @xmath247 } ) ) \\ & = \frac{1}{2 } \log 2 \pi e ( v+r h^2 ) - \frac{{\varepsilon}^2}{2 }   \left(\frac{r h^2}{v+rh^2}\right)^{3 } + o({\varepsilon}^2 ) + o(\delta)\end{aligned}\ ] ] @xmath248 and because of @xmath249 therefore , collecting all terms , we find that for @xmath250 and @xmath251 defined in and , expression reduces to @xmath252 where @xmath253 is equal to ( and is the mu - rate obtained with gaussian inputs ) . hence , if for some distribution of @xmath15 and some @xmath24 , we have that @xmath254 when @xmath255 and @xmath231 is optimal for @xmath136 , we can take @xmath74 and @xmath187 small enough in order to make strictly larger than @xmath253 .
we have shown how , if verified , inequality leads to counter - examples of the gaussian optimality , but with similar expansions , we would also get counter - examples if the following inequality holds for any power @xmath119 instead of 3 , as long as @xmath128 .
let us summarize what we obtained : let @xmath231 be optimal for @xmath136 , which means that holds if there is only one maxima ( not at the boarder ) .
then , non - gaussian codes along hermite s strictly outperforms gaussian codes , if , for some @xmath128 , holds .
if the maxima is unique , this becomes @xmath256 where @xmath257 so we want the jensen gap of @xmath258 for the power @xmath119 to be small enough compared to the jensen gap of @xmath259 .
we now give an example of a fading distribution for which the above conditions can be verified .
let @xmath15 be binary , taking values 1 and 10 with probability half and let @xmath260 .
let @xmath261 , then for any values of @xmath244 , the maximizer of is at @xmath262 , cf .
figure [ maxval ] , which corresponds in this case to the unique value of @xmath231 for which is satisfied .     for @xmath261 , @xmath260 , @xmath263 and @xmath15 binary @xmath264 .
maxima at @xmath262.,width=287,height=249 ]    hence if @xmath244 is larger than this value of @xmath231 , there is a corresponding fading bc for which the best gaussian code splits the power on @xmath22 and @xmath229 with @xmath262 to achieve the best mu - rate with @xmath261 . to fit the counter - examples with the choice of hermite perturbations made previously ,
we pick @xmath242 . finally , for these values of @xmath136 and @xmath231 , can be verified for @xmath265 , cf .
figure [ ineq ] , and the corresponding hermite code ( along @xmath266 ) strictly outperforms any gaussian codes .    , for @xmath261 , @xmath260 , @xmath267 and @xmath15 binary @xmath264 , positive at @xmath262.,width=287,height=249 ]
note that we can consider other non - gaussian encoders , such as when @xmath22 and @xmath229 are independent with @xmath22 gaussian and @xmath229 non - gaussian along hermite s .
then , we get the following condition .
if for @xmath128 and @xmath231 optimal for @xmath136 , we have @xmath268 , \label{weak}\end{aligned}\ ] ] then gaussian encoders are not optimal .
notice that previous inequality is stronger than the one in for fixed values of the parameters .
yet , it can still be verified for valid values of the parameters and there are also codes with @xmath22 gaussian and @xmath229 non - gaussian that outperform gaussian codes for some degraded fading bcs .
+ _ proof of theorem [ calc ] : _ + let @xmath269 and let @xmath37 and @xmath38 be respectively distributed as @xmath270)$ ] and @xmath271)$ ] , where @xmath272 .
we have @xmath273 where @xmath274 are independent gaussian 0-mean and @xmath14-variance random variables .
hence , we need to evaluate the contribution of each divergence appearing in previous expression , in order to know if the perturbations are improving on the gaussian distributions .
let us first analyze @xmath275 .
the density of @xmath276 is given by @xmath277 ) \star g_{a^2 p}(1-{\varepsilon}[h_k-\delta h_{4k } ] ) \star g_1 ,   \end{aligned}\ ] ] which , from theorem [ hermite2 ] , is equal to @xmath278 \\   & -\left [ \left(\frac{a^2p}{p+a^2 p + 1 } \right)^{\frac{k}{2 } } h_k -\delta \left(\frac{a^2p}{p+a^2 p + 1 } \right)^{2k } h_{4k } \right ]   \\ & - { \varepsilon}l \ } ) \end{aligned}\ ] ] where @xmath279 \star g_{a^2 p } [ h_k-\delta h_{4k } ]   \star g_1}{g_{p+a^2p+1}}.\ ] ] note that each direction in each line of the bracket @xmath280 above , including @xmath72 , satisfy and .
using lemma [ new ] , we have @xmath281 \star g_{a^2 p+1 } [ \left(\frac{a^2 p}{a^2 p+1 } \right)^{\frac{k}{2}}h_k- \left(\frac{a^2 p}{a^2 p+1 } \right)^{2k } \delta h_{4k } ]   } { g_{p+a^2p+1 } } \notag \\ & = c_1 h_{2k}^{[p+a^2 p + 1 ] } + c_2 h_{5k}^{[p+a^2 p + 1 ] } + c_3 h_{8k}^{[p+a^2 p + 1 ] } , \label{bug }    \ ] ] where @xmath282 are constants .
therefore , the density of @xmath276 is a gaussian @xmath283 perturbed along the direction @xmath116 in the order @xmath74 and several @xmath284 with @xmath285 in the order @xmath203 ( and other directions but that have a @xmath187 order ) .
so we can use lemma [ corr ] and write @xmath286 using lemma [ appro ] , we have @xmath287 \\   & -\left [ \left(\frac{a^2p}{p+a^2 p + 1 } \right)^{\frac{k}{2 } } h_k -\delta \left(\frac{a^2p}{p+a^2 p + 1 } \right)^{2k } h_{4k } \right ] \|^2\\ & =   \frac{{\varepsilon}^2}{2 } ( 1 - a^k)^2 \left(\frac{p}{p+a^2 p + 1 } \right)^k + \frac{{\varepsilon}^2}{2 } o(\delta).\end{aligned}\ ] ] hence @xmath288 similarly , we get @xmath289 and @xmath290 + \frac{{\varepsilon}^2}{2 } o(\delta ) .\end{aligned}\ ] ] finally , we have @xmath291 and @xmath292 + { \varepsilon}^2 o(\delta).\end{aligned}\ ] ] hence , if for some @xmath293 we have @xmath294 we can improve on the iid gaussian distributions @xmath65 by using the respective hermite perturbations .
now , we could have started with @xmath37 and @xmath38 distributed as @xmath295 and @xmath296 , where @xmath142 is defined in . with similar expansions , we would then get that we can improve on the gaussian distributions if for some some @xmath297 and @xmath272 we have @xmath298 ( b_k^2+c_k^2 ) \\    & - 4   \left(\frac { a p }
{ p+a^2p+1}\right)^k     b_k c_k > 0.\end{aligned}\ ] ] but the quadratic function @xmath299 with @xmath225 , can be made positive if and only if @xmath300 , and is made so by taking @xmath301 .
hence , the initial choice we made about @xmath37 and @xmath38 is optimal .
moreover , note that for this distribution of @xmath37 and @xmath38 , we could have actually chosen @xmath302 as well . because , even if lemma [ corr ] tells us that we must use correction terms , these correction terms will cancel out when we consider the sum - rate , since @xmath303 and since the correction is in @xmath74 .
there is however another problem when using @xmath302 , which is that @xmath304 has a larger second moment than @xmath14 .
however , if we use a scheme of block length 2 , we can compensate this excess on the first channel use with the second channel use , and because of the symmetry , we can achieve the desired rate .
but this is allowed only with synchronization .
we could also have used perturbations that are mixtures of hermite s , such as @xmath305 .
we would then get mixtures of previous equations as our condition .
but in the current problem this will not be helpful .
finally , perturbing iid gaussian inputs in a independent but non i.d .
way , i.e. , to perturb different components in different hermite directions , can not improve on our scheme , from previous arguments .
the only option which is not investigated here ( but in a work in progress ) , is to perturb iid gaussian inputs in a non independent manner . finally , if we work with @xmath175 , the proof sees the following modification . in
, we now have a term in @xmath114 .
however , even if this term is in the order @xmath203 , we can no longer neglect it , since from lemma [ corr ] , a @xmath306 term in the direction comes out as a @xmath307 term in the entropy .
hence , we do not get the above condition for @xmath175 , but the one obtained by replacing @xmath308 with @xmath309 , and the condition for positivity can never be fulfilled .
+ _ proof of proposition [ f2 ] : _ + from theorem [ calc ] , we know that when treating interference as noise and when @xmath310 , it is better to use encoders drawn from the 2 dimensional distributions @xmath311 and @xmath312 , where @xmath313 , @xmath314 , @xmath315 and @xmath316 , as opposed to using gaussian distributions . but perturbations in @xmath114 are changing the second moment of the input distribution . hence , this scheme is mimicing a time - sharing in our local setting .
moreover , a direct computation also allows to show that , constraining each user to use gaussian inputs of arbitrarily block length @xmath55 , with arbitrary covariances having a trace bounded by @xmath317 , the optimal covariances are @xmath318 if @xmath319 , and otherwise , are given by a time - sharing scheme ( cf .
definition [ ts ] for the definition of a gaussian time - sharing scheme and covariance matrices ) .
+ _ proof of proposition [ b2 ] : _ + note that when using blind time - sharing , no matter what the delay in the asynchronization of each user is , the users are interfering in @xmath320 channel uses and have each a non - intefering channel in @xmath320 channel uses ( the rest of the @xmath320 channel uses are not used by any users ) . hence ,
if the receiver have the knowledge of the asynchronization delay , the following sum - rate can be achieved : @xmath321 . and
if the delay is unknown to the receivers , the previous sum - rate can surely not be improved on .
+ _ disproof of conjecture [ sl ] : _ + this proof uses similar steps as previous proofs . using lemma
, we express @xmath49 as @xmath322 we then pick @xmath323 and assume that @xmath119 is even for now .
we then have @xmath324 and @xmath325 similarly , @xmath326 and @xmath327 finally , @xmath328 and @xmath329 therefore , is given by @xmath330 and if @xmath331 for some @xmath119 even and greater than 4 , we have a counter example to the strong conjecture .
note that , using the same trick as in previous proofs , that is , perturbing along @xmath142 instead of @xmath116 , we get that if holds for any @xmath128 , we have a counter example to the strong conjecture .
defining @xmath332 , is equivalent to @xmath333 as shown in figure [ slplot ] , this can indeed happen .
an interesting observation is that the range where holds is broader when @xmath334 is larger , i.e. , when snr is smaller . indeed , when @xmath335 , which corresponds to dropping the additive noise @xmath336 , we do not get a counter - example to the conjecture .
but in the presence of gaussian noise , the conjecture does not hold for some distributions of @xmath337 .
the conjecture had been numerically checked with binary inputs at low snr , and in this regime , it could not be disproved . with the hint described above
, we checked numerically the conjecture with binary inputs at high snr , and there we found counter - examples .
we have developed a technique to analyze codes drawn from non - gaussian ensembles using the hermite polynomials .
if the performance of non - gaussian inputs is usually hard to analyze , we showed how with this tool , it reduces to the analysis of analytic power series .
this allowed us to show that gaussian inputs are in general not optimal for degraded fading gaussian bc , although they might still be optimal for many fading distributions .
for the ic problem , we found that in the asynchronous setting and when treating interference as noise , using non - gaussian code ensembles ( @xmath166 perturbations ) can strictly improve on using gaussian ones , when the interference coefficient is above a given threshold , which significantly improves on the existing threshold ( cf .
we have also recovered the threshold of the moderate regime by using @xmath114 perturbations in the synch setting , showing that this global threshold is reflected in our local setting .
we also met mysteriously in our local setting the other global threshold found in @xcite , below which treating interference as noise with iid gaussian inputs is optimal .
it is worth noting that this two global thresholds ( moderate regime and noisy interference ) are recovered with our tool from a common analytic function .
we hope to understand this better with a work in progress .
the hermite technique provides not only counter - examples to the optimality of gaussian inputs but it also gives insight on the competitive situations in gaussian network problems .
for example , in the fading bc problem , the hermite technique gives a condition on what kind of fading distributions and degradedness ( values of @xmath24 ) non - gaussian inputs must be used .
it also points out that the perturbation in @xmath166 are most effective when carried in an opposite manner for the two users , so as to make the distribution of @xmath2 close to gaussian .
finally , in a different context , local results could be `` lifted '' to corresponding global results in @xcite .
there , the localization is made with respect to the channels and not the input distribution , yet , it would be interesting to compare the local with the global behavior for the current problem too . the fact that we have observed some global results locally , as mentioned previously ,
gives hope for possible local to global extensions .
a work in progress aims to use our tool beyond the local setting , in particular , by analyzing all sub - gaussian distributions .
moreover , there are interesting connections between the results developed in this paper and the properties of the ornstein - uhlenbeck process .
indeed , some of these properties have already been used in @xcite to solve the long standing entropy monotonicity conjecture , and we are currently investigating these relations from closer .
the authors would like to thank tie liu and shlomo shamai for pointing out problems relevant to the application of the proposed tool , daniel stroock for discussions on hermite polynomials , and emre telatar for stimulating discussions and helpful comments .    1 e. abbe and l. zheng , `` linear universal decoding : a local to global geometric approach '' , _ submitted ieee trans .
inform . theory _
annapureddy and v.v .
veeravalli , `` gaussian interference network : sum capacity in the low interference regime and new outer bounds on the capacity region '' , _ submitted ieee trans .
inform . theory _
, 2008 .
m. costa , `` on the gaussian interference channel , '' ieee trans .
theory , vol .
it-31 , pp .
607 - 615 , sept .
1985 . t. m. cover , `` comments on broadcast channel '' , _ ieee trans . on inform .
44 , n. 6 , pp .
2524 - 2530 , oct . 1998
. t. m. cover and j. a. thomas , `` elements of information theory '' , john wiley & sons , new york , ny , 1991 .
stam , `` some inequalities satisfied by the quantities of information of fisher and shannon '' , _ information and control _ , 2 : 101 - 112 , 1959 .
d. w. stroock , `` a concise introduction to the theory of integration '' , third edition , birkhuser , 1999 .


from theoretical and experimental observations it is known that reaction - diffusion systems play an important role in the study of generic spatiotemporal behavior of non - equilibrium systems .
usually such models admit main contributions related to both local dynamics ( chemical reactions type of birth - and - death processes ) and mass transport .
novel experimental methods , such as field ion microscopy , scanning tunneling microscopy can be used to monitor chemical reactions on the metal surfaces with atomic resolution .    in adsorption - desorption processes
when material can be deposited from the gaseous phase such experimental methods allow one to investigate formation of clusters or islands of adsorbed molecules / atoms @xcite . such islands can have linear size of nanometer range @xcite .
in refs.@xcite it was experimentally shown that nano - patterns on solid surface and nano - islands in adsorbed mono - atomic layers can be organized .
it was found that lateral interactions between adsorbate play an important role in processes of pattern formation on metal surfaces .
the adsorbate presence can modify the local crystallographic structures of the substrate s surface layer producing long range interactions between adsorbed atoms and their clusters ( see for example , refs.@xcite ) .
it was observed experimentally that nanometer - sized vacancy islands can be organized in a perfect triangular lattice when a single monolayer of ag was exposed on ru(0001 ) surface at room temperature @xcite .
the origin of such self - organization of vacancy islands is an elastic interaction between edges of neighboring islands @xcite .
such systems can be widely used in electronic devices in industrial applications . despite theoretical studies of nano - patterns
were widely reported ( see refs.@xcite and citations therein ) no detailed descriptions of elastic effects leading to formation vacancy / adsorbate islands lattice were presented .
it is known that such kind of processes can be studied by multiscale modeling including : _ ab - initio _ calculations allowing one to set parameters of primitive cell at zero temperature and some energetic parameters ; molecular dynamics simulations monitoring microscopic processes of interactions of atoms ; monte - carlo approach studying physics on bigger time scales .
unfortunately the corresponding computational time remains excessive and presented methods can be used to describe the system only on micron scales . to study pattern - forming transitions one can use mesoscopic approach to modelize self - organization processes on large sizes .
it can be generalized by taking into account corresponding fluctuating sources representing an influence of a microscopic level or fluctuations of environment @xcite .
properties of pattern formation in systems of adsorption - desorption type with dissipative dynamics were studied previously @xcite .
pattern formation in stochastic systems with internal multiplicative noise obeying fluctuation - dissipation relation was discussed in refs.@xcite .
it was shown that internal multiplicative noise can sustain stationary patterns of nano - size range ; the corresponding patterning processes are well related to non - equilibrium transitions induced by internal multiplicative noise ( entropy driven phase transitions @xcite ) .
recently modifications of this mesoscopic approach were proposed , where in the framework of continuous model one can describe processes at microscopic and diffusion scales .
this approach is known in literature as modified phase field crystals method @xcite .
it can be used to study crystal growth , solidification processes , motion of defects , phase transitions in crystalline systems ( having special kind of order ) .
its main advantage is in introduction of a crystalline order with elastic interactions leading to formation of crystal lattice and of inertia effects related to a finite propagation speed @xcite .
inertial contribution is essential at early stages in the system dynamics . on the diffusion time scale
one can observe patterning with a formation of crystalline lattice regarding to the hook s law .
formally , the elastic effects are modelized by expansion of a pair interatomic potential in the fourier space up to the fourth order in wave - numbers @xcite . as a result
the spatial interactions are described by the swift - hohenberg operator @xmath0 , @xmath1 , @xmath2 is a constant related to period of spatial patterns .
influence of both internal and external stochastic sources onto pattern - forming transitions in crystalline systems were discussed in refs.@xcite .
it was shown that in such kind of systems due to finite propagation speed pattern selection processes are possible ; the external noise can principally change topology of patterns realized in a pure deterministic case .    for class of systems with adsorption - desorption and diffusion
processes the typical deterministic equation , which will interest us is of the form @xmath3 where @xmath4 is the local coverage at surface defined as the quotient between the number of adsorbed molecules / atoms in a cell of the surface and the fixed number of available sites in each cell , @xmath5 .
the term @xmath6 stands for local dynamics and describes birth - and - death or adsorption - desorption processes ; the flux @xmath7 represents the mass transport
.    most of the theoretical studies deal with the standard fick law @xmath8 , @xmath9 is the diffusion constant .
it is known that at @xmath10 the ordinary diffusion equation @xmath11 has the unrealistic feature of infinitely fast ( infinite ) propagation .
to avoid such unphysical situation the diffusion flux can be generalized by taking into account memory ( non - equilibrium ) effects @xmath12 described by the memory kernel @xmath13 .
in such a case in the limit @xmath14 one gets the finite propagation speed @xmath15 . at @xmath16
the assymptotic @xmath17 leads to the fick law @xmath8 with infinite propagation .
as far as real systems ( molecules , atoms ) have finite propagation speed one should use eq.([eqj ] ) or the equivalent maxwell - cattaneo equation : @xmath18 @xcite . equations ( [ eqx],[eqj ] ) can be combined into one equation of the form @xmath19 where prime denotes derivative with respect to the argument . as was pointed out in ref.@xcite
this equation has some restrictions : ( i ) it typically does not preserve positivity of the solution @xmath20 ; ( ii ) the damping coefficient must be positive , i.e. , @xmath21 .    in this study we are aimed to describe pattern formation processes in system given by eq.([hypx ] ) generalized by introducing multiplicative noise satisfying fluctuation - dissipation relation .
we use suppositions of the phase field theory to construct effective model describing formation of nano - patterns with crystalline order and pattern selection processes .
it will be shown that formation of stationary patterns in stochastic systems is well related to noise - induced effects , namely , non - equilibrium transitions .
we shall argue that noise correlated in time can induce re - entrant ordering processes , where the stable or unstable patterns are realized in a window of control parameter values . studying behavior of islands
size as clusters of dense phase and vacancy islands of diluted phase we obtain its time asymptotics .
analytical results are verified by computer simulations .
the paper is organized as follows . in section 2
we describe the model of our system .
section 3 is devoted to study stability of patterns at early stages . in section 4
we derive the effective fokker - planck equation to study noise induced effects in section 5 and formation of stationary patterns in section 6 . in section 7 we discuss evolution of islands size .
we prove our analytical results by computer simulation in section 8 .
finally , conclusions are presented in last section .
we consider a model where only one class of particles is possible .
following refs.@xcite one assumes that the particles can be adsrobed , desorbed , can diffuse and interact among themselves .
therefore , we introduce the scalar field describing dynamics of the local coverage at surface @xmath22 $ ] .
the reaction term incorporates adsorption and desorption terms as follows @xmath23 . here
@xmath24 and @xmath25 are adsorption and desorption rates , respectively ; @xmath26 is the partial pressure of the gaseous phase ; @xmath27 is the interaction potential .
the total flux @xmath7 is a sum of both ordinary diffusion flux @xmath28 and flow of adsorbate @xmath29 . here
the multiplier @xmath30 denotes that the flux is only possible to the @xmath31 free sites .
using the phase field crystals methodology @xcite for the potential @xmath32 one assumes @xmath33 where @xmath34 , @xmath35 and @xmath36 . to estimate @xmath37 , @xmath38 and
@xmath39 one can suppose that there is maximal intensity @xmath40 and characteristic radius @xmath41 @xcite , therefore , one has @xmath42 , @xmath43 , @xmath44 , where the diffusion length is @xmath45 . for @xmath41
one has estimation @xmath46 , whereas @xmath47 @xcite . ] .
the obtained expression can be rewritten in a more simplest form @xmath48 , where @xmath49 is the swift - hohenberg spatial coupling operator where @xmath50 ; here we take @xmath51 that results to relations : @xmath52 , @xmath53 .
next , after substitution of all above terms we get a system of two equations kind of ( [ eqx ] , [ eqj ] ) with @xmath54 where the total diffusion flux @xmath7 satisfies the following relaxation equation @xmath55;\quad g(x)=\frac{1}{4\varepsilon}x(1-x).\ ] ] the equivalent equation for the coverage takes the form @xmath56,\ ] ] where @xmath57 . in our stochastic analysis we are based on this equation generalized by introduction of a multiplicative noise in _ ad hoc _ form using the fluctuation dissipation relation .
therefore , the corresponding langevin dynamics is described by the equation @xmath58 here @xmath59 $ ] ; @xmath60 in order to guarantee that the fluctuation - dissipation relation is satisfied @xcite ; @xmath61 represents langevin force . from the mathematical viewpoint
a generalized langevin equation defined with a memory kernel has a noise fulfilling the fluctuation - dissipation relation with the noise time - correlator as the memory kernel has .
therefore , in our case the memory kernel is reduced to the delta - function and hence , the process @xmath61 is a white noise in time ( with a zero correlation time , @xmath62 ) .
it is known that mathematically defined white noise ( with delta - correlation function ) has no physical meaning @xcite , for example : every next value for @xmath61 does not depend on it previous value ; a spread in values of the noise @xmath61 and a total power of the noise are infinite .
therefore , one needs to consider physical fluctuation source having small but finite correlation time @xmath63 ( see , for example @xcite ) .
hence , in our study we consider a generalized model of the gaussian noise with zero mean and a correlator @xmath64 where @xmath65 is the noise intensity .
hence , a set of main functions defining the system behavior is reduced to the reaction term @xmath6 and the field - dependent diffusion coefficient @xmath66 .
let us study fluctuation effects at early stages .
the stability analysis for the system can be done considering a behavior of the averaged quantity @xmath67 in the vicinity of the state @xmath68 related to unique minimum position of a bare potential @xmath69 .
averaging eq.([eq8 ] ) one gets @xmath70 \right>+\left < g(x)\zeta(\mathbf{r},t)\right>.\ ] ] making use of the novikov theorem @xcite , for the noise correlator we have @xmath71 .
therefore , in the vicinity of the state @xmath68 for the fourier component @xmath72 of the quantity @xmath73 the dispersion relation takes the form @xmath74
^ 2)-\frac{\tau_j\sigma^2}{2}f'''(x_0).\ ] ] it leads to the following dependencies for dispersion curves : @xmath75 ^ 2)-f'(x_0)}{\tau_j}+\frac{\sigma^2}{2}f'''(x_0)-\frac{\gamma^2(x_0)}{4\tau^2_j}\right]^{1/2}.\ ] ] one can see that @xmath76 can have real and imaginary parts , i.e. @xmath77 . the component @xmath78 is responsible for oscillatory solutions , whereas @xmath79 describes a stability of the solution @xmath72 .
analysis of both @xmath78 and @xmath79 allows us to set a threshold for a wave - number where oscillatory solutions are possible , from one hand , and , from another one , to define a wave - number for the first unstable solution . from the obtained dispersion relations it follows that at @xmath80 satisfying equation @xmath81 ^ 2)=\frac{1}{d_0}\left[f'(x_0)-\frac{\sigma^2\tau_j}{2}f'''(x_0)+\frac{\gamma^2(x_0)}{4\tau_j}\right]\ ] ] two branches of the dispersion relation degenerate .
the unstable mode appears at @xmath82 obtained from the equation
@xmath83 ^ 2)=\frac{1}{d_0}\left[f'(x_0)-\frac{\sigma^2\tau_j}{2}f'''(x_0)\right].\ ] ]    from eq.([sag ] ) it follows that a sign of the effective control parameter @xmath84 defines stability of the system .
moreover , as far as the sign of @xmath85 depends on the reaction term derivatives in the vicinity of the state @xmath68 , one can conclude that noise - induced effects in the linear regime are possible at @xmath86 and @xmath87 .    to make a quantitative analysis we consider firstly stability of the deterministic system setting @xmath88 in the vicinity of the unique minimum position @xmath68 of @xmath89 .
as far as @xmath2 increases with the lateral interaction energy , next , we take small values for @xmath37 and study two cases related to small and large values for the adsorption rate .
the corresponding dependencies of @xmath90 and @xmath91 are shown in fig.[det_x_small ] .
a ) for deterministic and stochastic systems ( thin and thick lines , respectively ) at @xmath92 , @xmath93 , @xmath94 : a ) @xmath95 , @xmath96 ( thin and thick lines relate to @xmath88 and @xmath97 , respectively ) ; b ) @xmath98 , @xmath99 ( thin and thick lines relate to @xmath88 and @xmath100 , respectively)[det_x_small],title="fig:",width=302 ] b ) for deterministic and stochastic systems ( thin and thick lines , respectively ) at @xmath92 , @xmath93 , @xmath94 : a ) @xmath95 , @xmath96 ( thin and thick lines relate to @xmath88 and @xmath97 , respectively ) ; b ) @xmath98 , @xmath99 ( thin and thick lines relate to @xmath88 and @xmath100 , respectively)[det_x_small],title="fig:",width=302 ]    from fig.[det_x_small]a obtained in the vicinity of diluted phase characterized by @xmath101 modes with wave - numbers @xmath102 manifest decaying oscillations ( here @xmath103 , @xmath104 ) . inside both intervals
@xmath105 $ ] and @xmath106 $ ] the corresponding modes are stable without oscillations ( here @xmath107 , @xmath108 ) . in the domain
@xmath109 unstable modes characterized by @xmath110 and @xmath107 appear . in the vicinity of the dense phase with @xmath111 ( see fig.[det_x_small]b ) there are no oscillatory solutions due to @xmath107 , but unstable modes can be realized .
here at small wave - numbers one has @xmath112 , whereas at @xmath113 all modes are unstable in the linear regime with @xmath110 in expansion ( [ expansionu ] ) one can find that in the vicinity of the diluted phase with above set of parameters @xmath37 and @xmath114 all modes manifest decaying oscillations , there are no unstable modes , here @xmath115 , @xmath110 . in the vicinity of the dense phase one
has @xmath116 , whereas @xmath110 at @xmath117 . in the stochastic case ( @xmath118 ) in the vicinity of the diluted phase the noise action results in a decrease of the number of oscillating modes shifting @xmath119 toward @xmath120 ; the interval for unstable modes @xmath121 $ ] enlarges ( see thick lines in fig.[det_x_small]a ) . in the vicinity of the dense phase the noise shifts position of the critical wave - number @xmath122 toward small @xmath123 increasing the number of unstable modes ( see thick lines in fig.[det_x_small]b ) .    in figure [ alpha(k )
] we plot stability diagram @xmath124 at fixed @xmath37 for both deterministic and stochastic cases .     at @xmath88 ( thin lines ) and @xmath97 ( thick lines ) .
other parameters are : @xmath92 , @xmath93 , @xmath94 .
domains of unstable modes and stable oscillating modes are bounded by solid and dashed lines , respectively .
[ alpha(k)],width=302 ]    here above solid lines ( thin for the noiseless case and thick for the stochastic one ) the corresponding modes with @xmath125 $ ] are unstable .
below dashed lines the stable modes manifest oscillating behavior . the first unstable mode @xmath126 is realized when equation @xmath127 has unique nonzero solution @xmath126 . at small lateral interaction strength of the adsorbate and small adsroption rate only stable modes without oscillations and stable oscillating modes can be formed . above the threshold @xmath128 unstable modes
appear leading to spatial self - organization of the adsorbate . when the noise intensity increases above pattern forming processes can be observed at lower value for the adsorption rate ; the corresponding patterns are characterized by large wave - lengths then in the noiseless case .
the thick dotted line denotes the dependence of the wave - number @xmath129 related to a position of the maximum for @xmath130 . at @xmath131
the corresponding value @xmath129 defines the most unstable mode giving a period of the spatial structure .
the most unstable mode at @xmath131 is characterized by @xmath132 it follows that @xmath129 depends only on the adsorption rate and the lateral interaction energy .
the noise action leads to a change in the position of the point @xmath126 lying in the line @xmath133 .
dependencies of the period of patterns at different values for lateral interaction energy are shown in fig.[km_alpha_epsilon ] as dashed lines @xmath134 starting at @xmath135 for @xmath97 . here choosing a fixed value for @xmath37 from dependencies
@xmath136 one can find both the threshold value @xmath137 and the corresponding @xmath126 . above the obtained @xmath126 the related period of patterns @xmath138 can be found at fixed system parameters .
it follows that as the noise intensity increases the unstable modes appear at lower values for both @xmath114 and @xmath37 . at elevated @xmath37
the selected patterns are characterized by small period .
( see left ordinate axis ) as dashed lines starting at @xmath135 shown for @xmath97 and threshold values as dependencies @xmath136 ( see right ordinate axis ) at @xmath88 and @xmath97 ( thin and thick solid lines , respectively ) .
other parameters are : @xmath93 , @xmath94 .
[ km_alpha_epsilon],width=302 ]    as far as the first statistical moment illustrates oscillatory behavior one should expect that processes of pattern selection can be observed in such system . considering properties of pattern selection
we need to find dynamical equation for the structure function @xmath139 as a fourier transform of the two point correlation function @xmath140 and to study its oscillatory behavior at small times . to that end
we obtain the linearized evolution equation for the fourier components @xmath141 and @xmath142 and compute @xmath143 . using novikov s theorem , and holding only leading terms in it ,
one gets a dynamical equation @xmath144 ^ 2)\right)s(\mathbf{k},t)+2\sigma^2\gamma(x_0 ) .
\end{split}\ ] ] its analytical solution can be found in the form @xmath145 , where @xmath146 ^ 2)-f'(x_0)\right)}{\tau_j}+\sigma^2f'''(x_0)-\frac{\gamma^2(x_0)}{4\tau^2_j}}.\ ] ] as in the previous case @xmath147 is responsible for stability of the system , whereas @xmath148 relates to pattern selection processes .
unstable modes are defined according to eq.([kc ] ) .
the quantity @xmath149 is determined as solution of the equation @xmath81 ^ 2)=\frac{1}{d_0}\left[f'(x_0)-\frac{\sigma^2\tau_j}{2}f'''(x_0)+\frac{\gamma^2(x_0)}{8\tau_j}\right].\ ] ]    an evolution of the structure function at early stages is shown in fig.[s(kt)]a .
a ) , 2.4 , 3 ( the time is increased from top to bottom ) in the vicinity of the state @xmath96 at @xmath92 , @xmath95 , @xmath93 , @xmath150 .
( b ) the structure function at @xmath151 and different values for @xmath37 , @xmath114 and @xmath65 .
[ s(kt)],title="fig:",width=188 ]  b ) , 2.4 , 3 ( the time is increased from top to bottom ) in the vicinity of the state @xmath96 at @xmath92 , @xmath95 , @xmath93 , @xmath150 .
( b ) the structure function at @xmath151 and different values for @xmath37 , @xmath114 and @xmath65 .
[ s(kt)],title="fig:",width=302 ]    it is seen that during the system evolution the main peak related to wave - number of most unstable mode is formed ; its height is increased . at small time interval except main unstable mode having largest peak in @xmath152 additional peaks at @xmath153 correspond to selecting modes @xcite .
as far as these modes are stable during the system evolution they disappear .
therefore , due to hyperbolic transport pattern selection processes can be realized in the adsorption - desorption systems . in fig.[s(kt)]b
we plot the structure function at early stages at different values for lateral interaction energy and noise intensity . as diagram in fig.[km_alpha_epsilon ] shows an increase in both @xmath37 and @xmath114 shifts the main peak position of @xmath154 toward large values for the wave - number .
the noise action does not shift the main peak position .
the peak spreads due to fluctuation effects resulting in more diffuse interfaces between two phases .
moreover , as the lateral interaction energy increases the pattern selection processes occur faster than at small @xmath37 .
the same is observed for large adsorption rate and noise intensity .
to consider statistical properties of the system at large time scales or in stationary case one can use the corresponding fokker - planck equation . to obtain it we exploit the standard technique and use the equivalent system of two differential equations ( [ eqx],[jtot ] ) written for both generalized coordinate @xmath155 and momentum @xmath156 on a grid of the mesh size @xmath157 : @xmath158 here @xmath159 enumerates cells .
then , the corresponding probability density function is given by the average of the density functional @xmath160 over noise : @xmath161 , where @xmath162 means average over initial conditions . to obtain an equation for the macroscopic density functional @xmath163 we exploit the conventional device and proceed from the continuity equation : @xmath164\overline\rho=0.\ ] ] inserting @xmath165 and @xmath166 from eq .
( [ eqpx ] ) , we obtain @xmath167 where the operators @xmath168 and @xmath169 are defined as follows : @xmath170    within the interaction representation , the density functional reads @xmath171 , that allows us to rewrite eq .
( [ 5_ce7 ] ) as @xmath172 the well - known cumulant expansion method serves as a standard and effective device to solve such a stochastic equation @xcite . neglecting terms of the order @xmath173
, we get the kinetic equation for the averaged quantity @xmath174 in the form @xmath175\langle\wp\rangle(t).\ ] ] integrating the right hand side , one gets @xmath176 within the original representation , the equation for the probability density reads @xmath177\right\}p.\ ] ] because the physical time is much larger than the correlation scale , one can replace the upper limit of the integration with @xmath178 .
to proceed we use the procedure proposed in ref.@xcite to obtain the effective fokker - planck equation for the stochastic model . then , expanding exponents , we obtain for the perturbation expansion @xmath179 where collision operator @xmath180 is defined as follows : @xmath181 here @xmath182 $ ] ( @xmath183 $ ] is the commutator ) ; moments of the noise temporal correlation function @xmath184 are : @xmath185 .    to perform the following calculations
we shall restrict ourselves to considering overdamped systems where the variation scales @xmath186 , @xmath157 , @xmath187 , @xmath188 , @xmath189 , @xmath190 , @xmath191 , and @xmath192 of the time @xmath193 , the coordinate @xmath194 , the quantity @xmath195 , the velocity @xmath196 , the damping coefficient @xmath197 , the force @xmath6 and the noise amplitude @xmath198 , respectively , obey the following conditions : @xmath199 these conditions means a hierarchy of the damping and the deterministic / stochastic forces characterized by relations @xmath200    as a result , the dimensionless system of equations ( [ eqpx ] ) takes the form @xmath201 . \end{split}\ ] ] respectively , the fokker
planck equation ( [ fp ] ) reads @xmath202 where the operator @xmath203 has the components @xmath204 the collision operator is defined through the expressions type of ( [ cop ] ) : @xmath205 after suppressing the factor @xmath206 , the collision operator written with accuracy up to the first order in @xmath207 takes on the explicit form up to second order in @xmath208 : @xmath209\right\}. \end{split}\ ] ] note that in further consideration we are interested in the behavior of the particular distribution @xmath210 , not the total one , @xmath211 .
the reduced distribution can be obtained according to the moments introduced as follows @xmath212,\ ] ] where the integration is provided over set  @xmath213 .
then , performing corresponding manipulations with eq .
( [ fp ] ) , we arrive at the recursive relations for the moments @xmath214 @xcite : @xmath215\\ & = n(n-1)\left(\mathcal{m}^{(0)}-\gamma_i \mathcal{m}^{(1 ) } \right)g_{i}^2{\mathcal p}^{(n-2)}-\epsilon n \mathcal{m}^{(1)}\left[g_{i}^2\frac{\partial{\mathcal p}^{(n-1)}}{\partial x_i}+n g_{i}\left(\frac{\partial g_{i}}{\partial x_i}\right){\mathcal p}^{(n-1)}\right ] . \end{split}\ ] ] here @xmath216 , the first moment @xmath217 can be considered as a flux of the probability density @xmath218 , i.e. @xmath219 . indeed , taking the zeroth moment
, we obtain the expected continuity equation @xmath220 the first moment calculation leads to @xmath221.\ ] ] to evaluate the second moment @xmath222 we put @xmath223 and take into account only terms of zeroth order in @xmath208 : @xmath224    combining all obtained terms we arrive at the system of two equations of the form @xmath225 . \end{split}\ ] ] in the stationary case with no probability density current one gets the solution of the quasi - gibbs form @xmath226\propto \exp\left(-\frac{\mathcal{u}_{ef}[x]}{\mathcal{m}^{(0)}}\right),\ ] ] where the effective energy functional is @xmath227=-\int{\rm d}{\mathbf{r}}\left(\int{\rm d}x ' \varphi(x';\nabla)+\frac{\mathcal{m}^{(1)}}{2}\gamma(x)\right).\ ] ] when we know the effective functional @xmath228 $ ] or its general construction in quadratures it can be used to describe possible phase transitions in the system or to analyze possibility of pattern formation .
let us consider non - equilibrium ( noise - induced ) transitions in a homogeneous system , setting @xmath229 .
the analysis of the effective potential @xmath230 results in noise - induced scenario for modality change of the stationary distribution function @xmath231 .
the corresponding diagrams for noise - induced transitions are shown in fig.[nit ] .
a ) b )    dashed lines in fig.[nit]a correspond to coexistence of two phases : equivalence of two minima of the function @xmath232 ( thin line ) and @xmath233 ( thick line ) .
the dash - dotted lines starting at the point @xmath234 relate to the system parameters where the minimum of @xmath235 at @xmath236 equals @xmath237 , where @xmath238 is the lower natural boundary of the diffusion process @xmath239 . from the mathematical viewpoint at such choice of the system parameters one has a unique minimum at @xmath236 , but due to existence of boundary @xmath240 one has an `` effective minimum '' of @xmath235 at @xmath241 .
this minimum relates to the probability density of the adsorbate absence .
therefore , the dash - dotted lines mean an equivalence of both the absence of the adsorbate and the existence of the dense state with @xmath242 .
the dotted line denotes values for @xmath37 and @xmath114 at which @xmath243 takes zeroth values and belongs to the boundary @xmath240 . at large noise intensity ( @xmath244 )
the dash - dot - dot line starting at the point @xmath245 denotes equivalence of two `` effective '' minima of @xmath235 at @xmath241 and at @xmath246 , where @xmath247 is the upper boundary of the diffusion process .
therefore , the dash - dot - dot line corresponds to equivalence of both states absence of the adsorbate and pure dense state . from fig.[nit]a
it follows that the colored noise results in a decrease of the critical values for the system parameters @xmath114 and @xmath37 ( cf .
cases of white and colored noise influence when @xmath248 and @xmath249 , respectively ) . therefore , if the system is subjected to the colored noise , then non - equilibrium unimodal / bimodal transitions are possible at lower lateral interaction energies @xmath37 and adsorption rates @xmath114 .
bifurcation diagram illustrating stationary values @xmath243 versus the adsorption rate @xmath114 is shown in fig.[nit]b . here for the solid line quantities @xmath250 and @xmath251 belong to spinodals ( solid lines in fig.[nit]a ) ; @xmath252 corresponds to the coexistence line ; @xmath253 ( see dashed line ) relates to the coexistence line of two states @xmath242 and @xmath254 .
the extreme of @xmath228 $ ] relates to the stationary patterns @xmath255 which can be computed setting the first variation of @xmath228 $ ] with respect to @xmath195 equal to zero , i.e. @xmath256=0 $ ] . from eq.([equx ] ) it follows that stationary patterns can be obtained as solutions of the problem @xmath257-\frac{c_0\tau_j\mathcal{m}^{(1)}}{2}f''(x)=0.\ ] ] in order to study stability of possible solution in the vicinity of the fixed point @xmath243 related to extrema positions of the homogeneous distribution we assume @xmath258 , where @xmath39 can have real and imaginary parts responsible for spatial modulations and stability of the solution , respectively .
after linearization of eq.([stpat ] ) one has @xmath259 ^ 2)=f'(x_m)-\frac{\tau_j\tau_c\sigma^2}{2}f'''(x_m),\ ] ] where we put @xmath260 , @xmath261 .    according to obtained values for @xmath39 from eq.([eq_kapa ] )
we have computed a stability diagram shown in fig.[kapa]a for both stochastic system with @xmath262 and deterministic one .
it is seen that in the case of the noiseless system there is unique line of a hyperbolic form dividing domains of unstable @xmath263 an stable @xmath264 patterns ( see thin line ) . at large @xmath37
it tends to the axis @xmath265 .
it follows that increasing the lateral interaction energy or the adsorption rate independently one can induce formation of stable stationary patterns .
if colored fluctuations ( @xmath249 ) are not negligibly small , then the picture of stationary pattern formation can be changed crucially . here
the thick curve dividing domains of stable and unstable stationary solutions at large @xmath37 tends to the dependence @xmath136 denoting formation of the first nonzero @xmath243 ( line ( ob ) in fig.[nit]a ) .
it follows that as lateral interaction energy increases at fixed small values for adsorption rate one can move from unstable domain into domain of stable stationary patterns .
at further increase in @xmath37 we move again into domain of unstable solutions .
therefore , a some kind of re - entrant pattern - forming transition is realized .
starting from domain of unstable stationary patterns at elevated values for @xmath114 we move into the domain of stable stationary solutions .
continuing an increase in @xmath37 we get a re - entrantce , where unstable stationary patterns can be observed .
the corresponding part of the stability diagram denoting re - entrant formation of unstable stationary patterns ( a `` beak '' ) is a part of phase diagram describing noise - induced transitions ( see fig.[nit]a ) . at further increase in @xmath37 one
can move again from domain of stable stationary patterns into domain of unstable stationary patterns . finally , at large @xmath114 the re - entrant ordering ( formation of stable stationary patterns ) can be observed .
one dimensional stationary patterns are shown in fig.[kapa]b where @xmath266 is the maximal period for one from three patterns obtained at @xmath114 and @xmath37 related to coordinates of points @xmath267 , @xmath234 and @xmath245 in fig.[kapa]a .
it is seen that at small @xmath114 one should expect a formation of islands of adsorbate .
such situation is possible for deposition of molecules ( polymers ) . at large @xmath114 relevant for deposition of atoms on metal surface islands of diluted phase ( vacancy islands )
can be organized .
\a )   and @xmath268 denote domains of unstable and stable patterns .
( b ) profiles of stationary patterns related to points @xmath269 , @xmath270 and @xmath271 in plot ( a ) at @xmath272 , @xmath262.[kapa],title="fig:",width=302 ] b )   and @xmath268 denote domains of unstable and stable patterns .
( b ) profiles of stationary patterns related to points @xmath269 , @xmath270 and @xmath271 in plot ( a ) at @xmath272 , @xmath262.[kapa],title="fig:",width=302 ]
let us consider evolution of islands of adsorbate or vacancies ( inclusions ) having radius @xmath273 taking into account that the most probable state @xmath243 is reached and the system fluctuates around it .
we assume @xmath274 where @xmath275 plays a role of thermodynamics concentration value , in other words a supersaturation is small .
let us introduce a quantity @xmath276 measuring a deviation of the inclusion form from a sphere of the radius @xmath273 . using the standard approach @xcite
one can assume the concentration in the vicinity of a spherical inclusion as follows : @xmath277 , @xmath278 is determined trough the surface tension @xmath279 , atomic volume @xmath280 and temperature @xmath281 .
then , the diffusion flux @xmath282 in the vicinity of the inclusion coincides with velocity of its radius change : @xmath283 . using the standard formalism of the coalescence theory @xcite
, we can replace the derivative operator @xmath284 in the flux @xmath285 by @xmath286 . as a result
we obtain an evolution equation for the radius of inclusion in the form @xmath287,\ ] ] where @xmath288 , @xmath289 , @xmath290 plays the role of a supersaturation , @xmath291 is the critical radius of the inclusion , @xmath292 ; next we drop prime for the time variable .
equation ( [ rho ] ) can be rewritten in the form @xmath293 , where a potential @xmath294 can be used in a further study .
at first let us consider the simplest case of @xmath295 where the critical radius @xmath296 remains constant .
in the stationary case we put @xmath297 in eq.([rho ] ) .
such assumption relates to autonomous regime of islands formations .
as far as @xmath298 the critical radius is large but finite .
one of the stationary solution having physical meaning is @xmath299 , where @xmath300 is the thermodynamic limit of the system size .
there is another solution @xmath301 responsible for formation of stationary patterns characterized by @xmath302 .
it can be computed from the condition @xmath303 .
the corresponding dependencies @xmath304 at @xmath305 at different @xmath37 and @xmath306 are shown in fig.[rho_0t]a .
its values are bounded by stability criteria @xmath307 .
it follows that @xmath301 relate to large stable values for @xmath243 we arrive at classical results of the coalescence theory @xcite . ] . using the potential @xmath294 one can find that @xmath294 has two asymptotics : @xmath308 , @xmath309 .
therefore , the quantity @xmath310 corresponds to a maximum of @xmath311 , whereas @xmath301 defines its minimum position .
hence , at @xmath312 all previously organized inclusions at early stages start to grow ; contrary , at @xmath313 the corresponding patterns relax to spherical inclusions having fixed radius @xmath301 . in order to illustrate the corresponding time dependencies of the inclusion radius we solve numerically eq.([rho ] ) ( see fig.[rho_0t]b ) .
a ) at different values for @xmath37 and @xmath306 ( plot ( a ) ) and typical dependencies @xmath314 at different initial conditions and at @xmath315 , @xmath316 , @xmath244 ( plot ( b ) ) .
other parameters are : @xmath272 , @xmath94 , @xmath317 .
[ rho_0t],title="fig:",width=302 ] b ) at different values for @xmath37 and @xmath306 ( plot ( a ) ) and typical dependencies @xmath314 at different initial conditions and at @xmath315 , @xmath316 , @xmath244 ( plot ( b ) ) .
other parameters are : @xmath272 , @xmath94 , @xmath317 .
[ rho_0t],title="fig:",width=302 ]    as far as @xmath301 has unique dependence versus the system parameters , one can solve eq.([rho ] ) analytically rewriting it for the quantity @xmath318 in the form @xmath319 , where @xmath320 , time derivative are taken with respect to @xmath321 ( next we drop prime for the time variable ) . dividing variables , from the formal solution of eq.([rho ] )
following asymptotics can be found : in the vicinity of @xmath322 there are exponentially increasing deviations from the critical radius ; in the vicinity of @xmath323 the exponential decaying @xmath324 is observed .
let us analyze the system behavior considering time dependence of the supersaturation @xmath325 .
hence there is a time dependent variable @xmath326 . next it will be convenient to move to new variables : @xmath327 , @xmath328 . as a new time variable
we introduce @xmath329 .
it yields the derivative @xmath330 , where @xmath331 . using eq.([rho ] ) and taking into account that there is unique value @xmath301
, one has @xmath332 .
therefore , an equation for the variable @xmath333 takes the form @xmath334 let us estimate time asymptotics for both @xmath314 and @xmath326 . to this end
we consider the limit @xmath335 .
it allows one to put @xmath336 or @xmath337 .
as far as the term in the left hand side has no time dependent coefficient and @xmath338 it means that all coefficients in the right hand side should be constants at @xmath339 .
when we assume @xmath340 the coefficient @xmath341 should take constant value too .
unfortunately , from the definition of @xmath342 we arrive that @xmath326 is an increasing function in time leading to non - physical situation where all coefficients from the right hand side of eq.([equ ] ) are constants , not time - dependent functions .
therefore , this assumption fails .
let us consider another case , setting @xmath343 at @xmath335 .
therefore , from the definition of @xmath342 one can find that @xmath344 and @xmath345 .
therefore , the term @xmath346 can be neglected at @xmath335 . using the definition of the variable @xmath333
one can set that the radius of inclusion increases in time according to the lifshitz - slyozov law @xmath347 @xcite .
in order to verify our analytical results we performed numerical simulations by solving the system ( [ eqpx ] ) in @xmath348-dimensional quadratic lattice of the linear size @xmath349 with periodic boundary conditions with the mesh size @xmath350 and the time step @xmath351 .
the noise term was modeled as the ornshtein - uhlenbeck process @xcite . as initial conditions
we have used @xmath352 , @xmath353 .
all measured quantities were averaged over 15 independent runs .
typical patterns organized during the system evolution are shown in fig.[evol_pat]a at @xmath354 .
it is seen that at small times patterns of different sizes are formed . after the system selects patterns with a fixed size . at large time intervals spherical patterns
are organized having preferably hexagonal symmetry as the swift - hohenberg operator predicts due to elastic interactions @xcite . in fig.[evol_pat]b
we plot the corresponding spherically averaged structure function at different times at relevant values for the wave - numbers to show the pattern selecting processes .
it is seen that except the main peak related to period of main patterns there are well defined peaks at small wave - numbers ( at elevated @xmath123 the corresponding peaks have very small intensity , not shown here ) . during the system evolution
its intensity decreases and only main peak is realized , its position is shifted toward @xmath355 defining the period of stationary patterns . in insertion
we have shown structure functions for two cases : @xmath354 and @xmath356 . at @xmath356
the period of patterns is smaller as the above stability analysis shows .
insertion in plot ( b ) shows structure function at @xmath357 and different values for adsorption rate .
other parameters are : @xmath272 , @xmath94 , @xmath262 .
[ evol_pat],title="fig:",width=226 ] b ) .
insertion in plot ( b ) shows structure function at @xmath357 and different values for adsorption rate .
other parameters are : @xmath272 , @xmath94 , @xmath262 .
[ evol_pat],title="fig:",width=302 ]    considering quasi - stationary limit ( time intervals where @xmath358 , @xmath359 and @xmath360 are time - independent quantities , @xmath361 ) we can compute radius of inclusions by measuring surfaces of islands ( number of points belonging to islands having spherical form ) .
evolution of the number of inclusion having fixed radius @xmath273 at different @xmath114 are shown in fig.[n_evol ] .
it is seen that during the system evolution a distribution of islands becomes more compact and it has one the most probable value for @xmath273 .
a ) versus the island size @xmath273 at different times .
plots ( a ) and ( b ) relate to @xmath354 and @xmath362 other parameters are : @xmath272 , @xmath94 , @xmath262 .
[ n_evol],title="fig:",width=302 ] b ) versus the island size @xmath273 at different times .
plots ( a ) and ( b ) relate to @xmath354 and @xmath362 other parameters are : @xmath272 , @xmath94 , @xmath262 .
[ n_evol],title="fig:",width=302 ]    after averaging this radius over independent runs we obtain the dependencies of the optimal characteristic length @xmath363 shown in fig.[rho_exp ] , where @xmath364 means average over experiments . on the top we present snapshots of patterns realized at different values for @xmath114 at fixed @xmath37 .
it is well seen that graphs @xmath365 repeat analytical dependencies @xmath304 in fig.[rho_0t]a at small and large @xmath114 where spherical islands forming hexagonal crystal structure are possible . at intermediate values for @xmath114 strip patterns
are realized and no spherical islands are possible .
therefore , in fig.[rho_exp ] there are gaps related to domains for @xmath114 where no spherical islands are realized .
fitting procedure allows us to find that at small @xmath114 where inclusions of the dense phase are possible the exponential law is realized , @xmath366 , where @xmath267 and @xmath234 are fitting parameters . at large @xmath114
when islands of vacancies are realized one has exponential decaying , @xmath367 .
therefore , an introduction of the swift - hohenberg operator into the model for the system dynamics allows one to describe formation of lattice with different kind of inclusions : at small @xmath114 one has molecular / atom inclusions forming hexagonal lattice , at large @xmath114 one has a situation when deposition of atoms on metal surface can be described with formation of lattice of vacancy islands .     at different @xmath37 and @xmath306 .
stationary patterns observed at @xmath368 at different @xmath114 are shown on the top .
other parameters are : @xmath272 , @xmath94 .
[ rho_exp],width=302 ]    to estimate the spatial size of patterns one can consider typical limits of @xmath369 and @xmath370 .
for the system of deposition of al on tin ( 100 ) the estimation in ref.@xcite shows that at room temperature with the lattice constant @xmath371 and pair interaction energy @xmath372 with lattice coordination number @xmath373 the size of patterns with @xmath374 is around @xmath375 m .
therefore , patterns such as inclusions of dense / diluted phase are of nano - scales and are _ nano - patterns _ @xcite .
patterns of diluted phase usually are realized as lattice of vacancy islands ( for example in absorbed monoatomic layer of ag deposited on ru(0001 ) ref.@xcite ) .
we have studied pattern formation processes in stochastic reaction - diffusion systems with the relaxation flux obeying maxwell - cattaneo equation .
the novelty of our approach is the usage of the phase field crystals formalism with special geometry of realized patterns and introduction of multiplicative noise obeying fluctuation - dissipation relation .
pattern formation is discussed on the model system describing adsorption - desorption processes . considering pattern formation at early stages of the system evolution
we have analyzed temporal stability of patterns and have shown that pattern selection processes are realized as a result of the diffusion flux relaxation to its stationary value .    to study stationary patterns
we have obtained the effective fokker - planck equation and its stationary solution .
it was shown that multiplicative noise can induce non - equilibrium transitions at low energies for the adsorbate interactions .
using stability analysis we have computed diagram illustrating domains of the system parameters where stationary patterns can be realized .
it was shown that introduced multiplicative noise can induce re - entrant pattern - forming transitions where stable / unstable stationary patterns exist inside fixed intervals for lateral interaction energy of the adsorbate .
stationary patterns organized during the system evolution form the crystalline structure of hexagonal symmetry .
it is explained by introduction of the swift - hohenberg operator into the expansion for the interaction potential for the adsorbate .
considering quasi - stationary limit we have discussed evolution of radius of the dense phase inclusions into diluted one or _
vice versa_. using this alternative approach related to coalescence theory it was shown that all patterns are described by the characteristic radius @xmath301 that less then the critical one @xmath296 . in the vicinity of @xmath296 spatial perturbations
grow exponentially , whereas in the vicinity of @xmath301 such perturbations exponentially decay .
we have found that both critical and characteristic radii have time asymptotics related to the lifshits - slyozov law .