Abstract
To better understand the pattern phase transition of both physical and biological systems, we investigate a two-dimensional spin particle lattice system using statistical mechanics methods together with XY model governed by Hamiltonian equations of motion. By tweaking the coupling strength and the intensity of the generalization field, we observe phase transitions among four patterns of spin particles, i.e., vortex, ferromagnet, worm and anti-ferromagnet. In addition, we analyze the effect of space boundaries on the formations of vortex and worm. Considering the inherent dynamics of individual particles, we revealed the forming mechanism of such phase transitions, which provides a new perspective for understanding the emergence of phase transition of spin particles systems.
1 Introduction
Various types of phase transition are appealing features of both physical and biological particle systems, which have been extensively studied during recent decades. In general, pattern phase transition means a system composed of a large number of interacting particles undergoes a transition from one pattern phase to another responding to the change of one or more external parameters. There are abundant pattern phase transitions in natural, biological, chemical and physical systems, such as bird flocks [1–4], insects [5–7], bacterial colonies [8–10], fish schools and shoals [11–13], groups of mammals and crowds [14–16], crystals and superfluids [17, 18], etc. Researchers in the fields of physics, control engineering, system science, artificial intelligence and computer science recently have become more and more interested in the appealing pattern phase transition for such abundant particle systems.
Revealing the mechanism of phase transition is fairly helpful to understand complex physical and biological collective behaviors. The various pattern phase inspired by self-propelled particle systems are highly dependent on the interaction between particles, but the interaction mechanism of the phase transition is unknown, and one of the difficulties is the lack of the ability to predict the final equilibrium state and its stability only through the interaction between self-driven particles. In biological, physical and engineering multi-agent systems, researchers have discovered a large number of phase transitions between different modes. Among these investigations, Carrillo et al [19] researched the formation conditions of single and double rotating mills. Single-vortex configurations can be observed under generally random initial conditions, however, double-vortex configurations (i.e., half of the particles rotating clockwise and the other half counterclockwise) are only possible under carefully chosen initial conditions. Birnir [20] found the circling and flocking solutions of ordinary differential equations derived from the Vicsek model. Building on the work of Gautrais et al. [21], Calovi et al. [22] showed a switched phase diagram between the circular (milling) patterns and migratory patterns of fish schools by varying the weights of two control terms in the model. Cheng et al. [23] discovered pattern transitions among gaseous, liquid, and crystalline flocks by slightly varying the zero crossing slope of the inter-particle interaction function.
Analogous to kinetic systems, study of structure and phase transitions like lattices has drawn increasing attentions in recent years. The well-known statistical-mechanical models of particles like XY model and Ising model have intrinsic phase transitions for regular and irregular lattices [24]. From the perspective of statistical physics, the XY model with Hamiltonian equations of motion is probably one of the most concise models which have continuous degrees of freedom. In the field of condensed matter physics, it is a model which realized the phase transition of superfluid , High- cuprate and Josephson junction arrays [25, 26]. The system undergoes the well-known Berezinsky-Kosterlitz-Thouless transition in two-dimensional spaces [27], whereas experiencing a second-order phase transition in three-dimensional spaces [28]. Nowadays, the study of physical particle systems has shed some lights onto both physical and biological collective motion investigation [29]. Interestingly, in recent years, researchers have found that the phase transitions that exist in natural biological or physical systems have similarities, and it is possible to learn from each other in the theoretical study of phase transition mechanisms. For instance, by slightly varying the vision range in a group of freely moving particles, Cheng et al. [30] discovered various pattern phase transitions between crystalline, liquid, gaseous, and mill-liquid coexistence states. Referring to a lattice-gas model for superfluid helium, Attanasi et al. [4] established a new model nourished by conservation-law and spontaneous symmetry breaking principle that nicely explains the universal fast synchronization behaviors in starling flocks.
With the pioneer efforts devoted to the collective motions on XY model, quite a few fascinating collective motion patterns are observed with the XY model, like flocking and torus. In this letter, we seek to refine the kinetic description of XY model which helps reveal the phase transition principles of XY model with Hamiltonian equations of motion. Following the dynamics of the XY model, we have meticulously described the role of model parameters in governing the state of the system so as to figure out the actual physical meaning of the parameters. By tweaking model parameters, a series of phase transitions emerge among vortex, ferromagnet, worm and anti-ferromagnet patterns.
2 Methods and results
The XY model has been widely applied to a large volume of physical systems like superfluids, nematic liquid crystals, electron nematics, planar magnets, and among others [25–28]. In this letter, we study an XY model considering an external field effects with Hamiltonian equations of motion and focus on the statistical properties from the perspective of particles dynamics. In the presence of a -fold generalization field [31], we consider the XY model with spin in the direction on a square lattice with particles.where is the intensity of a -fold generalization field and denotes the coupling strength within nearest neighbor. When and , the first term denotes a uniaxial and -fold generalization field, respectively [31]. is a fixed constant representing the direction of the -fold generalization field. The and label the lattice sites on a 2D square lattice, denotes a nearest-neighbor set of the particle , which is one unit away from particle .
Consider the -th particle’s spin, i.e., , the spin becomes rotator, and hence the XY model turns into a system of coupled rotators [32, 33]. By virtue of the process of taking the negative gradient for the Hamiltonian system in [34], the first derivative for the -th object is obtained by taking negative gradient for system (1). Considering the effects of external noises, the spin particles obey the stochastic ordinary equation below,where , , is the external Gaussian white noise, and is the corresponding magnitude.
Let , the phases of spin are quantified by three order parameters, i.e., the consensus state order of spin directionthe parallel state order quantifying the parallel tendencyand the vortex state order quantifying the vortex tendencywhere is the set of positions with vortices, , “” denotes cross product, and is the number of a set .
Apparently, is obtained by taking the norm of the sum of all particle directions and then averaging it, which is crucial for quantifying the consistency of spin directions and obtaining the alignment degree of the entire system. That is, indicates that the spin particles reach direction consensus state. is calculated by averaging the magnitudes of the outer products of the directions of any two adjacent particles, which is essential for representing the parallel state of adjacent particles and obtaining the parallel trend of the entire system. That is, indicates that the spin particles reach direction parallel state. Moreover, means that vortex state occurs, is derived by summing the directions of four adjacent particles and then averaging them, which is crucial for judging the vortex state of these four particles. Therefore, is vital for measuring the vortex state sequence and reflecting vortex formation by quantifying local rotational arrangements. That is, means that spin particles reach standard vortex state. Therefore, these parameters are chosen because they provide a straightforward method to differentiate between the various patterns emerging from our system. In order to show the order parameter characteristics of the four phases more clearly, the schematic diagram Figure 1 is exhibited. and mean the standard vortex (Figure 1A) and ferromagnet phase (Figure 1B), respectively. It is worth noting that we define a parallel state order parameter to distinguish worm and anti-ferromagnetic phases by combining the consensus and vortex state order parameters together. More precisely, and implies the worm phase (Figure 1C), whereas and represent the anti-ferromagnet phase (Figure 1D).
FIGURE 1
In the numerical simulations, particles are randomly initialized as . The parameters are picked as: , , and . In order to highlight the phase transition, we demonstrate the spin motional phases in the space spanned by the two parameters and .The parameters are presented on a logarithmic scale to encompass a wide range of values, this approach allows us to explore the behavior of the system across multiple orders of magnitude. This is important for understanding the fundamental properties and scalability of the system. With different combinations of and , three phase transitions emerge from vortex (Figures 2D–H) to ferromagnet (Figure 2A), to worm (Figure 2B), and then to anti-ferromagnet (Figure 2C). More precisely, in the ferromagnet phase, the spin particles form parallel arrays and direction consensus, whereas in the anti-ferromagnet phase, the spin particles form anti-parallel pattern. In the worm phase, spin particles form a curved line arrangement configuration pattern like a worm. By contrast, in the vortex pattern, there are four adjacent spin particles that form an ordered rotation, crossover, relative arrangement, etc. In the vortex pattern, multiple vortex modes occur simultaneously.
FIGURE 2
Figure 3 demonstrates the phase transitions among vortex, ferromagnet, worm and anti-ferromagnet patterns. At the top-left corner of the - space, vortex pattern emerges by four adjacent spin particles, where the parameter , implying the appearance of vortices. Appealingly, five distinct vortex modes have been observed in this phase, namely, rotational vortex, positive vortex, cross vortex, relative vortex and anti-vortex, as illustrated in Figures 2D–H. With increasing parameter , that is, the upper middle of Figure 3, spin particles finally reach direction synchronization, which naturally forms a ferromagnet phase, i.e., . In the middle diagonal stripe of Figure 3, group direction no longer reaches synchronization (i.e., ), the parameter and implying the emergence of worm phase. By keeping increasing the parameter , that is, the bottom-right corner of Figure 3, anti-ferromagnet phase (i.e., and ) is yielded. In order to show the above four phase transition processes more vividly, we provide numerical simulation videos in the supplementary material and the link http://imds.aia.hust.edu.cn/info/1247/2403.htm.
FIGURE 3
In order to investigate more deeply into the principles governing the pattern phase transitions, we present three heat maps of , and in Figures 4–6 along increasing parameters and , respectively. Each point is an average over 50 independent runs, with each run lasting running steps guaranteeing the attainment of steady states.Note that the choice of stepsis a conservative estimate to ensure that steady states are attained uniformly across the entire parameter range.
FIGURE 4
FIGURE 5
FIGURE 6
Significantly, the order parameter implies ferromagnet. It is observed from Figure 4 that the red region corresponds to the ferromagnetic phase . It is noteworthy that the value presented here is an average over 50 independent runs, thus there is a very small probability that there will be a non-ferromagnetic phase (i.e., ). Consequently, the red region does not exactly equal 1, especially near the boundaries of pattern phase transition where the probability of a non-ferromagnetic phase increases, leading to the average value less than 1. Therefore, the pattern phase discussed here is a ferromagnetic phase in a probabilistic sense. The order parameter means that the spin particles reach parallel state, from Figure 5, the parallel state in the blue region can be observed, which serves as an important basis for identifying the worm phase and the anti-ferromagnet phase . Notably, in the bottom-right corners of Figure 5, the spin motion of particles in the system is governed by the 2-fold generalization field term, i.e., . However, there are two stable equilibria, namely, and , so an anti-ferromagnet effect will be formed. The order parameter implies vortex phase. It is observed from Figure 6 that the red region corresponds to the vortex phase. Meanwhile, in the middle-upper diagonal striped region of Figure 6, the system undergoes a phase transition from a vortex to a ferromagnet . In conclusion, by combining Figures 4–6, a sequential phase transition process can be observed, including the vortex phase , the ferromagnetic phase , the worm phase , and finally the anti-ferromagnetic phase .
To investigate more deeply into the phase transition mechanism, we conduct the effects of the side length parameter on vortex-worm phase transitions. Along ascending from 5 to 20, as shown in Figure 7A, the vortex phase begins to appear. Upon reaching , the occurrence frequency of vortex phase exceeds that of ferromagnet phase, where occurrence frequency refers to the number of times a specific phase (e.g., vortex or worm) is observed across multiple simulations (same below). Until , only the vortex phase exists, whose number increases rapidly with rising (see Figures 7B, C), where vortex number quantifies the number of vortices formed in the system during a single simulation, and the error bars representing the standard deviation across 200 independent runs (same below). Analogously, in Figure 8A, when , the occurrence frequency of the worm phase exceeds that of the parallel patterns (i.e., ferromagnet and anti-ferromagnet phases). Upon reaching only the worm phase exists. Figure 8B shows the evolution of the parallel order parameter value along increasing , which can quantitatively identifies the worm phase. A plausible explanation is given below. From the spin dynamics XY model (2), it can be observed that the model has many stable equilibria, including spin particles aligned in a parallel or vertical direction, relative to each other. However, according to the characteristics of the dynamic model (2), spin particles on the four corners and four sides of a square lattice have only two and three neighbors, respectively, which forms a boundary effect. Specifically, in the scenario of random initial values, the spin particles on the boundary and their neighbors will be prone to reach the consensus of heading directions, which drives the inner side spin particles to reach the same direction. However, with increasing number of the spin particles, the proportion of spin particles on the boundary decreases, and hence the influence of the boundary effect on the inner side particles will be weakened. In such a situation, vortex and worm patterns can be formed.
FIGURE 7
FIGURE 8
So far, the phases and their transitions discussed above are noise-free, i.e., . Next, we conduct an investigation on the effect of on the formation of vortex and worm phase in model (2). As shown in Figure 9A, with the increase of , the number of vortices gradually decreases until the vortex phase is totally ruined by noise. It is observed that vortex phase can be maintained for . In the blue curves of Figure 9B, it is shown that with moderate increase of external noise, a new phase appeared, that is, when , would approach 1, but could not exactly reach 1 due to the influence of the external noise, i.e., , which is called quasi ferromagnet phase. As shown in Figure 9B, when , the quasi ferromagnet phase appears, and the variance of in the worm phase grows larger.
FIGURE 9
3 Conclusion
This letter investigates the pattern phase transition mechanism for a class of the XY model with Hamiltonian equations of motion. Three phase transitions among four patterns, i.e., vortex, ferromagnet, worm and anti-ferromagnet, are revealed simply by tweaking two parameters and . Specifically, the occurrence frequency of vortex phase exceeds that of ferromagnet phase for a sufficiently large space. This study is expected to provide an insight for understanding the emergence of single vortex and tight vortices pairs in pattern phase transition of spin particle groups. The observed phase transition may shed some lights onto the self-assembly dynamics analysis of magnets, electron nematics, and quantum gases.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: http://imds.aia.hust.edu.cn/info/1247/2403.htm.
Author contributions
YW: Writing–original draft, Writing–review and editing. JY: Writing–original draft. BX: Writing–review and editing. YZ: Writing–review and editing. DC: Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China under Grant. 62273053, Fundamental Research Funds for the Central Universities of China (BLX202128), China Postdoctoral Science Foundation under Grant No. 2024M754222, and the Natural Science Foundation of Shanxi Province under Grant No. 2024JC-YBQN-0663.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
FlackANagyMFiedlerWCouzinIDWikelskiM. From local collective behavior to global migratory patterns in white storks. Science (2018) 360:911–4. 10.1126/science.aap7781
2.
ChenDVicsekTLiuXZhouTZhangH-T. Switching hierarchical leadership mechanism in homing flight of pigeon flocks. Europhysics Lett (2016) 114:60008. 10.1209/0295-5075/114/60008
3.
LukemanRLiY-XEdelstein-KeshetL. Inferring individual rules from collective behavior. Proc Natl Acad Sci (2010) 107:12576–80. 10.1073/pnas.1001763107
4.
AttanasiACavagnaACastelloLDGiardinaIGrigeraTSJelićAet alInformation transfer and behavioural inertia in starling flocks. Nat Phys (2014) 10:691–6. 10.1038/nphys3035
5.
BuhlJSumpterDCouzinIHaleJJDesplandEMillerERet alFrom disorder to order in marching locusts. Science (2006) 312:1402–6. 10.1126/science.1125142
6.
FranksNRPrattSCMallonEBBrittonNFSumpterDJT. Information flow, opinion polling and collective intelligence in house–hunting social insects. Philosophical Trans R Soc Lond Ser B: Biol Sci (2002) 357:1567–83. 10.1098/rstb.2002.1066
7.
BazaziSBuhlJHaleJJAnsteyMLSwordGASimpsonSJet alCollective motion and cannibalism in locust migratory bands. Curr Biol (2008) 18:735–9. 10.1016/j.cub.2008.04.035
8.
CuratoloAZhouNZhaoYLiuCDaerrATailleurJet alCooperative pattern formation in multi-component bacterial systems through reciprocal motility regulation. Nat Phys (2020) 16:1152–7. 10.1038/s41567-020-0964-z
9.
SzaboBSzollosiGJGönciBJurányiZSelmecziDVicsekT. Phase transition in the collective migration of tissue cells: experiment and model. Phys Rev E—Statistical, Nonlinear, Soft Matter Phys (2006) 74:061908. 10.1103/PhysRevE.74.061908
10.
DeisboeckTSCouzinID. Collective behavior in cancer cell populations. Bioessays (2009) 31:190–7. 10.1002/bies.200800084
11.
WardAJWSumpterDJTCouzinIDHartPJBKrauseJ. Quorum decision-making facilitates information transfer in fish shoals. Proc Natl Acad Sci (2008) 105:6948–53. 10.1073/pnas.0710344105
12.
HeinAMGilMATwomeyCRCouzinIDLevinSA. Conserved behavioral circuits govern high-speed decision-making in wild fish shoals. Proc Natl Acad Sci (2018) 115:12224–8. 10.1073/pnas.1809140115
13.
HoareDCouzinIDGodinJ-GJKrauseJ. Context-dependent group size choice in fish. Anim Behav (2004) 67:155–64. 10.1016/j.anbehav.2003.04.004
14.
HelbingDFarkasIVicsekT. Simulating dynamical features of escape panic. Nature (2000) 407:487–90. 10.1038/35035023
15.
FischhoffISundaresanSRCordingleyJLarkinHMSellierMJRubensteinDI. Social relationships and reproductive state influence leadership roles in movements of plains zebra, equus burchellii. Anim Behav (2007) 73:825–31. 10.1016/j.anbehav.2006.10.012
16.
SarovaRSpinkaMPanamáJSimecekP (2010). Graded leadership by dominant animals in a herd of female beef cattle on pasture. Anim Behav79, 1037–45. 10.1016/j.anbehav.2010.01.019
17.
LiuHLvJChenQ. Ising-like phase transition in the fully frustrated xyz model with weak disorder. Europhysics Lett (2008) 84:66004. 10.1209/0295-5075/84/66004
18.
JesariewDIlczyszynMMPietraszkoA. The crystal structure and the phase transitions of pyridinium trifluoromethanesulfonate. Mater Res Express (2014) 1:015705. 10.1088/2053-1591/1/1/015705
19.
CarrilloJAD’OrsognaMPanferovV. Double milling in self-propelled swarms from kinetic theory. Kinetic Relat Models (2009) 2:363–78. 10.3934/krm.2009.2.363
20.
BirnirB. An ode model of the motion of pelagic fish. J Stat Phys (2007) 128:535–68. 10.1007/s10955-007-9292-2
21.
GautraisJGinelliFFournierRBlancoSSoriaMChatéHet alDeciphering interactions in moving animal groups. PLoS Comput Biol (2012) 8:e1002678. 10.1371/journal.pcbi.1002678
22.
CaloviDSLopezUNgoSSireCChatéHTheraulazG. Swarming, schooling, milling: phase diagram of a data-driven fish school model. New J Phys (2014) 16:015026. 10.1088/1367-2630/16/1/015026
23.
ChengZ-XSunBWangS-JGaoYZhangYMZhouHXet alNuclear factor-κb–dependent epithelial to mesenchymal transition induced by HIF-1α activation in pancreatic cancer cells under hypoxic Conditionsα activation in pancreatic cancer cells under hypoxic conditions. PLoS One (2011) 6:e23752. 10.1371/journal.pone.0023752
24.
BaekSKMinnhagenPKimBJ. Phase transition of xyz model in heptagonal lattice. Europhysics Lett (2007) 79:26002. 10.1209/0295-5075/79/26002
25.
de SouzaLCde SouzaAJFLyraML. Hamiltonian short-time critical dynamics of the three-dimensional xyz model. Phys Rev E (2019) 99:052104. 10.1103/PhysRevE.99.052104
26.
MonKTeitelS. Phase coherence and nonequilibrium behavior in josephs junction arrays. Phys Rev Lett (1989) 62:673–6. 10.1103/PhysRevLett.62.673
27.
KosterlitzJM. The critical properties of the two-dimensional xyz model. J Phys C: Solid State Phys (1974) 7:1046–60. 10.1088/0022-3719/7/6/005
28.
LiY-HTeitelS. Finite-size scaling study of the three-dimensional classical xyz model. Phys Rev B (1989) 40:9122–5. 10.1103/PhysRevB.40.9122
29.
VicsekTZafeirisA. Collective motion. Phys Rep (2012) 517:71–140. 10.1016/j.physrep.2012.03.004
30.
ChengZChenZVicsekTChenDZhangHT. Pattern phase transitions of self-propelled particles: gases, crystals, liquids, and mills. New J Phys (2016) 18:103005. 10.1088/1367-2630/18/10/103005
31.
MinchauBJPelcovitsRA. Two-dimensional xyz model in a random uniaxial field. Phys Rev B (1985) 32:3081–7. 10.1103/PhysRevB.32.3081
32.
AsadAZhengB. Non-equilibrium critical dynamics of the two-dimensional xyz model with Hamiltonian equations of motion. J Phys A: Math Theor (2007) 40:9957–68. 10.1088/1751-8113/40/33/001
33.
LeonciniXVergaADRuffoS. Hamiltonian dynamics and the phase transition of the xyz model. Phys Rev E (1998) 57:6377–89. 10.1103/PhysRevE.57.6377
34.
MatteoPMarconiUMBMaggiC. Effective equilibrium picture in the xyz model with exponentially correlated noise. Phys Rev E (2018) 97:022605. 10.1103/PhysRevE.97.022605
35.
BasakSDahmenKCarlsonE. Period multiplication cascade at the order-by-disorder transition in uniaxial random field xyz magnets. Nat Commun (2020) 11:4665. 10.1038/s41467-020-18270-6
Summary
Keywords
collective motions, pattern phase transition, spin particle, Hamiltonian equations, complex networks
Citation
Wu Y, Yan J, Xu B, Zheng Y and Chen D (2024) Pattern phase transition of spin particle lattice system. Front. Phys. 12:1472564. doi: 10.3389/fphy.2024.1472564
Received
29 July 2024
Accepted
21 October 2024
Published
08 November 2024
Volume
12 - 2024
Edited by
Sanja Janicevic, University of Kragujevac, Serbia
Reviewed by
Tapas Bar, Catalan Institute of Nanoscience and Nanotechnology (CIN2), Spain
Diana Thongjaomayum, Tezpur University, India
Updates
Copyright
© 2024 Wu, Yan, Xu, Zheng and Chen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yili Zheng, zhengyili@bjfu.edu.cn; Duxin Chen, chendx@seu.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.