Phase transitions of two spin-1/2 Baxter–Wu layers coupled with Ising-type interactions

Using a Monte Carlo simulation and the single histogram reweighting technique,we study the critical behaviors and phase transitions of the Baxter-Wu(BW)model on a two-layer triangular lattice with Ising-type interlayer couplings.Via the finite-size analysis,we obtain the transition temperatures and critical exponents at repulsive and attractive interlayer couplings.The data for the repulsive interlayer coupling suggest continuous transitions,and the critical behaviors are the same as those of the 2D BW model,belonging to the four-state Potts universality class.The reduced energy cumulants and the histograms reveal that attractive coupling leads to weak firstorder phase transitions.The pseudocritical exponents with the existence of the interlayer couplings indicate that the first-order transition is very close to the critical point of the 2D standard BW model.     

A new cluster algorithm for the Baxter-Wu model

We propose a new cluster algorithm for the Baxter-Wu model that significantly reduces critical slowing down. We examine the behavior of the created clusters as we vary the temperature and then specify dynamic exponents. Comparison is made with the Metropolis algorithm and with the other existing cluster algorithm.     

Phase transitions and critical phenomena in a three-dimensional site-diluted Potts model

The phase transitions and critical phenomena in the three-dimensional (3D) site-diluted q-state Potts models on a simple cubic lattice are explored. We systematically study the phase transitions of the models for q=3 and q=4 on the basis of Wolff high-effective algorithm by the Monte–Carlo (MC) method. The calculations are carried out for systems with periodic boundary conditions and spin concentrations p=1.00–0.65. It is shown that introducing of weak disorder (p∼0.95) into the system is sufficient to change the first order phase transition into a second order one for the 3D 3-state Potts model, while for the 3D 4-state Potts model, such a phase transformation occurs when introducing strong disorder (p∼0.65). Results for 3D pure 3-state and 4-state Potts models (p=1.00) agree with conclusions of mean field theory. The static critical exponents of the specific heat α, susceptibility γ, magnetization β, and the exponent of the correlation radius ν are calculated for the samples on the basis of finite-size scaling theory.     

Phase diagrams of the mixed spin-1 and spin-1/2 Ising-Heisenberg model on the diamond-like decorated Bethe lattice

The ground-state and finite-temperature behavior of the mixed spin-1 and spin-1/2 Ising-Heisenberg model on the diamond-like decorated Bethe lattice is investigated within the framework of two rigorous methods: the decoration-iteration transformation and exact recursion relations. The model under consideration describes a hybrid classical-quantum system consisting of the Ising and Heisenberg spins, which interact among themselves either through the Ising or XXZ Heisenberg nearest-neighbor interaction. Both sublattice magnetizations of the Ising and Heisenberg spins are exactly calculated with the aim to examine phase diagrams, thermal variations of the total and sublattice magnetizations. The finite-temperature phase diagrams form continuous (second-order) phase transition lines only, which exhibit a small reentrant region if the diamond-like decorated Bethe lattice with a sufficiently high coordination number is considered.     

Dynamic phase transitions in the kinetic mixed spin-1/2 and spin-5/2 Ising model under a time-dependent oscillating magnetic field

We calculate the dynamic phase transition (DPT) temperatures and present the dynamic phase diagrams in the kinetic mixed spin-1/2 and spin-5/2 Ising model under the presence of a time-dependent oscillating external magnetic field. We employ the Glauber transition rates to construct the set of mean-field dynamic equations. The time variation of the average magnetizations and the thermal behavior of the dynamic magnetizations are investigated, extensively. The nature (continuous or discontinuous) of the transitions is characterized by studying the thermal behaviors of the dynamic magnetizations. The DPT points are obtained and the phase diagrams are presented in two different planes. Phase diagrams contain four fundamental phases and three coexistence or mixed phases, which strongly depend on interaction parameters. The phase diagrams are discussed and a comparison is made with the results of the other mixed spin Ising systems.     

Thermal and magnetic properties of a ferrimagnetic nanoparticle with spin-3/2 core and spin-1 shell structure

Intensive Monte Carlo simulations based on standard Metropolis algorithm have been applied to investigate the phase diagrams of a ferrimagnetic cubic nanoparticle (nanocube) with a spin-3/2 core surrounded by a spin-1 shell layer with antiferromagnetic interface coupling. It has been shown that occupation of sites of particle core by spin-3/2 plays an important role on the shape of the phase diagrams. In addition, effects of the crystal field interaction as well as exchange interactions on the thermal and magnetic properties of the system have also been discussed in detail, and some interesting features have been observed for the temperature dependence of total magnetization curves of particle.     

