pith. sign in

arxiv: 2605.26668 · v1 · pith:XQJE56OBnew · submitted 2026-05-26 · ❄️ cond-mat.stat-mech · nlin.PS

Spatiotemporal Structures of Parametrically Driven Nonlinear Lattices

Yu Funami , Kazushi Aoyama This is my paper

Pith reviewed 2026-07-01 16:12 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech nlin.PS
keywords Frenkel-Kontorova modelparametric drivingspatiotemporal orderingsubharmonic resonancewavenumber selectionnonlinear latticelinear stability analysis
0
0 comments X

The pith

Parametric driving of the Frenkel-Kontorova lattice produces subharmonic spatiotemporal order whose wavenumber rises continuously then jumps discontinuously to π.

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

The paper studies parametric driving applied to the one-dimensional Frenkel-Kontorova model. Numerical results show that the drive creates ordered patterns that oscillate at half the drive frequency, with a selected wavenumber k that depends on drive strength P_ac and frequency ω_ex. As P_ac increases, k grows gradually before jumping to π. An extended linear stability analysis identifies the continuous-k states as unstable modes and the π state as stable, with both selections arising from the combined action of parametric amplification and collective fluctuations generated by the lattice nonlinearity.

Core claim

In the one-dimensional Frenkel-Kontorova model, parametric vibration induces spatiotemporal ordering characterized by the subharmonic frequency ω_ex/2 and wavenumber k which takes nontrivial values depending on the vibration strength P_ac and frequency ω_ex. With increasing P_ac, k gradually increases and discontinuously jumps up to k=π. Based on an extended linear stability analysis, the former k resonance corresponds to the unstable mode and the latter π resonance to the stable mode, made possible by the interplay between the parametric amplification and the collective fluctuation developing through the nonlinear effect.

What carries the argument

Extended linear stability analysis that distinguishes unstable and stable modes selected by the interplay of parametric amplification with nonlinear collective fluctuations in the driven lattice.

If this is right

  • The wavenumber k increases continuously with drive strength before a discontinuous jump to π occurs.
  • The lower-k resonance corresponds to an unstable mode while the π resonance corresponds to a stable mode.
  • Wavenumber selection requires both parametric amplification and the collective fluctuations produced by nonlinearity.
  • The patterns remain locked at the subharmonic frequency ω_ex/2 throughout the range of drive strengths examined.

Where Pith is reading between the lines

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

  • The same mechanism may produce analogous wavenumber jumps in other one-dimensional nonlinear lattices under parametric drive.
  • The discontinuous jump to π suggests a possible first-order-like transition in pattern selection that could be tested by varying system size or temperature.
  • Similar ordering might appear in experimental realizations such as driven mechanical chains or Josephson-junction arrays.

Load-bearing premise

The extended linear stability analysis is sufficient to identify which observed patterns are unstable versus stable and to establish that the wavenumber selection arises specifically from the interplay of parametric amplification with nonlinear collective fluctuations rather than from other mechanisms.

What would settle it

A simulation or experiment in which the observed wavenumber values, the location of the jump to π, or the stability of the selected modes fails to match the predictions of the extended linear stability analysis.

Figures

