Pith. sign in

REVIEW 4 major objections 5 minor 61 references

A harmonic crystal's out-of-plane ripples thermalize in two power-law stages, while frequencies and defects proliferate together.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-13 12:02 UTC

load-bearing objection Useful MD phenomenology on a clamped harmonic lattice, but the supplied text is truncated so the power-law and concurrency claims cannot be checked; still worth a referee if the full results exist. the 4 major comments →

arxiv 2604.03913 v1 submitted 2026-04-05 cond-mat.stat-mech cond-mat.mtrl-scicond-mat.softphysics.class-phphysics.comp-ph

A molecular dynamics simulation of thermalization of crystalline lattice with harmonic interaction

Zhenwei Yao This is my paper
classification cond-mat.stat-mech cond-mat.mtrl-scicond-mat.softphysics.class-phphysics.comp-ph
keywords thermalizationharmonic latticemolecular dynamicstopological defectsout-of-plane fluctuationspower-law proliferationgeometric nonlinearitybroken up-down symmetry
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper asks how a many-body system forgets its initial state and reaches thermal equilibrium, using the simplest possible crystal: a circular patch of particles linked by harmonic springs with a fixed rim, started with random velocities. Because the lattice can buckle out of plane, the geometry itself supplies the nonlinearity needed for thermalization. Molecular-dynamics trajectories show that the transverse and longitudinal parts of the velocity relax at different rates, that the number of dominant frequencies grows as a power law, and that this growth occurs together with a rapid rise in topological defects. The persistent out-of-plane deformations themselves follow two successive power laws with fractional exponents once the up-down symmetry is broken. The work therefore supplies a concrete, atomistic picture of how geometric nonlinearity drives the dynamical adaptation of a many-body system to a sudden disturbance.

Core claim

In a clamped harmonic triangular lattice evolved under Hamiltonian dynamics from random initial velocities, thermalization proceeds by distinct relaxation of transverse versus longitudinal velocity components, by a power-law proliferation of dominant frequencies that coincides with the rapid appearance of topological defects, and by two-stage out-of-plane fluctuations whose fractional power laws are tied to the spontaneous breaking of up-down symmetry.

What carries the argument

The drum-like harmonic triangular lattice with anchored circular boundary: purely harmonic springs whose geometric nonlinearity (out-of-plane buckling) breaks integrability and thereby permits thermalization under Verlet Hamiltonian dynamics.

Load-bearing premise

The claim rests on the premise that geometric nonlinearity alone, in an otherwise purely harmonic lattice of finite size with fixed rim, is enough to produce genuine thermalization rather than quasi-periodic motion or finite-size artifacts.

What would settle it

A controlled comparison, at several system sizes, of the same initial-velocity ensemble under purely planar (strictly two-dimensional) constraints versus free three-dimensional motion: if the two-stage fractional power laws, concurrent defect proliferation, and loss of initial-state memory disappear when out-of-plane motion is forbidden, the geometric-nonlinearity mechanism is falsified.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript studies thermalization in a classical many-body system via molecular-dynamics simulation of a circular, boundary-anchored triangular lattice with purely harmonic nearest-neighbor interactions, evolved under Verlet Hamiltonian dynamics from random initial velocities. From atomic-level trajectories the author reports: (i) distinct relaxation rates of transverse versus longitudinal velocity components; (ii) a power-law proliferation of dominant frequencies; (iii) concurrent rapid growth of those frequencies and topological defects (disclinations); and (iv) two-stage, fractional-exponent power-law fluctuations of persistent out-of-plane deformations linked to broken up-down symmetry. The work is framed as bridging nonlinear dynamics and statistical mechanics, and as a model for thermal fluctuations and mechanical instabilities in 2D crystalline systems.

Significance. Thermalization of nearly integrable or weakly nonlinear lattices remains a central open problem (FPU and related lines). A geometrically nonlinear but energetically harmonic lattice is a clean setting in which to isolate the role of configuration-space nonlinearity. If the reported power laws, concurrency of frequency and defect proliferation, and two-stage out-of-plane scaling survive controlled checks, the paper would supply a concrete microscopic dynamical picture of how a crystalline sheet loses memory of its initial state and develops topological disorder under Hamiltonian evolution. The multi-channel analysis (velocity components, spectral content, defects, and out-of-plane height) is a genuine strength relative to purely spectral or purely configurational studies. The connection to 2D melting and thermally driven mechanical instabilities is well motivated.

