discrete nonlinear models play a very important role in many physical applications @xcite . an important class of these models consists of a few completely integrable lattice equations , such as the toda lattice @xcite , the ablowitz - ladik lattice @xcite , and the integrable sine - gordon lattice @xcite .
the fact that these lattices possess exact soliton solutions demonstrates that , in principle , discreteness of the host media does not preclude the propagation of localized coherent structures .
moreover , the mobility of soliton - like excitations in discrete media is a key issue in many physical contexts ; for instance , the mobility of dislocations , a kind of topological solitons , is of paramount importance in the physics of plastic deformation of metals and other crystalline bodies @xcite .
a prototypical class of discrete models , relevant to a variety of applications , consists of the klein - gordon dynamical lattices @xcite .
one of the main representatives of this family of models is the so - called @xmath0 model @xcite , which features a cubic nonlinearity .
this simple power law nonlinearity renders this model a ripe testbed for studying the existence and stability of nonlinear solutions , and comparing their properties in continua and lattices .    in the @xmath1-dimensional continuum framework ( and in the absence of spatially dependent external potentials ) , a solution can be shifted arbitrarily along @xmath2 by any @xmath3 ( @xmath2 is the spatial coordinate and @xmath4 ) , due to the existing translational invariance . on the other hand , in the discrete system
, translational invariance is generically lost and equilibrium static solutions exist only for a discrete rather than for a continuum set of @xmath3 @xcite .
some of these equilibrium solutions correspond to energy maxima and are unstable , while others , corresponding to energy minima , are stable .
the difference between such maxima and minima of the energy is typically referred to as the peierls - nabarro barrier ( pnb ) .
it is of particular interest to develop discretizations that do not feature such barriers . in such cases
, one might expect that the ensuing models would be more faithful representations of their continuum counterparts , regarding both symmetry properties and travelling solution features ( even though there are some notes of caution that should be made ; see e.g. the discussion of @xcite ) .    in that vein , recently , a number of non - integrable discrete klein - gordon equations free of the peierls - nabarro barrier ( pnb - free ) has been systematically constructed .
the first set of such models which were , by construction , hamiltonian ones , was obtained by speight and co - workers @xcite using the bogomolnyi argument @xcite .
a second successful attempt led to the construction of momentum - conserving discretizations @xcite .
it was then demonstrated , surprisingly , that the pnb - free models of that kind conserving both energy and classical momentum do not exist @xcite .
new pnb - free @xmath0 lattices were derived by barashenkov _
. however , we note that for the lattices derived in @xcite , only the kink - type solution has been considered , while it is well - known that the continuum klein - gordon equation can support a number of other solutions @xcite .    on the other hand , a general approach to the construction of the pnb - free lattices
was recently reported in @xcite .
this approach is based on the use of the _
discretized first integral _ ( dfi ) of the corresponding static field equation , and the integration constant that enters
dfi generates the one - parameter solution space .
the dfi is in fact a two - point nonlinear algebraic equation from which the exact static solutions of the three - point pnb - free discrete models can be found .    in this work
we systematically use the dfi approach and derive the hamiltonian and non - hamiltonian pnb - free klein - gordon lattices supporting the one - parameter space of solutions , generalizing the lattices offered in @xcite .
we then focus on the @xmath0 field and discuss in detail the energy - conserving pnb - free model proposed very recently in @xcite , as well as the momentum - conserving model of @xcite ( the two are identical in the static case ) .
for these models we describe the full solution space of the underlying static problem that also includes the solutions derived earlier in @xcite .
our setting is the following : we consider the hamiltonian of the klein - gordon field , @xmath5 , with the kinetic and potential energy functionals respectively defined as , @xmath6 @xmath7dx\,,\ ] ] where @xmath8 is the unknown field and @xmath9 is a given potential function .
the corresponding equation of motion is @xmath10 where @xmath11
.    equation ( [ kg ] ) will be discretized on the lattice @xmath12 , where @xmath13 , and @xmath14 is the lattice spacing .
we would like to construct a nearest - neighbor discrete analog to eq .
( [ kg ] ) of the form @xmath15 such that in the continuum limit @xmath16 we have @xmath17 _ and _ that the solution to the three - point static problem corresponding to eq .
( [ kleingordondiscrete1 ] ) , @xmath18 can be found from a reduced two - point problem @xmath19 if this reduction is achieved , then the exact static solutions can be constructed upon solving the algebraic eq .
( [ twopointstatic ] ) iteratively , starting from arbitrary admissible value of @xmath20 or @xmath21 .
arbitrariness in the choice of the initial condition implies that the static solution can be placed anywhere with respect to the lattice and , for that reason , such lattices are called translationally invariant .    discretization eq . ( [ kleingordondiscrete1 ] ) may result in a non - hamiltonian model due to the non - potential nature of the background forces . in the latter
setting , the absence of an energy functional renders ambiguous the definition of the pn barrier , hence we clarify this point in what follows .
suppose that we have two equilibrium solutions , @xmath22 and @xmath23 ( often , in hamiltonian models with pnb , the second one corresponds to a linearly unstable energy maximum while the first one to a linearly stable energy minimum ) .
the work done by the inter - particle and background forces to move the @xmath24th particle from the position @xmath22 to the position @xmath23 is @xmath25 and the total work performed to `` transform '' the solution @xmath22 to @xmath23 is @xmath26 .
the pn barrier is defined to be equal to @xmath27 . for hamiltonian models
this definition is equivalent to the classical one because @xmath27 is equal to the energy difference between the second ( possibly higher energy unstable ) configuration and the first ( possibly lower energy stable ) configuration . for non - hamiltonian models
, @xmath28 will depend on the path connecting the initial and final configurations of particles . to calculate @xmath27 one
, therefore , has to specify such a path .
while the height of the pn barrier in a non - hamiltonian lattice depends on the path , the pnb - free , non - hamiltonian lattice can be unambiguously defined as the translationally invariant lattice where the quasi - static transformation between the two configurations of interest can be done continuously through the set of _ equilibrium _ configurations . along the path through the equilibrium configurations , forces acting on particles are zero and thus , @xmath29 for any @xmath24 which , in turn , results in @xmath30 .
this notion of the pnb - free lattice is also applicable to the hamiltonian lattices and will serve as our definition of the pnb - free models ; notice , however , that we do not overlook the mathematical subtleties involved in this definition , including the question of whether @xmath27 is zero along all realizable paths connecting @xmath22 and @xmath23 in the non - hamiltonian models , among equilibrium configurations , as is the case for their hamiltonian siblings . these and related questions , including an appropriate modified definition of the relevant quantities for models such as those of the discrete nonlinear schrdinger type , will be left for a future publication .
our scope in what follows is to generalize the approach developed in @xcite to show how to construct _ all possible _ exact static solutions of the models of @xcite and @xcite ( see also @xcite ) , including those derived in @xcite .
this will also lead us to introduce a number of solutions ( both localized and extended ones ) that , to the best of our knowledge , have not been discussed / analyzed previously , such as the ones that will be termed `` inverted '' ( see below ) .
we will also discuss the stability of certain solutions among the obtained ones , for each of the models of interest ( since their stability properties are different ) .
we also note in passing that while our presentation will be geared towards the @xmath0 models , our results regarding the existence of solutions can equally well be applied to discrete equations of the nonlinear schrdinger ( nls ) type @xcite , such as e.g. the ablowitz - ladik model .
in particular , let us consider equations of the form : @xmath31 ( where the overbar denotes complex conjugate ) with @xmath32 and @xmath33 . then ,
looking for standing wave solutions of the form @xmath34 , with @xmath21 real , one would obtain the klein - gordon static problem for the standing wave spatial profile @xmath21 .
hence all the discussion given below for the existence of such solutions can be appropriately translated in the existence of standing waves of the corresponding discrete nls models .
the reader should be cautioned however that the stability properties in the latter context may differ [ a relevant example will be discussed later in the text ] .
our presentation will be structured as follows . in section [ sec :
pnpfreekg ] , we develop the dfi approach and derive pnb - free klein - gordon models supporting one - parameter space of static solutions . in section [ sec : phi4 ] , we present some pnb - free @xmath0 lattices and describe some of their basic properties . in section [ exactsolutionsdiscrete ] , we discuss the details of the construction of the general exact static solutions of the models both in the localized ( hyperbolic function ) and in the extended ( general elliptic function ) form , and for both signs of the nonlinearity prefactor @xmath35 .
@xmath36 corresponds to the so - called defocusing case , while @xmath37 corresponds to the focusing case in the standard terminology of nls equations .
examples of the solutions are then given in section [ analysis ] . in section [ linearstability ] , we analyze the stability of the obtained solutions . in section [ dynamics ] , slow kink dynamics in the pnb - free models is compared numerically with that in the classical @xmath0 model . in section [ continuumsolutions ]
we give a complete list of bounded and unbounded exact solutions to the continuum @xmath0 field . in section [ conclusions ] , we summarize our findings and present our conclusions .
following the lines of the dfi approach of ref .
@xcite , we start from the first integral of the static eq .
( [ kg ] ) , @xmath38 where @xmath39 is the integration constant .
the first integral can also be taken in modified forms @xcite , e.g. , as @xmath40 = 0 \,,\end{aligned}\ ] ] or as @xmath41 = 0 \,,\end{aligned}\ ] ] where @xmath42 and @xmath43 are continuous functions and @xmath44 .
note that @xmath45 and @xmath46 are equivalent only if @xmath47 .
we will consider the case of @xmath48 , @xmath47 , i.e. , the unchanged form of the first integral , eq .
( [ firstintu ] ) , together with the case of @xmath48 , @xmath49 , for which we have the two possibilities @xmath50 and @xmath51 of which only the first one will be discussed below .
we then construct the dfis corresponding to eq .
( [ firstintu ] ) and eq .
( [ firstintv ] ) , which are respectively given by @xmath52 @xmath53 where we demand that @xmath54 in the continuum limit ( @xmath55 ) .
first we construct a pnb - free klein - gordon lattice using the unchanged form of the first integral , eq .
( [ firstintu ] ) , and corresponding dfi , eq .
( [ twopointmapu ] ) . calculating @xmath56 and multiplying the result by @xmath57
, we find @xmath58 discretizing the left - hand side of eq .
( [ d1 ] ) we obtain the lattice klein - gordon equation @xmath59 whose static solutions can be found from the two - point problem , eq .
( [ twopointmapu ] ) , and thus , the lattice is pnb - free .    the lattice of this type was first derived in @xcite where it was also demonstrated that it conserves the momentum @xmath60    the integration constant @xmath39 that appears in eq .
( [ twopointmapu ] ) cancels out in eq .
( [ d1discr ] ) .
this means that all kinds of static solutions derived from eq .
( [ twopointmapu ] ) for different @xmath39 ( kink solution corresponds to @xmath61 ) will be the static solutions to one and the same klein - gordon lattice , eq .
( [ d1discr ] ) , since it does not depend on @xmath39 .    the right - hand side of eq .
( [ d1discr ] ) becomes non - singular when the discretization of the potential @xmath9 is a polynomial function having symmetry @xmath62 .
the most general expression of this type was given in @xcite .
to construct hamiltonian pnb - free lattices we discretize not the equation of motion , eq . ( [ kg ] ) , but the hamiltonian , eq .
( [ kinen ] ) and eq .
( [ poten ] ) .
we now use the modified first integral in the form of eq .
( [ firstintv ] ) with the upper sign and rewrite the potential energy functional , eq .
( [ poten ] ) , as follows , @xmath63 ^ 2     + 2\phi_x\sqrt{2v(\phi)-c } \right)dx\,,\ ] ] where we omitted the constant term .
discretizing the kinetic energy , eq .
( [ kinen ] ) , and the potential energy , eq .
( [ potenfi ] ) , we obtain the discrete hamiltonian @xmath64 ^ 2 \nonumber \\     & + & 2\frac{\phi_n-\phi_{n-1}}{h}\sqrt{2v(\phi_{n-1},\phi_n)-c}\big\}\,.\end{aligned}\ ] ] if the background potential is discretized as suggested in @xcite , @xmath65 then the last term of the hamiltonian eq . (
[ hamiltoniandiscr ] ) reduces to @xmath66 $ ] and it disappears in the telescopic summation . with the choice of eq .
( [ pnpfreecondition ] ) , the equations of motion derived from eq .
( [ hamiltoniandiscr ] ) are @xmath67 where , according to eq .
( [ twopointmapv ] ) and eq .
( [ pnpfreecondition ] ) , @xmath68 this lattice conserves the hamiltonian ( total energy ) @xmath69 ^ 2\big\}\,.\end{aligned}\ ] ] obviously , static solutions to the lattice eq .
( [ hamiltpnpfree ] ) can be found from the two - point dfi @xmath70 and it is also clear that the hamiltonian eq .
( [ hamiltoniandiscr1 ] ) is pnb - free .
the energy - conserving lattice involves the integration constant @xmath39 ( through the function @xmath71 ) and this is different from what we had for the momentum - conserving model . the special case of @xmath61 yields the original model by speight @xcite , which supports a kink solution . models with @xmath72 describe solutions different from the kink solution .
suppose that we use the same function @xmath73 in the dfis eq .
( [ twopointmapu ] ) and eq .
( [ twopointmapv ] ) to construct different @xmath74 terms .
then equations @xmath75 and @xmath76 are equivalent . linear combination of those terms can be used to write the following pnb - free klein - gordon model @xmath77 where the constant coefficients satisfy the continuity constraint @xmath78 assuming that @xmath79 can depend on @xmath14 .
static solutions of eq .
( [ kleingordondiscrall ] ) can be found iteratively from eq .
( [ twopointmapu ] ) or eq .
( [ twopointmapv ] ) .    the model of eq .
( [ kleingordondiscrall ] ) can be generalized in a number of ways .
for example , one can append terms which disappear in the continuum limit and ones that vanish upon substituting @xmath75 .
let us call such terms as @xmath80-terms .
furthermore , any term @xmath81 can be modified by multiplying by a continuous function @xmath82 , which never vanishes and whose continuum limit is unity . such multiplication will not change either the continuum limit , or the static solutions of the model .    particularly we will study the pnb - free model of the form @xmath83,\end{aligned}\ ] ] which is the lattice eq .
( [ d1discr ] ) modified by a multiplier @xmath84 and augmented with the @xmath80-term with arbitrary coefficient @xmath85 .
this @xmath80-term was constructed to fit to the @xmath86 invariant offered in @xcite for the discretization of the @xmath0 field .
finally we note that expressions similar to eq .
( [ d1 ] ) can be derived from derivatives @xmath87 and @xmath88 and they can produce new pnb - free models .
to give examples of klein - gordon lattices , we will discretize the well - known @xmath0 field with the potential @xmath89 where the parameter @xmath90 .
the corresponding equation of motion is @xmath91 the above , so - called , @xmath0 equation supports moving periodic solutions that can be expressed in terms of the jacobi elliptic functions .
the latter can be reduced to localized hyperbolic function solutions ( when the elliptic modulus @xmath92 )
. bounded solutions of this sort were discussed in the context of structural phase transitions @xcite and were used as the starting point for derivation of exact solutions to a discrete @xmath0 model @xcite . in section [ continuumsolutions ]
we complete the list of the solutions presented in @xcite by providing also the set of unbounded solutions of the present model .
the simplest discretization of eq .
( [ kleingordon ] ) is of the form @xmath93 and it possesses a pn barrier .
however , as mentioned above , one can construct pnb - free discrete klein - gordon models by discretizing the nonlinear term @xmath94 typically over three neighboring points , @xmath95 , in contrast to the classical discretization @xmath96 .
the three - point discretization may be physically meaningful in some settings @xcite , but is also interesting from the more fundamental point of view of developing pnb - free discretizations and obtaining analytically ( or semi - analytically ) explicit waveforms of their solutions .
we discretize the potential of eq .
( [ potential ] ) as follows @xmath97 then , the dfis , eqs .
( [ twopointmapu ] ) and ( [ twopointmapv ] ) , become respectively , @xmath98 @xmath99 it has been demonstrated that at @xmath61 eq .
( [ firstintdiscrete ] ) and eq .
( [ twopointmapvkg ] ) support kink solutions @xcite . for @xmath100 they support solutions different from kink .
let us denote @xmath101 substituting eq .
( [ firstintdiscrete ] ) and eq .
( [ twopointmapvkg ] ) into eq .
( [ bop ] ) with @xmath102 we arrive at the following pnb - free @xmath0 lattice @xmath103 equation ( [ barashenkov1 ] ) gives a one - parameter @xmath104 family of pnb - free models .
for the special case of @xmath61 and @xmath36 we get one of the lattices derived in @xcite , @xmath105 , and this lattice supports the kink solution .
lattices with @xmath106 support solutions different from the kink .    for @xmath107
we get from eq .
( [ barashenkov1 ] ) the @xmath39-independent model @xmath108 this non - hamiltonian pnb - free @xmath0 model conserves the momentum eq .
( [ mom ] ) and it will be referred to as the momentum - conserving ( mc ) model .
substituting eq .
( [ firstintdiscrete ] ) into eq .
( [ bop ] ) with @xmath109 and @xmath107 we obtain another @xmath39-independent pnb - free model discovered in @xcite , @xmath110 this model will be called the energy - conserving ( ec ) model because it possesses the hamiltonian @xcite @xmath111 , \label{calh}\end{aligned}\ ] ] where the potential @xmath112 is given by @xmath113    ) , @xmath112 , defined by eq .
( [ hamsaxena ] ) , for @xmath114 and @xmath115 ( solid line ) and @xmath116 ( dashed line ) . for @xmath117
the potential is smooth and it has one minimum at @xmath118 and two maxima at @xmath119 . for @xmath120
the potential has two minima at @xmath119 and a maximum at @xmath118 ; as @xmath121 the potential @xmath122 . ]    in fig .
[ figure1 ] we plot the potential @xmath112 for @xmath114 ( i.e. , @xmath123 ) in both cases @xmath36 ( solid line ) and @xmath37 ( dashed line ) . for @xmath117 the potential is smooth and it has one minimum at @xmath118 and two maxima at @xmath119 . for @xmath120
the potential has two minima at @xmath119 and a maximum at @xmath118 ; note that in the limit @xmath121 , the potential @xmath122 .
it is not possible to plot an analog of fig .
[ figure1 ] for the mc model since the relevant background forces are of many - body type ( i.e. , involve nearest neighbors ) and are non - potential .
pnb - free models given by eq .
( [ barashenkov1 ] ) , eq .
( [ physicad ] ) , and eq .
( [ saxena ] ) have exactly same static solutions derivable from dfi eq .
( [ firstintdiscrete ] ) .
to find _ all _ static solutions to the pnb - free models of eq .
( [ barashenkov1 ] ) , eq .
( [ physicad ] ) , and eq .
( [ saxena ] ) we solve the dfi of eq .
( [ firstintdiscrete ] ) : @xmath124 where @xmath21 and @xmath20 can be interchanged due to the symmetry of the equation .
starting from any admissible `` initial '' value @xmath125 , by iterating eq .
( [ pulse ] ) and its counterpart written as @xmath126 , one can construct recurrently the static solution to both the mc model of eq .
( [ physicad ] ) and the ec model of eq .
( [ saxena ] ) , or to a linear combination thereof .
arbitrariness in the choice of @xmath125 implies the absence of pnb in these models , which has been also demonstrated in @xcite .    as can be seen from eq .
( [ pulse ] ) , once the values of @xmath127 and @xmath128 are fixed , there are certain restrictions on the choice of the values of @xmath125 . in particular , inadmissible initial values are those for which the denominator becomes zero , i.e. , @xmath129 for @xmath120 .
an exceptional case is that of @xmath130 , @xmath131 when an arbitrary sequence of @xmath132 is a solution of eq .
( [ firstintdiscrete ] ) .
inadmissible values of @xmath125 are also ones for which @xmath133 .
the condition @xmath134 leads to a biquadratic algebraic equation determining the borders of admissible region ; the roots of this equation are : @xmath135 let us introduce the following notations for these roots : @xmath136     in the nonlinear map eq .
( [ pulse ] ) for different values of @xmath127 at @xmath137 .
there are three inadmissible regions marked with `` no '' and two inadmissible lines @xmath138 ( horizontal dashed lines ) . on these lines
the on - site potential @xmath112 diverges ( see fig . [ figure1 ] ) .
@xmath139 ( @xmath140 ) designate different branches of the borders of the admissible regions [ see eq .
( [ borders ] ) and eq .
( [ rootnotation ] ) ] .
roots @xmath141 and @xmath142 merge at @xmath131 and at @xmath143 .
roots @xmath142 and @xmath144 merge at @xmath145 .
vertical dashed lines separate regions with different jacobi elliptic function solutions . vertical dotted line , situated ( for the chosen parameters ) at @xmath146 ,
divides the region of the @xmath147 solution into two portions corresponding to two roots of the first equation in eq .
( [ p1m1p1d ] ) . ]    the admissible regions for the values of @xmath148 of the nonlinear map eq .
( [ pulse ] ) are shown in fig .
[ figure2 ] for different values of @xmath127 at @xmath137 .
the corresponding result for @xmath149 is presented in fig .
[ figure3 ] .
these graphs present a road map for constructing the various possible solutions of the above models .    the symmetry of eq .
( [ firstintdiscrete ] ) suggests that the topology of the admissible regions is such that once started from an admissible value of @xmath148 , one can not leave the admissible region iterating eq .
( [ pulse ] ) , so that the static solution will surely be constructed for the whole chain .
this is so because eq .
( [ pulse ] ) serves for calculating both back and forth points of the map , and if one is admissible , the other one is also admissible .    equation ( [ pulse ] ) possesses two roots , which means that for an admissible initial condition one can construct two different solutions , e.g. , a kink and an antikink . when iterating , to keep moving along the same solution
, one must take @xmath150 different from @xmath151 [ if the roots of eq .
( [ pulse ] ) are different ] .
indeed , setting in the three - point static problem , eq .
( [ physicad ] ) , @xmath152 , we find @xmath153 comparing this with eq .
( [ pulse ] ) , it is readily seen that the equilibrium in the three - point equation in the case of @xmath152 can be achieved only if @xmath134 and the two roots coincide . as mentioned above ,
the latter requirement is equivalent to the condition that @xmath21 is on the border of the admissible region .    as it can be deduced from eq .
( [ borders ] ) , the topology of the admissible regions presented in fig .
[ figure2 ] and fig .
[ figure3 ] does not significantly change among different but positive and among different but negative values of @xmath128 respectively ; of course , one has to exclude the particular case of @xmath130 and also the continuum limit ( see sec .
[ continuumsolutions ] ) . here
we will not discuss in detail the case of extremely high discreteness @xmath154 , even though the analysis of this case does not present any additional difficulties .    to summarize , inside the admissible region , @xmath155 ( shown in figs .
[ figure2 ] and [ figure3 ] ) , the nonlinear map of eq .
( [ pulse ] ) generates static solutions to the mc model of eq .
( [ physicad ] ) and to the ec model of eq .
( [ saxena ] ) .
but for @xmath149 . roots @xmath141 and @xmath142 merge at @xmath131 , and roots @xmath142 and @xmath144 merge at @xmath145 .
the @xmath156 solution , eq .
( [ n0n0p1d ] ) , is defined for @xmath157 .
the @xmath158 solution , eq .
( [ n0p1n0d ] ) , is defined for @xmath159 , and this region is divided into two parts , @xmath160 and @xmath161 , each corresponding to a particular root of the first equation in eq .
( [ n0p1n0d ] ) .
the border between these two regions is shown by the dotted line situated ( for the chosen parameters ) at @xmath162 . ]
static solutions to the discrete pnb - free models of eq .
( [ physicad ] ) and eq .
( [ saxena ] ) have been reported @xcite in the form of the jacobi elliptic functions , @xmath163 and @xmath156 @xcite .
below , we will derive such solutions .
the general form of the solutions is @xmath164 where @xmath165 is the modulus of the jacobi elliptic functions , @xmath166 and @xmath167 are the parameters of the solution , and @xmath3 is the arbitrary initial position . finally the integers @xmath168 specify a particular form of the solution .    in the limit of @xmath92 , eq . ( [ ellipticdiscrete ] ) reduces to the hyperbolic function form @xmath169 and , in the limit of @xmath170 , to the trigonometric function form @xmath171    substituting the ansatz eq .
( [ ellipticdiscrete ] ) into eq .
( [ firstintdiscrete ] ) and equalizing the coefficients in front of similar terms we find that it can be satisfied for a limited number of combinations of integer powers @xmath168 .
solving the equations for the coefficients we relate the parameters of the solution eq .
( [ ellipticdiscrete ] ) , @xmath167 and @xmath166 , to @xmath128 and @xmath172 and also find the relation between @xmath39 , that enters eq .
( [ firstintdiscrete ] ) , and @xmath172 .    for some of the sets @xmath173 ,
e.g. , for @xmath174 and for @xmath175 , we obtain imaginary amplitude @xmath166 in the whole range of parameters .
essentially different , real amplitude solutions described by eq .
( [ ellipticdiscrete ] ) , have the following form and are characterized by the following parameters :    the @xmath176 solution , @xmath177 , @xmath178 the @xmath179 solution , @xmath180 , @xmath181 the @xmath156 solution , @xmath182 , @xmath183 the @xmath184 solution , @xmath185 , @xmath186    the @xmath187 solution , @xmath188 , @xmath189    the @xmath147 solution , @xmath190 , @xmath191    it is worth making the following remarks :    the solutions shown in eq .
( [ n0p1n0d ] ) and eq .
( [ n0n0p1d ] ) have real amplitudes for @xmath192 while the others for @xmath193 .
the solutions should be interpreted in the following form . for a given @xmath128
, one can find @xmath167 by solving the first equation in eqs .
( [ p1n0n0d]-[p1m1p1d ] ) , and then @xmath166 from the second one . substituting these values in eq .
( [ ellipticdiscrete ] ) results in the static solutions of the original discrete model .
the expressions for @xmath127 in eqs .
( [ p1n0n0d]-[p1m1p1d ] ) link the elliptic jacobi function solutions and the solution in the form of the nonlinear map , eq .
( [ pulse ] ) . as for the other free parameter of the solutions eq . ( [ p1n0n0d]-[p1m1p1d ] ) , the arbitrary shift @xmath3 , its counterpart in the nonlinear map , eq .
( [ pulse ] ) , is effectively the initial value @xmath125 .
the solutions of eqs .
( [ p1n0n0d]-[p1m1p1d ] ) can be expressed in a number of other forms using the well - known identities for the jacobi elliptic functions @xcite .
for example , shifting the argument by a quarter period , one can transform the @xmath176 solution to the form of @xmath194 , or , applying the ascending landen transformation , to the form of @xmath195 .
mathematically , these three expressions look as different members of eq .
( [ ellipticdiscrete ] ) , but physically they are indistinguishable .
in what follows , we analyze the static solutions derivable from eq .
( [ pulse ] ) discussing their relation to eqs .
( [ p1n0n0d])-([p1m1p1d ] ) . since the topology of the admissible regions is different for positive and negative @xmath128 , these two cases are studied separately .      in the @xmath176 solution ,
( [ p1n0n0d ] ) , when the elliptic modulus increases from its smallest value @xmath170 to its largest value @xmath92 , the amplitude @xmath166 increases from @xmath196 to @xmath197 , since @xmath198 [ see eq .
( [ borders ] ) and eq .
( [ rootnotation ] ) ] . as a result ,
the parameter @xmath127 monotonically decreases from @xmath199 to @xmath196 .
thus , the @xmath176 solution is defined in the portion of the @xmath200-plane , @xmath201 and @xmath202 ( see fig . [ figure2 ] ) .
the @xmath203 solution , eq .
( [ m1n0n0d ] ) , is defined for @xmath165 , and it is complementary to the @xmath176 solution since it is also valid for @xmath204 , but for @xmath205 ( see fig .
[ figure2 ] ) .    the @xmath206 solution , eq .
( [ n0m1n0d ] ) , is defined for @xmath165 and it occupies the region @xmath207 , @xmath205 ( see fig .
[ figure2 ] ) .
the @xmath147 solution , eq .
( [ p1m1p1d ] ) , is defined for unlimited @xmath125 in the range @xmath208 ( see fig . [ figure2 ] ) .
this solution is only valid for @xmath209 , where , for fixed @xmath35 , @xmath210 is an increasing function of @xmath14 and @xmath211 .
for @xmath212 the first expression in eq .
( [ p1m1p1d ] ) does not have solutions for @xmath167 .
for @xmath209 , the equation has two roots , @xmath213 . for the limiting value , @xmath214
, one can find the corresponding amplitude @xmath215 from the second expression of eq .
( [ p1m1p1d ] ) , and then @xmath216 , from the last expression . for the case of @xmath137 presented in fig .
[ figure2 ] , we find @xmath217 and @xmath146 .    it should be noticed that we have not found a solution of the form of eq .
( [ ellipticdiscrete ] ) valid in the range of @xmath218 ( portion marked with the question mark in fig .
[ figure2 ] ) .
it is likely that static solutions in this range can not be expressed in terms of the jacobi elliptic functions because they do not survive in the continuum limit ( see sec .
[ continuumsolutions ] ) .
however , the solution can easily be constructed from the nonlinear map in eq .
( [ pulse ] ) .
let us now discuss further several particular examples of the above solutions .
first of all , in the limit @xmath219 [ see eq .
( [ hyperbolicdiscrete ] ) ] , the @xmath176 solution , eq .
( [ p1n0n0d ] ) , reduces to the kink solution @xcite , @xmath220 , \label{kinkcooper}\end{aligned}\ ] ] while the @xmath203 solution , eq .
( [ m1n0n0d ] ) , to the solution called hereafter the `` inverted '' kink , @xmath221}. \label{invertedkink}\end{aligned}\ ] ] in eqs .
( [ kinkcooper ] ) and ( [ invertedkink ] ) , @xmath3 is the ( arbitrary ) position of the solution and @xmath222 .
this limiting case corresponds to @xmath131 , for which a heteroclinic connection is possible between the fixed points @xmath223 and @xmath224 ( see fig .
[ figure2 ] ) , giving rise to the kink or the inverted kink . in this case , eq .
( [ pulse ] ) assumes the following simple form @xmath225 where one can choose either the upper or the lower signs .
the kink , eq .
( [ kinkcooper ] ) , and the inverted kink , eq .
( [ invertedkink ] ) , can be derived from eq .
( [ kink ] ) taking initial values from @xmath226 and @xmath227 , respectively .    in fig .
[ figure4 ] we show ( a ) the kink and ( b ) the inverted kink solutions taking @xmath228 and @xmath229 , respectively , for @xmath137 .    : ( a ) kink and ( b ) inverted kink at @xmath137 .
the solutions can be found from the nonlinear map eq .
( [ kink ] ) with the initial conditions @xmath228 and @xmath229 for ( a ) and ( b ) , respectively .
the kink and the inverted kink are also given by eqs .
( [ kinkcooper ] ) and ( [ invertedkink ] ) respectively , for @xmath230 .
horizontal dashed lines , @xmath231 , show positions of singular points of the potential @xmath112 ( see fig . [ figure1 ] ) . ]
close to 0 , at @xmath137 , obtained from the map of eq .
( [ pulse ] ) : ( a ) @xmath232 , @xmath233 and ( b ) @xmath232 , @xmath234 .
the solution in ( a ) is the @xmath235 solution , eq .
( [ p1n0n0d ] ) , and the solution in ( b ) is the @xmath147 solution , eq .
( [ p1m1p1d ] ) .
horizontal dashed lines , @xmath231 , show positions of singular points of the potential @xmath112 ( see fig .
[ figure1 ] ) . ]    in fig .
[ figure5 ] , and for @xmath137 , we present two examples of solutions for positive and negative @xmath127 close to @xmath196 . in particular , for @xmath236 , taking initial value @xmath232 , we obtain from the map of eq .
( [ pulse ] ) the solution presented in panel ( a ) . in fact , it is the @xmath176 solution close to the hyperbolic function limit having the form of a periodic chain of kinks and anti - kinks . on the other hand , for @xmath234 and @xmath232 , we obtain from the map of eq .
( [ pulse ] ) the solution shown in panel ( b ) .
this is the @xmath147 solution close to the hyperbolic function limit and has a form of a chain of kinks and inverted anti - kinks .    in fig .
[ figure6 ] , and for @xmath137 , we present the @xmath237 solution , eq .
( [ m1n0n0d ] ) , ( a ) close to the hyperbolic and ( b ) at the trigonometric limits .
the solution shown in panel ( a ) is a chain of inverted kinks and inverted anti - kinks .
these solutions are obtained from the nonlinear map eq .
( [ pulse ] ) setting @xmath238 , @xmath239 for ( a ) , and @xmath240 , @xmath241 for ( b ) .    in summary
, we have shown that the well - established ( hyperbolic and elliptic function ) solutions of the model correspond to a very narrow region of the two - parameter ( @xmath242 ) admissible space .
a natural question is what is the typical profile outcome stemming from the use of other pairs of @xmath125 and @xmath127 in eq .
( [ pulse ] ) .
generically , upon testing the different regions of the admissible regime we have observed that arbitrary choices may lead to seemingly erratic solutions with very large amplitudes . a different sign choice in the right hand side of eq .
( [ pulse ] ) may , however , lead to a periodically locked tail structure .
a simple example of such a solution is given below .
solution , eq .
( [ m1n0n0d ] ) , ( a ) close to the hyperbolic and ( b ) at the trigonometric limits .
we set @xmath137 .
the solution in ( a ) consists of a chain of inverted kinks and inverted anti - kinks . to obtain these solutions from the nonlinear map eq .
( [ pulse ] ) we set @xmath238 , @xmath239 for ( a ) , and @xmath240 , @xmath241 for ( b ) . ]      for negative @xmath128 , we have the @xmath179 solution , eq .
( [ n0p1n0d ] ) , and the @xmath156 solution , eq .
( [ n0n0p1d ] ) .
let us first start from the latter : when @xmath172 increases from @xmath196 to @xmath197 in eq .
( [ n0n0p1d ] ) , the parameter @xmath39 decreases monotonically from @xmath196 to @xmath199 , i.e. , this solution is defined in the `` dn '' portion of fig .
[ figure3 ] .    on the other hand , the @xmath179 solution of eq .
( [ n0p1n0d ] ) is only valid for @xmath209 , where @xmath210 is an increasing function of @xmath14 and @xmath211 .
for @xmath212 the second expression of eq .
( [ n0p1n0d ] ) does not have solutions for @xmath167 .
for @xmath209 , the equation has two roots , @xmath213 . for the limiting value , @xmath214
, one can find the corresponding amplitude @xmath215 from the second expression of eq .
( [ n0p1n0d ] ) and then @xmath216 from the last expression . when @xmath172 increases from @xmath214 to @xmath197 in eq .
( [ n0p1n0d ] ) , the parameter @xmath127 of the nonlinear map eq .
( [ pulse ] ) corresponding to the root @xmath160 ( @xmath161 ) increases from @xmath216 ( decreases from @xmath216 ) to @xmath199 ( to @xmath244 ) .
thus , the @xmath158 solution occupies the rest of the admissible region for the case of @xmath243 ( see fig .
[ figure3 ] ) . in the case of @xmath149 presented in fig .
[ figure3 ] , we find @xmath245 and @xmath162 .    both @xmath179 and @xmath156 solutions , in the limit @xmath246 ( @xmath247 ) , reduce to a homoclinic to @xmath196 pulse solution ( see also figs .
[ figure1 ] and [ figure3 ] where this possibility is illustrated ) .
this solution has the form : @xmath248 , \label{pulsecooper}\end{aligned}\ ] ] where @xmath249 this is illustrated by fig .
[ figure7 ] where , taking @xmath149 , we show the solution obtained from the map eq .
( [ pulse ] ) at ( a ) @xmath250 ( @xmath156 solution ) and ( b ) @xmath251 ( @xmath179 solution ) , for initial value of @xmath252 .
the figure clearly illustrates the two limits ( to the left and to the right of @xmath145 ) and their correspondence to pairs of pulse  pulse solutions and ones of pulse  anti - pulse solutions , as one enters the two different regimes `` dn '' and `` cn '' of fig . [ figure3 ] .
solutions of this special form can be expressed as @xmath253 with constant @xmath166 and @xmath254 . for the sake of simplicity , we set @xmath255 and substitute the ansatz into eq .
( [ physicad ] ) or eq .
( [ saxena ] ) to find the zigzag solution @xmath256 to obtain this solution from the map eq .
( [ pulse ] ) , one has to set @xmath257 , which is the general condition for getting @xmath258 . substituting here @xmath166 from eq .
( [ zigzag ] ) , we get @xmath259 . for the case of @xmath137 presented in fig .
[ figure2 ] , we have @xmath260 .    the zigzag solution is an exceptional one , as @xmath261 for any @xmath24 .
however , as shown above , this is only possible when eq .
( [ pulse ] ) has multiple roots , which means that the two - point static problem , eq .
( [ firstintdiscrete ] ) , is factorized .
this solution does not have a counterpart in the continuum limit .     at ( a ) @xmath250 and ( b ) @xmath262 , for initial value of @xmath252 . since @xmath127 is close to @xmath263 , corresponding to @xmath92 , the dn solution in ( a ) and the cn solution in ( b ) look like chains of separated pulses given by eq .
( [ pulsecooper ] ) . ]
the static solutions of the discrete models eq . ( [ physicad ] ) and eq . ( [ saxena ] )
are exactly the same , but the dynamical properties of the two models are different . as the corresponding stability analysis takes into account the dynamical form of each model , it will be carried out separately for each of the two .
let us first consider the mc model of eq .
( [ physicad ] ) . introducing the ansatz @xmath264 ( where @xmath265 is an equilibrium solution and @xmath266 is a small perturbation ) , we linearize eq .
( [ physicad ] ) with respect to @xmath267 and obtain the following equation : @xmath268 for the small - amplitude phonons , @xmath269 , with frequency @xmath270 and wave number @xmath271 , around the uniform steady states @xmath272 ( @xmath120 ) , eq . ( [ linphysicad ] )
is reduced to the following dispersion relation : @xmath273 while the spectrum of the vacuum solution for @xmath117 , @xmath274 , is @xmath275    for an arbitrary stationary solution @xmath265 , stability is inferred in the mc model if the eigenvalue problem obtained from eq .
( [ linphysicad ] ) by replacing @xmath276 with @xmath277 has only non - negative solutions @xmath278 .
recalling that the mc model is a non - hamiltonian one , the resulting eigenvalue problem is a non - self - adjoint one involving a non - symmetric matrix .
similarly , we obtain analogous expressions for the ec model of eq .
( [ saxena ] ) .
the linearized equation reads @xmath279 ^ 2 } \varepsilon_{n } , \label{linsaxena}\end{aligned}\ ] ] and the corresponding dispersion relation for the linear phonon modes around @xmath280 has the form : @xmath281 while that of vacuum solution @xmath118 ( for @xmath117 ) is @xmath282    the stationary solution @xmath265 is stable in the ec model of eq .
( [ saxena ] ) if the self - adjoint eigenvalue problem obtained from eq .
( [ linsaxena ] ) by replacing @xmath276 with @xmath277 has only non - negative solutions @xmath278 .
we note that all stable and unstable static solutions of the mc and the ec models , except for the zigzag solution eq .
( [ zigzag ] ) , possess a zero - frequency mode .
this is a consequence of the effective translational invariance of the discrete pnb - free models ; this is related also to the freedom of selecting the free parameter @xmath3 in the corresponding solution expressions .
our aim here is not to study the whole bunch of solutions in the whole range of parameters , but rather to demonstrate the existence of stable solutions and also provide some examples where solutions , being stable in one model , may be unstable in another . the results of stability analysis for some characteristic solutions are summarized in table [ table ] .
first , we consider the kink and inverted kink solutions given in eq .
( [ kinkcooper ] ) and eq .
( [ invertedkink ] ) respectively : in a numerical experiment , these solutions were placed in the middle of a chain of @xmath283 particles with clamped boundary conditions . we have found that for the chosen set of parameters , the kink is stable in both mc and ec models , while _ the inverted kink is stable in the ec and unstable in the mc model_.    in fig .
[ figure8 ] we show the spectra of the kink and the inverted kink for their different positions with respect to the lattice , @xmath3 . in the figure , the horizontal lines show the borders of the phonon band of the vacuum @xmath119 and the dots show the frequencies of the kink s internal modes lying outside of the band . in all cases , the spectra are shown as functions of @xmath3 , even though , admittedly , the kink in the mc model presented in ( a ) demonstrates very weak sensitivity of its spectrum to variations of @xmath3 .
let us now consider the periodic solutions depicted in fig .
[ figure5](b ) ( @xmath147 ) and in fig .
[ figure6](a ) ( @xmath203 ) . for these classes of solutions we used periodic boundary conditions and
the length of the lattice was commensurate with the period of the solutions containing a number of periods .
close to the hyperbolic function limit , these solutions are stable in the ec model but unstable in the mc model , following the stability of the building blocks ( inverted kinks ) associated with these structures . on the other hand ,
the @xmath203 solution becomes unstable in both models at the trigonometric limit depicted in fig .
[ figure6](b ) .
nonetheless , we find it remarkable that some of these `` apparently non - smooth '' solutions ( or maybe more correctly , solutions apparently containing a non - smooth continuum limit ) may be potentially stable in the discrete setting
. it would be wortwhile to analyze the stability of such structures in more detail and as a function of their elliptic modulus @xmath172 .
next , we consider the zigzag solution , eq .
( [ zigzag ] ) , which exists for @xmath284 .
it can be shown that this solution can be stable in the ec model : indeed , substituting eq .
( [ zigzag ] ) into the relevant eigenvalue problem , we obtain the dispersion relation for the linear phonon modes , @xmath285 $ ] . the solution is stable when @xmath278 are non - negative , i.e. , when @xmath286 . combining this with the existence condition we find that the zigzag solution exists and it is stable in the ec model for @xmath287 .
note that this exceptional solution does not possess the zero frequency mode , as it can be seen from the above dispersion relation .
the amplitude of the solution , @xmath288 , is always greater than 1 .
the zigzag solution is always unstable in the mc model , a result that can similarly be demonstrated .
, @xmath137 ( @xmath114 ) .
the inverted kink in mc model is unstable for the chosen parameters .
horizontal lines show the borders of the phonon bands . ]
the @xmath176 , @xmath179 , and @xmath156 solutions close to the hyperbolic limit were found to be generically unstable in both the mc and ec models .
the instability of the pulse solution of eq .
( [ pulsecooper ] ) can be analytically justified in the present setting .
it can be easily inferred from the invariance of the solution with respect to @xmath3 that the eigenvector leading to the zero eigenvalue of fig .
[ figure8 ] is proportional to @xmath289 ( where @xmath290 is the pulse solution profile ) .
this is an anti - symmetric eigenvector , given the pulse s symmetric ( around its center ) nature .
but then , from sturm - liouville theory for discrete operators @xcite , there should be an eigenvalue such that @xmath291 with a symmetric ( i.e. , even ) eigenvector , resulting in the instability of the relevant pulse solution .
our numerical simulations fully support this conclusion .
the pulse solution was constructed by iterating eq .
( [ pulse ] ) for @xmath243 , @xmath145 , and @xmath292 , or directly from eq .
( [ pulsecooper ] ) .
we found that for @xmath37 the pulse solution is unstable in both models and over a wide range of the discreteness parameter @xmath14 , including the case of rather small @xmath293 .
we have also confirmed that the instability mode is similar to the pulse profile ( i.e. , of even parity ) .    both the cn and dn solutions ( see ,
e.g. , fig .
[ figure7 ] ) are also found to be unstable for the different parameter values that we have used in both the mc and ec models .
this is rather natural to expect given that their `` building blocks '' , namely the pulse - like solutions are each unstable ( hence their concatenation will only add to the unstable eigendirections ) .
finally , as far as the sn solutions are concerned , two basic instability modes were revealed .
one resembles the instability of the pulse solution , when the whole structure , being waved around the top of the background potential , tends to slide down to one of the potential wells .
the other instability modes can be understood if one regards the sn solution as a set of alternating kinks and anti - kinks [ see fig . [ figure5](a ) ] interacting with each other by means of overlapping tails .
such a system is unstable because the kinks tend to annihilate with the anti - kinks .
it is interesting to note that these stability results suggest a disparity with the earlier numerical findings of @xcite where the pulse - like solution was generically reported to be stable and analogous statements were made for the cn- and dn - type solutions .    at this point
, it is worth discussing a connection with nls - type models mentioned in the introduction . in the case of such models , the aforementioned unstable eigendirection ( with an eigenvector similar in profile to the pulse itself ) is `` prohibited '' by the additional conservation law of the @xmath294 norm , resulting in a separate neutral eigendirection with zero eigenfrequency . as a consequence
, it is an interesting twist that the instability reported for such pulse - like solutions in scalar @xmath0 models would be absent in their discrete nls counterparts .
.results of stability analysis [ table ] [ cols="<,^,^ " , ]
in @xcite , we have compared the spectra and the long - term dynamics of the kinks in the two pnb - free models , eq .
( [ physicad ] ) and eq .
( [ saxena ] ) , with that of the classical discrete model , eq .
( [ kleingordondiscrete ] ) ( see fig . 1 of that paper ) .
the case of very small kink velocities was not studied there . here
we would like to focus on this case to demonstrate the qualitatively different behavior for slowly moving kinks in the pnb - free and classical models .    to boost the kink in the pnb - free models eq .
( [ physicad ] ) and eq .
( [ saxena ] ) we used the following dynamical solution corresponding to the multiple eigenvalue @xmath295 , @xmath296 , where @xmath265 is a static kink solution , @xmath267 is the normalized translational eigenmode and @xmath297 is the amplitude playing the role of kink velocity . this construction yields a more accurate approximate solution for very small @xmath297 when linearized equations eq .
( [ linphysicad ] ) and eq .
( [ linsaxena ] ) are accurate .
increasing @xmath297 leads to a decrease in the accuracy of the dynamical solution used for boosting .
on the other hand , we note in passing that there are techniques developed recently for klein - gordon @xcite and even nonlinear schrdinger lattices @xcite that also allow the construction of finite speed , numerically exact travelling solutions in these classes of models .    in fig .
[ figure9 ] we present the lowest frequency normalized eigenmodes for the inter - site kinks in the three models , for the case of @xmath114 , @xmath36 ( @xmath137 ) . both pnb - free models , eq .
( [ physicad ] ) and eq .
( [ saxena ] ) , have the same translational eigenmodes with @xmath298 , shown by dots connected with solid lines .
one can show by straightforward algebra that the pnb - free model modified by a non - vanishing multiplier @xmath299 as described in sec .
[ sec : generalizations ] has the same translational eigenmode as the original model . assuming that the two models are respectively : @xmath300 where the function @xmath301 , then the steady state solutions @xmath265 will satisfy @xmath302 , for all @xmath24 , while the corresponding linearization equations ( using @xmath303 ) are respectively of the form @xmath304 . \label{extra4}\end{aligned}\ ] ] hence , when solving the corresponding eigenvalue problem ( again substituting @xmath305 by @xmath306 ) , for @xmath298 , the eigenvalue problems become identical , both satisfying @xmath307 hence the coincidence of the corresponding eigenvectors .
we are extremely grateful to an anonymous referee for this remark and its proof .
we do note , in passing , also that while this justifies the coincidence of the zero eigenfrequency modes of the pnb - free models in fig .
[ figure9 ] , it is also in tune with the results of fig . [ figure8 ] for nonzero eigenfrequencies @xmath270 .
the latter are not identical between the different models , as indicated by the left and middle panels of that figure .
this is due to the differences between the corresponding eigenvalue problems of eqs .
( [ extra3 ] ) and ( [ extra4 ] ) , when @xmath308 .    for the classical model , eq .
( [ kleingordondiscrete ] ) , the lowest - frequency mode has @xmath309 , and it is shown in fig .
[ figure9 ] by open circles and dashed lines .
actually , this mode is not a translational mode ( since , strictly speaking , there is no translational invariance ) but we use it to boost the kink .
one can say that this will _ become _ the translational mode for this model in the continuum limit .
we define the kink center of mass as @xmath310    the evolution of the kink coordinates is shown in fig .
[ figure10 ] .
kinks were boosted with two different amplitudes of the normalized lowest - frequency eigenmodes , @xmath311 and @xmath312 .
results for the pnb - free models practically coincide , as depicted by the dashed and solid lines for the mc and ec models , respectively .
it is readily seen that kinks in the pnb - free lattices propagate with roughly constant velocities .
the oscillatory trajectories in fig .
[ figure10 ] correspond to the kink in the classical discrete @xmath0 model .
the faster kink , boosted with @xmath312 , propagates along the lattice but its velocity gradually decreases . the kinetic energy of the translational motion is partly lost to the excitation of the kink internal mode with @xmath313 , lying below the phonon frequency band .
higher harmonics of this mode interact with the phonon spectrum producing radiation , which also slows the kink down .
an even more dramatic difference is observed for the slower kinks , boosted with @xmath311 . here , the classical kink can not overcome the pn barrier and can not propagate , oscillating near the stable inter - site configuration . on the other hand ,
the kinks in the pnb - free models are not trapped by the lattice and propagate due to the absence of the pn barrier .
alternatively , one can say that such waveforms can be accelerated by arbitrarily small external fields .    ) and eq .
( [ saxena ] ) , have the same translational eigenmode with @xmath298 , shown by dots connected with solid lines ( see also the relevant discussion in the text ) . for the classical model , eq .
( [ kleingordondiscrete ] ) , the lowest - frequency mode has @xmath314 ( shown by open circles and dashed lines ) .
results for @xmath114 , @xmath36 ( @xmath137 ) . ]    , as a function of time for kinks propagating in the pnb - free mc ( dashed line ) and ec ( solid line ) models , as well as in the classical discrete @xmath0 model ( oscillatory line ) .
the kinks were boosted with two different amplitudes of the normalized lowest - frequency eigenmodes , @xmath311 and @xmath312 .
the faster classical kink is able to propagate , while the slower one is not , since it can not overcome the pn barrier .
kinks in pnb - free models are not trapped by the lattice and can therefore propagate . ]
in the continuum limit , @xmath55 , the borders of the admissible region , eq .
( [ borders ] ) , become @xmath315     field for ( a ) @xmath316 and ( b ) @xmath37 obtained as the continuum limits ( @xmath55 ) of those presented in fig .
[ figure2 ] and fig .
[ figure3 ] , respectively . in each panel there is one inadmissible region marked with `` no '' . ]    in fig .
[ figure11 ] we plot the admissible regions for ( a ) @xmath36 and ( b ) @xmath37 .
the topology of the admissible regions for the @xmath0 field is simpler than the one pertaining to the discrete models . in the continuum limit there
exists a sole inadmissible region for both cases @xmath90 , while three inadmissible regions exist for the discrete models in the case @xmath193 .
one more simplification is that the domains of the @xmath147 and @xmath179 solutions do not split into two parts since the smaller root @xmath160 disappears in the continuum limit .
particularly we note that the region marked with the question mark in fig .
[ figure2 ] , @xmath317 , disappears in the continuum limit .
that might be the reason why we failed to find a jacobi elliptic function expression for the discrete static solutions in this case .
static solutions eqs .
( [ p1n0n0d]-[p1m1p1d ] ) obtained for the discrete models have their continuum counterparts as the traveling solutions to the @xmath0 field , eq .
( [ kleingordon ] ) .
the general form of the solutions is @xmath318 where @xmath165 is the modulus of the jacobi elliptic functions , @xmath166 and @xmath167 are the parameters of the solution , @xmath3 is the arbitrary initial position and @xmath319 is the velocity of the solution .
the integers @xmath168 once again specify the particular form of the solution .
continuum analogues of eqs .
( [ p1n0n0d])-([p1m1p1d ] ) have the following form and are characterized by the following parameters :    the @xmath176 solution , @xmath177 , @xmath320    the @xmath179 solution , @xmath180 , @xmath321    the @xmath156 solution , @xmath182 , @xmath322    the @xmath203 solution , @xmath185 , @xmath323    the @xmath206 solution , @xmath188 , @xmath324    the @xmath147 solution , @xmath190 , @xmath325    the above six solutions can be rewritten in a great variety of forms using the properties of the jacobi elliptic functions @xcite .
however , we believe that they may constitute the full list of physically different solutions to the continuum @xmath0 equation since they fill the whole two - parameter space , @xmath326 , obtained as the continuum limit of corresponding space of the discrete models .
all solutions are conveniently parameterized by a single parameter @xmath327 for @xmath120 , and @xmath328 for @xmath117 , as it is presented in fig .
[ figure11 ] .    the solutions in eqs .
( [ p1n0n0])-([n0n0p1 ] ) are bounded while the other ones are unbounded .
the solutions in eq .
( [ n0p1n0 ] ) and eq .
( [ n0n0p1 ] ) are defined for @xmath329 while the others for @xmath120 .
the solutions in eq .
( [ n0p1n0 ] ) and eq .
( [ p1m1p1 ] ) are valid for @xmath330 , the solution in eq .
( [ n0m1n0 ] ) is valid for @xmath331 , while the other solutions for @xmath165 .
in the present paper we have shown that the reduction of the static problem of a discrete klein - gordon ( and by extension of the standing wave problem of a discrete nls ) equation to a two - point problem is a powerful tool for obtaining all possible static solutions of the corresponding model .
we have applied this general idea to a momentum conserving @xcite and an energy conserving @xcite discretization of the @xmath0 model analyzing the full two - parameter plane of solutions and giving a natural parametrization for it . in particular , we have examined the admissible regions of the field value at a given point and of the constant entering the two - point function pertinent to the model .
we have specifically illustrated how to use different choices of these two parameters to obtain not only the well - known hyperbolic function solutions and the established elliptic function solutions , but also novel classes of solutions including the inverted ( non - monotonic ) kink solutions presented herein and multi - kink generalizations thereof .
we performed such computations both for the attractive ( focusing ) and for the repulsive ( defocusing ) types of nonlinearity .
the presented methodology has the significant advantage over earlier work that by introducing the integration constant in the discretization of the first integral of the static part of the pde , it allows one to construct the _ full _ family of the static solutions of the corresponding model .
earlier work had implicitly set this additional free parameter to @xmath196 ; this choice was sufficient for obtaining the important hyperbolic function solutions , but the present formalism systematically illustrates how to generalize the latter .
we have also derived the continuum analogues of the discrete solutions and found that they fill the whole space of parameters obtained in the limit of @xmath332 from the parameter space of the discrete models .
in addition , we have examined some of the key stability features of the solutions obtained in the various models , illustrating the different stability properties of the models considered herein ( even if their static solutions are identical ) .
we have obtained interesting stability properties , including , e.g. , the counter - intuitive stability of the inverted kink and of some of its periodic generalizations for the ec model .
the present study generates a variety of interesting questions . for instance
, it would be relevant to examine whether the solutions presented herein ( including the inverted ones ) have counterparts in the `` standard '' discrete @xmath0 model and to analyze their respective stability .
it would also be relevant to examine how the stability ( and existence since some of them may not survive the continuum limit ) of such solutions depend on the lattice spacing @xmath14 .
finally , these and related ( e.g. stability ) questions in the context of discrete nls lattices ( see @xcite ) would also be equally or even more ( given the multitude of relevant applications of the latter model ) interesting to answer .
such studies are currently in progress and will be reported elsewhere .
the authors would like to thank d.  e.  pelinovsky for a number of stimulating and insightful discussions .
we are also grateful to an anonymous referee for offering a proof that the mc and ec models have the same translational eigenmodes .
pgk gratefully acknowledges support from nsf - dms-0204585 , nsf - dms-0505663 and nsf - career .