Scaling properties of the Baxter-Wu model

We use a recently developed cluster algorithm for the Baxter-Wu model to study characteristic features of its behavior. Magnetic scaling properties of the model are investigated for second-order phase transitions. We improve significantly the accuracy of Monte Carlo simulation results and we provide new results for the behavior of the Binder cumulants for temperatures below, above and at the critical point.     

Effective-field and Monte Carlo studies of mixed spin-2 and spin-1/2 Ising diamond chain

The magnetic properties of a mixed spin-2 and spin-1/2 ferromagnetic diamond chain are studied by effective-field theory and Monte Carlo(MC) simulation based on the Ising model.The temperature dependences of magnetization,magnetic susceptibility,internal energy,and specific heat are studied,respectively.The exchange interaction dependences of magnetization and the critical temperature are calculated by MC simulation.The changes of magnetization depending on the field increasing and then the field decreasing under steady-static conditions are also given.     

Phase transitions and universality in nonequilibrium steady states of stochastic Ising models

We present results of direct computer simulations and of Monte Carlo renormalization group (MCRG) studies of the nonequilibrium steady states of a spin system with competing dynamics and of the voter model. The MCRG method, previously used only for equilibrium systems, appears to give useful information also for these nonequilibrium systems. The critical exponents are found to be of Ising type for the competing dynamics model at its second-order phase transitions, and of mean-field type for the voter model (consistent with known results for the latter).     

超薄膜磁性材料上XY模型的临界现象

采用经典XY模型，阐明了三角格子上层状超薄膜磁性材料的相变和临界现象.并用Monte-Carlo方法对该模型的内部能量、比热、chirality等热力学量进行了计算.使用finite-size scaling 分析法对临界现象的性质进行了论述.通过上述分析和计算，发现该模型在反强磁性层与强磁性层内的chirality在一定范围内随温度的变化急剧增加，这是一种新的chirality相变；而在另一范围内存在Kosterliz-Thouless相变和通常的chirality相变.这种新的chirality相变的 关键词： 层状XY模型 临界现象 Monte-Carlo模拟 磁性薄膜     

自旋-1/2量子罗盘链的量子相与相变

基于矩阵乘积态表述的无限时间演化块算法,研究了具有x,y,z三个自旋方向的轨道自由度和轨道序竞争的量子罗盘自旋链模型.为了刻画该模型的量子相和相变,计算了基态能量、局域序参量、弦关联序参量、临界指数、冯诺依曼熵、有限纠缠标度和中心荷.结果表明:该量子基态相图由条纹反铁磁相、反铁磁相、单调奇数Haldane相和振荡奇数Haldane相构成.从条纹反铁磁相到反铁磁相,以及从单调奇数Haldane相到振荡奇数Haldane相发生了非连续相变;从振荡奇数Haldane相到条纹反铁磁相,以及从反铁磁相到单调奇数Haldane相发生了连续相变;连续相变线和非连续相变线的交点是多临界点.此外,连续相变点处的临界指数β=1/8和中心荷c=1/2表明连续相变的普适类属于Ising类.由此揭示了该模型量子基态相图的本性,对今后研究更高自旋以及更为复杂轨道序竞争的量子罗盘链模型的量子相与相变具有一定借鉴与参考意义.     

Phase transitions in the Ising model with transverse field

We show the existence of a phase transition in the Ising model with transverse field for dimensionsv 2 provided the transverse term is sufficiently small. This is done by proving long-range order occurs using the reflection positivity of the Hamiltonian and localization of eigenvectors.     

Phase transitions in a driven lattice gas in two planes

We report on a Monte Carlo study of ordering in a nonequilibrium system. The system is a lattice gas that comprises two equal, parallel square lattices with stochastic particle-conserving irreversible dynamics. The particles are driven along a principal direction under the competition of the heat bath and a large, constant external electric field. There is attraction only between particles on nearest-neighbor sites within the same lattice. Particles may jump from one plane to the other; therefore, density fluctuations have an extra mechanism to decay and build up. It helps to obtain the steady-state accurately. Spatial correlations decay with distance according to a power law at high enough temperature, as for the ordinary two-dimensional case. We find two kinds of nonequilibrium phase transitions. The first one has a critical point for half occupation of the lattice, and seems to be related to the anisotropic phase transition reported before for the plane. This transition becomes discontinuous for low enough density. The difference of density between the planes changes discontinuously for any density at a lower temperature. This seems to correspond to a phase transition that does not have a counterpart in equilibrium nor in the two-dimensional nonequilibrium case.     

