pith. sign in

arxiv: 2511.14519 · v1 · pith:LDMMWDGZnew · submitted 2025-11-18 · 🪐 quant-ph · nlin.SI

Perturbative nonlinear J-matrix method of scattering in two dimensions

T. J. Taiwo , A. D. Alhaidari , U. Al Khawaja This is my paper

Pith reviewed 2026-05-17 20:51 UTC · model grok-4.3

classification 🪐 quant-ph nlin.SI
keywords nonlinear Schrödinger equationJ-matrix methodtwo-dimensional scatteringcircular symmetrybifurcationperturbative methodorthogonal polynomialsGauss quadrature
0
0 comments X

The pith

A perturbative nonlinear J-matrix method yields the scattering matrix for the two-dimensional nonlinear Schrödinger equation with circular symmetry and detects energy-dependent bifurcations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a perturbative formulation that extends the J-matrix method to treat nonlinear scattering in two dimensions. It applies this to the time-independent nonlinear Schrödinger equation with circular symmetry for general odd-powered nonlinearities. Linearization of products of orthogonal polynomials combined with Gauss quadrature enables computation of the scattering matrix for cubic and quintic cases. At specific energies the solutions bifurcate, producing two stable branches as a direct signature of the nonlinearity. If the method holds, it supplies a practical route to scattering data in nonlinear quantum problems where exact analytic solutions are unavailable.

Core claim

We obtain the scattering matrix for the time-independent nonlinear Schrödinger equation in two dimensions with circular symmetry using a perturbative nonlinear extension of the J-matrix method. The formulation relies on the linearization of products of orthogonal polynomials and on the tools of the J-matrix method. Gauss quadrature integral approximation is used in the numerical implementation. Results are given for the cubic and quintic nonlinearities, and at certain values of the energy bifurcation occurs with two stable solutions.

What carries the argument

Perturbative linearization of products of orthogonal polynomials inside the J-matrix basis, with Gauss quadrature for the resulting integrals.

If this is right

  • The scattering matrix is obtained for a general ψ^{2n+1} nonlinearity.
  • Explicit results follow for the cubic (ψ^{3}) and quintic (ψ^{5}) cases.
  • Bifurcation into two stable solutions appears at certain discrete energy values.
  • The method preserves the J-matrix treatment of asymptotic boundary conditions while adding nonlinearity perturbatively.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same linearization step could be checked by comparing bifurcation thresholds against direct numerical solutions of the nonlinear wave equation at increasing nonlinearity strength.
  • The observed bifurcation may connect to multistability phenomena in other two-dimensional nonlinear wave systems such as optical or condensate models.
  • Extending the perturbative order or testing higher odd nonlinearities would show whether the bifurcation pattern persists or changes character.

Load-bearing premise

Linearizing products of orthogonal polynomials under the perturbative treatment stays accurate enough to capture the essential nonlinear effects without introducing uncontrolled errors in the scattering matrix or the bifurcation points.

What would settle it

An independent numerical integration of the nonlinear Schrödinger equation for a chosen circular potential and nonlinearity strength that produces bifurcation energies or scattering-matrix elements differing from those computed by the perturbative J-matrix method.

read the original abstract

We introduce a perturbative formulation for a nonlinear extension of the J-matrix method of scattering in two dimensions. That is, we obtain the scattering matrix for the time-independent nonlinear Schr\"odinger equation in two dimensions with circular symmetry. The formulation relies on the linearization of products of orthogonal polynomials and on the utilization of the tools of the J-matrix method. Gauss quadrature integral approximation is instrumental in the numerical implementation of the approach. We present the theory for a general \psi ^{2n + 1} nonlinearity, where n is a natural number, and obtain results for the cubic and quintic nonlinearities, \psi ^3 and \psi ^5. At certain value(s) of the energy, we observe the occurrence of bifurcation with two stable solutions. This curious and interesting phenomenon is a clear signature and manifestation of the underlying nonlinearity.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper introduces a perturbative nonlinear extension of the J-matrix method for computing the scattering matrix of the time-independent nonlinear Schrödinger equation in two dimensions with circular symmetry. The approach linearizes products of orthogonal polynomials in the nonlinear term (for general ψ^{2n+1}, with explicit results for cubic ψ^3 and quintic ψ^5 nonlinearities), employs Gauss quadrature for the resulting integrals, and reports the appearance of bifurcations yielding two stable solutions at specific energies.

Significance. If the perturbative linearization and quadrature steps prove accurate, the work would provide a practical extension of established J-matrix techniques to nonlinear scattering, potentially useful for exploring bifurcation phenomena in 2D NLSE problems. The reliance on standard orthogonal-polynomial properties and the J-matrix framework is a natural starting point, but the lack of verification limits the strength of the contribution.