major comments (4)
  1. The central modeling premise (Introduction) is that geometric nonlinearity of an otherwise harmonic triangular lattice with clamped circular boundary is sufficient to break integrability and produce genuine thermalization (loss of initial-state information and exploration of permissible states). The available text does not supply an operational thermalization criterion—e.g., approach to equipartition across modes, decay of velocity or mode autocorrelations, Lyapunov spectra, or comparison against a small-amplitude (near-integrable) control. Without such diagnostics, the reported power laws and concurrency could equally be finite-size quasi-periodic or transient phenomena. A clear definition and at least one quantitative diagnostic of thermalization are load-bearing for all main claims.
  2. The strongest quantitative claims—power-law proliferation of dominant frequencies, concurrency with topological-defect growth, and two-stage fractional-exponent out-of-plane fluctuations—are stated in the Abstract and Introduction but are not supported by methods parameters, fitting windows, error bars, system-size dependence, or number of independent runs in the supplied manuscript body (which jumps from the Introduction to Fig. 7 and the reference list). Robustness against spectral-extraction thresholds, Verlet time step, initial-velocity amplitude, and lattice size N must be shown; otherwise the fractional exponents and concurrency remain unverifiable.
  3. Fig. 7 and the surrounding discussion introduce four-, five-, seven-, and eight-fold disclinations as markers of topological disorder concurrent with frequency proliferation. For a triangular lattice the natural topological charges are 5 and 7; the appearance and dynamical role of 4- and 8-fold defects need a precise identification protocol (e.g., Voronoi/Delaunay criteria under large out-of-plane displacement) and a demonstration that they are not analysis artifacts of projecting a strongly buckled sheet. The claimed concurrency with frequency proliferation requires a quantitative time-series comparison (not a single snapshot).
  4. The two-stage out-of-plane fluctuation law is attributed to broken up-down symmetry. The manuscript should clarify whether the symmetry breaking is spontaneous (and how it is measured, e.g., by a global height moment) or is seeded by the initial velocity draw, and whether the two fractional exponents are stable under reversal of the initial out-of-plane bias. Without that, the association of the two stages with broken up-down symmetry remains interpretive rather than demonstrated.
minor comments (5)
  1. The manuscript body as provided is truncated after the Introduction (page content jumps to Fig. 7 and references). Ensure the full Methods, Results, and figure set are present and consistently numbered in the submission.
  2. Introduction: several self-citations to the author’s prior Lennard-Jones and three-body disturbance papers are appropriate as background, but a short explicit contrast with those anharmonic models would help the reader see what is new in the purely harmonic geometric-nonlinearity setting.
  3. Fig. 7 caption and color legend for defect types should state the coordination-number algorithm and whether out-of-plane neighbors are included in the coordination count.
  4. Notation for transverse/longitudinal velocity components and for the ‘dominant frequency’ threshold should be defined once, early, and used consistently.
  5. References include standard FPU, KAM, and 2D-melting sources; adding a brief pointer to modern numerical thermalization diagnostics in classical lattices would help place the power-law claims in context.

Circularity Check

0 steps flagged

No circularity found: pure observational MD study whose reported power laws and concurrency claims are outputs of Verlet trajectories, not reductions of fitted inputs or self-definitional constructions.

full rationale

The manuscript (even in its incomplete form) presents a classical molecular-dynamics investigation of a harmonic triangular lattice under clamped circular boundary conditions, initialized with random velocities and integrated by the Verlet algorithm. All central claims—distinct transverse/longitudinal velocity relaxation, power-law proliferation of dominant frequencies, concurrent growth of topological defects, and two-stage fractional-exponent out-of-plane fluctuations—are empirical observations extracted from the simulated trajectories. No equation, fitting procedure, or uniqueness theorem is used to define a quantity that is later re-presented as a prediction. The few self-citations ([28], [29]) appear only as background illustrations of nonlinear dynamics in related Lennard-Jones systems; they are not load-bearing for the harmonic-lattice results. Geometric nonlinearity is invoked as a modeling premise that breaks integrability, but that premise is an input assumption, not a circular derivation of the reported exponents. Consequently the derivation chain contains no self-definitional loops, no fitted-input-as-prediction steps, and no uniqueness imported from the author’s prior work. Score 0 is therefore required by the analyzer rules.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 0 invented entities

