pith. sign in

arxiv: 2605.16125 · v1 · pith:TCPIKRQ2new · submitted 2026-05-15 · ❄️ cond-mat.stat-mech · nlin.CD· physics.flu-dyn

Staggering domino-like blast front motion in a one-dimensional cold gas

Pith reviewed 2026-05-19 18:42 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech nlin.CDphysics.flu-dyn
keywords one-dimensional gaselastic collisionsblast frontshock waveexact solutionmass ratiodomino motionalternating masses
0
0 comments X

The pith

For an infinite family of special mass ratios, a one-dimensional alternating particle chain with equidistant spacing admits an exact solution where the blast front advances ballistically because only one triplet of particles moves at any时刻.

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

The paper studies elastic collisions in a one-dimensional chain of alternating heavy and light particles started by kicking the leftmost particle. For a mass ratio of 2 and random initial positions the motion follows standard hydrodynamic predictions, including a decelerating shock front. For an infinite set of discrete mass ratios and exactly regular initial positions the dynamics instead collapse to a precise staggered pattern. In this pattern only three consecutive particles ever move while the rest remain at rest, so the front travels at constant average speed equal to the initial velocity. Explicit formulas for the allowed mass ratios, velocities, and positions follow directly from the construction.

Core claim

For an infinite family of mass ratios M_k the dynamics admit an exact solution in which at each moment only a single triplet m, μ, m is in motion, all other particles are at rest, and the shock front moves ballistically with average velocity equal to the initial one. The solution yields explicit formulas for the sequence M_k together with the time-dependent velocities and positions of every particle.

What carries the argument

The staggering triplet mechanism in which collisions propagate sequentially through isolated groups of three particles while all others remain motionless.

Load-bearing premise

The initial positions must be exactly equidistant and the mass ratio must be tuned precisely to one of the discrete values M_k; any deviation from these exact conditions destroys the staggered-triplet regime.

What would settle it

A direct numerical integration with mass ratio exactly equal to one computed M_k and perfectly equidistant initial positions should show only one triplet moving at any time and a front advancing at constant average speed; introducing even tiny random displacements in the initial positions should immediately restore the usual decelerating hydrodynamic front.

Figures

Figures reproduced from arXiv: 2605.16125 by Krzysztof Pilorz, Taras Holovatch, Yurij Holovatch, Yuri Kozitsky.