major comments (2)
  1. [Numerical implementation and results sections] The central numerical results on bifurcations rest on the perturbative linearization of polynomial products and Gauss quadrature without any reported error analysis, convergence tests with basis size, or comparison to independent benchmarks (e.g., direct numerical integration of the NLSE or known linear limits). This directly affects the reliability of the claimed bifurcation energies and stability of the two solutions.
  2. [Formulation of the perturbative nonlinear J-matrix method] The linearization step for the nonlinear term (e.g., products arising from ψ^3 or ψ^5) is applied perturbatively, yet no estimate is given for the regime where this approximation remains controlled relative to the linear terms; at the finite amplitudes where bifurcations occur, truncation errors could shift or artifactually create the reported bifurcation points.
minor comments (2)
  1. [Results on bifurcations] Clarify in the text how 'stable solutions' are identified numerically (e.g., via eigenvalue analysis of the linearized operator or iteration convergence).
  2. [Theory section] The abstract states results for general ψ^{2n+1} but the main text focuses on n=1 and n=2; a brief outline of the general case would improve completeness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We agree that additional numerical validation and an assessment of the perturbative regime would strengthen the presentation. We address the major comments point by point below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Numerical implementation and results sections] The central numerical results on bifurcations rest on the perturbative linearization of polynomial products and Gauss quadrature without any reported error analysis, convergence tests with basis size, or comparison to independent benchmarks (e.g., direct numerical integration of the NLSE or known linear limits). This directly affects the reliability of the claimed bifurcation energies and stability of the two solutions.

    Authors: We acknowledge that the current manuscript does not present explicit error analysis, convergence tests with basis size, or direct comparisons to independent numerical solutions of the NLSE. In the revised version we will add a dedicated subsection on numerical convergence, showing how the scattering-matrix elements and bifurcation energies stabilize as the basis dimension and quadrature order are increased. We will also include the linear limit (vanishing nonlinearity strength) as a benchmark, recovering the known results of the standard J-matrix method. While a full direct numerical integration of the two-dimensional NLSE is computationally demanding, we will provide a limited comparison for weak nonlinearities to support the perturbative results. These additions will improve the reliability assessment of the reported bifurcation points. revision: yes

  2. Referee: [Formulation of the perturbative nonlinear J-matrix method] The linearization step for the nonlinear term (e.g., products arising from ψ^3 or ψ^5) is applied perturbatively, yet no estimate is given for the regime where this approximation remains controlled relative to the linear terms; at the finite amplitudes where bifurcations occur, truncation errors could shift or artifactually create the reported bifurcation points.

    Authors: The linearization of the nonlinear term is performed within the J-matrix framework by expressing products of radial functions in the chosen orthogonal-polynomial basis. We agree that an explicit estimate of the validity range is missing. In the revision we will include a brief error analysis that compares the magnitude of the retained linear terms to the neglected higher-order contributions in the nonlinear potential, thereby delineating the amplitude regime where the approximation is controlled. We will also evaluate the wave-function amplitudes at the observed bifurcation energies to confirm that they remain consistent with the perturbative ordering. This discussion will clarify the conditions under which the reported bifurcations are expected to be accurate. revision: yes

Circularity Check

0 steps flagged

Minor self-citation of J-matrix foundations; central perturbative extension remains independent

full rationale

The paper constructs a perturbative nonlinear extension of the J-matrix method by linearizing products of orthogonal polynomials and applying Gauss quadrature to the NLSE nonlinearity. The scattering matrix and observed bifurcations follow directly from this construction applied to the time-independent NLSE with circular symmetry. No step reduces a claimed prediction or result to a fitted parameter or self-citation by definition. Prior J-matrix work by co-author Alhaidari is cited for the linear foundation, but this is standard and does not bear the load of the nonlinear results or bifurcation findings, which are generated within the present perturbative framework.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of linearizing products of orthogonal polynomials for the perturbative treatment and on the accuracy of Gauss quadrature for the resulting integrals; both are standard numerical techniques rather than new postulates.

axioms (2)
  • domain assumption Products of orthogonal polynomials can be linearized to sufficient accuracy for perturbative treatment of the nonlinear term
    Invoked to convert the nonlinear Schrödinger equation into a form compatible with the J-matrix expansion.
  • domain assumption Gauss quadrature provides a sufficiently accurate approximation to the integrals arising after linearization
    Used for numerical implementation of the scattering-matrix calculation.