The central claims rest on classical Hamiltonian MD of a harmonic lattice whose only nonlinearity is geometric. Free parameters are the usual simulation knobs (initial velocity scale, system size, time step, analysis thresholds). Axioms are standard classical mechanics plus the modeling choice that geometric nonlinearity alone thermalizes the system. No new particles or forces are invented; topological defects are standard lattice defects identified in the trajectories.

free parameters (4)
  • initial velocity distribution scale / random kick amplitude
    Sets the energy injected into the dissipationless system; controls how far the lattice is driven from the flat ground state and therefore the observed defect and frequency proliferation rates.
  • lattice size N and circular radius (number of particles)
    Finite-size cutoff for long-wavelength out-of-plane modes and defect statistics; power-law regimes can depend on N.
  • Verlet time step and total integration time
    Numerical accuracy and the temporal window over which power laws and two-stage fluctuation regimes are fitted.
  • thresholds for identifying 'dominant frequencies' and topological defects
    Analysis cutoffs that define the observables whose proliferation and concurrency are claimed; not fixed by the Hamiltonian alone.
axioms (4)
  • domain assumption Particles interact via a purely harmonic pair potential; nonlinearity arises only from geometry (large displacements / out-of-plane motion).
    Stated in the Introduction as the model that 'breaks the integrability of the system' while remaining simple; all thermalization claims rest on this.
  • domain assumption The system is a dissipationless Hamiltonian system evolved by Verlet integration from random initial velocities with anchored circular boundary.
    Introduction; defines the dynamical ensemble in which relaxation rates and power laws are measured.
  • standard math Classical Newtonian mechanics and standard continuum/lattice notions of topological defects (disclinations of coordination 4,5,7,8) apply.
    Background used to interpret Fig. 7 and defect proliferation claims.
  • ad hoc to paper Geometric nonlinearity is sufficient for the system to lose initial-state information and explore permissible states (thermalization).
    Load-bearing modeling premise in the Introduction; not independently proved for this drum geometry in the accessible text.

pith-pipeline@v1.1.0-grok45 · 10204 in / 3070 out tokens · 36319 ms · 2026-07-13T12:02:44.986728+00:00 · methodology

0 comments
read the original abstract

Understanding the realization of thermal equilibrium through the thermalization process in a many-body system is a fundamental and complex scientific question, bridging thermodynamics and classical dynamics and connecting to a host of physical phenomena, such as mechanical instabilities in a thermal environment. In this work, based on the harmonic lattice model, we investigate the thermalization process in both velocity and coordinate spaces, by examining microscopic dynamics on the atomic level. We show the distinct relaxation rates of the transverse and longitudinal components of the velocity, reveal the power law governing the nonlinear proliferation of dominant frequencies, and observe the concurrent rapid proliferations of frequencies and topological defects. We also show that the lattice system's persistent out-of-plane deformations exhibit two-stage fluctuation behaviors, characterized by distinct power laws of fractional exponents and associated with the broken up-down symmetry. This work demonstrates the rich dynamics underlying the thermalization process, and advances our understanding on the dynamical adaptations of many-body systems to external disturbances.

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

