首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 13 毫秒
1.
In this paper we introduce a new Monte Carlo procedure based on the Markov property. This procedure is particularly well suited to massively parallel computation. We illustrate the method on the critical phenomena of the well known one-dimensional Ising model. In the course of this work we found that the autocorrelation time for the Metropolis Monte Carlo algorithm is closely given by the square of the correlation length. We find speedup factors of the order of 1 million for the method as implemented on the CM2 relative to a serial machine. Our procedure gives error estimates which are quite consistent with the observed deviations from the analytically known exact results.  相似文献   

2.
The mixed spin-1/2 and spin-1 Ising chain with both longitudinal and transverse single-ion anisotropies is solved exactly by means of a mapping to the spin-1/2 Ising chains with alternating transverse fields and the Jordan-Wigner transformation. Within this scheme, the thermodynamic quantities of this model are rigorously determined by a recursion formula derived for the partition function based on the reduced spin-1/2 transverse Ising model. The corresponding thermodynamic properties are calculated and discussed.  相似文献   

3.
The magnetization of mixed spin-1/2, 1 Ising diamond chain with antiferromagnetic spin exchange interaction has been investigated by using the Monte Carlo simulation. Magnetization behaviors with different antiferromagnetic interactions are simulated under magnetic field at low temperature. The behaviors are studied in two different cases. In one case the absolute value of the interaction between two spin-1/2 atoms is the largest and in the other case it is not the largest. Multiple plateaus of magnetization M, such as plateaus at M=±1, M=0 and M=±0.566 are obtained. The plateaus at M=±0.566 correspond to two metastable states. The M=0 plateau may disappear at some specific values of spin exchange interactions. The spin configurations were checked when the plateaus appear. The dimmer and trimmer states are found.  相似文献   

4.
Ising models in high dimensions are used to compare high-temperature series expansions with Monte Carlo simulations. Simulations of the magnetization on four-, six-, and seven-dimensional hypercubic lattices give consistent values of the critical temperature from both equilibrium and nonequilibrium data ford=6 and 7. We tabulate 15 terms of series expansions for the susceptibility for generald and giveJ/k B T c =0.092295 (3) and 0.077706 (2) ford=6 and 7. In contrast to five dimensions, where earlier series found nonanalytic scaling corrections, for d=6 and 7 the leading scaling correction may be analytic inT-T c . In most cases these expansions gave more accurate results than these simulations.  相似文献   

5.
The magnetization of a cylindrical Ising nano-tube is investigated by the use of the effective field theory with correlations. The effects of the crystal field couplings at the surface shell to the order parameters, susceptibility, internal energy, specific heat and free energy are investigated. Some characteristic phenomena are examined in the thermal variations, depending on crystal field term. Moreover, tricritical and critical points are found on the (D/J,kT/J) plane, where D/J and kT/J are reduced crystal-field and temperature, respectively.  相似文献   

6.
The phase diagrams in the mixed spin-3/2 and spin-2 Ising system with two alternative layers on a honeycomb lattice are investigated and discussed by the use of the effective-field theory with correlations. The interaction of the nearest-neighbour spins of each layer is taken to be positive (ferromagnetic interaction) and the interaction of the adjacent spins of the nearest-neighbour layers is considered to be either positive or negative (ferromagnetic or anti-ferromagnetic interaction). The temperature dependence of the layer magnetizations of the system is examined to characterize the nature (continuous or discontinuous) of the phase transitions and obtain the phase transition temperatures. The system exhibits both second-and first-order phase transitions besides triple point (T P ), critical end point (E), multicritical point (A), isolated critical point (C) and reentrant behaviour depending on the interaction parameters. We have also studied the temperature dependence of the total magnetization to find the compensation points, as well as to determine the type of behaviour, and N-type behaviour in N′eel classification nomenclature existing in the system. The phase diagrams are constructed in eight different planes and it is found that the system also presents the compensation phenomena depending on the sign of the bilinear exchange interactions.  相似文献   

