Physical Review Research
- Welcome
- Recent
- Subjects
- Accepted
- Collections
- Authors
- Referees
- Search
- About
- Scope
- Editorial Team
- Open Access
Origin of symmetry breaking in the grasshopper model
David Llamas, Jaron Kent-Dobias, Kun Chen, Adrian Kent, and Olga Goulko
Phys. Rev. Research 6, 023235 – Published 3 June 2024
- Article
- References
- No Citing Articles
- Supplemental Material
PDFHTMLExport Citation
Abstract
The planar grasshopper problem, originally introduced by Goulko and Kent [Proc. R. Soc. A 473, 20170494 (2017)], is a striking example of a model with long-range isotropic interactions whose ground states break rotational symmetry. In this paper we analyze and explain the nature of this symmetry breaking with emphasis on the importance of dimensionality. Interestingly, rotational symmetry is recovered in three dimensions for small jumps, which correspond to the nonisotropic cogwheel regime of the two-dimensional problem. We discuss simplified models that reproduce the symmetry properties of the original system in dimensions. For the full grasshopper model in two dimensions we obtain quantitative predictions for optimal perturbations of the disk. Our analytical results are confirmed by numerical simulations.
5 More
- Received 20 November 2023
- Accepted 18 March 2024
DOI:https://doi.org/10.1103/PhysRevResearch.6.023235
![Origin of symmetry breaking in the grasshopper model (8) Origin of symmetry breaking in the grasshopper model (8)](https://i0.wp.com/cdn.journals.aps.org/files/icons/creativecommons.png)
Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
Published by the American Physical Society
Physics Subject Headings (PhySH)
- Research Areas
Classical statistical mechanicsMathematical physicsOptimization problemsPattern formationPatterns in complex systemsQuantum foundationsQuantum information theory
- Techniques
Ising modelPerturbative methodsSimulated annealing
Interdisciplinary PhysicsQuantum Information, Science & TechnologyStatistical Physics & Thermodynamics
Authors & Affiliations
David Llamas1, Jaron Kent-Dobias2,*, Kun Chen3, Adrian Kent4,5,†, and Olga Goulko1,‡
- 1Department of Physics, University of Massachusetts Boston, Boston, Massachusetts 02125, USA
- 2Istituto Nazionale di Fisica Nucleare, Sezione di Roma I, 00185 Rome, Italy
- 3Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, USA
- 4Centre for Mathematical Sciences, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom
- 5Perimeter Institute for Theoretical Physics, 31 Caroline Street North, Waterloo, Ontario N2L 2Y5, Canada
- *jaron.kent-dobias@roma1.infn.it
- †apak@cam.ac.uk
- ‡olga.goulko@umb.edu
Article Text
Click to Expand
Supplemental Material
Click to Expand
References
Click to Expand
Issue
Vol. 6, Iss. 2 — June - August 2024
Subject Areas
- Interdisciplinary Physics
- Quantum Physics
- Statistical Physics
![Origin of symmetry breaking in the grasshopper model (12) Origin of symmetry breaking in the grasshopper model (12)](https://i0.wp.com/cdn.journals.aps.org/development/journals/images/author-services-placard.png)
Images
Figure 1
The stability of flat half-space in dimensions to plane-wave perturbations of wavenumber . All values are negative except at (which corresponds to translation of the interface and is not a probability-conserving perturbation) and for that are multiples of in the case . These values are zero, meaning that there is a marginal stability to perturbations with wavenumber commensurate with in two dimensions but not an instability. For larger dimension of space, the result becomes increasingly insensitive to the commensurability of and .
Figure 2
The coefficient of stability for the disk to small perturbations of -fold symmetry. A large dot is drawn at the smallest value of where each curve first becomes positive. The black bar shows the point in at which the transition to disconnected shapes occurs as measured in [1].
Figure 3
The scaled coefficient for increasing on the disk as a function of , along with the result for the semi-infinite plane () from (14). As is increased, the disk result asymptotically approaches that of the semi-infinite plane. Since the finite- curves tend to the zeros of the limit curve from above, the marginally stable points of the half-plane are destabilized at finite .
Figure 4
First and second most unstable modes at a given for . The blue and yellow lines show the location of the first and second peak in , respectively. The location of the first peak is extremely similar to (3.4) of [1] but not identical. If the cutoff in is increased, more bands appear at higher . Black markers denote the corresponding numbers of cogs for the optimal solutions found through numerical simulation in Ref.[1] (the same data are shown in the left panel of Fig.5 in Ref.[1]). The numerical results from [1] are in very close agreement with the current prediction.
Figure 5
(Left) Grasshopper probability for lawns with fixed cog number for two values of , obtained numerically for discrete lawns with cells. The maxima corresponding to the two leading unstable modes are clearly visible ( and for ; and for ) and are marked with thin vertical lines. These are indeed (local) maxima of the grasshopper problem. Horizontal-dashed lines denote the corresponding disk probabilities given by Eq.(3). (Right) The same probability rescaled by the squared amplitude of the cogs (symbols with error bars). Solid lines show from (25), which corresponds roughly with the finite- numerical data.
Figure 6
Study of discretization effects. (Left) The exact continuous probability functional for the solid 3-ball of unit volume given by Eq.(4) (red solid line) compared with the corresponding discrete (black dots) as function of the grasshopper jump distance . For the 3-ball configuration is the optimal lawn shape. (Right) Relative deviation of for the solid 3-ball configuration from the corresponding as function of the lattice spacing for two representative values of the grasshopper jump: (blue dots and line) and (red dots and line). For the highest resolutions considered () the discretization error is well below . Lines are to guide the eye.
Figure 7
Study of discretization effects. (Left) The exact continuous probability functional for the 3-shell where the inner radius is selected as (solid blue line) compared with the corresponding discrete (black dots) as function of the grasshopper jump distance . The corresponding probability for the solid 3-ball (red line) is also shown for comparison. For the 3-shell has a higher success probability than the 3-ball. (Right) Relative deviation of for the 3-shell configuration from the corresponding as function of the lattice spacing for two representative values of the grasshopper jump: (blue dots and line) and (red dots and line). The inner radius for each is , as before. For the highest resolutions considered () the discretization error is below . Lines are to guide the eye.
Figure 8
(Left) Cross sectionof the optimal configuration for found numerically for a system with spins. The configuration has the shape of a solid 3-ball. This configuration was found to be optimal for all . (Right) Cross sectionof the optimal configuration for found numerically for a system with spins. If the jump length exceeds the configurations remain isotropic (for ) but develop a spherical hole in the center; the radius of the hole grows with increasing . Note that the outer radius of the configuration is slightly larger than to ensure that it has unit volume.
Figure 9
(Left) The optimal inner radius of the 3-shell vs grasshopper jump. The numerically found inner radii (black dots) match very well the analytical result (red solid line). The value (blue-dashed lines) is shown for comparison. (Right) The corresponding optimal grasshopper probabilities for the optimal inner radius (red-solid line for analytical value and black dots for the numerical value) and for (blue-dashed lines). The isotropic 3-shell ceases to be optimal for jumps exceeding a critical value of approximately .
Figure 10
Histograms of the radial coordinates of the configuration boundary points for, from left to right, (isotropic 3-shell), (eightfold perturbation), and (sixfold perturbation). Results shown were obtained for systems with spins. The vertical-red lines denote the theoretical values for the optimal (for the respective value of ) inner and outer radii of the corresponding 3-shell configurations. The shaded regions mark the interval around the optimal radii.
Figure 11
Maximal values found numerically for the discrete grasshopper lawn probability . Results shown were obtained for systems with spins. Vertical lines denote boundaries between the different regimes. From left to right these are: isotropic solid ball, isotropic shells, eightfold shell perturbations, sixfold shell perturbations, ring with caps, nested crescents. Insets show examples of representative configurations. 3d animations displaying in full the features summarized here are given within the Supplemental Material [14].
Figure 12
The points give the difference in probability between the perturbed disk and the disk as a function of perturbation size for and . They were computed using numeric integration on the full expression. The solid line gives . The results agree at small , as expected.