Non-kirchhoff vertices and nonlinear schrodinger ground states on graphs

: We review some recent results and announce some new ones on the problem of the existence of ground states for the Nonlinear Schrödinger Equation on graphs endowed with vertices where the matching condition, instead of being free (or Kirchhoff’s), is non-trivially interacting. This category includes Dirac’s delta conditions, delta prime, Fülöp-Tsutsui, and others.


Introduction
In this paper we aim at presenting some new trends in the study of the Nonlinear Schrödinger equation on metric graphs.Specifically, we focus on cases where non-Kirchhoff's conditions are imposed at the vertices.
Roughly speaking, metric graphs, also called networks, are a simple mathematical scheme used to build up models of branched structures, as they are made of one-dimensional segments, named edges or arcs, that can intersect in points called vertices.The only geometric information available for such structures is the length of the edges, whereas no curvature, angles, nor embedding in higher dimensional settings are considered.Therefore, metric graphs are suitable to provide simple models for the propagation of signals of various origin: acoustic, electric and also matter waves (for an exhaustive physical introduction see for instance [50]).
As it is intuitive, in order to be defined, a wave dynamics on a network requires an evolution equation, describing the change of the profile in time inside the edges, together with a matching condition at the vertices, that rules the transmission and the reflection of the signals when crossing a junction.In the context of quantum graphs, namely metric graphs in which the ruling equation is the linear Schrödinger, the problem of finding all admissible transmission and reflection rates, maps to the issue of finding all possible self-adjoint extensions of the restriction of the Laplacian to functions that vanish in a neighbourhood of every vertex.This problem was exhaustively treated in a seminal paper by Kostrykin and Schrader [47].In fact, self-adjointness is the translation in the language of operators of the consevation of the total probability, that is a fundamental and inescapable requirement for every genuine quantum theory.In the context of quantum graphs, self-adjointness, and therefore conservation of probability, is encoded exactly in the law ruling the transmission and the reflection at the vertices.More generally, the task of finding all self-adjoint extensions of a symmetric operator is the key ingredient of the definition of point interactions, i.e. potentials located at a single point in space.For a collection of fundamental results we refer the reader to the treatise [19].On the other hand, for the application of the theory of self-adjoint extensions to graphs, the basic reference is the monograph [21].Among all the possible matching conditions, the most used are those named after Kirchhoff as a reminescence of Kirchhoff's law for linear circuits.Such conditions prescribe that the sum of the derivatives of the wave function ingoing into every vertex equals zero (for the mathematical formulation of such conditions see (6)).In the particular case of a vertex attached to one edge only, Kirchhoff's condition reduces to Neumann's, while in the case of two edges only concurring to the same vertex, Kirchhoff's condition restores the requirements of continuity and differentiability at the point occupied by the vertex.Furthermore, Kirchhoff's conditions naturally arise when dealing with the search for ground states, namely with the minimization of the energy functional under, as a unique constraint, the value of the mass.Finally, when imposing Kirchhoff's conditions, transmission and reflection coefficients are independent of the momentum.For all these reasons, Kirchhoff's conditions have been assumed as the most natural, and hence the most widely studied.However, it is not clear whether in relevant physical contexts such conditions prove able to fit the dynamics well.On the contrary, Fülöp, Tsutsui and Cheon ([38, 55]) suggested that some other conditions could be more satisfactory from the point of view of invariance laws.Furthermore, the presence of non-trivial, localized interactions near junctions make quite unlikely that Kirchhoff's conditions would supply a predictive model.Therefore, quite recently a programme for the study of nonlinear dynamics on graphs including non-Kirchhoff's conditions has started.Here we report on the first results and perspectives of this project with a particular emphasis on the existence and stability of ground states for the Nonlinear Schrödinger equation on graphs, endowed with non-Kirchhoff's conditions at the vertices.
In the literature, metric graphs are graphs G = (V, E), where V is the set of the vertices and each edge e ∈ E is identified either with a bounded and closed interval I e = [0, l e ] or with a halfline I e = [0, +∞).
Owing to the metric structure, on this kind of graphs one can define real or complexvalued functions and the related ordinary functional spaces, like Lebesgue and Sobolev spaces, as and Consequently, a function ψ on a metric graph can be seen as a collection of functions ψ = (ψ e ) e∈E and the variable x on the graph runs through the collection of all variables x e defined on each edge.
As already specified, dynamics on metric graphs can provide a good approximation for the evolution of systems located on ramified structures, i.e. systems for which locally there is a privileged direction for the propagation of signals, since the dimensions transverse to that of propagation are negligible compared to the longitudinal one.Such structures are often referred to as quasi one-dimensional.This line of investigation has recently undergone a dramatic boost due to its closeness to the problem of the evolution of Bose-Einstein Condensates (BEC).Let us just briefly recall that a Bose-Einstein Condensate is a system composed by a large number of identical bosons (often alkali atoms) whose spatial confinement is usually realized by magneto-optical traps.Bose in 1924 [23] and Einstein in 1925 [35] predicted that, under a critical value of the temperature, the state of the whole system collapses into a non-classical state in which each particle acquires the same wave function, called wave function of the condensate, that solves the following variational problem min where Ω is the trap in which the particles are confined, N is the number of the particles of the system and finally E GP is the Gross-Pitaevskii functional defined as In equation (1), α is the scattering length of the two-body interaction between the particles in the condensate.
Provided that it exists, a solution u ω to the variational problem (1) must satisfy the Euler-Lagrange equation where ω arises as a Lagrange multiplier and depends on N , and it is immediately seen that the function ψ(t, x) = e iωt u ω (x) is a solution, in particular a standing wave, to the Gross-Pitaevskii equation The interaction between atoms in a BEC is usually repulsive, so that the sign of α is in general positive, but it is nowadays possible to tune such interaction through a mechanism called Feshbach resonance [24], so that it become possible to create collapsing condensates, by making α to become negative.This fact makes interesting to study Gross-Pitaevskii equation with a focusing nonlinearity.Concerning the trap Ω, it is generically intended to be a smooth region of the three-dimensional space.It is then possible to arrange experimentally elongated traps, and also branched traps by the insertion of suitable junctions.In this cases, the trap remains genuinely three-dimensional, but it is commonly accepted, even though a general and rigorous proof is still lacking, that, as a result of a suitable shrinking limit, Gross-Pitaevskii equation for a three dimensional system can be approximated by a focusing nonlinear Schrödinger equation on a metric graph like where the graph can be understood as a quasi one-dimensional skeleton of the original elongated and branched trap.In (2), H is a self-adjoint extension of the Laplace operator where Ie is the open interval corresponding to the edge e, and the space C ∞ 0 Ie collects all smooth and compactly supported functions on Ie .A prolific research line is focused on the search for standing waves of (2), i.e. solutions of (2) of the form ψ(t, x) = e iωt ϕ ω (x), where ω ∈ R and ϕ ω solves the stationary equation A significant part of the literature has dealt with stationary solutions for the focusing nonlinear Schrödinger equation when Kirchhoff's boundary conditions are imposed at the vertices, namely coupled with the boundary conditions where e ≻ v means that e is an incoming edge of the vertex v and dϕe dxe (v) stands for ϕ ′ e (l e ) if x e = l e at v and for ϕ ′ e (0) if x e = 0 at v. In this context, two main variational approaches have been used to find solutions to this equation, both consisting in minimizing a proper functional under some additional constraints.In the first approach one looks for critical points of the energy functional on the manifold where p > 2 and µ is a fixed parameter usually called the mass.In the second case, instead, one is interested in critical points of the action functional under a constraint called Nehari's constraint, namely where . This constraint is considered as natural since it hosts all stationary points of the action functional (8).In the following, we will refer to global minimizers of (7) or (8) as ground states, regardless of the functional they minimize.The difference between the two approaches reflects on the parameter ω ∈ R in the equation (4).In fact, it can be an unknown of the problem and be interpreted as a Lagrange multiplier like in the former approach or it can be given, like in the latter.
The forerunner of the study of nonlinear evolution on metric graphs is considered to be a paper by Ali Mehmeti [20] that dates back to 1984, but it is in the last decade that the study of the NLSE with Kirchhoff's conditions has developed.Several papers [16,17,15] study the existence of ground states for the energy functional under the mass constraint.They analyse the problem for metric graphs with a finite number of vertices and at least one halfline, distinguishing between cases p ∈ (2, 6) (known as subcritical) and p = 6 (called the critical case).In the first case, the authors show how the topology of the graph can affect the existence of ground states, or, on the other hand, when it depends on the interplay between the metric features of the graph and the mass µ.In the critical case, instead, the mass assumes a crucial role for the existence of ground states.In particular, differently from the case of the real line, the authors prove that ground states can exist not only for a critical value of the mass, but for a whole interval of masses.In addition, the existence of ground states for periodic graphs has been studied in [10,11,32,9], while the problem on infinite metric trees has been approached and partially solved in [33].Finally, existence of ground and/or bound states for the NLS on graphs with a nonlinearity concentrated on a subgraph has been variously explored, for instance in [51,34,54].
On the other hand, the search for ground states for the NLS on the line with delta or delta-prime interactions, that can be thought of as a graph with a vertex and two infinite edges with a non-Kirchhoff's condition, is older than the extended specific research on graphs [37,12].Point interactions on the real line have been previously studied in the time-dependent setting by Caudrelier et al. using integrability tools [28,29].We stress however that these results are bounded to the case of the cubic nonlinearity, i.e. p = 4, that corresponds to the integrable case [56].Other non-free boundary conditions on graphs have been extensively investigated in connection with the integrability features.A breakthrough result, due to Matrasulov and coworkers, is the discovery of a class of non-reflecting matching conditions that make the cubic NLS on graphs inherit the integrability from the corresponding one-dimensional system ( [53,52]) and, at least for star graphs, make possible to restore the structure and the methods typical of integrable systems, like Lax pairs and inverse scattering: this was accomplished in [27], whose results extend to non-integrable boundary conditions too.Other milestone results on graphs with the same non-reflecting (hence non-Kirchhoff) conditions have been obtained by Pelinovsky and collaborators [44,45,46] in a series of works where the spectral stability of special solutions (like half-solitons or shifted states) was investigated.The main purpose of this review is to present some of the results in the literature or still in preparation which involve nonlinear Schrödinger equations on metric graphs when non-Kirchhoff's conditions are imposed at the vertices.The physical motivations for the introduction of such conditions at the vertices rely on the necessity to represent some inhomogeneity or defect in the medium in which the dynamics occurs.Rigorous studies of NLSE in presence of impurities described by point interactions have been given along several lines, with a special consideration for the so-called delta interaction.As already specified, delta interactions have been the oldest non-Kirchhoff's conditions to be studied and mathematically they are described as conditions localized at the vertices v ∈ V, involving both the value of the function and its derivative.Specifically, they are defined by with α ∈ R, and they are obtained as the result of a proper self-adjoint extension of the Laplace operator in (3).In constrast to Kirchhoff's conditions, up to now the study of delta conditions has been confined on simple graphs with a single vertex, namely the real line R or more generally star graphs S N with N halflines.This is the starting point for possible future studies on more general graphs, a step that turns out to be highly non-trivial, considering that for general power nonlinearities it is not possible to directly make recourse to abstract methods used e.g. in [30].
Through the study of delta interactions one is naturally led to consider two different classes of non-Kirchhoff's conditions: the first one is the class of linear non-Kirchhoff's conditions, which ensure self-adjointness of the operator defined on the metric graph and include delta prime, dipole and Fülöp-Tsutsui's conditions, while the second one is the class of nonlinear delta interactions, obtained replacing the real number α with a nonlinear function of ϕ.In the latter case, two nonlinearities coexist: the standard one, given by the p-th power of the L p norm, and a pointwise one.The review is organized as follows: in Section 2 we will introduce the linear non-Kirchhoff's conditions, collecting some old results and presenting other ones more recent.In Section 3 some new results concerning the nonlinear delta conditions are shown.

