We present a family of graphical representations for the O(N)spin model,where N≥1 represents the spin dimension,and N=1,2,3 corresponds to the Ising,XY and Heisenberg models,respectively.With an integer parameter 0≤...We present a family of graphical representations for the O(N)spin model,where N≥1 represents the spin dimension,and N=1,2,3 corresponds to the Ising,XY and Heisenberg models,respectively.With an integer parameter 0≤ℓ≤N/2,each configuration is the coupling of ℓ copies of subgraphs consisting of directed flows and N−2ℓ copies of subgraphs constructed by undirected loops,which we call the XY and Ising subgraphs,respectively.On each lattice site,the XY subgraphs satisfy the Kirchhoff flow-conservation law and the Ising subgraphs obey the Eulerian bond condition.Then,we formulate worm-type algorithms and simulate the O(N)model on the simple-cubic lattice for N from 2 to 6 at all possibleℓ.It is observed that the worm algorithm has much higher efficiency than the Metropolis method,and,for a given N,the efficiency is an increasing function ofℓ.Besides Monte Carlo simulations,we expect that these graphical representations would provide a convenient basis for the study of the O(N)spin model by other state-of-the-art methods like the tensor network renormalization.展开更多
To investigate the band structure is one of the key approaches to study the fundamental properties of a novel material.We report here the precision band mapping of a 2-dimensional(2D) spin-orbit(SO) coupling in an opt...To investigate the band structure is one of the key approaches to study the fundamental properties of a novel material.We report here the precision band mapping of a 2-dimensional(2D) spin-orbit(SO) coupling in an optical lattice.By applying the microwave spin-injection spectroscopy, the band structure and spin-polarization distribution are achieved simultaneously.The band topology is also addressed with observing the band gap close and re-open at the Dirac points.Furthermore, the lattice depth and the Raman coupling strength are precisely calibrated with relative errors in the order of 10^(-3).Our approach could also be applied for exploring the exotic topological phases with even higher dimensional system.展开更多
We simulate the bond and site percolation models on several three-dimensional lattices, including the diamond, body-centered cubic, and face-centered cubic lattices. As on the simple-cubic lattice [Phys. Rev. E, 2013,...We simulate the bond and site percolation models on several three-dimensional lattices, including the diamond, body-centered cubic, and face-centered cubic lattices. As on the simple-cubic lattice [Phys. Rev. E, 2013, 87(5): 052107], it is observed that in comparison with dimensionless ratios based on cluster-size distribution, certain wrapping probabilities exhibit weaker finite-size corrections and are more sensitive to the deviation from percolation threshold Pc, and thus provide a powerful means for determining Pc. We analyze the numerical data of the wrapping probabilities simultaneously such that universal parameters are shared by the aforementioned models, and thus significantly improved estimates of Pc are obtained.展开更多
We formulate an irreversible Markov chain Monte Carlo algorithm for the self-avoiding walk (SAW), which violates the detailed balance condition and satisfies tile balance condition. Its performance improves signific...We formulate an irreversible Markov chain Monte Carlo algorithm for the self-avoiding walk (SAW), which violates the detailed balance condition and satisfies tile balance condition. Its performance improves significantly compared to that of the Berretti-Sokal algorithm, which is a variant of the Metropolis Hastings method. The gained efficiency increases with spatial dimension (D), from approximately 10 times in 2D to approximately 40 times in 5D. We simulate the SAW on a 5D hypercubic lattice with periodic boundary conditions, for a linear system with a size up to L = 128, and confirm that as for the 5D Ising model, the finite-size scaling of the SAW is governed by renormalized exponents, υ^* = 2/d and γ/υ^* = d/2. The critical point is determined, which is approximately 8 times more precise than the best available estimate.展开更多
基金supported by the National Natural Science Foundation of China(under Grant No.12275263)the Innovation Program for Quantum Science and Technology(under Grant No.2021ZD0301900)the Natural Science Foundation of Fujian Province of China:2023J02032.
文摘We present a family of graphical representations for the O(N)spin model,where N≥1 represents the spin dimension,and N=1,2,3 corresponds to the Ising,XY and Heisenberg models,respectively.With an integer parameter 0≤ℓ≤N/2,each configuration is the coupling of ℓ copies of subgraphs consisting of directed flows and N−2ℓ copies of subgraphs constructed by undirected loops,which we call the XY and Ising subgraphs,respectively.On each lattice site,the XY subgraphs satisfy the Kirchhoff flow-conservation law and the Ising subgraphs obey the Eulerian bond condition.Then,we formulate worm-type algorithms and simulate the O(N)model on the simple-cubic lattice for N from 2 to 6 at all possibleℓ.It is observed that the worm algorithm has much higher efficiency than the Metropolis method,and,for a given N,the efficiency is an increasing function ofℓ.Besides Monte Carlo simulations,we expect that these graphical representations would provide a convenient basis for the study of the O(N)spin model by other state-of-the-art methods like the tensor network renormalization.
基金supported by the National Key R&D Program of China (2016YFA0301601 and 2016YFA0301604)the National Natural Science Foundation of China (11674301, 11761161003, and 11625522)the Thousand-Young-Talent Program of China
文摘To investigate the band structure is one of the key approaches to study the fundamental properties of a novel material.We report here the precision band mapping of a 2-dimensional(2D) spin-orbit(SO) coupling in an optical lattice.By applying the microwave spin-injection spectroscopy, the band structure and spin-polarization distribution are achieved simultaneously.The band topology is also addressed with observing the band gap close and re-open at the Dirac points.Furthermore, the lattice depth and the Raman coupling strength are precisely calibrated with relative errors in the order of 10^(-3).Our approach could also be applied for exploring the exotic topological phases with even higher dimensional system.
文摘We simulate the bond and site percolation models on several three-dimensional lattices, including the diamond, body-centered cubic, and face-centered cubic lattices. As on the simple-cubic lattice [Phys. Rev. E, 2013, 87(5): 052107], it is observed that in comparison with dimensionless ratios based on cluster-size distribution, certain wrapping probabilities exhibit weaker finite-size corrections and are more sensitive to the deviation from percolation threshold Pc, and thus provide a powerful means for determining Pc. We analyze the numerical data of the wrapping probabilities simultaneously such that universal parameters are shared by the aforementioned models, and thus significantly improved estimates of Pc are obtained.
基金Acknowledgements This work was supported by the National Natural Science Foundation of China under Grant Nos. 11275185 and 11625522, and the Open Project Program of State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, China (No. Y5KF191CJ1). Y. Deng acknowledges the Ministry of Education (of China) for the Fundamental Research Funds for the Central Universities under Grant No. 2340000034.
文摘We formulate an irreversible Markov chain Monte Carlo algorithm for the self-avoiding walk (SAW), which violates the detailed balance condition and satisfies tile balance condition. Its performance improves significantly compared to that of the Berretti-Sokal algorithm, which is a variant of the Metropolis Hastings method. The gained efficiency increases with spatial dimension (D), from approximately 10 times in 2D to approximately 40 times in 5D. We simulate the SAW on a 5D hypercubic lattice with periodic boundary conditions, for a linear system with a size up to L = 128, and confirm that as for the 5D Ising model, the finite-size scaling of the SAW is governed by renormalized exponents, υ^* = 2/d and γ/υ^* = d/2. The critical point is determined, which is approximately 8 times more precise than the best available estimate.
基金supported by the National Natural Science Foundation of China (11874340)the National Key R&D Program of China (2018YFA0306501)+2 种基金the CASthe Anhui Initiative in Quantum Information Technologiesthe Fundamental Research Funds for the Central Universities (WK2340000081)