Quantum phase transitions in matrix product states of one-dimensional spin-1/2 chains

For the matrix product system of a one-dimensional spin-1/2 chain, we present a new model of quantum2 phase transitions and find that in the thermodynamic limit, both sides of the critical point are respectively described by phases ｜Ψa 〉=｜1··· 1 representing all particles spin up and ｜Ψb 〉=｜0··· 0 representing all particles spin down, while the phase transition point is an isolated intermediate-coupling point where√ the two phases coexist equally, which is2 described by the so-called N-qubit maximally entangled GHZ state ｜Ψpt =√2/2（｜1··· 1 ＋｜0··· 0）. At the critical point,2the physical quantities including the entanglement are not discontinuous and the matrix product system has longrange correlation and N-qubit maximal entanglement. We believe that our work is helpful for having a comprehensive understanding of quantum phase transitions in matrix product states of one-dimensional spin chains and of potential directive significance to the preparation and control of one-dimensional spin lattice models with stable coherence and N-qubit maximal entanglement.     

Dynamic phase transitions in the kinetic spin-1 Blume-Capel model on the Bethe lattice

The stationary states of the kinetic spin-1 Blume-Capel (BC) model on the Bethe lattice are analyzed in detail in terms of recursion relations. The model is described using a Glauber-type stochastic dynamics in the presence of a time-dependent oscillating external magnetic field (h) and crystal field (D) interactions. The dynamic order parameter, the hysteresis loop area and the dynamic correlation are calculated. It is found that the magnetization oscillates around nonzero values at low temperatures (T) for the ferromagnetic (F) phase while it only oscillates around zero values at high temperatures for the paramagnetic (P) phase. There are regions of the phase space where the two solutions coexist. The dynamic phase diagrams are obtained on the (kT/J,h/J) and (kT/J,D/J) planes for the coordination number q=4. In addition to second-order and first-order phase transitions, dynamical tricritical points and triple points are also observed.     

Novel scaling behavior of the Ising model on curved surfaces

We demonstrate the nontrivial scaling behavior of Ising models defined on (i) a donut-shaped surface and (ii) a curved surface with a constant negative curvature. By performing Monte Carlo simulations, we find that the former model has two distinct critical temperatures at which both the specific heat C(T) and magnetic susceptibility χ(T) show sharp peaks. The critical exponents associated with the two critical temperatures are evaluated by the finite-size scaling analysis; the result reveals that the values of these exponents vary depending on the temperature range under consideration. In the case of the latter model, it is found that static and dynamic critical exponents deviate from those of the Ising model on a flat plane; this is a direct consequence of the constant negative curvature of the underlying surface.     

Study of cluster fluctuations in the two-dimensional q-state Potts model

The two-dimensional Potts model with 2 to 10 states is studied using a cluster algorithm to calculate fluctuations in cluster size as well as commonly used quantities like equilibrium averages and the histograms for energy and the order parameter. Results provide information about the variation of cluster sizes depending on the temperature and the number of states. They also give evidence for first-order transition when energy and the order parameter related measurables are inconclusive on small size lattices.     

自旋为1的双层平方晶格阻挫模型的相变

采用双时格林函数方法研究了自旋为1的双层平方晶格阻挫模型的相变行为.详细探讨了层间耦合相互作用J_c和单离子各向异性参数D对奈尔态(AF1)和共线态(AF2)之间相转换的影响.结果显示:只要参数J_c和D不同时为零,奈尔态和共线态在J₂=J₁/2(这里J₁和J₂分别描述的是系统自旋间最近邻和次近邻交换作用)时的相变温度相等,两个态共存.在低于相变点的温度范围内,AF1-AF2态之间可以发生相转换,其相变类型为一阶相变.当J₂≠J₁/2时,尽管AF1-AF2态有不同相变温度,但它们也可以共存.如果AF1(AF2)态的相变温度大,在低温,AF1(AF2)态更稳定;在高温,AF2(AF1)态更稳定;在中间温度范围内,AF1-AF2态之间也可以发生一阶相转换.