pith. sign in

arxiv: 2603.28609 · v1 · submitted 2026-03-30 · ❄️ cond-mat.soft · cond-mat.stat-mech

Rounded hard squares confined in a circle

Zhongtian Yuan , Yao Li This is my paper

Pith reviewed 2026-05-14 01:29 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords rounded hard squarescircular confinementtopological defectsdisclinationsMonte Carlo simulationsentropic effectstetratic orderstructural transitions
0
0 comments X

The pith

Rounded hard squares confined in a circle form a six-domain structure with six +1/4 disclinations and a central -1/2 disclination as roundness increases.

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

The paper uses NPT Monte Carlo simulations to study how circular confinement shapes the ordering of hard squares whose corners are rounded to different degrees. At low roundness the particles organize into one cross-shaped domain that shows tetratic order and four +1/4 disclinations at the boundary. When roundness is increased the system switches to a partitioned arrangement in which the particles divide into six separate domains separated by six +1/4 disclinations plus one -1/2 disclination at the center. A sympathetic reader would care because the result shows that modest changes in particle shape can select entirely different topological defect patterns through purely entropic effects.

Core claim

As roundness of the hard squares is increased under circular confinement the system undergoes a structural transition to a new partition structure. Particles self-assemble into six domains separated by six +1/4 disclinations together with a central -1/2 disclination. This replaces the integrated cross-shaped domain that exhibits tetratic order and four +1/4 disclinations at lower roundness. The transition is driven by the interplay between the circular boundary and the changing anisotropy of the particles.

What carries the argument

the roundness parameter that continuously tunes the particle corners from sharp to curved, thereby controlling the switch between the four-defect cross-shaped state and the six-plus-one defect partitioned state

Load-bearing premise

The NPT Monte Carlo simulations reach true equilibrium configurations rather than becoming trapped in long-lived metastable states for the rounded hard-square systems.

What would settle it

A simulation started from a completely random initial condition or run for significantly longer times that fails to produce the six-domain arrangement with exactly six +1/4 disclinations and one central -1/2 disclination would indicate the reported structure is not the equilibrium state.

Figures

Figures reproduced from arXiv: 2603.28609 by Yao Li, Zhongtian Yuan.

