Stationary Nonlinear Schr"odinger Equation on Simplest Graphs: Boundary conditions and exact solutions

We treat the stationary (cubic) nonlinear Schr"odinger equation (NSLE) on simplest graphs. Formulation of the problem and exact analytical solutions of NLSE are presented for star graphs consisting of three bonds. It is shown that the method can be extended for the case of arbitrary number of bonds of star graphs and for other simplest topologies such as tree and loop graphs. The case of repulsive and attractive nonlinearities are treated separately.

💡 Research Summary

The paper addresses the stationary (time‑independent) cubic nonlinear Schrödinger equation (NLSE) on the simplest one‑dimensional quantum graphs. Starting with the most elementary topology—a star graph composed of three finite bonds meeting at a single vertex—the authors formulate the NLSE on each bond as
‑ψ″_j + β_j|ψ_j|²ψ_j = λ²ψ_j, β_j > 0, j = 1,2,3,
and impose physically motivated boundary conditions at the external ends (Dirichlet ψ = 0) and at the central vertex. At the vertex they require continuity of the wave function (up to scaling factors A₂, A₃) and conservation of the probability current, which together constitute a δ′′‑type condition:

ψ₁(L₁) = A₂ψ₂(L₁) = A₃ψ₃(L₁),
(1/√β₁)ψ′₁(L₁) − (1/√β₂)ψ′₂(L₁) − (1/√β₃)ψ′₃(L₁) = 0.

The authors separate the phase by writing ψ_j(x)=e^{iγ_j}f_j(x) and show that the phases must be equal constants, reducing the problem to a real ordinary differential equation for f_j(x). The general solution of this ODE on a finite interval is expressed through Jacobi elliptic functions (sn, cn, dn). For the repulsive (β > 0) case the solution is taken as f_j(x)=B_j sn(α_j x+δ_j|k_j); for the attractive case (β < 0) the cn function is used. The parameters B_j, α_j, k_j, and the shift δ_j are linked by the NLSE itself (relations (9) in the paper) and by the normalization condition ∑_j∫|ψ_j|²dx = 1.

Applying the vertex conditions yields a coupled nonlinear algebraic system (Eqs. 10–12) involving β_j, bond lengths L_j, and the elliptic parameters. In the general situation this system must be solved numerically (e.g., Newton or Krylov iterations). However, the authors identify two analytically tractable families of solutions. In the first family the three bonds share identical elliptic modulus (k₁ = k₂ = k₃ = k) and the scaling parameters α_j are chosen such that the arguments of the sn functions at the vertex are integer multiples of the complete elliptic integral K(k). This reduces the whole system to a single scalar equation g(k)=0 (Eq. 13). By evaluating the limits g(0)=−1 and g(1)=+∞ and noting continuity on (0,1), the existence of at least one root is guaranteed. The second family allows opposite signs for the amplitudes on two bonds (σ₂=σ₃=−1) and imposes a linear relation among the β_j coefficients; again a single equation g(k)=0 (Eq. 14) emerges with the same limiting behavior, ensuring a solution.

The attractive‑nonlinearity case is treated analogously. The solution ansatz uses the cn function, the phase shifts δ_j are fixed by the Dirichlet conditions at the bond ends, and the same type of algebraic system (Eqs. 17–20) is derived. Two special parameter sets are exhibited, leading to scalar equations (23) and (25) whose solvability follows from the same continuity argument.

Having established exact solutions for the star graph, the authors then outline how the method extends to other elementary topologies: general star graphs with N bonds, tree graphs, and loop (ring) graphs. The key is to impose δ′′‑type conditions at every internal vertex while keeping ψ = 0 at all external termini. On each bond the same elliptic‑function ansatz applies, and the matching conditions at each vertex generate a larger but structurally similar algebraic system. The paper sketches the procedure for a simple tree (Fig. 2) and indicates that loop graphs require periodicity constraints in addition to the vertex conditions.

From a physical standpoint, the work provides the first explicit analytical stationary states of the NLSE on finite‑length quantum graphs with nontrivial topology. The results are directly relevant to Bose‑Einstein condensates confined in network‑shaped traps, to optical soliton propagation in waveguide lattices, and to models of energy transport in molecular chains (e.g., DNA). The existence of exact solutions demonstrates that nonlinear wave transmission, reflection, and resonance phenomena on graphs can be analyzed without resorting solely to numerical simulations. Moreover, the dependence of the solutions on the bond lengths, nonlinearity strengths β_j, and the choice of vertex coupling (δ′′ versus other types) offers a rich parameter space for engineering desired transport properties.

In conclusion, the paper makes a substantial contribution by (i) formulating the stationary NLSE on simple graphs with rigorous vertex boundary conditions, (ii) deriving exact elliptic‑function solutions for both repulsive and attractive nonlinearities, (iii) proving the existence of solutions for specific parameter families, and (iv) outlining a systematic extension to more complex graph topologies. Future work could address the time‑dependent NLSE on graphs, inclusion of external potentials, higher‑order nonlinearities, and the impact of disorder or loss, thereby broadening the applicability of the present analytical framework.