^{1}

^{1}

^{1}

^{*}

^{2}

After some reduction procedure made on the nonlinear evolution equation for nerve pulses, based on thermodynamic principles, new classic and non-classic traveling solutions have been obtained. We have studied this model for particular values in the parameter space, and obtained both the bell and compacton like solutions. These nonlinear traveling waves could be responsible for transmitting efficiently the necessary information along the axons. The non-classic structures named as compactons, due to their robust configuration, could be considered in some sense a more realistic type of nonlinear chargers of information. The last solutions do not have tails and as adiabatic waves could propagate along the nerve with constant velocity that could be equal, smaller or higher than the sound velocity.

One of key fundamental problems in biophysics is to understand how nature makes to carry information without significant distortion along distances between two considerable long separated centers. Now it is well established that the dynamic of ionic currents through voltage channels is responsible for the change of the membrane potential in nerve tissues. As is well known the first measurements on ionic currents were performed by Hodgkin and Huxley (HH) in the 50s [

Later, the Hodgkin and Huxley system was developed independently by Richard FitzHugh, and Jin-ichi Nagumo [3,4]. Based on the work of Balthazar Van Der Pol, FitzHugh proposed a simplified neuronal model of Hodgkin and Huxley. For its part, Nagumo suggested as analogous neuronal, a nonlinear electrical circuit, controlled by an equation system also similar to those of Van Der Pol currents. The proposed simplified analog of these authors, is called FitzHugh-Nagumo model. Being susceptible fairly complete analysis, the FHN system allows a qualitative understanding of the phenomenon of excitability, from the point of view of dynamic systems [

Recently, Heimburg and coworkers have developed a model for nerve pulses that support soliton like solutions [7-10]. The model is constructed considering the nerve axon as a one dimensional cylinder with lateral density excitations moving along the axes that is represented by the coordinate z. This model shows the appearance of lipid phase transition slightly below physiological temperatures. Given measured values of the compression modulus as a function of lateral density and frequency, soliton properties can be determined by the velocity of the traveling waves. That is, resuming we can say that this theory is based on the lipid transition from a fluid to a gel phase at slightly below of body temperatures. The effects of nonlinearity and dispersion as is common would be responsible for appearance of soliton like structures in nerve membrane in the gel state [

On the other hand, as is well known, classical solitons possess infinite tails that could implicate some distortions at long distances far from the excited zone. This constriction could also be applied to nerve solitons as well. However, according to several important results from experiments made on nerve axon, there is a suggestion that the nerve pulse need to be narrow and probably they could posses depreciable tails. Fortunately, several types of non-classic solutions without or depreciable tails were found in a wide class of nonlinear processes in physics.

We suppose that along the axon, not only well famous bell solitons could propagate, but also non-classical soliton like solutions named compactons could be excited and travel carrying needed information. Thus, our main interest in this contribution is to show the appearance of compactons along the nerve, that is, solitons without tails. This nonlinear robust entity will eventually fit more realistically to the features of that waves that are reported on nerve pulses in the literature. We use the model of Heimburg and coworkers [

The detailed discussion on methods and proposals for obtaining the nonlinear differential equation which is the subject of our analysis, could be found in the appropriated literature, see for example [7-10]. The theory is based on hydrodynamic properties of a density pulse in the presence of dispersion. The authors started with the classic wave propagation equation in the absence of dispersion along the quasi-one dimensional axon:

being the change of density in the membrane, is the density of the membrane at physiological condition slightly above of melting transition. This “field” is supposed to have dependence of coordinate x and time t. The excitations are moving along the coordinate x. Here c is the sound velocity with the value

, being the compressibility. Subsequently, it is added the additional term responsible for dispersive effects. This term appears by experimental evidence on the relation between frequencies and velocities. For getting available nonlinear real dispersion relation the next term after the second derivative with respect to the variable x would stay the corresponding fourth derivative [

Finally, one can conclude that the partial differential equation for propagating waves can be represented as

As can be easily seen this equation is similar if we restrict ourselves to take only the first two terms of the Taylor series in the rhs of Equation (2), to the well known Boussinesq like equation.

It is possible then to study the nerve pulse in the framework of a completely integrable system like Boussinesq one (with restrictions) as a very special case for which the solutions obtained there could fit some experimental results. The Boussinesq like Equation (4) describes in addition to the bright or regular solutions, also singular or anti-cusp travelling solitons [

Despite the great effort for research concerning these type of models, and results, there is a question that still is waiting its successful solution. Indeed these theories predict the appearance of solitary waves with broad extension i.e. broad excitation region that contradicts experiment results. The experiments made with nerve pulses show the existence of narrow localized region of excitation. Thus, there is a big possibility that these pulses are in fact solitons but without tails. As is well known the classical solitons possess tails that could eventually interact at long distances that subsequently would degenerate the whole pulse to charge the needed information. So, the effort in this sense will be assumed to find solitons without tails i.e. solitary waves named regular or singular compactons.

As is well known, such solution with finite span is currently found in many natural processes. Thus, in this direction Rosenau and Hyman [

For the case m = 2, n = 2 It is possible to found

where the solution is valid only inside the interval, and vanishes outside of this interval. These solutions including interact elastically like hard spheres do not possess infinite tails. The invariance under the transformation and produces the anticompactons.

After this important discovery, compactons have found diverse applications in natural science, we cite here a few of them, for example in viscous films [

Let us investigate different kind of solutions that can be obtained by making suitable transformations of parameter values on the Equation (3). By making the standard transformation for considering the existence of traveling waves along the nerve, we introduce the new independent variable. Thus, the Equation (3) can be transformed to the next one

here C is a constant that is obtained after integration. This equation is very similar to a common equation that represents the Newton second law of mechanics. After subsequent standard transformation one can obtain the next equation

with being and the parameters that appear in the Taylor expansion and the constant of integration.

Now let us study this equation by taking in mind the trivial boundary condition, that means at long distances from the main excited zone the perturbation pulse vanishes and also its first derivative. Thus, we have for this boundary condition the standard requirement

By applying this restriction the constant of integration satisfies:. It is important to mention that this value does not depend on any experiment conditions, the only restriction is that far from the excited zone of excitation the perturbation wave disappears.

First, let us consider the case:, when the nonlinear wave will move with the sound velocity along the axis z. Analyzing the possible consequences of this reduction, one can find that at the maximum amplitude value of solitary waves the right hand side of the Equation (8) for trivial boundary condition has three solutions, i.e. we face with an algebraic cubic equation. At the simplest case when we can obtain two solutions for the U difference density. So, in both cases the maximum amplitude (when and) of possible solitons is not equal to zero and the solution could exist as we demonstrate in the next subsection. Thus, we perform analytic studies for two sub cases, indeed.

1) Sub-case:. In this simple case by integrating the Equation (8) we obtain the following solution

By avoiding singular behaviors, the possible solution in this case is achieved by assuming the negative value of the parameter q. If this is satisfied, we can obtain a solitary wave that resemble the well known bell like soliton, see

off with a power-law distribution, i.e., algebraically.

2) Sub-case:. For this sub-case it was also possible to find a particular physical interesting compact solution when the parameters of the equation of motion (8) satisfy the condition

For all compact type solutions we will follow the next strategy. Compacton solutions will be constructed by patching a compact portion of a periodic solution that is zero at both ends to a solution that vanishes outside the compact region to give a weak solution to the equation. We see that the condition for a weak solution for Equation (3) in the case of traveling waves with independent variable is that the jump across the boundary of the equation of motion at is

The integration of the Equation (8) yields the following compact soliton like solution

with the new variable and the parameter.

This solution survive only in the segment mentioned above and exactly vanishes outside of this region and satisfies the trivial boundary condition. The points are that points where the solutions and its first derivative vanish and satisfy the condition of solvability (12). Nextit is assumed that or. The

Now we can show other type of solutions that support the thermodynamic model [

with the renamed parameters

Without lost of generality we can choice. In order to solve the differential Equation (14) it seems natural to introduce the following ansatz: Soliton like structures are available if the parameters of the strong model satisfy the following relation:

Then, after integration we again obtain a perfect compact like solution.

When this is satisfied, one trivially can obtain the solution for the density displacements that takes the form of compact soliton like solution

with. As can be observed this specific solution exists by avoiding singularities i.e. in the particular case only for negative parameter values

The

This expression depends on the soliton velocity by virtue of Equation (15). Outside of this dependence between the parameters the solutions probably will not admit compact support. The

of the Equation (19). This expression gives to us the wide spectrum for using the parameter values.

By deploying the constraint (18) we obtain the parameter region of available velocities. Thus, when the velocity corresponds to subsonic waves, i.e.

If the values of we obtain supersonic soliton waves with the next segment of values

One of the implications of these relations is that it would travel along the axon robust solitonic waves with speed greater than the sound velocity.

The energy of solutions obtained for the ordinary nonlinear differential equation i.e. for traveling solutions can be calculated by using the energy density, which in some sense could be considered as energy density of a mechanical analog particle. The integration with respect to the variable in the interval for compactons can be directly evaluated using the energy density. From the Equation (8) one can notice that the kinetic energy is equal to the potential one for the “mechanical analog particle”. This leads to find the total energy as

with the “potential” piece of energy density

Thus for the first algebraic soliton that moves with the same sound velocity the energy in dependence of the main parameters p and q is expressed as (where h = 1)

The energy for the second compacton solution that moves with the sound velocity is

with.

Finally the last compacton solution has the following expression for its energy.

.

This total energy depends on the velocity of the nonclassical solution. For avoiding singular solutions, when the values, traveling compact solution could emerge only when the velocities satisfy the relation, that is in complete accordance with the previous analysis made on the parameter space (18) for obtaining subsonic and also supersonic compact solutions. Here the parameters A and are determined by the Equations (15) and (16).

We have discussed in this contribution the appearance of classical and non-classical soliton like pulses along the axon based on the pioneering work of [

Authors are indebted to Professor V. G. Makhankov for constant support and discussions. We acknowledge the efforts of SIEA-UAEMEX for supporting the fundamental research and also the financial aid of Secretaría de Educación Pública de México (SEP) under the project FEO1/2012 103.5/12/2126 for developing research scientific groups.