Many phenomena show that in a favorable circumstance an agent still has an updating possibility, and in an unfavor- able circumstance an agent also has a possibility of holding its own state and reselecting its neighb...Many phenomena show that in a favorable circumstance an agent still has an updating possibility, and in an unfavor- able circumstance an agent also has a possibility of holding its own state and reselecting its neighbors. To describe this kind of phenomena an Ising model on evolution networks was presented and used for consensus formation and separation of opinion groups in human population. In this model the state-holding probability p and selection-rewiring probability q were introduced. The influence of this mixed dynamics of spin flips and network rewiring on the ordering behavior of the model was investigated, p hinders ordering of opinion networks and q accelerates the dynamical process of networks. Influence of q on the ordering and separating stems from its effect on average path length of networks.展开更多
An overview of the mathematical structure of the three-dimensional(3D) Ising model is given from the points of view of topology,algebra,and geometry.By analyzing the relationships among transfer matrices of the 3D I...An overview of the mathematical structure of the three-dimensional(3D) Ising model is given from the points of view of topology,algebra,and geometry.By analyzing the relationships among transfer matrices of the 3D Ising model,Reidemeister moves in the knot theory,Yang-Baxter and tetrahedron equations,the following facts are illustrated for the 3D Ising model.1) The complex quaternion basis constructed for the 3D Ising model naturally represents the rotation in a(3+1)-dimensional space-time as a relativistic quantum statistical mechanics model,which is consistent with the 4-fold integrand of the partition function obtained by taking the time average.2) A unitary transformation with a matrix that is a spin representation in 2 n·l·o-space corresponds to a rotation in 2n·l·o-space,which serves to smooth all the crossings in the transfer matrices and contributes the non-trivial topological part of the partition function of the 3D Ising model.3) A tetrahedron relationship would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model,and its existence is guaranteed by the Jordan algebra and the Jordan-von Neumann-Wigner procedures.4) The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases φx,φy,and φz.The relationship with quantum field and gauge theories and the physical significance of the weight factors are discussed in detail.The conjectured exact solution is compared with numerical results,and the singularities at/near infinite temperature are inspected.The analyticity in β=1/(kBT) of both the hard-core and the Ising models has been proved only for β〉0,not for β=0.Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.展开更多
The review paper by Zhang Zhi-Dong (Zhang Z D 2013 Chin. Phys. B 22 030513, arXiv:1305.2956) contains many errors and is based on several earlier works that are equally wrong.
We use a new updated algorithm scheme to investigate the critical behaviour of the two-dimensional ferromagnetic Ising model on a triangular lattice with the nearest neighbour interactions. The transition is examined ...We use a new updated algorithm scheme to investigate the critical behaviour of the two-dimensional ferromagnetic Ising model on a triangular lattice with the nearest neighbour interactions. The transition is examined by generating accurate data for lattices with L=8, 10, 12, 15, 20, 25, 30, 40 and 50. The updated spin algorithm we employ has the advantages of both a Metropolis algorithm and a single-update method. Our study indicates that the transition is continuous at Тc=3.6403(2). A convincing finite-size scaling analysis of the model yields ν=0.9995(21), β/ν=0.12400(17), γ/v=1.75223(22), γ^1/ν=1.7555(22), α/ν=0.00077(420) (scaling) and α/ν=0.0010(42) (hyperscaling). The present scheme yields more accurate estimates for all the critical exponents than the Monte Carlo method, and our estimates are shown to be in excellent agreement with their predicted values.展开更多
We use the quantum renormalization-group(QRG) method to study the entanglement and quantum phase transition(QPT) in the one-dimensional spin-1/2 Heisenberg-Ising model [Lieb E,Schultz T and Mattis D 1961 Ann.Phys....We use the quantum renormalization-group(QRG) method to study the entanglement and quantum phase transition(QPT) in the one-dimensional spin-1/2 Heisenberg-Ising model [Lieb E,Schultz T and Mattis D 1961 Ann.Phys.(N.Y.) 16 407].We find the quantum phase boundary of this model by investigating the evolution of concurrence in terms of QRG iterations.We also investigate the scaling behavior of the system close to the quantum critical point,which shows that the minimum value of the first derivative of concurrence and the position of the minimum scale with an exponent of the system size.Also,the first derivative of concurrence between two blocks diverges at the quantum critical point,which is directly associated with the divergence of the correlation length.展开更多
Recently, Shiet al. [2008 Phys. Left. A 372 5922] have studied the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field and presented the dynamic phase diagrams by using an e...Recently, Shiet al. [2008 Phys. Left. A 372 5922] have studied the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field and presented the dynamic phase diagrams by using an effective-field theory (EFT) and a mean-field theory (MFT). The MFT results are in conflict with those of the earlier work of Tome and de Oliveira, [1990 Phys. Rev. A 41 4251]. We calculate the dynamic phase diagrams and find that our results are similar to those of the earlier work of Tome and de Oliveira; hence the dynamic phase diagrams calculated by Shiet al. are incomplete within both theories, except the low values of frequencies for the MFT calculation. We also investigate the influence of external field frequency (w) and static external field amplitude (h0) for both MFT and EFT calculations. We find that the behaviour of the system strongly depends on the values of w and h0.展开更多
In order to calculate 3-dimensional Ising model,we develop a method to build a much smaller transfer matrix containing the largest eigenvalue from the original 2^(N_1N_2) × 2^(N_1N_2) matrix V_1. Firstly,the tran...In order to calculate 3-dimensional Ising model,we develop a method to build a much smaller transfer matrix containing the largest eigenvalue from the original 2^(N_1N_2) × 2^(N_1N_2) matrix V_1. Firstly,the transfer matrix V_1 is written as the linear combination of several basic vectors. Secondly,we divide the basic vectors into several subgroups. The multiplication of a basic vector and V_1 can be written as the linear combination of basic vectors from the same subgroup. Finally,we use a new transfer matrix V_2 to describe the relationship between basic vectors of the same subgroup.V_2 is much smaller than the original transfer matrix and contains the largest eigenvalue of V_1.We use this method to calculate the specific heat per atom Cpaand the magnetic momentum per atom Mpa. The results show that there exists a pair of temperature and magnetic field intensity where the specific heat gets to its maximum value. When N_1N_2 increases,the maximum value of specific heat becomes larger.展开更多
An elastic Ising model for a one-dimensional diatomic spin chain is proposed to explain the ferroelectricity induced by the collinear magnetic order with a low-excited energy state. A statistical theory based on this ...An elastic Ising model for a one-dimensional diatomic spin chain is proposed to explain the ferroelectricity induced by the collinear magnetic order with a low-excited energy state. A statistical theory based on this model is developed to calculate the electrical and magnetic properties of Ca3CoMnO6, a typical quasi-one-dimensional diatomic spin chain system. The calculated ferroelectric polarization and dielectric susceptibility show a good agreement with recently reported data on Ca3Co2-xMnxO6 (x ≈ 0.96) (Phys. Rev. Lett. 100 047601 (2008)), although the predicted magnetic susceptibility does not coincide well with experiment. We also address the rationality and deficiency of this model by including a first-order correction which improves the consistency between the model and experiment.展开更多
We use the Monte Carlo method to study an antiferromagnetical Ising spin system on a centred honeycomb lattice, which is composed of two kinds of 1/2 spin particles A and B. There exist two different bond energies JA-...We use the Monte Carlo method to study an antiferromagnetical Ising spin system on a centred honeycomb lattice, which is composed of two kinds of 1/2 spin particles A and B. There exist two different bond energies JA-A and JA--B in this lattice. Our study is focused on how the ratio of JA-B to JA--A influences the critical behaviour of this system by analysing the physical quantities, such as the energy, the order parameter, the specific heat, susceptibility, etc each as a function of temperature for a given ratio of JA-B to JA-A. Using these results together with the finite-size scaling method, we obtain a phase diagram for the ratio JA-B / JA--A. This work is helpful for studying the phase transition problem of crystals composed of compounds.展开更多
This paper proposes a scheme based on the Potts and Ising models for simulating polarization switching of polycrystalline ferroelectrics using the Monte Carlo method. The polycrystalline texture with different average...This paper proposes a scheme based on the Potts and Ising models for simulating polarization switching of polycrystalline ferroelectrics using the Monte Carlo method. The polycrystalline texture with different average grain size is produced from the Potts model. Then Ising model is implemented in the polycrystalline texture to produce the domain pattern and hysteresis loop. The domain patterns and hysteresis loops have been obtained for polycrystalline texture with different average grain size. From the results of domain pattern evolution process under an applied electric field using this scheme, an extended domain, which covers more than one grain with polarization aligned roughly in the same direction, has been observed during the polarization reversal. This scheme can well reproduce the basic properties of polycrystalline ferroelectrics and is a valuable tool for exploring the physical properties of polycrystalline ferroelectrics.展开更多
We present a study of the dynamic behavior of a two-sublattice spin-5/2 Ising model with bilinear and crystal-field interactions in the presence of a time-dependent oscillating external magnetic field on alternating l...We present a study of the dynamic behavior of a two-sublattice spin-5/2 Ising model with bilinear and crystal-field interactions in the presence of a time-dependent oscillating external magnetic field on alternating layers of a hexagonal lattice by using the Glauber-type stochastic dynamics.The lattice is formed by alternate layers of spins σ=5/2 and S=5/2.We employ the Glauber transition rates to construct the mean-field dynamic equations.First,we investigate the time variations of the average sublattice magnetizations to find the phases in the system and then the thermal behavior of the dynamic sublattice magnetizations to characterize the nature(first-or second-order) of the phase transitions and to obtain the dynamic phase transition(DPT) points.We also study the thermal behavior of the dynamic total magnetization to find the dynamic compensation temperature and to determine the type of the dynamic compensation behavior.We present the dynamic phase diagrams,including the dynamic compensation temperatures,in nine different planes.The phase diagrams contain seven different fundamental phases,thirteen different mixed phases,in which the binary and ternary combination of fundamental phases and the compensation temperature or the L-type behavior strongly depend on the interaction parameters.展开更多
It was first proposed that an extended chain Ising (ECI) model contains the Ising chain model, single spin doublewell potentials and a pure phonon heat bath of a specific energy exchange with the spins. The extensio...It was first proposed that an extended chain Ising (ECI) model contains the Ising chain model, single spin doublewell potentials and a pure phonon heat bath of a specific energy exchange with the spins. The extension method is easy to apply to high dimensional cases. Then the single spin-flip probability (rate) of the ECI model is deduced based on the Boltzmann principle and general statistical principles of independent events and the model is simplified to an extended chain Glauber-Ising (ECGI) model. Moreover, the relaxation dynamics of the ECGI model were simulated by the Monte Carlo method and a comparison with the predictions of the special chain Clauber-Ising (SCGI) model was presented. It was found that the results of the two models are consistent with each other when the Ising chain length is large enough and temperature is relative low, which is the most valuable case of the model applications. These show that the ECI model will provide a firm physical base for the widely used single spin-flip rate proposed by Glauber and a possible route to obtain the single spin-flip rate of other form and even the multi-spin-flip rate.展开更多
The Blume-Capel model in the presence of external magnetic field H has been simulated using a cellular automaton algorithm improved from the Creutz cellular automaton in three-dimension lattice. The field critical exp...The Blume-Capel model in the presence of external magnetic field H has been simulated using a cellular automaton algorithm improved from the Creutz cellular automaton in three-dimension lattice. The field critical exponent 5 is estimated using the power law relations and the finite size scaling functions for the magnetization and the susceptibility in the range -0.1≤ h = H/J ≤0. The estimated value of the field critical exponent 5 is in good agreement with the universal value (δ = 5) in three dimensions. The simulations are carried out on a simple cubic lattice under periodic boundary conditions.展开更多
In this paper it is shown that the thermodynamic limit of the partition function of the statistical models under consideration on a one-dimensional lattice with an arbitrary finite number of interacting neighbors is e...In this paper it is shown that the thermodynamic limit of the partition function of the statistical models under consideration on a one-dimensional lattice with an arbitrary finite number of interacting neighbors is expressed in terms of the principal eigenvalue of a matrix of finite size. The high sparseness of these matrices for any number of interactions makes it possible to perform an effective numerical analysis of the macro characteristics of these models.展开更多
Further to the investigation of the critical properties of the Potts model with <em>q</em> = 3 and 8 states in one dimension (1D) on directed small-world networks reported by Aquino and Lima, which present...Further to the investigation of the critical properties of the Potts model with <em>q</em> = 3 and 8 states in one dimension (1D) on directed small-world networks reported by Aquino and Lima, which presents, in fact, a second-order phase transition with a new set of critical exponents, in addition to what was reported in Sumour and Lima in studying Ising model on non-local directed small-world for several values of probability 0 < <em>P</em> < 1. In this paper the behavior of two models discussed previously, will be re-examined to study differences between their behavior on directed small-world networks for networks of different values of probability <em>P</em> = 0.1, 0.2, 0.3, 0.4 and 0.5 with different lattice sizes <em>L</em> = 10, 20, 30, 40, and 50 to compare between the important physical variables between Ising and Potts models on the directed small-world networks. We found in our paper that is a phase transitions in both Ising and Potts models depending essentially on the probability <em>P</em>.展开更多
The phase diagrams and the other crtical properties of S-3/2 Ising model in random longitudinal and transverse fields(RLIM) are dicussed with the approximate scheme combined by mean-field renormalization group theory(...The phase diagrams and the other crtical properties of S-3/2 Ising model in random longitudinal and transverse fields(RLIM) are dicussed with the approximate scheme combined by mean-field renormalization group theory(MFRG) and the discretized path-integral representation(DPIR).展开更多
From the mathematical point of view, the modeling of epidemics—in other words, the spread of infectious diseases transmitted from individual to individual—is very similar to the modeling of the magnetic systems stud...From the mathematical point of view, the modeling of epidemics—in other words, the spread of infectious diseases transmitted from individual to individual—is very similar to the modeling of the magnetic systems studied by statistical physics. In this work, we use this analogy between mathematical epidemiology and statistical physics to study the classical mathematical model of epidemiology SI (Susceptible-Infected) approached through the Ising-Glauber model, in which individuals would be represented by atoms with spins -1 (susceptible) and 1 (infected). A Monte Carlo computational simulation was also performed for the Ising-Glauber model in a square network, where each network point represents an individual and the down and up spins represent susceptible and infected individuals.展开更多
The Isin g model is used to investigate (0001), (10^-10), (10^-11) interfaces, etc. , for the hcp lattice. The system studied here have the nearest-neighbor (NN) and next-nearest-neighbor (NNN) pair interactions. The ...The Isin g model is used to investigate (0001), (10^-10), (10^-11) interfaces, etc. , for the hcp lattice. The system studied here have the nearest-neighbor (NN) and next-nearest-neighbor (NNN) pair interactions. The interface energy, interface phase transition at zero temperature and roughening temperature are used to analyse the properties of these interfaces. As a special case of the hcp crystals, we give the equilibrium shape of the He crystal at T=OK.展开更多
Quantum computing is a field with increasing relevance as quantum hardware improves and more applications of quantum computing are discovered. In this paper, we demonstrate the feasibility of modeling Ising Model Hami...Quantum computing is a field with increasing relevance as quantum hardware improves and more applications of quantum computing are discovered. In this paper, we demonstrate the feasibility of modeling Ising Model Hamiltonians on the IBM quantum computer. We developed quantum circuits to simulate these systems more efficiently for both closed and open boundary Ising models, with and without perturbations. We tested these various geometries of systems in both 1-D and 2-D space to mimic two real systems: magnetic materials and biological neural networks (BNNs). Our quantum model is more efficient than classical computers, which can struggle to simulate large, complex systems of particles.展开更多
This study demonstrates that beyond standard model (BSM) cosmic fundamental interactions—weak, strong, and electromagnetic forces—can be unified through a common basis of representation. This unification allows for ...This study demonstrates that beyond standard model (BSM) cosmic fundamental interactions—weak, strong, and electromagnetic forces—can be unified through a common basis of representation. This unification allows for the derivation of the fine structure constant with running points of α(t) ≈ 1/(136.9038) at high energy scales, based on electroweak interactions. Through the application of the Ising model, the running point of the elementary charge e at high energy scales is determined, and Coulomb’s law is actually derived from the Yukawa potential. Theoretically, based on S. Weinberg’s electroweak interaction theory, this study unifies the strong and electromagnetic forces by representing them with rYuka, and further advances the reconstruction of the SU(3)C×SU(1)L×U(1)EMframework on the basis of electroweak interaction concepts. In fact, the cosmic fundamental forces can interchange at the mass gap, defined as the Yukawa turning phase at rYuka ≃1.9404 fm, with the SU(3)Diag structural constant fijk on glueballs calculated, estimating a spectrum mass gap of ∆0 > 0.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.11304123)the Scientific Research Foundation of Jianghan University(Grant No.2010014)
文摘Many phenomena show that in a favorable circumstance an agent still has an updating possibility, and in an unfavor- able circumstance an agent also has a possibility of holding its own state and reselecting its neighbors. To describe this kind of phenomena an Ising model on evolution networks was presented and used for consensus formation and separation of opinion groups in human population. In this model the state-holding probability p and selection-rewiring probability q were introduced. The influence of this mixed dynamics of spin flips and network rewiring on the ordering behavior of the model was investigated, p hinders ordering of opinion networks and q accelerates the dynamical process of networks. Influence of q on the ordering and separating stems from its effect on average path length of networks.
基金Project supported by the National Natural Science Foundation of China (Grant No. 50831006)
文摘An overview of the mathematical structure of the three-dimensional(3D) Ising model is given from the points of view of topology,algebra,and geometry.By analyzing the relationships among transfer matrices of the 3D Ising model,Reidemeister moves in the knot theory,Yang-Baxter and tetrahedron equations,the following facts are illustrated for the 3D Ising model.1) The complex quaternion basis constructed for the 3D Ising model naturally represents the rotation in a(3+1)-dimensional space-time as a relativistic quantum statistical mechanics model,which is consistent with the 4-fold integrand of the partition function obtained by taking the time average.2) A unitary transformation with a matrix that is a spin representation in 2 n·l·o-space corresponds to a rotation in 2n·l·o-space,which serves to smooth all the crossings in the transfer matrices and contributes the non-trivial topological part of the partition function of the 3D Ising model.3) A tetrahedron relationship would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model,and its existence is guaranteed by the Jordan algebra and the Jordan-von Neumann-Wigner procedures.4) The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases φx,φy,and φz.The relationship with quantum field and gauge theories and the physical significance of the weight factors are discussed in detail.The conjectured exact solution is compared with numerical results,and the singularities at/near infinite temperature are inspected.The analyticity in β=1/(kBT) of both the hard-core and the Ising models has been proved only for β〉0,not for β=0.Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.
文摘The review paper by Zhang Zhi-Dong (Zhang Z D 2013 Chin. Phys. B 22 030513, arXiv:1305.2956) contains many errors and is based on several earlier works that are equally wrong.
基金Project supported partially by Guangdong Natural Science Foundation (GDNSF) of China (Grant No 07300793)One of authors(Loan Mushtaq) was partially supported by the Guangdong Ministry of Education,China
文摘We use a new updated algorithm scheme to investigate the critical behaviour of the two-dimensional ferromagnetic Ising model on a triangular lattice with the nearest neighbour interactions. The transition is examined by generating accurate data for lattices with L=8, 10, 12, 15, 20, 25, 30, 40 and 50. The updated spin algorithm we employ has the advantages of both a Metropolis algorithm and a single-update method. Our study indicates that the transition is continuous at Тc=3.6403(2). A convincing finite-size scaling analysis of the model yields ν=0.9995(21), β/ν=0.12400(17), γ/v=1.75223(22), γ^1/ν=1.7555(22), α/ν=0.00077(420) (scaling) and α/ν=0.0010(42) (hyperscaling). The present scheme yields more accurate estimates for all the critical exponents than the Monte Carlo method, and our estimates are shown to be in excellent agreement with their predicted values.
文摘We use the quantum renormalization-group(QRG) method to study the entanglement and quantum phase transition(QPT) in the one-dimensional spin-1/2 Heisenberg-Ising model [Lieb E,Schultz T and Mattis D 1961 Ann.Phys.(N.Y.) 16 407].We find the quantum phase boundary of this model by investigating the evolution of concurrence in terms of QRG iterations.We also investigate the scaling behavior of the system close to the quantum critical point,which shows that the minimum value of the first derivative of concurrence and the position of the minimum scale with an exponent of the system size.Also,the first derivative of concurrence between two blocks diverges at the quantum critical point,which is directly associated with the divergence of the correlation length.
基金Project supported by the Scientific and Technological Research Council of Turkey (TBTAK) (Grant No. 107T533)the Erciyes University Research Funds (Grant Nos. FBA-06-01 and FBD-08-593)
文摘Recently, Shiet al. [2008 Phys. Left. A 372 5922] have studied the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field and presented the dynamic phase diagrams by using an effective-field theory (EFT) and a mean-field theory (MFT). The MFT results are in conflict with those of the earlier work of Tome and de Oliveira, [1990 Phys. Rev. A 41 4251]. We calculate the dynamic phase diagrams and find that our results are similar to those of the earlier work of Tome and de Oliveira; hence the dynamic phase diagrams calculated by Shiet al. are incomplete within both theories, except the low values of frequencies for the MFT calculation. We also investigate the influence of external field frequency (w) and static external field amplitude (h0) for both MFT and EFT calculations. We find that the behaviour of the system strongly depends on the values of w and h0.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51271070)the Natural Science Foundation of Hebei Province(Grant No.E2014202008,E2016202406)+2 种基金the National Science and Technology Program-International Collaborative Research Project(Grant No.2010DFA51920)the Key Program from the science and Technology Research for Colleges and Universities in Hebei Province(Grant No.ZH2011202)the Key Lab.for New Type of Functional Materials in Hebei Province
文摘In order to calculate 3-dimensional Ising model,we develop a method to build a much smaller transfer matrix containing the largest eigenvalue from the original 2^(N_1N_2) × 2^(N_1N_2) matrix V_1. Firstly,the transfer matrix V_1 is written as the linear combination of several basic vectors. Secondly,we divide the basic vectors into several subgroups. The multiplication of a basic vector and V_1 can be written as the linear combination of basic vectors from the same subgroup. Finally,we use a new transfer matrix V_2 to describe the relationship between basic vectors of the same subgroup.V_2 is much smaller than the original transfer matrix and contains the largest eigenvalue of V_1.We use this method to calculate the specific heat per atom Cpaand the magnetic momentum per atom Mpa. The results show that there exists a pair of temperature and magnetic field intensity where the specific heat gets to its maximum value. When N_1N_2 increases,the maximum value of specific heat becomes larger.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 50832002,10674061 and 10874075)the National Key Projects for Basic Research of China (Grant Nos 2006CB921802 and 2009CB623303)
文摘An elastic Ising model for a one-dimensional diatomic spin chain is proposed to explain the ferroelectricity induced by the collinear magnetic order with a low-excited energy state. A statistical theory based on this model is developed to calculate the electrical and magnetic properties of Ca3CoMnO6, a typical quasi-one-dimensional diatomic spin chain system. The calculated ferroelectric polarization and dielectric susceptibility show a good agreement with recently reported data on Ca3Co2-xMnxO6 (x ≈ 0.96) (Phys. Rev. Lett. 100 047601 (2008)), although the predicted magnetic susceptibility does not coincide well with experiment. We also address the rationality and deficiency of this model by including a first-order correction which improves the consistency between the model and experiment.
基金supported by the National Natural Science Foundation of China(Grant No 10571091)the National Science and Technology Supporting Program of China(Grant No 2006BAD11A07)
文摘We use the Monte Carlo method to study an antiferromagnetical Ising spin system on a centred honeycomb lattice, which is composed of two kinds of 1/2 spin particles A and B. There exist two different bond energies JA-A and JA--B in this lattice. Our study is focused on how the ratio of JA-B to JA--A influences the critical behaviour of this system by analysing the physical quantities, such as the energy, the order parameter, the specific heat, susceptibility, etc each as a function of temperature for a given ratio of JA-B to JA-A. Using these results together with the finite-size scaling method, we obtain a phase diagram for the ratio JA-B / JA--A. This work is helpful for studying the phase transition problem of crystals composed of compounds.
基金Project supported by National Natural Science Foundation of China (Grant No 10474057)
文摘This paper proposes a scheme based on the Potts and Ising models for simulating polarization switching of polycrystalline ferroelectrics using the Monte Carlo method. The polycrystalline texture with different average grain size is produced from the Potts model. Then Ising model is implemented in the polycrystalline texture to produce the domain pattern and hysteresis loop. The domain patterns and hysteresis loops have been obtained for polycrystalline texture with different average grain size. From the results of domain pattern evolution process under an applied electric field using this scheme, an extended domain, which covers more than one grain with polarization aligned roughly in the same direction, has been observed during the polarization reversal. This scheme can well reproduce the basic properties of polycrystalline ferroelectrics and is a valuable tool for exploring the physical properties of polycrystalline ferroelectrics.
文摘We present a study of the dynamic behavior of a two-sublattice spin-5/2 Ising model with bilinear and crystal-field interactions in the presence of a time-dependent oscillating external magnetic field on alternating layers of a hexagonal lattice by using the Glauber-type stochastic dynamics.The lattice is formed by alternate layers of spins σ=5/2 and S=5/2.We employ the Glauber transition rates to construct the mean-field dynamic equations.First,we investigate the time variations of the average sublattice magnetizations to find the phases in the system and then the thermal behavior of the dynamic sublattice magnetizations to characterize the nature(first-or second-order) of the phase transitions and to obtain the dynamic phase transition(DPT) points.We also study the thermal behavior of the dynamic total magnetization to find the dynamic compensation temperature and to determine the type of the dynamic compensation behavior.We present the dynamic phase diagrams,including the dynamic compensation temperatures,in nine different planes.The phase diagrams contain seven different fundamental phases,thirteen different mixed phases,in which the binary and ternary combination of fundamental phases and the compensation temperature or the L-type behavior strongly depend on the interaction parameters.
基金supported by the National Natural Science Foundation of China (Grant Nos.10774064 and 30860076)Xinjiang High-Tech Development Foundation (Grant No.200916126)the Key Natural Science Foundation of Xinjiang Science-Technology Department (Grant Nos.200821104 and 200821184)
文摘It was first proposed that an extended chain Ising (ECI) model contains the Ising chain model, single spin doublewell potentials and a pure phonon heat bath of a specific energy exchange with the spins. The extension method is easy to apply to high dimensional cases. Then the single spin-flip probability (rate) of the ECI model is deduced based on the Boltzmann principle and general statistical principles of independent events and the model is simplified to an extended chain Glauber-Ising (ECGI) model. Moreover, the relaxation dynamics of the ECGI model were simulated by the Monte Carlo method and a comparison with the predictions of the special chain Clauber-Ising (SCGI) model was presented. It was found that the results of the two models are consistent with each other when the Ising chain length is large enough and temperature is relative low, which is the most valuable case of the model applications. These show that the ECI model will provide a firm physical base for the widely used single spin-flip rate proposed by Glauber and a possible route to obtain the single spin-flip rate of other form and even the multi-spin-flip rate.
文摘The Blume-Capel model in the presence of external magnetic field H has been simulated using a cellular automaton algorithm improved from the Creutz cellular automaton in three-dimension lattice. The field critical exponent 5 is estimated using the power law relations and the finite size scaling functions for the magnetization and the susceptibility in the range -0.1≤ h = H/J ≤0. The estimated value of the field critical exponent 5 is in good agreement with the universal value (δ = 5) in three dimensions. The simulations are carried out on a simple cubic lattice under periodic boundary conditions.
文摘In this paper it is shown that the thermodynamic limit of the partition function of the statistical models under consideration on a one-dimensional lattice with an arbitrary finite number of interacting neighbors is expressed in terms of the principal eigenvalue of a matrix of finite size. The high sparseness of these matrices for any number of interactions makes it possible to perform an effective numerical analysis of the macro characteristics of these models.
文摘Further to the investigation of the critical properties of the Potts model with <em>q</em> = 3 and 8 states in one dimension (1D) on directed small-world networks reported by Aquino and Lima, which presents, in fact, a second-order phase transition with a new set of critical exponents, in addition to what was reported in Sumour and Lima in studying Ising model on non-local directed small-world for several values of probability 0 < <em>P</em> < 1. In this paper the behavior of two models discussed previously, will be re-examined to study differences between their behavior on directed small-world networks for networks of different values of probability <em>P</em> = 0.1, 0.2, 0.3, 0.4 and 0.5 with different lattice sizes <em>L</em> = 10, 20, 30, 40, and 50 to compare between the important physical variables between Ising and Potts models on the directed small-world networks. We found in our paper that is a phase transitions in both Ising and Potts models depending essentially on the probability <em>P</em>.
文摘The phase diagrams and the other crtical properties of S-3/2 Ising model in random longitudinal and transverse fields(RLIM) are dicussed with the approximate scheme combined by mean-field renormalization group theory(MFRG) and the discretized path-integral representation(DPIR).
文摘From the mathematical point of view, the modeling of epidemics—in other words, the spread of infectious diseases transmitted from individual to individual—is very similar to the modeling of the magnetic systems studied by statistical physics. In this work, we use this analogy between mathematical epidemiology and statistical physics to study the classical mathematical model of epidemiology SI (Susceptible-Infected) approached through the Ising-Glauber model, in which individuals would be represented by atoms with spins -1 (susceptible) and 1 (infected). A Monte Carlo computational simulation was also performed for the Ising-Glauber model in a square network, where each network point represents an individual and the down and up spins represent susceptible and infected individuals.
文摘The Isin g model is used to investigate (0001), (10^-10), (10^-11) interfaces, etc. , for the hcp lattice. The system studied here have the nearest-neighbor (NN) and next-nearest-neighbor (NNN) pair interactions. The interface energy, interface phase transition at zero temperature and roughening temperature are used to analyse the properties of these interfaces. As a special case of the hcp crystals, we give the equilibrium shape of the He crystal at T=OK.
文摘Quantum computing is a field with increasing relevance as quantum hardware improves and more applications of quantum computing are discovered. In this paper, we demonstrate the feasibility of modeling Ising Model Hamiltonians on the IBM quantum computer. We developed quantum circuits to simulate these systems more efficiently for both closed and open boundary Ising models, with and without perturbations. We tested these various geometries of systems in both 1-D and 2-D space to mimic two real systems: magnetic materials and biological neural networks (BNNs). Our quantum model is more efficient than classical computers, which can struggle to simulate large, complex systems of particles.
文摘This study demonstrates that beyond standard model (BSM) cosmic fundamental interactions—weak, strong, and electromagnetic forces—can be unified through a common basis of representation. This unification allows for the derivation of the fine structure constant with running points of α(t) ≈ 1/(136.9038) at high energy scales, based on electroweak interactions. Through the application of the Ising model, the running point of the elementary charge e at high energy scales is determined, and Coulomb’s law is actually derived from the Yukawa potential. Theoretically, based on S. Weinberg’s electroweak interaction theory, this study unifies the strong and electromagnetic forces by representing them with rYuka, and further advances the reconstruction of the SU(3)C×SU(1)L×U(1)EMframework on the basis of electroweak interaction concepts. In fact, the cosmic fundamental forces can interchange at the mass gap, defined as the Yukawa turning phase at rYuka ≃1.9404 fm, with the SU(3)Diag structural constant fijk on glueballs calculated, estimating a spectrum mass gap of ∆0 > 0.