Linear non-Kirchhoff's condition
As outlined in the Introduction, linear non-Kirchhoff's conditions are a family of conditions imposed at the vertices of a metric graph G in such a way that the operator H in (2) turns out to be self-adjoint.Among the first ones who studied non-Kirchhoff's conditions, in the integrable cubic case there were Caudrelier et al. [29], who presented a family of point interactions that preserves the quantum integrability and Goodman et al. [39] and Holmer et al. [42], who introduced delta interactions and started the study of the existence and the stability of solutions of the NLSE.As for the problem with Kirchhoff's conditions, two main approaches have been carried on: the first one concerns the minimization of the energy under the mass constraint, while the second one is based on the minimization of the action functional restricted to the Nehari's manifold.

Minimization of the energy under the mass constraint
This first approach has been used in [13], where the authors study existence of 1D ground states and their orbital stability when a point interaction is placed at the origin of the real line and the standard nonlinearity is subcritical, i.e. when 2 < p < 6.In particular, they are interested in three different conditions at the origin: • attractive delta conditions, i.e.
We recall that these boundary conditions are induced by three of the possible selfadjoint extensions of the 1D Laplacian on the real line and up to now they are the most studied in the stationary setting.In this context, the authors of [13] prove an abstract theorem that revisits the concentration-compactness method by Cazenave and Lions [26] and which is suitable to treat all these three inhomogeneities.Applying this general result to all these three cases, it is straightforward that ground states exist and they are orbitally stable: moreover, thanks to the one dimensional structure, it is possible to compute explicitly solutions.These results do not reveal any substantial novelty moving from the Kirchhoff's case to these non-Kirchhoff's conditions: indeed, a ground state exists for every mass also in the Kirchhoff's case.
The relevance of point interactions becomes clearer passing from the real line R to star graphs S N .Indeed, in the paper [3], the concentration-compactness method was adapted to the case of an attractive delta interaction localized at the vertex of a star graph and prove that there exists a threshold value of the mass µ * such that, under this value, the ground state exists, is symmetric, decreasing on each halfline and orbitally stable: among all the stationary states, this is called the N tail state (see Figure 2).In particular, the result is valid both in the subcritical and in the critical case.
It is important to notice that, in the Kirchhoff's case, ground states do not exist for any value of the mass and for this reason the delta interaction is crucial to gain the existence of ground states for small masses.Moreover, the N tail state turns out to be orbitally stable for any value of the mass µ > 0 (see [4]) : this means that the orbital stability stands even when the N tail state is not a ground state but only a local minimizer of the energy functional.