Figure 1
Figure 1. Figure 1: Average of local bond angular orders (Ψ¯ 4 and Ψ¯ 6) as functions of roundness ζ. The state diagram is segmented into four distinct regions, with representative snapshots illus￾trating the corresponding structures at ζ = 0.0, 0.41, 0.6, and 0.8, respectively. The particles sketched above the diagram depict rounded-corner hard-squares with roundness values of ζ = 0.1, 0.4, 0.6, and 0.8. into a square lattic… view at source ↗
Figure 2
Figure 2. Figure 2: Topological defects under square and partition [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Column shifts in square structure. Snapshots for [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Detailed structures of the square and partition states in the confined RCHS system. Each panel contains two [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: Large-scale simulation results for the partition [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: Order parameters as functions of roundness [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
read the original abstract

Packing under confinement could generate rich ordered structures through entropic effects, which is a fundamental problem in condensed matter, biophysics and material science. The influence of confinement to the anisotropic hard particles--particularly regarding the emergence of topological defect structures--remains poorly understood. Recent studies have shown that granular rods confined within circular boundaries can cluster into square like super-particles, forming four disclinations. In this study, we employ Monte Carlo simulations in the NPT ensemble to investigate how circular confinement influences the ordered structures of rounded-corner hard-squares with varying roundness. At low roundness, the system forms an integrated cross-shaped domain with tetratic order and four +1/4 disclinations in the corners, along with some column shifts. As roundness increases, we found a new partition structure, where particles self-assemble into six domains separated by six +1/4 disclinations and a central -1/2 disclination. Our findings reveal that the interplay between confinement geometry and colloid shape can drive entropy governed structural transitions, offering new insights for the design of topological metamaterials.

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 reports NPT Monte Carlo simulations of rounded hard squares confined in a circle. At low roundness the particles form an integrated cross-shaped domain with tetratic order and four +1/4 disclinations; at higher roundness a new six-domain structure appears, separated by six +1/4 disclinations and containing a central -1/2 disclination.

Significance. If the reported structures are equilibrium states, the work shows that modest changes in particle roundness can drive a topological transition between distinct defect patterns under circular confinement. This supplies a concrete example of entropy-driven self-assembly in anisotropic hard particles and may guide the design of metamaterials that exploit controlled disclination networks.

major comments (1)
  1. [Methods and Results] The central claim of an equilibrium six-domain structure rests on NPT Monte Carlo trajectories whose equilibration is not demonstrated. No system sizes, run lengths, acceptance rates, order-parameter autocorrelation times, or results from independent initial conditions (random vs. lattice) are reported, leaving open the possibility that the six-domain arrangement is a long-lived metastable state rather than the global free-energy minimum.
minor comments (2)
  1. [Abstract] The abstract refers to 'some column shifts' at low roundness; define this term and show the relevant local order parameter or snapshot in the main text.
  2. [Results] Clarify how the disclination charges (+1/4, -1/2) are identified and counted from the particle orientations; state the precise criterion used.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive evaluation of the significance of our work and for the constructive comment on the need for more rigorous documentation of the simulation protocols. We address this point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: The central claim of an equilibrium six-domain structure rests on NPT Monte Carlo trajectories whose equilibration is not demonstrated. No system sizes, run lengths, acceptance rates, order-parameter autocorrelation times, or results from independent initial conditions (random vs. lattice) are reported, leaving open the possibility that the six-domain arrangement is a long-lived metastable state rather than the global free-energy minimum.

    Authors: We agree that the original manuscript lacked sufficient technical details to fully demonstrate equilibration. In the revised version we will add a dedicated Methods subsection specifying the system sizes (N = 256 particles for the primary results, with checks at N = 512), total run lengths (10^7 equilibration steps followed by 5×10^7 production steps), acceptance rates (maintained at 0.25–0.35 via adaptive maximum displacements), and order-parameter autocorrelation times (tetratic order parameter decorrelates within ~2×10^5 steps). We will also include new figures showing the time series of the global tetratic order parameter and defect positions from at least five independent runs started from both random and pre-ordered lattice configurations; these trajectories converge to the same six-domain structure with a central −1/2 disclination at higher roundness. Thermodynamic integration estimates of the free energy will be added to confirm that the six-domain partition is the global minimum relative to the cross-shaped tetratic state. revision: yes

Circularity Check

0 steps flagged

No circularity: direct simulation outputs with no reduction to fitted inputs or self-citations

full rationale

The paper reports structural observations obtained from NPT Monte Carlo simulations of rounded hard squares in circular confinement. No analytical derivation chain, self-referential equations, fitted parameters renamed as predictions, or load-bearing self-citations appear in the provided text. The six-domain structure with +1/4 and -1/2 disclinations is presented as a direct simulation result rather than a quantity derived from or equivalent to the simulation inputs by construction. This is the expected non-finding for a purely computational study whose claims rest on observed configurations.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard hard-particle model (purely repulsive, entropy-driven) and the assumption that NPT Monte Carlo adequately samples the equilibrium ensemble for these confined systems.

free parameters (1)
  • roundness parameter
    Degree of corner rounding is varied as the control parameter; specific numerical values are not given in the abstract.
axioms (2)
  • domain assumption Hard-particle interactions are purely repulsive and entropic with no attractive forces
    Standard modeling choice for hard squares; invoked implicitly throughout the simulation description.
  • domain assumption NPT Monte Carlo sampling reaches equilibrium configurations for the confined system
    Required for interpreting the observed structures as stable phases rather than simulation artifacts.

pith-pipeline@v0.9.0 · 5483 in / 1393 out tokens · 45296 ms · 2026-05-14T01:29:54.942648+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

40 extracted references · 40 canonical work pages · 1 internal anchor

  1. [1]

    Rounded hard squares confined in a circle

    shows that particles form a hexagonal rotator crys- tal (RX) in which particles self-assemble in a triangular lattice, but their orientations are randomly rotated at ∗ liyao@nankai.edu.cn the center of the mess. Recent studies [7, 17] have sum- marized the phase behavior and melting transitions of hard regular polygons by MC simulation. Apart from the reg...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  2. [2]

    In this region,¯Ψ4 de- creases significantly, while¯Ψ6 rises to approximately0.8, since both hexatic and rhombic domains exhibit higher ¯Ψ6

    This structure is characterized by the coexistence of colloids forming a square domain with colloids forming either hexatic or rhombic domains. In this region,¯Ψ4 de- creases significantly, while¯Ψ6 rises to approximately0.8, since both hexatic and rhombic domains exhibit higher ¯Ψ6. The final green region corresponds to the RHX phase. Forζ >0.7, particle...

    work page doi:10.5281/zenodo.17971145

  3. [3]

    J. M. Kosterlitz and D. J. Thouless, Ordering, metasta- bility and phase transitions in two-dimensional systems, J. Phys. C6, 1181 (1973)

    work page 1973

  4. [4]

    B. I. Halperin and D. R. Nelson, Theory of Two- Dimensional Melting, Phys. Rev. Lett.41, 121 (1978)

    work page 1978

  5. [5]

    A. P. Young, Melting and the vector Coulomb gas in two dimensions, Phys. Rev. B19, 1855 (1979)

    work page 1979

  6. [6]

    C. H. Mak, Large-scale simulations of the two- dimensional melting of hard disks, Phys. Rev. E73, 065104 (2006)

    work page 2006

  7. [7]

    Engel, J

    M. Engel, J. A. Anderson, S. C. Glotzer, M. Isobe, E. P. Bernard, and W. Krauth, Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods, Phys. Rev. E87, 042134 (2013)

    work page 2013

  8. [8]

    Y. Han, N. Y. Ha, A. M. Alsayed, and A. G. Yodh, Melt- ing of two-dimensional tunable-diameter colloidal crys- tals, Phys. Rev. E77, 041406 (2008)

    work page 2008

  9. [9]

    J. A. Anderson, J. Antonaglia, J. A. Millan, M. Engel, and S. C. Glotzer, Shape and Symmetry Determine Two- Dimensional Melting Transitions of Hard Regular Poly- gons, Phys. Rev. X7, 021001 (2017)

    work page 2017

  10. [10]

    E. P. Bernard and W. Krauth, Two-Step Melting in Two Dimensions: First-Order Liquid-Hexatic Transition, Phys. Rev. Lett.107, 155704 (2011)

    work page 2011

  11. [11]

    Wojciechowski and D

    K. Wojciechowski and D. Frenkel, Tetratic phase in the planar hard square system, Comput. Methods Sci. Tech- nol.10, 235 (2004)

    work page 2004

  12. [12]

    Löffler, L

    R. Löffler, L. Siedentop, and P. Keim, Tetratic phase in 2D crystals of squares, Soft Matter21, 2026 (2025)

    work page 2026

  13. [13]

    Walsh and N

    L. Walsh and N. Menon, Ordering and dynamics of vi- brated hard squares, J. Stat. Mech.2016, 083302 (2016)

    work page 2016

  14. [14]

    van Anders, D

    G. van Anders, D. Klotsa, N. K. Ahmed, M. Engel, and S.C.Glotzer,Understandingshapeentropythroughlocal dense packing, Proc. Natl. Acad. Sci.111, E4812 (2014)

    work page 2014

  15. [15]

    K. Zhao, R. Bruinsma, and T. G. Mason, Local chiral symmetry breaking in triatic liquid crystals, Nat. Com- mun.3, 801 (2012)

    work page 2012

  16. [16]

    Donev, J

    A. Donev, J. Burton, F. H. Stillinger, and S. Torquato, Tetratic order in the phase behavior of a hard-rectangle system, Phys. Rev. B73, 054109 (2006)

    work page 2006

  17. [17]

    A. P. Gantapara, W. Qi, and M. Dijkstra, A novel chi- ral phase of achiral hard triangles and an entropy-driven demixing of enantiomers, Soft Matter11, 8684 (2015)

    work page 2015

  18. [18]

    Schilling, S

    T. Schilling, S. Pronk, B. Mulder, and D. Frenkel, Monte Carlo study of hard pentagons, Phys. Rev. E71, 036138 (2005)

    work page 2005

  19. [19]

    W. Shen, J. Antonaglia, J. A. Anderson, M. En- gel, G. van Anders, and S. C. Glotzer, Symmetries in hard polygon systems determine plastic colloidal crys- tal mesophases in two dimensions, Soft Matter15, 2571 (2019)

    work page 2019

  20. [20]

    M. A. Bates and D. Frenkel, Phase behavior of two- dimensional hard rod fluids, J. Chem. Phys.112, 10034 (2000)

    work page 2000

  21. [21]

    Bolhuis and D

    P. Bolhuis and D. Frenkel, Tracing the phase boundaries of hard spherocylinders, J. Chem. Phys.106, 666 (1997)

    work page 1997

  22. [22]

    K. Zhao, R. Bruinsma, and T. G. Mason, Entropic crys- tal–crystal transitions of Brownian squares, Proc. Natl. Acad. Sci.108, 2684 (2011)

    work page 2011

  23. [23]

    Avendaño and F

    C. Avendaño and F. A. Escobedo, Phase behavior of rounded hard-squares, Soft Matter8, 4675 (2012)

    work page 2012

  24. [24]

    Z. Hou, K. Zhao, Y. Zong, and T. G. Mason, Phase be- havior of two-dimensional Brownian systems of corner- rounded hexagons, Phys. Rev. Mater.3, 015601 (2019)

    work page 2019

  25. [25]

    Z.-K. Lin, Q. Wang, Y. Liu, H. Xue, B. Zhang, Y. Chong, and J.-H. Jiang, Topological phenomena at defects in acoustic, photonic and solid-state lattices, Nat. Rev. Phys.5, 483 (2023)

    work page 2023

  26. [26]

    J. Lu, W. Deng, X. Huang, M. Ke, and Z. Liu, Non- Hermitian Topological Phononic Metamaterials, Adv. Mater.37, 10.1002/adma.202307998 (2025)

    work page doi:10.1002/adma.202307998 2025

  27. [27]

    Zhang, Y

    T. Zhang, Y. Jiang, Z. Song, H. Huang, Y. He, Z. Fang, H. Weng, and C. Fang, Catalogue of topological elec- tronic materials, Nature566, 475 (2019)

    work page 2019

  28. [28]

    Fu, Topological Crystalline Insulators, Phys

    L. Fu, Topological Crystalline Insulators, Phys. Rev. Lett.106, 106802 (2011)

    work page 2011

  29. [29]

    Slager, A

    R.-J. Slager, A. Mesaros, V. Juričić, and J. Zaanen, The space group classification of topological band-insulators, Nat. Phys.9, 98 (2013)

    work page 2013

  30. [30]

    Y. Liu, S. Leung, F.-F. Li, Z.-K. Lin, X. Tao, Y. Poo, and J.-H. Jiang, Bulk–disclination correspondence in topolog- ical crystalline insulators, Nature589, 381 (2021)

    work page 2021

  31. [31]

    C. W. Peterson, T. Li, W. Jiang, T. L. Hughes, and G. Bahl, Trapped fractional charges at bulk defects in topological insulators, Nature589, 376 (2021)

    work page 2021

  32. [32]

    S. Liu, A. Vishwanath, and E. Khalaf, Shift Insulators: Rotation-Protected Two-Dimensional Topological Crys- talline Insulators, Phys. Rev. X9, 031003 (2019)

    work page 2019

  33. [33]

    K. W. Desmond and E. R. Weeks, Random close packing of disks and spheres in confined geometries, Phys. Rev. E80, 051305 (2009)

    work page 2009

  34. [34]

    Fortini and M

    A. Fortini and M. Dijkstra, Phase behaviour of hard spheres confined between parallel hard plates: manip- ulation of colloidal crystal structures by confinement, J. Phys. Condens. Matter18, L371 (2006)

    work page 2006

  35. [35]

    Geigenfeind, S

    T. Geigenfeind, S. Rosenzweig, M. Schmidt, and D. de Las Heras, Confinement of two-dimensional rods in slit pores and square cavities., J. Chem. Phys.142, 174701 (2015)

    work page 2015

  36. [36]

    Y. Li, H. Miao, H. Ma, and J. Z. Y. Chen, Topological defects of tetratic liquid-crystal order on a soft spherical surface, Soft Matter9, 11461 (2013)

    work page 2013

  37. [37]

    G. Zhu, L. Gao, Z. Xu, X. Dai, X. Zhang, and L.-T. Yan, Entropy-Driven Unconventional Crystallization of Spher- ical Colloidal Nanocrystals Confined in Wide Cylinders, Nano Lett.21, 8439 (2021)

    work page 2021

  38. [38]

    González-Pinto, F

    M. González-Pinto, F. Borondo, Y. Martínez-Ratón, and E. Velasco, Clustering in vibrated monolayers of granular rods, Soft Matter13, 2571 (2017). 8

    work page 2017

  39. [39]

    See supplemental material

    work page

  40. [40]

    D. Wang, M. Hermes, R. Kotni, Y. Wu, N. Tasios, Y. Liu, B. de Nijs, E. B. van der Wee, C. B. Murray, M. Dijk- stra, and A. van Blaaderen, Interplay between spherical confinement and particle shape on the self-assembly of rounded cubes, Nat. Commun.9, 2228 (2018)

    work page 2018