pith-pipeline@v0.9.0 · 5451 in / 1368 out tokens · 50248 ms · 2026-05-17T20:51:17.204558+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [1]

    G. P. Agrawal, Nonlinear Fiber Optics, 6th ed. (San Diego: Academic Press, 2019)

    work page 2019

  2. [2]

    S. B. Pope, Turbulent Flows (Cambridge: Cambridge University Press, 2000)

    work page 2000

  3. [3]

    Kono and M

    M. Kono and M. M. Skoric, Nonlinear Physics of Plasmas (Berlin: Springer, 2010)

    work page 2010

  4. [4]

    R. C. Hilborn, Chaos and Nonlinear Dynamics: An Introduction for Scientists and Engineers, 2nd ed. (Oxford: Oxford University Press, 2000)

    work page 2000

  5. [5]

    P. G. Kevrekidis, D . J. Frantzeskakis, and R . Carretero-González (Eds.) Emergent Nonlinear Phenomena in Bose -Einstein Condensates: Theory and Experiment (Berlin: Springer, 2008)

    work page 2008

  6. [6]

    Carballido-Landeira and B

    J. Carballido-Landeira and B. Escribano (Eds.) Biological Systems: Nonlinear Dynamics Approach, SEMA SIMAI, Vol. 20 (Cham: Springer, 2019)

    work page 2019

  7. [7]

    S. H. Strogatz, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering, 2nd ed. (Boulder, CO: Westview Press, 2015)

    work page 2015

  8. [8]

    Zheng, Nonlinear Evolution Equations (Boca Raton, FL: Chapman and Hall/CRC, 2004)

    S. Zheng, Nonlinear Evolution Equations (Boca Raton, FL: Chapman and Hall/CRC, 2004)

    work page 2004

  9. [9]

    Al Sakkaf and U

    L. Al Sakkaf and U . Al Khawaja, Superposition Principle and Composite Solutions to Coupled Nonlinear Schrödinger Equations, Math. Meth. Appl. Sci. 43 (2020) 10841

    work page 2020

  10. [10]

    C. J. Pethick and H. Smith , Bose–Einstein Condensation in Dilute Gases . 2nd ed. (Cambridge: Cambridge University Press, 2008)

    work page 2008

  11. [11]

    Kharif, E

    C. Kharif, E. Pelinovsky, and A. Slunyaev, Rogue Waves in the Ocean (Berlin: Springer, 2009)

    work page 2009

  12. [12]

    V. E. Zakharov and A. B. Shabat , Exact Theory of Two-Dimensional Self-Focusing and One-Dimensional Self-Modulation of Waves in Nonlinear Media, JETP 34 (1972) 62

    work page 1972

  13. [13]

    Al Khawaja and L

    U. Al Khawaja and L . Al Sakkaf , Handbook of Exact Solutions to the Nonlinear Schrödinger Equations, 2nd ed. (Bristol: IOP Publishing, 2024). −25−

    work page 2024

  14. [14]

    Almualla, L

    Z. Almualla, L . Al Sakkaf, and U . Al Khawaja , Resonant Soliton Scattering by Topological States of a Finite Array of Reflectionless Potential Wells , Phys. Lett. A 552 (2025) 130646

    work page 2025

  15. [15]

    A. D. Alhaidari, E. J. Heller, H. A. Yamani, and M. S. Abdelmonem (Eds.), The J-Matrix Method: Developments and Applications (Springer, Dordrecht, Netherlands, 2008)

    work page 2008

  16. [16]

    A. D. Alhaidari and T. J. Taiwo, Nonlinear extension of the J-matrix method of scattering: A toy model, Arab. J. Math. 14 (2025) https://doi.org/10.1007/s40065-025-00535-x

    work page doi:10.1007/s40065-025-00535-x 2025

  17. [17]

    E. J. Heller and H. A. Yamani, New L2 approach to quantum scattering: Theory , Phys. Rev. A 9 (1974) 1201

    work page 1974

  18. [18]

    E. J. Heller and H. A. Yamani, J-matrix method: Application to S-wave electron-hydrogen scattering, Phys. Rev. A 9 (1974) 1209

    work page 1974

  19. [19]

    H. A. Yamani and L. Fishman, J-matrix method: Extensions to arbitrary angular momentum and to Coulomb scattering, J. Math. Phys. 16 (1975) 410

    work page 1975

  20. [20]

    M. E. H. Ismail and E. Koelink, The J-matrix method, Adv. Appl. Math. 46 (2011) 379

    work page 2011

  21. [21]

    A. D. Alhaidari, Exact and simple formulas for the linearization coefficients of products of orthogonal polynomials and physical application, J. Comput. Appl. Math. 436 (2024) 115368

    work page 2024

  22. [22]

    M. E. H. Ismail, Classical and Quantum orthogonal polynomials in one variable (Cambridge University press, 2009)

    work page 2009

  23. [23]

    A. D. Alhaidari, Gauss Quadrature for Integrals and Sums, Int. J. Pure Appl. Math. Res. 3 (2023) 1

    work page 2023

  24. [24]

    H. M. Srivastava, H. M. Mavromatis and R. S. Alassar, Remarks on Some Associated Laguerre Integral Results, Appl. Math. Lett. 16 (2003) 1131

    work page 2003

  25. [25]

    Y. K. Ho, The method of complex coordinate rotation and its applications to atomic collision processes, Phys. Rep. 99 (1983) 1

    work page 1983

  26. [26]

    I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products , 7th edition (Academic, New York, 1980)

    work page 1980

  27. [27]

    Magnus, F

    W. Magnus, F. Oberhettinger and R. P. Soni, Formulas and theorems for the special functions of mathematical physics, 3rd ed. (Springer, 1966)

    work page 1966

  28. [28]

    A. D. Alhaidari, Matrix representation of the resolvent operator in square -integrable basis and physical application, J. Appl. Math. Computation 8 (2024) 296 −26− Figure Captions Fig. 1: Plot of the linear scattering matrix, 02i ( )1 Ee − , as a function of the energy associated with the potential 2( ) 7.5 rV r r e −= for 0== m using the linear J -matri...

    work page 2024