Minimization of the action under the Nehari's constraint
As noticed in the introduction, there is a second way to approach the study of the NLS equation and it consists in minimizing a constrained action funcional.The usual aims in this setting are to identify stationary states, characterize ground states and show if standing waves are stable or not for any value of the power p.
We remark that in this section stability always stands for orbital stability and the outcomes concerning it are achieved taking advantage of the well known theory by Grillakis, Shatah and Strauss [40,41].
Fukuizumi et al. (see [37], [36] and [48]) investigate the previous problems on the real line with a delta type defect located at the origin, analysing both the attractive case (α > 0) and the repulsive one (α < 0).In particular they use the conditions defined in (10) with α ∈ R and study the existence and the stability of global minimizers of the action functional under the Nehari's constraint.A first difference between the results obtained in the attractive case and in the repulsive one concerns the functional space in which the minimization holds.In fact, if for α > 0 the existence of the ground state is proved in H 1 (R, C), for α < 0 the same result is shown only on the subspace of the even functions of H 1 (R, C), namely H 1 r .At the same time, also the stability of the ground state changes if the problem is set in the attractive case or the repulsive one, gaining more stability in the former, as depicted in the following theorems.
We note that the value ω = α 2 4 corresponds to the frequency of the linear ground state in the attractive case and in particular it represents the threshold after which we can observe the presence of stationary states for the NLSE with delta conditions at the origin.Another remark is that, while in Theorem 2.1 the stability and the instability outcomes hold in H 1 , in Theorem 2.2 only the instability results are valid in H 1 and for the stability ones authors restrict to H 1 r .
These results have been generalized on a star graphs S N in [2], where the search for stationary states has been still conducted both in the attractive and in the repulsive regime.However, the ground state has been identified with the N tail state and characterized as the minimizer of a constrained action functional only in presence of a strong attractive interaction α * .In addition, it has been proved that it is stable in the subcritical and critical regime.
A second family of linear non-Kirchhoff's conditions are the so-called delta prime conditions, introduced in the previous section and defined in (11).In [12], we can find a deep investigation about the existence and the orbital stability of ground states using the constrained action functional In this work, the authors prove that there exists a critical value ω * for which an interesting bifurcation result occurs.In particular it follows that Theorem 2.3 (Theorem 5.3, Proposition 6.11 and Theorem 6.13 in [12]).Let β > 0.
Then, there exists ω * = 4 , then there exists a unique ground state, which is odd and orbitally stable for any p ∈ (2,6].
• If ω ≥ ω * , then there exist two non-symmetric ground states that are stable if the power nonlinearity does not exceed a critical value p * > 6 and become unstable for p > p * .The branch of odd solutions continues to exist at any ω > ω * , but they become a family of orbitally unstable stationary states.
As it appears from the results, the delta prime conditions give rise to a much richer structure of the family of ground states, including a pitchfork bifurcation with symmetry breaking.In fact, for frequency higher than ω * , the ground states display no symmetry, making not possible to reduce the problem to the halfline (contrarily to what happens in the case of a delta, where all ground states are even functions).This higher level of complexity of the whole picture arises from the fact that the energy space is larger, including functions with arbitrary jumps, and no relationship between the positive and negative halflines.Such a connection is restored by the interacting term of the energy.More recently, an other type of conditions characterized by a discontinuity has been studied on the real line.They arise from a particular self-adjoint extension of the 1D Laplacian and they are called Fülöp-Tsutsui's conditions as in [31].They are defined as ϕ(0 where τ ∈ R\{0, 1} and v > 0. Roughly speaking they can be seen as weighted delta conditions that allow discontinuities at the origin.Some studies have been conducted on these conditions [31], but to the knowledge of the authors, up to now no investigations concerning the existence and the stability of the ground states have been done.To fill the gap and give a more complete review, in the following we present some results to appear [14] obtained studying the minimization problem for the action functional The following result proves the existence of ground states for the previous constrained functional Even if the proof is quite standard and exploits Banach-Alaoglu's theorem and Brezis-Lieb's lemma in order to obtain a convergence results for the minimizing sequences, a crucial role is played by the following result that allows us to study an equivalent problem.
In particular, the two minimization problems are equivalent.
Even though these conditions prescribe a discontinuity at the origin, they share the same qualitative behaviour of delta interactions for what concerns the orbital stability of ground states.In particular, since we consider an attractive interaction (v > 0), we get the following result, analogous to Theorem 2.1 valid for classical delta conditions.More precisely: then the ground state is stable for any ω ∈ v 2 (τ 2 +1) 2 , +∞ .
To conclude, we want to remark that, while the proof of the stability result for p ∈ (2, 6] has been conducted relying on the Grillakis-Shatah-Strauss' theory, the stability conjecture for p > 6 leans on both numerical simulations and the asymptotical analysis of the behaviour of the L 2 norm of the ground state for ω ∈