Figure 1
Figure 1. Figure 1: FIG. 1. Snapshots of the particle motion for [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Time dependence of [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Time dependence of the Shannon entropy [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Mass [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
read the original abstract

One-dimensional alternating particle systems are widely used to study interconnections between the hydrodynamics of blast waves in a gas-like medium and the Newtonian dynamics of its corpuscular constituents. We study the model in which point particles with masses $m,\mu, m,\mu,\dots, (m\geq\mu)$ are distributed on the positive half-line $\mathbb{R}_{+}$. Their dynamics are initiated by giving a positive velocity to the leftmost particle; in its course, the particles undergo elastic collisions. For this model with $m/\mu=2$, it has previously been established that the dynamics that start from random initial positions are consistent with predictions based on Euler's hydrodynamic equation. In particular, they have the following properties: (i) the position of the rightmost particle (shock front) evolves as $t^\delta$ with $\delta<1$; (ii) recoiled particles behind the front enter the negative half-axis; (iii) particles with locations $x\leq0$ move ballistically and eventually take over the total energy of the system. In this paper, we present numerical and analytical results for the dynamics of this model with nonrandom (typically equidistant) initial positions and various values of $m/\mu$. For $m/\mu=2$ and equidistant initial positions, our results qualitatively agree with those just mentioned. At the same time, we found an infinite family of numbers $\{\mathcal{M}_k\}$ such that, for $m/\mu=\mathcal{M}_k$, the hydrodynamic behavior mentioned changes drastically to the following. At each moment, only a single triplet $m,\mu, m$ is in motion, whereas all other particles are at rest. As a result, the shock front moves ballistically with an average velocity equal to the initial one. Such a `staggering domino-like' picture is obtained as an exact solution, which yields, in particular, explicit formulas for $\mathcal{M}_k$ and the particle velocities and positions.

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 examines one-dimensional elastic collisions in an alternating-mass system (m, μ, m, μ, ...) on the positive half-line, initiated by imparting positive velocity to the leftmost particle. For generic mass ratios and random initial positions the dynamics are consistent with Euler hydrodynamics (sub-ballistic shock front, recoils into x<0, energy takeover by leftward particles). The central claim is that for a countable family of special ratios M_k the dynamics admit an exact closed-form solution in which motion is confined at all times to a single isolated triplet m-μ-m while every other particle remains at rest; the resulting shock front therefore advances ballistically with the same average speed as the initial impulse. Explicit algebraic expressions for the sequence M_k together with the successive velocities and positions are derived from the two-body collision map.

Significance. If the construction is correct, the result supplies a rare, fully explicit solvable case inside a model whose generic behavior is hydrodynamic and non-integrable. It demonstrates that the staggered-triplet regime is an exact, non-perturbative solution for discrete mass ratios under perfectly equidistant initial data, thereby furnishing a concrete benchmark against which numerical schemes and hydrodynamic approximations can be tested. The explicit formulas also make the dependence of the front velocity on the mass ratio immediately falsifiable.

major comments (2)
  1. [derivation of M_k and propagation argument] The manuscript asserts that the staggered pattern propagates indefinitely without out-of-sequence collisions. This claim is load-bearing for the exact-solution statement; an explicit inductive step or closed-form verification that the arrival time of the next triplet always precedes any possible crossing collision should be supplied (e.g., in the section deriving the recursion for M_k).
  2. [explicit formulas for M_k] The explicit formula for M_k is obtained by solving the finite collision map at each step. It is not immediately clear whether the resulting algebraic expression for general k remains positive and greater than unity for all k, which is required for physical masses; a short table or asymptotic analysis of M_k for large k would confirm this.
minor comments (2)
  1. [numerical results for m/μ = 2] The abstract states that for m/μ = 2 and equidistant positions the results 'qualitatively agree' with the hydrodynamic picture; a brief quantitative comparison (e.g., measured exponent δ versus the expected hydrodynamic value) would strengthen the contrast with the M_k cases.
  2. [velocity and position formulas] Notation for the successive velocities after each triplet collision is introduced without a compact summary table; adding such a table would improve readability of the exact solution.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and the recommendation of minor revision. The two major comments identify points where additional explicit verification will improve clarity; we address each below and will incorporate the requested material in the revised manuscript.

read point-by-point responses
  1. Referee: [derivation of M_k and propagation argument] The manuscript asserts that the staggered pattern propagates indefinitely without out-of-sequence collisions. This claim is load-bearing for the exact-solution statement; an explicit inductive step or closed-form verification that the arrival time of the next triplet always precedes any possible crossing collision should be supplied (e.g., in the section deriving the recursion for M_k).

    Authors: We agree that an explicit inductive argument strengthens the presentation. In the revised version we will insert a new subsection immediately after the recursion for M_k. The base case (k=1) is verified by direct computation of the three-particle collision sequence. For the inductive step we show that if the pattern holds up to triplet k with mass ratio M_k, then the arrival time of the (k+1) triplet at the next stationary particle is strictly earlier than any possible crossing collision; the ordering follows from the velocity inequalities preserved by the two-body map and the recursive definition of M_{k+1}. revision: yes

  2. Referee: [explicit formulas for M_k] The explicit formula for M_k is obtained by solving the finite collision map at each step. It is not immediately clear whether the resulting algebraic expression for general k remains positive and greater than unity for all k, which is required for physical masses; a short table or asymptotic analysis of M_k for large k would confirm this.

    Authors: We thank the referee for highlighting this point. Although the recursive construction ensures M_k > 1 by design, we will add both a short table listing the first ten explicit values of M_k (all >1) and a brief asymptotic analysis showing that M_k = 1 + O(1/k) as k → ∞, thereby confirming positivity and physical admissibility for every k. revision: yes

Circularity Check

0 steps flagged

No significant circularity; exact solution constructed directly from collision rules

full rationale

The paper derives an exact staggered-triplet solution for a countable family of mass ratios M_k by iteratively solving the finite elastic collision map under equidistant initial positions. Explicit algebraic formulas for M_k, velocities, and positions are obtained step-by-step from the Newtonian rules without parameter fitting, self-referential predictions, or load-bearing self-citations. The construction is shown to propagate consistently without generating out-of-sequence collisions, making the derivation self-contained against the stated initial conditions and hypotheses. No reduction of the central claim to its own inputs occurs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on elastic two-body collision rules, conservation of momentum and energy, and the assumption of perfectly equidistant initial positions. No additional free parameters or invented entities are introduced beyond the discrete mass ratios that are solved for.

axioms (2)
  • domain assumption Particles interact only via instantaneous elastic collisions; no other forces act between collisions.
    Standard assumption for the cold-gas model stated in the abstract.
  • domain assumption Initial positions are exactly equidistant (or otherwise deterministic and non-random).
    Explicitly contrasted with the random-initial-position case studied previously.

pith-pipeline@v0.9.0 · 5918 in / 1297 out tokens · 26529 ms · 2026-05-19T18:42:32.318419+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

52 extracted references · 52 canonical work pages

  1. [1]

    The exponents corresponding to other observables, obtained in a similar way, are also listed in Table I

    using scaling arguments applied to the corresponding hydrodynam ic equation and then confirmed by molecular (Newtonian) dynamics calculations with random initial positions of the particles. The exponents corresponding to other observables, obtained in a similar way, are also listed in Table I. Their physical significance is discussed below. In this work, we...

  2. [2]

    Note that H(t) = 0 for all t > 0 for the degenerate domino dynamics

    to secure equality ∑ l pl(t) = 1, which is standard in this case. Note that H(t) = 0 for all t > 0 for the degenerate domino dynamics. Hence, the higher the value of H(t), the more uniform the distribution of the initial energy between the particles. IV. MOLECULAR DYNAMICS SIMULATIONS Here, we present and discuss the numerical calculation of the afor emen...

  3. [3]

    By combining the latter with (

    and ( 6) H(t) = − 2 m ∑ l≥ 0 mlu2 l (t) log |ul(t)|+ 2 m ∑ l:odd mlu2 l (t) ∼ tα . By combining the latter with (

  4. [4]

    Note that for m = 2, α is much smaller than γ, which means that the increase of the latter sum is much slower than that of ∑ l≥ 0 ml|ul(t)|

    we obtain ∑ l≥ 0 mlu2 l (t) = m, ∑ l≥ 0 ml|ul(t)| ≥ P x≥ 0(t) ∼ tγ , (11) ∑ l≥ 0 mlu2 l (t) log |ul(t)| ∼ tα . Note that for m = 2, α is much smaller than γ, which means that the increase of the latter sum is much slower than that of ∑ l≥ 0 ml|ul(t)|. 10 3 10 4 10 5 /u1D635 5 × 10 0 6 × 10 0 7 × 10 0  FIG. 3. Time dependence of the Shannon entropy H(t); ...

  5. [5]

    Then A := A1A0 =    θ 1 − θ 0 − θ(1 + θ) θ2 1 + θ 1 − θ2 − θ(1 − θ) θ    (16) 8 describes the complete round of collisions

    and ( 12), the velocities after the next round are    v0,k +1 v1,k +1 v2,k +1    =    1 0 0 0 − θ 1 + θ 0 1 − θ θ       θ 1 − θ 0 1 + θ − θ 0 0 0 1       v0,k v1,k v2,k    , which we abbreviate to the following formula Vk+1 = A1A0Vk, V k :=    v0,k v1,k v2,k    , (15) where A0 and A1 describe collisions 0 ↔ 1 and 1 ↔ 2, respe...

  6. [6]

    Then, in the Appendix B, we solve the recursion relations in ( 32) – ( 34)

    – ( 23). Then, in the Appendix B, we solve the recursion relations in ( 32) – ( 34). The first line of ( 13) allows one to exclude v2,k without losing the linearity of the recursion relation in ( 15). To this end, we introduce w0,k = v0,k , w 1,k = v1,k θ + 1 , (A1) and then get v2,k = 1 − w0,k − (1 − θ)w1,k . (A2) Then by ( 16) and ( A1), ( A2), we obtain...

  7. [7]

    (A25) By ( A23), for such α , v0,k = v1,k = 0, v 2,k = 1, and v0,s > 0, for all 2 ≤ s ≤ k − 1, (A26) cf

    (A24) Now, for a fixed k ≥ 2, let us find α ∈ [π/ 3, π/ 2) (hence m, see ( A24)) such that the following holds sin(2k + 1)α = 0, sin(2s − 1)α sin(2s + 1)α < 0, for all 2 ≤ s ≤ k − 1. (A25) By ( A23), for such α , v0,k = v1,k = 0, v 2,k = 1, and v0,s > 0, for all 2 ≤ s ≤ k − 1, (A26) cf. item ( A) of the staggering domino scenario. By the first condition in (...

  8. [8]

    Spohn, Kinetic equations from Hamiltonian dynamics: Markovian limits, Rev

    H. Spohn, Kinetic equations from Hamiltonian dynamics: Markovian limits, Rev. Mod. Phys. 52, 569 (1980)

  9. [9]

    J. L. Lebowitz, Microscopic origins of macroscopic behavior (2021), arXiv:2105.03470 [math-ph]

  10. [10]

    A. N. Gorban, Hilbert’s sixth problem: the endless road t o rigour, Philosophical Transactions of the Royal Society A: Mathema tical, Physical

  11. [11]

    Antal, P

    T. Antal, P. L. Krapivsky, and S. Redner, Exciting hard sp heres, Phys. Rev. E 78, 030301 (2008)

  12. [12]

    Chakraborti, A

    S. Chakraborti, A. Dhar, and P. L. Krapivsky, A splash in a one-dimensional cold gas, SciPost Phys. 13, 074 (2022)

  13. [13]

    Holovatch, Y

    T. Holovatch, Y. Kozitsky, K. Pilorz, and Y. Holovatch, B reakdown of hydrodynamics in a one-dimensional cold gas, Phys. Rev. E 112, L052101 (2025)

  14. [14]

    G. I. Taylor, The formation of a blast wave by a very intens e explosion I. Theoretical discussion, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 201, 159 (1950)

  15. [15]

    Sakurai, Blast Wave Theory, in Basic Developments in Fluid Dynamics , Vol

    A. Sakurai, Blast Wave Theory, in Basic Developments in Fluid Dynamics , Vol. I, edited by M. Holt (Academic Press, New York and London, 1965) Chap. 4, p. 101

  16. [16]

    Sedov, Similarity and Dimensional Methods in Mechanics (Academic Press, New York and London, 1959) p

    L. Sedov, Similarity and Dimensional Methods in Mechanics (Academic Press, New York and London, 1959) p. 363

  17. [17]

    J. v. Neumann, Collected Works, Volume VI: Theory of Games, Astrophysics, Hydrodynamics and Meteorology (Perga- mon Press, Oxford, 1963) p. 548

  18. [18]

    Y. B. Zel’dovich and Y. P. Raizer, Physics of Shock Waves and High-Temperature Hydrodynamic P henomena, edited by W. D. Hayes and R. F. Probstein (Dover Publications, Inc., Mi neola, New York, 2002) p. 916

  19. [19]

    Landau and E

    L. Landau and E. Lifshitz, Course of Theoretical Physics, Volume VI: Fluid Mechanics (Academic Press, New York and London, 1987) p. 539

  20. [20]

    G. I. Taylor, The formation of a blast wave by a very inten se explosion. - II. The atomic explosion of 1945, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 201, 175 (1950)

  21. [21]

    M. J. Edwards, A. J. MacKinnon, J. Zweiback, K. Shigemor i, D. Ryutov, A. M. Rubenchik, K. A. Keilty, E. Liang, B. A. Remington, and T. Ditmire, Investigation of u ltrafast laser-driven radiative blast waves, Phys. Rev. Lett. 87, 085004 (2001)

  22. [22]

    A. D. Edens, T. Ditmire, J. F. Hansen, M. J. Edwards, R. G. Adams, P. Rambo, L. Ruggles, I. C. Smith, and J. L. Porter, Study of high Mach number laser driven blast waves, Physics of Plasmas 11, 4968 (2004)

  23. [23]

    A. S. Moore, D. R. Symes, and R. A. Smith, Tailored blast w ave formation: Developing experiments pertinent to labora tory astrophysics, Physics of Plasmas 12, 052707 (2005)

  24. [24]

    J. P. Ostriker and C. F. McKee, Astrophysical blastwave s, Rev. Mod. Phys. 60, 1 (1988)

  25. [25]

    Tang and R

    X. Tang and R. A. Chevalier, Shock evolution in non-radi ative supernova remnants, Monthly Notices of the Royal Astronomical Society 465, 3793 (2017)

  26. [26]

    Y. Du, H. Li, and L. P. Kadanoff, Breakdown of hydrodynami cs in a one-dimensional system of inelastic particles, Phys. Rev. Lett. 74, 1268 (1995)

  27. [27]

    Jabeen, R

    Z. Jabeen, R. Rajesh, and P. Ray, Universal scaling dyna mics in a perturbed granular gas, Europhysics Letters 89, 34001 (2010)

  28. [28]

    Barbier, D

    M. Barbier, D. Villamaina, and E. Trizac, Blast dynamic s in a dissipative gas, Phys. Rev. Lett. 115, 214301 (2015)

  29. [29]

    J. Joy, S. Pathak, and R. Rajesh, Shock propagation foll owing an intense explosion: Comparison between hydrodynam ics and simulations, Jounal of Statistical Physics 182, 34 (2021)

  30. [30]

    Joy and R

    J. Joy and R. Rajesh, Shock propagation in the hard spher e gas in two dimensions: Comparison between simulations and hydrodynamics, Jounal of Statistical Physics 184, 3 (2021)

  31. [31]

    Kumar and R

    A. Kumar and R. Rajesh, Shock propagation following an i ntense explosion in an inhomogeneous gas: Core scaling and hydrodynamics, Journal of Statistical Physics 192, 17, 18 (2025)

  32. [32]

    Kumar, A

    U. Kumar, A. Dhar, and P. L. Krapivsky, Impurity dynamics in a zero-temperature gas (2025), arXiv:2505.02225 [cond-mat.stat-mech]

  33. [33]

    Y. G. Sinai, Dynamical systems with elastic reflections , Russian Mathematical Surveys 25, 137 (1970)

  34. [34]

    MacKay and J

    R. MacKay and J. Meiss, eds., Hamiltonian Dynamical Systems (Adam Hilger, Bristol and Philadelphia, 1987)

  35. [35]

    Tabachnikov, Geometry and Billiards , Student Mathematical Library, Vol

    S. Tabachnikov, Geometry and Billiards , Student Mathematical Library, Vol. 30 (Americal Mathemat ical Society, 2005) p. 178

  36. [36]

    Bunimovich, D

    L. Bunimovich, D. Burago, N. Chernov, E. Cohen, C. Dettm ann, S. Dorfman, J.R. Ferleger, R. Hirschl, A. Kononenko, J. Lebowitz, C. Liverani, T. Murphy, J. Piasecki, H. Posch, N . Sim´ anyi, Y. Sinai, D. Sz´ asz, T. T´ el, H. Beijeren, R. Zon, J. Vollmer, and L. Young, Hard Ball Systems and the Lorentz Gas , edited by D. Sz´ asz, Encyclopaedia of Mathemati...

  37. [37]

    Grassberger, W

    P. Grassberger, W. Nadler, and L. Yang, Heat conduction and entropy production in a one-dimensional hard-particle gas, Phys. Rev. Lett. 89, 180601 (2002)

  38. [38]

    P. I. Hurtado and P. L. Garrido, A violation of universal ity in anomalous Fourier’s law, Scientific Reports 6, 38823, 10 (2016)

  39. [39]

    Lepri, R

    S. Lepri, R. Livi, and A. Politi, Too close to integrable : Crossover from normal to anomalous heat diffusion, Phys. Rev. Lett. 125, 040604 (2020) . 20

  40. [40]

    Grassberger, D

    P. Grassberger, D. Dhar, and P. K. Mohanty, Many univers ality classes in an interface model restricted to non-negat ive heights, Phys. Rev. E 107, 044112 (2023)

  41. [41]

    P. I. Hurtado, Breakdown of hydrodynamics in a simple on e-dimensional fluid, Phys. Rev. Lett. 96, 010601 (2006)

  42. [42]

    Ganapa, S

    S. Ganapa, S. Chakraborti, P. L. Krapivsky, and A. Dhar, Blast in the one-dimensional cold gas: Comparison of micros copic simulations with hydrodynamic predictions, Physics of Fluids 33, 087113 (2021)

  43. [43]

    Chakraborti, S

    S. Chakraborti, S. Ganapa, P. L. Krapivsky, and A. Dhar, Blast in a one-dimensional cold gas: From Newtonian dynamic s to hydrodynamics, Phys. Rev. Lett. 126, 244503 (2021)

  44. [44]

    Dhar, Heat conduction in a one-dimensional gas of ela stically colliding particles of unequal masses, Phys

    A. Dhar, Heat conduction in a one-dimensional gas of ela stically colliding particles of unequal masses, Phys. Rev. Lett. 86, 3554 (2001)

  45. [45]

    P. L. Garrido, P. I. Hurtado, and B. Nadrowski, Simple on e-dimensional model of heat conduction which obeys Fourier ’s law, Phys. Rev. Lett. 86, 5486 (2001)

  46. [46]

    Casati and T

    G. Casati and T. Prosen, Anomalous heat conduction in a o ne-dimensional ideal gas, Phys. Rev. E 67, 015203 (2003)

  47. [47]

    Cipriani, S

    P. Cipriani, S. Denisov, and A. Politi, From anomalous e nergy diffusion to L´ evy walks and heat conductivity in one- dimensional systems, Phys. Rev. Lett. 94, 244301 (2005)

  48. [48]

    Williams, Algorithm 232: Heapsort, Communications of the ACM 7, 347 (1964)

    J. Williams, Algorithm 232: Heapsort, Communications of the ACM 7, 347 (1964)

  49. [49]

    Holovatch, Y

    T. Holovatch, Y. Kozitsky, K. Pilorz, and Y. Holovatch, Effective and asymptotic scaling in a one-dimensional billi ard problem, Condens. Matter Phys. 28, 23401 (2025)

  50. [50]

    Baxter, Exactly Solved Models in Statistical Mechanics , Dover books on physics (Dover Publications, 2007)

    R. Baxter, Exactly Solved Models in Statistical Mechanics , Dover books on physics (Dover Publications, 2007)

  51. [51]

    Sheng and T

    W. Sheng and T. Zhang, The Riemann problem for the transp ortation equations in gas dynamics, Memoirs of the American Mathematical Society 137, 1 (1999)

  52. [52]

    Swamy, Generalized Fibonacci and Lucas polynomials , and their associated diagonal polynomials, The Fibonacci Quarterly 37, 213 (1999)

    M. Swamy, Generalized Fibonacci and Lucas polynomials , and their associated diagonal polynomials, The Fibonacci Quarterly 37, 213 (1999) . 21