61 extracted references

  1. [1]

    J. C. Maxwell, V. illustrations of the dynamical theory of gases.—part i. on the motions and collisions of perfectly elastic spheres, Philos. Mag.19, 19 (1860)

  2. [2]

    Boltzmann,Lectures On Gas Theory(University of California Press, Berkeley, 1964)

    L. Boltzmann,Lectures On Gas Theory(University of California Press, Berkeley, 1964)

    1964

  3. [3]

    Ehrenfest and T

    P. Ehrenfest and T. Ehrenfest,The Conceptual Founda- tions of The Statistical Approach in Mechanics(Courier Corporation, Massachusetts, 2002)

    2002

  4. [4]

    I. G. Sinai,Dynamical Systems II: Ergodic Theory with Applications to Dynamical Systems and Statistical Me- chanics(Springer, Berlin, 1989)

    1989

  5. [5]

    Prigogine,Non-equilibrium Statistical Mechanics (Dover Publications, New York, 2017)

    I. Prigogine,Non-equilibrium Statistical Mechanics (Dover Publications, New York, 2017)

    2017

  6. [6]

    H. S. Dumas,The KAM Story(World Scientific Publish- ing Company, Singapore, 2014)

    2014

  7. [7]

    H. B. Callen and T. A. Welton, Irreversibility and gen- eralized noise, Phys. Rev.83, 34 (1951)

    1951

  8. [8]

    Zheng, G

    Z. Zheng, G. Hu, and J. Zhang, Ergodicity in hard-ball systems and boltzmann’s entropy, Phys. Rev. E53, 3246 (1996)

    1996

  9. [9]

    Frenkel, Order through entropy, Nat

    D. Frenkel, Order through entropy, Nat. Mater.14, 9 (2015)

    2015

  10. [10]

    Z. Wang, W. Fu, Y. Zhang, and H. Zhao, Thermaliza- tion of two-and three-dimensional classical lattices, Phys. Rev. Lett.132, 217102 (2024)

    2024

  11. [11]

    Onsager, Reciprocal relations in irreversible processes

    L. Onsager, Reciprocal relations in irreversible processes. i., Phys. Rev.37, 405 (1931)

    1931

  12. [12]

    Kubo, The fluctuation-dissipation theorem, Rep

    R. Kubo, The fluctuation-dissipation theorem, Rep. Prog. Phys.29, 255 (1966)

    1966

  13. [13]

    Li and J

    T.-Y. Li and J. A. Yorke, Period three implies chaos, Am. Math. Mon.82, 985 (1975)

    1975

  14. [14]

    Scheck,Mechanics: From Newton ’s Laws to Determin- istic Chaos(Springer Science & Business Media, Berlin, 2010)

    F. Scheck,Mechanics: From Newton ’s Laws to Determin- istic Chaos(Springer Science & Business Media, Berlin, 2010)

    2010

  15. [15]

    Kaplan and L

    D. Kaplan and L. Glass,Understanding Nonlinear Dy- namics(Springer Science & Business Media, Berlin, 2012)

    2012

  16. [16]

    Trachenko and V

    K. Trachenko and V. Brazhkin, Collective modes and thermodynamics of the liquid state, Rep. Prog. Phys.79, 016502 (2015)

    2015

  17. [17]

    Saporta Katz and E

    O. Saporta Katz and E. Efrati, Self-driven fractional rotational diffusion of the harmonic three-mass system, Phys. Rev. Lett.122, 024102 (2019)

    2019

  18. [18]

    Saporta Katz and E

    O. Saporta Katz and E. Efrati, Regular regimes of the harmonic three-mass system, Phys. Rev. E101, 032211 (2020)

    2020

  19. [19]

    Yao, Collective dynamics and shattering of disturbed two-dimensional lennard–jones crystals, Eur

    Z. Yao, Collective dynamics and shattering of disturbed two-dimensional lennard–jones crystals, Eur. Phys. J. E 45, 88 (2022)

    2022

  20. [20]

    Yao, Non-linear dynamics and emergent statistical reg- ularity in classical lennard–jones three-body system upon disturbance, Eur

    Z. Yao, Non-linear dynamics and emergent statistical reg- ularity in classical lennard–jones three-body system upon disturbance, Eur. Phys. J. B96, 159 (2023)

    2023

  21. [21]

    L. P. Kadanoff, On two levels, Physics today39, 7 (1986)

    1986

  22. [22]

    L. P. Kadanoff,From Order to Chaos II(World Scien- tific, 1999)

    1999

  23. [23]

    Fermi, P

    E. Fermi, P. Pasta, S. Ulam, and M. Tsingou,Studies of the nonlinear problems, Tech. Rep. (Los Alamos Scientific Lab., 1955)

    1955

  24. [24]

    Dauxois, S

    T. Dauxois, S. Ruffo, E. Arimondo, and M. Wilkens, Dynamics and Thermodynamics of Systems with Long- Range Interactions(Springer, 2002)

    2002

  25. [25]

    Berman and F

    G. Berman and F. Izrailev, The fermi–pasta–ulam prob- lem: fifty years of progress, Chaos15(2005)

    2005

  26. [26]

    Gallavotti,The Fermi-Pasta-Ulam Problem: A Status Report(Springer Berlin, Heidelberg, 2007)

    G. Gallavotti,The Fermi-Pasta-Ulam Problem: A Status Report(Springer Berlin, Heidelberg, 2007)

    2007

  27. [27]

    Mulansky, K

    M. Mulansky, K. Ahnert, A. Pikovsky, and D. L. She- pelyansky, Dynamical thermalization of disordered non- linear lattices, Phys. Rev. E80, 056212 (2009)

    2009

  28. [28]

    A. C. Ribeiro-Teixeira, F. P. Benetti, R. Pakter, and Y. Levin, Ergodicity breaking and quasistationary states in systems with long-range interactions, Phys. Rev. E89, 022130 (2014)

    2014

  29. [29]

    Horn and H

    T. Horn and H. L¨ owen, How does a thermal binary crys- tal break under shear?, J. Chem. Phys.141(2014)

    2014

  30. [30]

    Koˇ smrlj and D

    A. Koˇ smrlj and D. R. Nelson, Response of thermalized ribbons to pulling and bending, Phys. Rev. B93, 125431 (2016). 11

    2016

  31. [31]

    Zaccone,Theory of Disordered Solids(Springer, New York, 2023)

    A. Zaccone,Theory of Disordered Solids(Springer, New York, 2023)

    2023

  32. [32]

    Jonay and T

    C. Jonay and T. Zhou, Physical theory of two-stage ther- malization, Phys. Rev. B110, L020306 (2024)

    2024

  33. [33]

    Vargas, T

    A. Vargas, T. Puccinelli, and J. R. Bordin, Order- disorder transitions and thermal pathways in frustrated 2d colloidal crystals, Braz. J. Phys.55, 281 (2025)

    2025

  34. [34]

    A. S. De Wijn and A. Fasolino, Relating chaos to deter- ministic diffusion of a molecule adsorbed on a surface, J. Phys. Condens. Matter21, 264002 (2009)

    2009

  35. [35]

    P. Keim, G. Maret, U. Herz, and H.-H. von Gr¨ unberg, Harmonic lattice behavior of two-dimensional colloidal crystals, Phys. Rev. Lett.92, 215504 (2004)

    2004

  36. [36]

    L. D. Landau and E. M. Lifshitz,Theory of Elasticity, 3rd edition(Butterworth-Heinemann, Oxford, UK, 1986)

    1986

  37. [37]

    Audoly and Y

    B. Audoly and Y. Pomeau,Elasticity and Geometry(Ox- ford University Press, Oxford, UK, 2010)

    2010

  38. [38]

    Halperin and D

    B. Halperin and D. R. Nelson, Theory of two-dimensional melting, Phys. Rev. Lett.41, 121 (1978)

    1978

  39. [39]

    K. J. Strandburg, Two-dimensional melting, Rev. Mod. Phys.60, 161 (1988)

    1988

  40. [40]

    D. R. Nelson,Defects and Geometry in Condensed Matter Physics(Cambridge University Press, Cambridge, 2002)

    2002

  41. [41]

    M. K. Blees, A. W. Barnard, P. A. Rose, S. P. Roberts, K. L. McGill, P. Y. Huang, A. R. Ruyack, J. W. Kevek, B. Kobrin, D. A. Muller, and P. L. McEuen, Graphene kirigami, Nature524, 204 (2015)

    2015

  42. [42]

    D. Wan, D. R. Nelson, and M. J. Bowick, Thermal stiffen- ing of clamped elastic ribbons, Phys. Rev. B96, 014106 (2017)

    2017

  43. [43]

    Le Doussal and L

    P. Le Doussal and L. Radzihovsky, Thermal buckling transition of crystalline membranes in a field, Phys. Rev. Lett.127, 015702 (2021)

    2021

  44. [44]

    P. Z. Hanakata, S. S. Bhabesh, M. J. Bowick, D. R. Nelson, and D. Yllanes, Thermal buckling and symme- try breaking in thin ribbons under compression, Extreme Mech. Lett.44, 101270 (2021)

    2021

  45. [45]

    Z. Chen, D. Wan, and M. J. Bowick, Spontaneous tilt of single-clamped thermal elastic sheets, Phys. Rev. Lett. 128, 028006 (2022)

    2022

  46. [46]

    Rapaport,The Art of Molecular Dynamics Simulation (Cambridge University Press, Cambridge, UK, 2004)

    D. Rapaport,The Art of Molecular Dynamics Simulation (Cambridge University Press, Cambridge, UK, 2004)

    2004

  47. [47]

    Acharya, S

    P. Acharya, S. Sengupta, B. Chakraborty, and K. Ramola, Athermal fluctuations in disordered crystals, Phys. Rev. Lett.124, 168004 (2020)

    2020

  48. [48]

    Maharana, Athermal fluctuations in three dimensional disordered crystals, J

    R. Maharana, Athermal fluctuations in three dimensional disordered crystals, J. Stat. Mech. Theory Exp.2022, 103201 (2022)

    2022

  49. [49]

    N. Xu, M. Wyart, A. J. Liu, and S. R. Nagel, Excess vibrational modes and the boson peak in model glasses, Phys. Rev. Lett.98, 175502 (2007)

    2007

  50. [50]

    M. L. Manning and A. J. Liu, Vibrational modes identify soft spots in a sheared disordered packing, Phys. Rev. Lett.107, 108302 (2011)

    2011

  51. [51]

    Z. W. Wu, Y. Chen, W.-H. Wang, W. Kob, and L. Xu, Topology of vibrational modes predicts plastic events in glasses, Nat. Commun.14, 2955 (2023)

    2023

  52. [52]

    K. W. Wojciechowski and J. Klos, On the minimum energy structure of soft, two-dimensional matter in a strong uniform field:gravity’s rainbow’revisited, J. Phys. A: Math. Gen.29, 3963 (1996)

    1996

  53. [53]

    Mughal and M

    A. Mughal and M. Moore, Topological defects in the crys- talline state of one-component plasmas of nonuniform density, Phys. Rev. E76, 011606 (2007)

    2007

  54. [54]

    Mughal and D

    A. Mughal and D. Weaire, Curvature in conformal map- pings of two-dimensional lattices and foam structure, Proc. R. Soc. London, Ser. A465, 219 (2009)

    2009

  55. [55]

    Yao and M

    Z. Yao and M. Olvera de la Cruz, Topological defects in flat geometry: the role of density inhomogeneity, Phys. Rev. Lett.111, 115503 (2013)

    2013

  56. [56]

    V. Soni, L. R. G´ omez, and W. T. Irvine, Emergent ge- ometry of inhomogeneous planar crystals, Phys. Rev. X 8, 011039 (2018)

    2018

  57. [57]

    F. C. Silva, R. M. Menezes, L. R. Cabral, and C. C. de Souza Silva, Formation and stability of conformal spi- rals in confined 2d crystals, J. Phys.: Condens. Matter 32, 505401 (2020)

    2020

  58. [58]

    Yao, Intrinsic conformal order revealed in geometri- cally confined long-range repulsive particles, Europhys

    Z. Yao, Intrinsic conformal order revealed in geometri- cally confined long-range repulsive particles, Europhys. Lett.146, 46003 (2024)

    2024

  59. [59]

    P. C. Hohenberg, Existence of long-range order in one and two dimensions, Phys. Rev.158, 383 (1967)

    1967

  60. [60]

    N. D. Mermin and H. Wagner, Absence of ferromag- netism or antiferromagnetism in one- or two-dimensional isotropic heisenberg models, Phys. Rev. Lett.17, 1133 (1966)

    1966

  61. [61]

    Coleman, There are no goldstone bosons in two dimen- sions, Commun

    S. Coleman, There are no goldstone bosons in two dimen- sions, Commun. Math. Phys.31, 259 (1973)

    1973