Nonlinear delta conditions
Recently, we started studying the nonlinear Schrödinger equation with attractive nonlinear delta interactions at the vertices of the graph.More precisely, starting from the case of the attractive linear delta interaction, it is natural to generalize (9) replacing the positive number α by |ϕ(v)| q−2 with q > 2, getting the condition Such a condition generalizes the model introduced and studied in [18], where the effects of nonlinear point interactions are treated, then extended to the three-dimensional setting in [7,8], and only recently to space dimension two in [25,5,6].Such models were introduced in order to collect several results coming from theoretical physics and to include concentrated nonlinearities in a new class of mathematically rigorous models.It is worth recalling the application of such models to resonant tunneling [43,49].An immediate remark is that, differently from all the conditions presented before, ( 14) is nonlinear and does not follow from any self-adjoint extension H of the Laplace operator (3).
As for the linear non-Kirchhoff's conditions, the problem of existence of ground states has been studied only when G = R (see [22]) or when G is a star graph with N halflines [1] and not for more general metric graphs yet.In particular, we have looked for global minimizers of the energy functional among all the continuous functions u ∈ H 1 (G) satisfying the mass constraint: it is immediate to show that ground states of (15), if they exist, are solutions of the stationary equation ( 5) on each edge, are continuous at the vertices by definition and fulfill (14).As anticipated in the Introduction, in (15) we have coexistence of two nonlinearities, the standard and the pointwise one.We denote by F p,q : [0, +∞) → [−∞, +∞) the function defined as Due to the presence of two nonlinearities p and q and the validity of Gagliardo-Nirenberg inequalities, if 2 < p < 6 and 2 < q < 4, then the functional F p,q is bounded from below on H 1 µ (G) for every µ > 0 and we are in the so-called subcritical case.When the boundedness from below of F p,q depends on the value of the mass µ, i.e. in the so-called critical cases, it is important to distinguish the single critical case in which p = 6 and 2 < q < 4 or 2 < p < 6 and q = 4 and the doubly critical case in which p = 6 and q = 4.