7.
In this expanded version of an earlier letter, we consider many computational details that were omitted for want of space. Ford = 2 Ising spins with up to 13 different short-range interactions, we construct the critical surface in the vicinity of (Onsager's) nearest-neighbor (nn) critical point by using the body of the available information on the solvable nn case. We then see if the Monte Carlo renormalization group flows generated from this point do indeed lie on this surface and quantify the errors if they do not.  相似文献   

8.
The alternating-bond mixed spin-1/2 and spin-1 Ising chain with both longitudinal and transverse single-ion anisotropies are solved exactly by means of a mapping of the spin-1/2 transverse Ising chain and the Jordan-Wigner transformation. The ground state quantities are strongly dependent on the model Hamiltonian parameters J1, J2, Dx and Dz. We obtain the quasi-particles' spectra Λk, the dimerization gap Δd, the minimal energy Δ0 for exciting a fermion quasi-particle, the minimal energy gap Δh for exciting a hole and the ground state energy Eg. The phase diagram of the ground state is also given. The results show that the alternating bond just quantitatively changes the ground state properties; no matter the nearest-neighbor exchange interactions J1 and J2 are equal or not, when Dz≥0 for any finite value of Dx, there is no quantum critical point and the ground state is always in a spin ordered phase.  相似文献   

9.
The dynamical responses of Ising metamagnet (layered antiferromagnet) in the presence of a sinusoidally oscillating magnetic field are studied by Monte Carlo simulation. The time average staggered magnetisation plays the role of dynamic order parameter. A dynamical phase transition was observed and a phase diagram was plotted in the plane formed by field amplitude and temperature. The dynamical phase boundary is observed to shrink inward as the relative antiferromagnetic strength decreases. The results are compared with that obtained from pure ferromagnetic system. The shape of dynamic phase boundary observed to be qualitatively similar to that obtained from previous meanfield calculations.  相似文献   

10.
An algorithm for the simulation of the 3-dimensional random field Ising model with a binary distribution of the random fields is presented. It uses multi-spin coding and simulates 64 physically different systems simultaneously. On one processor of a Cray YMP it reaches a speed of 184 million spin updates per second. For smaller field strength we present a version of the algorithm that can perform 242 million spin updates per second on the same machine.  相似文献   

11.
Extensive Monte Carlo simulations have been performed to analyze the dynamical behavior of the three-dimensional Ising model with local dynamics. We have studied the equilibrium correlation functions and the power spectral densities of odd and even observables. The exponential relaxation times have been calculated in the asymptotic one-exponential time region. We find that the critical exponentz=2.09 ±0.02 characterizes the algebraic divergence with lattice size for all observables. The influence of scaling corrections has been analyzed. We have determined integrated relaxation times as well. Their dynamical exponentz int agrees withz for correlations of the magnetization and its absolute value, but it is different for energy correlations. We have applied a scaling method to analyze the behavior of the correlation functions. This method verifies excellent scaling behavior and yields a dynamical exponentz scal which perfectly agrees withz.  相似文献   

12.
Monte Carlo simulations are performed for pure and site-diluted Ising ferro- and ferrimagnets on a simple cubic lattice with up to 403 sites and with impurity concentrationx. For the diluted ferromagnet (x=0.2) the exponent= 0.392±0.03 is definitely larger than the pure model value of=0.304±0.03. In contrast, for ferrimagnetic systems (x=0, 0.1, 0.2) the values appear to be independent ofx and within the error limits consistent with the value for the pure ferromagnet, possibly because the width of the asymptotic random critical regime (or of the crossover regime) is even smaller than in the case of ferromagnets.  相似文献   

13.
基于伊辛模型的单自旋反转蒙特卡洛算法,考虑了粒子间的最近邻以及次近邻相互作用,研究了无序 合金的磁化强度和磁熵变。首先,强调了粒子间的次近邻相关作用对体系的磁性和热力学性质的影响,明确了次近邻相互作用系数,证实了低温合金阻挫的存在;其次,研究了在相变温度处(不同磁场下)磁化强度随外加磁场(温度)的变化情况以及磁性粒子对磁化强度的贡献,发现反铁磁性粒子Mn在低温区对 合金的相变起了主要作用,而高温区体系的相变是由铁磁性粒子Fe贡献的;最后,分析了体系在相变温度处磁熵变数值随外加磁场的变化情况以及磁熵变在不同的外磁场下随温度的变化情况,当外加磁场h=0.14时,Mn粒子在冻结温度处的平均磁化强度为零,体系处于最无序的状态,对应的磁熵变 达到了正向最大值,极值的位置对应于体系的相变温度。  相似文献   

14.
Monte Carlo simulations are performed for pure and site-diluted Ising antiferromagnets on a simple cubic lattice with up to 403 sites and impurity concentrationx=0, 0.2. and 0.5. Forx=0.5 a cusp emerges for the temperature dependence of the uniform susceptibility at the critical temperature which is contrasted with the smooth behavior for a pure antiferromagnet, in agreement with the theoretical prediction of Fishman and Aharony.  相似文献   

15.
基于伊辛模型的单自旋反转蒙特卡洛算法,考虑了粒子间的最近邻以及次近邻相互作用,研究了无序Fe0.5Mn0.1Al0.4合金的磁化强度和磁熵变.首先,强调了粒子间的次近邻相关作用对体系的磁性和热力学性质的影响,明确了次近邻相互作用系数,证实了低温合金阻挫的存在;其次,研究了在相变温度处(不同磁场下)磁化强度随外加磁场(温度)的变化情况以及磁性粒子对磁化强度的贡献,发现反铁磁性粒子Mn在低温区对Fe0.5Mn0.1Al0.4合金的相变起了主要作用,而高温区体系的相变是由铁磁性粒子Fe贡献的;最后,分析了体系在相变温度处磁熵变数值随外加磁场的变化情况以及磁熵变在不同的外磁场下随温度的变化情况,当外加磁场H=0.14(a.u.)时,Mn粒子在冻结温度处的平均磁化强度为零,体系处于最无序的状态,对应的磁熵变ΔS(0.1,0.14)达到了正向最大值,极值的位置对应于体系的相变温度.  相似文献   

16.
王志军  董丽芳  尚勇 《物理学报》2005,54(2):880-885
采用蒙特卡罗方法,对源料气体为CH4/H2混合气的电子助进化学气相沉积(EACVD)中 的氢原子(H)、碳原子(C)以及CH基团的发射过程进行了模拟.研究了CH4浓度、反应室气压 和衬底偏压等工艺参数对发射光谱及成膜的影响.研究发现,CH基团可能是有利于金刚石薄 膜生长的活性基团,而碳原子不是;偏压的升高可提高电子平均温度及衬底表面附近氢原子 的相对浓度;通过氢原子谱线可测定电子平均温度并找到最佳成膜实验条件.该结果对EACVD 生长金刚石薄膜过程中实时监测电子平均温度,有效控制工艺条件,生长出高质量的金刚石 薄膜具有重要的意义. 关键词: 蒙特卡罗模拟 金刚石薄膜 发射光谱  相似文献   

17.
Ahmed Zaim  Mohamed Kerouad 《Physica A》2010,389(17):3435-3703
A Monte Carlo simulation has been used to study the magnetic properties and the critical behaviors of a single spherical nanoparticle, consisting of a ferromagnetic core of spins surrounded by a ferromagnetic shell of S=±1, 0 or , spins with antiferromagnetic interface coupling, located on a simple cubic lattice. A number of characteristic phenomena has been found. In particular, the effects of the shell coupling and the interface coupling on both the critical and compensation temperatures are investigated. We have found that, for appropriate values of the system parameters, two compensation temperatures may occur in the present system.  相似文献   

18.
We apply a new entropic scheme to study the critical behavior of the square-lattice Ising model with nearest- and next-nearest-neighbor antiferromagnetic interactions. Estimates of the present scheme are compared with those of the Metropolis algorithm. We consider interactions in the range where superantiferromagnetic (SAF) order appears at low temperatures. A recent prediction of a first-order transition along a certain range (0.5–1.2) of the interaction ratio (R=Jnnn/Jnn) is examined by generating accurate data for large lattices at a particular value of the ratio (R=1). Our study does not support a first-order transition and a convincing finite-size scaling analysis of the model is presented, yielding accurate estimates for all critical exponents for R=1. The magnetic exponents are found to obey “weak universality” in accordance with a previous conjecture.  相似文献   

19.
We present a new algorithm for Monte Carlo simulation of the Ising model. The usual serial architecture of a computer is exploited in a novel way, enabling parallel but independent calculations to be carried out on as many spins as there are bits in a computer word in each fundamental move. The algorithm enjoys a substantial increase in execution speed over more usual multispin coding algorithms. By its very nature, the algorithm constitutes a design for a special-purpose processor.  相似文献   

20.
Xiao-Chen Na 《中国物理 B》2022,31(8):87502-087502
Using an effective field theory with correlations, magnetic properties of an octahedral chain described by a mixed spin Ising model are investigated. Unique phenomena such as multiple hysteresis loops, saturation magnetization, and reverse flip of the magnetization plateaus occur when certain typical parameters are applied. These results may be helpful to further investigate the magnetic properties of one-dimensional systems and could potentially be utilized in the design of spin devices.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号