Figures reproduced from arXiv: 2605.26668 by Kazushi Aoyama, Yu Funami.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Schematic of the [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. ST structures (left panels) and the corresponding [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: Since in the right panels in Fig. 2, the background [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
read the original abstract

We theoretically investigate the effects of parametric driving on the one-dimensional Frenkel-Kontorova model, a nonlinear many-body lattice system. It is numerically found that a parametric vibration induces spatiotemporal ordering characterized by the subharmonic frequency $\omega_{\mathrm{ex}}/2$ and wavenumber $k$ which takes nontrivial values depending on the vibration strength $P_{\mathrm{ac}}$ and frequency $\omega_{\mathrm{ex}}$. With increasing $P_{\mathrm{ac}}$, $k$ gradually increases and discontinuously jumps up to $k=\pi$. Based on an extended linear stability analysis, we show that the former $k$ resonance (the latter $\pi$ resonance) corresponds to the unstable (stable) mode and that the $k$ resonance is made possible by the interplay between the parametric amplification and the collective fluctuation developing through the nonlinear effect.

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

1 major / 2 minor

Summary. The manuscript investigates parametric driving in the one-dimensional Frenkel-Kontorova model. Numerical simulations show that parametric vibration induces spatiotemporal ordering with subharmonic frequency ω_ex/2 and wavenumber k that depends on drive strength P_ac and frequency ω_ex, increasing gradually before jumping discontinuously to k=π. An extended linear stability analysis attributes the nontrivial k resonance to unstable modes enabled by the interplay of parametric amplification and nonlinear collective fluctuations, while the π resonance corresponds to stable modes.

Significance. If the mechanism attribution holds, the work would contribute to understanding pattern formation and mode selection in driven nonlinear lattices, relevant to statistical mechanics of many-body systems. The combination of numerics and stability analysis is a positive feature, though the extension of the analysis requires clarification.

major comments (1)
  1. [Abstract and extended linear stability analysis section] Abstract (and the section presenting the extended linear stability analysis): The central claim that the observed nontrivial k (and its jump to π) arises specifically from the interplay between parametric amplification and nonlinear collective fluctuations is load-bearing. The extended linear stability analysis remains linear (e.g., Floquet multipliers around a uniform or periodic background) and therefore cannot directly demonstrate that nonlinearity is required for wavenumber selection, as opposed to linear parametric resonance alone. An explicit comparison, such as results obtained after suppressing the nonlinear terms while retaining the parametric drive, is needed to support the mechanism.
minor comments (2)
  1. [Abstract and numerical methods] The abstract (and methods section) provides no information on system size, integration method, boundary conditions, or the precise manner in which the linear stability analysis was extended. These details are required to evaluate robustness against finite-size effects or numerical artifacts, particularly regarding the discontinuous jump in k.
  2. [Abstract] A brief statement of the explicit form of the parametric drive term added to the Frenkel-Kontorova Hamiltonian or equations of motion would improve clarity for readers unfamiliar with the specific implementation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the major comment below and are prepared to revise the manuscript to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [Abstract and extended linear stability analysis section] Abstract (and the section presenting the extended linear stability analysis): The central claim that the observed nontrivial k (and its jump to π) arises specifically from the interplay between parametric amplification and nonlinear collective fluctuations is load-bearing. The extended linear stability analysis remains linear (e.g., Floquet multipliers around a uniform or periodic background) and therefore cannot directly demonstrate that nonlinearity is required for wavenumber selection, as opposed to linear parametric resonance alone. An explicit comparison, such as results obtained after suppressing the nonlinear terms while retaining the parametric drive, is needed to support the mechanism.

    Authors: We agree that an explicit comparison with the purely linear case (nonlinear terms suppressed) is required to directly substantiate that the nontrivial wavenumber selection arises from the interplay with nonlinear collective fluctuations rather than linear parametric resonance alone. In the revised manuscript we will add both numerical simulations and stability analysis results for the linear case, showing that only the k=π resonance persists without nonlinearity. We will update the abstract and the extended linear stability analysis section to incorporate this comparison and clarify the mechanism. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper reports independent numerical simulations of the driven Frenkel-Kontorova lattice that discover the P_ac- and ω_ex-dependent wavenumber k and its discontinuous jump to π. An extended linear stability analysis is then applied to the same model to classify the observed k resonance as unstable and the π resonance as stable, with the attribution to parametric amplification plus nonlinear fluctuations offered as an interpretive conclusion rather than a definitional identity. No equation reduces to its own input by construction, no fitted parameter is relabeled as a prediction, and no load-bearing premise rests on a self-citation chain. The numerical observations and the subsequent linear analysis constitute separate steps whose outputs are not forced by the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are described in the abstract; the work relies on the standard Frenkel-Kontorova model and conventional numerical and stability methods.

pith-pipeline@v0.9.1-grok · 5669 in / 1227 out tokens · 51660 ms · 2026-07-01T16:12:29.077984+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

52 extracted references

  1. [1]

    Faraday, On a Peculiar Class of Acoustical Figures; and on Certain Forms Assumed by Groups of Particles upon Vibrating Elastic Surfaces, Phil

    M. Faraday, On a Peculiar Class of Acoustical Figures; and on Certain Forms Assumed by Groups of Particles upon Vibrating Elastic Surfaces, Phil. Trans. R. Soc. Lon- don 121, 299 (1831)

  2. [2]

    T. B. Benjamin and F. Ursell, The Stability of the Plane Free Surface of a Liquid in Vertical Periodic Motion, Proc. R. Soc. A 225, 505 (1954)

    1954

  3. [3]

    Porter, C

    J. Porter, C. Topaz, and M. Silber, Pattern Control via Multifrequency Parametric Forcing, Phys. Rev. Lett. 93, 034502 (2004)

    2004

  4. [4]

    Dinesh, J

    B. Dinesh, J. Livesay, I. B. Ignatius, and R. Narayanan, Pattern formation in Faraday instability—experimental validation of theoretical models, Phil. Trans. R. Soc. A. 381, 20220081 (2023)

    2023

  5. [5]

    Staliunas, C

    K. Staliunas, C. Hang, and V. V. Konotop, Parametric patterns in optical fiber ring nonlinear resonators, Phys. Rev. A 88, 023846 (2013)

    2013

  6. [6]

    J. J. García-Ripoll, V. M. Pérez-García, and P. Torres, Ex- tended Parametric Resonances in Nonlinear Schrödinger Systems, Phys. Rev. Lett. 83, 1715 (1999)

    1999

  7. [7]

    Kivshar, Dark optical solitons: physics and applica- tions, Physics Reports 298, 81 (1998)

    Y. Kivshar, Dark optical solitons: physics and applica- tions, Physics Reports 298, 81 (1998)

    1998

  8. [8]

    H. Fu, L. Feng, B. M. Anderson, L. W. Clark, J. Hu, J. W. Andrade, C. Chin, and K. Levin, Density Waves and Jet Emission Asymmetry in Bose Fireworks, Phys. Rev. Lett. 121, 243001 (2018)

    2018

  9. [9]

    J. H. V. Nguyen, M. C. Tsatsos, D. Luo, A. U. J. Lode, G. D. Telles, V. S. Bagnato, and R. G. Hulet, Parametric Excitation of a Bose-Einstein Condensate: From Faraday Waves to Granulation, Phys. Rev. X 9, 011052 (2019)

    2019

  10. [10]

    A. I. Nicolin, R. Carretero-González, and P. G. Kevrekidis, Faraday waves in Bose-Einstein condensates, Phys. Rev. A 76, 063609 (2007)

    2007

  11. [11]

    Engels, C

    P. Engels, C. Atherton, and M. A. Hoefer, Observation of Faraday Waves in a Bose-Einstein Condensate, Phys. Rev. Lett. 98, 095301 (2007)

    2007

  12. [12]

    Buks and M

    E. Buks and M. L. Roukes, Electrically tunable collective response in a coupled micromechanical array, J. Microelec- tromech. Syst. 11, 802 (2002)

    2002

  13. [13]

    Lifshitz and M

    R. Lifshitz and M. C. Cross, Response of parametri- cally driven nonlinear coupled oscillators with application to micromechanical and nanomechanical resonator arrays, Phys. Rev. B 67, 134302 (2003)

    2003

  14. [14]

    Lifshitz and M

    R. Lifshitz and M. C. Cross, Nonlinear dynamics of nanomechanical and micromechanical resonators, Rev. Nonlinear Dyn. Complexity 1, 1 (2008)

    2008

  15. [15]

    Kovacic, R

    I. Kovacic, R. Rand, and S. Mohamed Sah, Mathieu’s equation and its generalizations: Overview of stability charts and their features, Appl. Mech. Rev. 70, 020802 (2018)

    2018

  16. [16]

    Eichler and O

    A. Eichler and O. Zilberberg, Classical and Quantum Parametric Phenomena (Oxford University Press, New York, 2023)

    2023

  17. [17]

    M. C. Cross and P. C. Hohenberg, Pattern formation outside of equilibrium, Rev. Mod. Phys. 65, 851 (1993)

    1993

  18. [18]

    L. I. Reyes, L. M. Pérez, L. Pedraja-Rejas, P. Díaz, J. Mendoza, J. Bragard, M. G. Clerc, and D. Laroze, Char- acterization of Faraday patterns and spatiotemporal chaos in parametrically driven dissipative systems, Chaos, Soli- tons & Fractals 186, 115244 (2024)

    2024

  19. [19]

    Z. G. Nicolaou, D. J. Case, E. B. van der Wee, M. M. Driscoll, and A. E. Motter, Heterogeneity-stabilized homo- geneous states in driven media, Nat. Commun. 12, 4486 (2021)

    2021

  20. [20]

    Y. S. Kivshar and M. Peyrard, Modulational instabilities in discrete lattices, Phys. Rev. A 46, 3198 (1992)

    1992

  21. [21]

    Citro, E

    R. Citro, E. G. Dalla Torre, L. D’Alessio, A. Polkovnikov, M. Babadi, T. Oka, and E. Demler, Dynamical stability of a many-body Kapitza pendulum, Ann. Phys. 360, 694 (2015)

    2015

  22. [22]

    Huang, Gap solitons in damped and parametrically driven nonlinear diatomic lattices, Phys

    G. Huang, Gap solitons in damped and parametrically driven nonlinear diatomic lattices, Phys. Rev. E 49, 5893 (1994)

    1994

  23. [23]

    Denardo, W

    B. Denardo, W. Wright, S. Putterman, and A. Larraza, Observation of a kink soliton on the surface of a liquid, Phys. Rev. Lett. 64, 1518 (1990)

    1990

  24. [24]

    I. Bena, C. Van Den Broeck, R. Kawai, M. Copelli, and K. Lindenberg, Collective behavior of parametric oscilla- tors, Phys. Rev. E 65, 036611 (2002)

    2002

  25. [25]

    W. I. Newman, R. H. Rand, and A. L. Newman, Dynam- ics of a nonlinear parametrically excited partial differential equation, Chaos 9, 242 (1999)

    1999

  26. [26]

    Chen, Experimental observation of solitons in a 1D nonlinear lattice, Phys

    W. Chen, Experimental observation of solitons in a 1D nonlinear lattice, Phys. Rev. B 49, 15063 (1994)

    1994

  27. [27]

    Y. S. Kivshar, Class of localized structures in nonlinear lattices, Phys. Rev. B 46, 8652 (1992). 6

    1992

  28. [28]

    Huang, J

    G. Huang, J. Shen, and H. Quan, Noncutoff kinks in damped and parametrically driven nonlinear lattices, Phys. Rev. B 48, 16795 (1993)

    1993

  29. [29]

    Denardo and W

    B. Denardo and W. B. Wright, Structural properties of kinks and domain walls in nonlinear oscillatory lattices, Phys. Rev. E 52, 1094 (1995)

    1995

  30. [30]

    Denardo, B

    B. Denardo, B. Galvin, A. Greenfield, A. Larraza, S. Put- terman, and W. Wright, Observations of localized struc- tures in nonlinear lattices: Domain walls and kinks, Phys. Rev. Lett. 68, 1730 (1992)

    1992

  31. [31]

    W. Chen, W. Lin, and Y. Zhu, Onset instability of a parametrically excited pendulum array, Phys. Rev. E 75, 016606 (2007)

    2007

  32. [32]

    N. Y. Yao, C. Nayak, L. Balents, and M. P. Zaletel, Clas- sical discrete time crystals, Nat. Phys. 16, 438 (2020)

    2020

  33. [33]

    T. L. Heugel, M. Oscity, A. Eichler, O. Zilberberg, and R. Chitra, Classical Many-Body Time Crystals, Phys. Rev. Lett. 123, 124301 (2019)

    2019

  34. [34]

    Z. G. Nicolaou and A. E. Motter, Anharmonic classical time crystals: A coresonance pattern formation mecha- nism, Phys. Rev. Res. 3, 023106 (2021)

    2021

  35. [35]

    Yi-Thomas and J

    S. Yi-Thomas and J. D. Sau, Theory for Dissipative Time Crystals in Coupled Parametric Oscillators, Phys. Rev. Lett. 133, 266601 (2024)

    2024

  36. [36]

    D. V. Else, C. Monroe, C. Nayak, and N. Y. Yao, Discrete Time Crystals, Annu. Rev. Condens. Matter Phys. 11, 467 (2020)

    2020

  37. [37]

    M. P. Zaletel, M. Lukin, C. Monroe, C. Nayak, F. Wilczek, and N. Y. Yao, Colloquium : Quantum and clas- sical discrete time crystals, Rev. Mod. Phys. 95, 031001 (2023)

    2023

  38. [38]

    T. L. Heugel, A. Eichler, R. Chitra, and O. Zilberberg, The role of fluctuations in quantum and classical time crystals, SciPost Phys. Core 6, 053 (2023)

    2023

  39. [39]

    Machado, Q

    F. Machado, Q. Zhuang, N. Y. Yao, and M. P. Zaletel, Absolutely Stable Time Crystals at Finite Temperature, Phys. Rev. Lett. 131, 180402 (2023)

    2023

  40. [40]

    Sarkar, Anurag, J

    A. Sarkar, Anurag, J. A. Mondal, R. Singh, A. A. Makki, A. K. Rathi, R. J. T. Nicholl, S. Chakraborty, K. I. Bolotin, and S. Ghosh, Observation of tunable discrete- time-crystalline phases, Phys. Rev. Applied 25, L021001 (2026)

    2026

  41. [41]

    O. M. Braun and Y. S. Kivshar, Nonlinear dynamics of the Frenkel–Kontorova model, Phys. Rep. 306, 1 (1998)

    1998

  42. [42]

    Ameye, A

    O. Ameye, A. Eichler, and O. Zilberberg, Parametric in- stability landscape of coupled Kerr parametric oscillators, Phys. Rev. Res. 7, 033204 (2025)

    2025

  43. [43]

    Funami and K

    Y. Funami and K. Aoyama, Supplemental Material which includes Refs.[15, 27–32], where the dynamics of domain walls, the derivation of the dispersion relation Eq. (5), the parametric instability, and the method of multiple scales are presented

  44. [44]

    Creutz, L

    M. Creutz, L. Jacobs, and C. Rebbi, Monte Carlo study of Abelian lattice gauge theories, Phys. Rev. D 20, 1915 (1979)

    1915

  45. [45]

    Aoyama, M

    K. Aoyama, M. Gen, and H. Kawamura, Effects of spin- lattice coupling and a magnetic field in classical Heisen- berg antiferromagnets on the breathing pyrochlore lattice, Phys. Rev. B 104, 184411 (2021), Appendix A

    2021

  46. [46]

    B. C. Denardo, Observations of Nonpropagating Oscil- latory Solitons. PhD thesis, University of California, Los Angeles, (1990)

    1990

  47. [47]

    Funami and K

    Y. Funami and K. Aoyama, Fractal and subharmonic re- sponses driven by surface acoustic waves during charge density wave sliding, Phys. Rev. B 108, L100508 (2023)

    2023

  48. [48]

    Funami and K

    Y. Funami and K. Aoyama, Shapiro Steps and Surface Acoustic Waves in Charge Density Wave Dynamics, New Phys.: Sae Mulli 73, 1086 (2023)

    2023

  49. [49]

    M. V. Nikitin, V. Ya. Pokrovskii, D. A. Kai, and S. G. Zybtsev, On the Fundamental Difference between the Ef- fects of Electrical and Mechanical Vibrations on the Dy- namics of a Charge Density Wave, JETP Lett. 118, 861 (2023)

    2023

  50. [50]

    Fujiwara, T

    K. Fujiwara, T. Kawada, N. Nikaido, J. Park, N. Jiang, S. Takada, and Y. Niimi, Observation of Shapiro Steps in the Charge Density Wave State Induced by Strain on a Piezoelectric Substrate, Phys. Rev. Lett. 135, 256304 (2025)

    2025

  51. [51]

    S. A. Nikonov, S. G. Zybtsev, and V. Ya. Pokrovskii, RF wave mixing with sliding charge-density waves, Appl. Phys. Lett. 118, 253108 (2021)

    2021

  52. [52]

    Spatiotemporal Structures of Parametrically Driven Nonlinear Lattices

    Y. Funami and K. Aoyama, Effects of frequency mix- ing on Shapiro-step formations in sliding charge-density- waves, Appl. Phys. Lett. 125, 173102 (2024). 7 Supplemental Material for “Spatiotemporal Structures of Parametrically Driven Nonlinear Lattices” S1. DOMAIN W ALL AND ITS DYNAMICS In the simulations, we often encounter a domain wall (DW), a phase-mi...

    2024