Subcritical case
For what concerns the subcritical case, existence results on R and on S N are very different since sufficient conditions for the compactness of minimizing sequences drastically change.Indeed, one can prove that on R inf F p,q (u, R) < 0 ⇒ a ground state of F p,q (•, R) exists and consequently the following theorem holds.
As one can observe, there is no significant interaction between the two nonlinear terms of the energy since the existence result is the same when the delta interaction is linear.
On the other hand, if G = S N , then First, notice that, in order to have ground states, the energy F p,q (•, S N ) has to be smaller than the energy E(•, R) of the soliton and not only than 0.More precisely, it can be shown that the existence of a ground state of F p,q (•, S N ) is equivalent to the existence of a function u ∈ H 1 µ (S N ) such that F p,q (u, S N ) ≤ E(ϕ µ ).Existence and non-existence results are obtained taking advantage of this equivalence and the interplay between pointwise and standard nonlinearities becomes evident.In particular, if q < p 2 + 1, then existence of ground states holds for small masses and does not hold for large masses and this behaviour is similar to the one described in [3] in the case of a linear delta interaction at the origin: this suggests that when the nonlinear delta interaction is not too strong, then it has qualitatively the same effect as the linear delta on existence results.If instead q > p 2 , then ground states exist for large masses and do not exist for small masses.In both the cases just described, the passage from existence of ground states to non-existence or viceversa identifies a unique threshold value of the mass µ * which varies depending on the two powers p and q and on the number of halflines N .When the two nonlinearities are in a perfect balance, that is the case in which q = p 2 + 1, existence and non-existence of ground states depend only on the number of halflines of the star graph and not on the value of the mass.These results are summarized in the following theorem, which is in preparation [1].

Critical cases
Existence of ground states in critical cases has been studied only when G = R.The first theorem deals with the cases in which only one power is critical.Theorem 3.3 (Theorem 1.4 in [22]).Let µ > 0.
(i) If p = 6 and 2 < q < 4, then there exists a unique positive ground state at mass µ if and only if µ < √ 3 2 π, and (ii) If 2 < p < 6 and q = 4, then there exists a unique positive ground state at mass µ if and only if µ < 2 These particular regimes show the interplay between a subcritical and a critical power nonlinearity: indeed, while the ground state level moves on from 0 to −∞ in correspondence of a threshold value of the mass as usual in critical cases, the presence of the subcritical power ensures existence of ground states for all the masses under the critical mass, highlighting an important difference with what one expects in critical cases.
The last result concerns the doubly critical case, where simultaneously p = 6 and q = 4.Here we recover the typical structure of a purely critical setting, with the ground state energy level lifting from 0 to −∞ when exceeding a critical value of the mass and solutions existing only at the threshold.A quite remarkable feature due to the interaction between the two nonlinearities is given by the fact that the critical mass (19) is lower than the critical masses √ 3 2 π and 2 for the standard and pointwise nonlinearity.

Conclusion and perspectives
We reviewed on recent developments and announced new results on the the existence of Ground States for the Nonlinear Schrödinger Equation on infinite star graphs, stressing possible applications to the description of the behaviour of Bose-Einstein condensates.From the physical point of view we are still at a preliminary step, in which the correct conditions are to be singled out, but there are hints that the correct conditions could not be free, or Kirchhoff's.On the other hand, such a topic has already shown a very rich potential from a mathematical point of view, underlying how the interplay between a point interaction and a standard power nonlinearity can give rise to phenomena that are new, in the sense that they do not occur neither in the standard Nonlinear Schrödinger equation, nor in the linear one with point interactions.
The most spectacular among them is the occurrence of bifurcations, for which a non-trivial energy space is needed, a requirement that is fulfilled by specific (in particular, discontinuous) boundary conditions at the vertex, e.g.delta prime or Fülöp-Tsutsui.Further phenomena to be still completely investigated arise when the vertex interaction is nonlinear too, so that there is competition between the two nonlinearities and possibly coexistence of two criticalities.In our results it is shown how this competition is ruled, and in which cases one nonlinearity prevails on the others.
As it is natural, a prosecution of the research will be concerned with the deepening of the resarch in this direction.Another one, completely to be explored, will deal with considering more-dimensional analogues of the graphs, the so-called hybrids, with components of different dimensions to be glued together.Preliminary discussions and formal computations made on this, seem to show that the phenomenology given by the nonlinear hybrids could be still richer than the one exhibited by nonlinear graphs.