Critical behavior of the specific heat for pure and site-diluted simple cubic Ising systems

The exponent of the specific heatC is determined for the pure and the site-diluted simple cubic Ising model (concentrationx=0, 0.2, 0.4 of nonmagnetic sites) by a finite-size scaling analysis of the peak value C_max(L) for systems of linear dimensionsL=8, 16, 32, and 64. The C_max values are obtained by the Ferrenberg-Swendsen algorithm, using Monte Carlo data from a fully-vectorized multi-spin coding program. We obtain =0.11 for x=0 and a crossover to a negative value upon dilution, with =–0.029(4) both forx=0.2 andx=0.4.     

Monte Carlo study of the critical behavior of pure and site-diluted Ising ferro-and ferrimagnets

Monte Carlo simulations are performed for pure and site-diluted Ising ferro- and ferrimagnets on a simple cubic lattice with up to 40³ 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.     

Interface motion in a two-dimensional Ising model with a field

We determine by Monte Carlo simulations the width of an interface between the stable phase and the metastable phase in a two-dimensional Ising model with a magnetic field, in the case of nonconversed order parameter (Glauber dynamics). At zero temperature, the width increases ast with–1/3, as predicted by earlier theories. As temperature increases, the value of the effective exponent that we measure decreases toward the value 1/4, which is the value in the absence of magnetic field.     

Critical wetting in the square Ising model with a boundary field

The Ising square lattice with nearest-neighbor exchangeJ>0 and a free surface at which a boundary magnetic fieldH ₁ acts has a second-order wetting transition. We study the surface excess magnetization and the susceptibility ofL×M lattices by Monte Carlo simulation and probe the critical behavior of this wetting transition, applying finite-size scaling methods. For the cases studied, the results are not consistent with the presumably exactly known values of the critical exponents, because the asymptotic critical region has not yet been reached. Implication of our results for critical wetting in three dimensions and for the application of the present model to adsorbed wetting layers at surface steps are briefly discussed.Alexander von Humboldt-Fellow     

Wang-Landau study of the critical behavior of the bimodal 3D random field Ising model

We apply the Wang-Landau method to the study of the critical behavior of the three-dimensional random field Ising model with a bimodal probability distribution. For high values of the random field intensity we find that the energy probability distribution at the transition temperature is double peaked, suggesting that the phase transition is of first order. On the other hand, the transition looks continuous for low values of the field intensity. In spite of the large sample to sample fluctuations observed, the double peak in the probability distribution is always present for high fields.     

Monte Carlo evidence for a cusp in the uniform susceptibility of the site-diluted simple cubic Ising antiferromagnet

Monte Carlo simulations are performed for pure and site-diluted Ising antiferromagnets on a simple cubic lattice with up to 40³ 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.     

Critical behaviour of the ferromagnetic Ising model on a triangular lattice

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 T_c= 3.6403({2}). A convincing finite-size scaling analysis of the model yields υ=0.9995(21), β / υ = 0.12400({17}), γ / υ = 1.75223(22), γ '/υ=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.     

Monte Carlo study of the antiferromagnetical Ising model on a centred honeycomb lattice

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 J_A-A and J_A-B in this lattice. Our study is focused on how the ratio of J_A-B to J_A-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 J_A-B to J_A-A. Using these results together with the finite-size scaling method, we obtain a phase diagram for the ratio J_A-B / J_A-A. This work is helpful for studying the phase transition problem of crystals composed of compounds.     

Fast vectorized algorithm for the Monte Carlo simulation of the random field Ising model

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.     

Critical dynamics of the three-dimensional Ising model: A Monte Carlo study

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.     

Ising model for the alloys according to the effective field theory

FeAl alloys in their disordered structural phase have been investigated through an Ising model where besides exchange interactions between nearest-neighbors Fe atoms, a superexchange interaction mediated by Al atoms is also taken into account. The model has been approximately treated according to the effective field theory. Although the phase diagram, as a function of Al concentration, is similar to the one previously obtained from Bogoliubov variational approach for the free energy, a different behavior for the superexchange interaction is achieved, which can also be physically accepted for this system.     

Computer investigations of the three-dimensional Ising model

We have been studying the three-dimensional Ising model using some finite-size scaling ideas. The simulation is done by a fast microcanonical method. Here we present our results for the critical exponents and.     

Surface Critical Behavior of Two-Dimensional Dilute Ising Models

Ising models with nearest neighbor ferromagnetic random couplings on a square lattice with a (1, 1) surface are studied, using Monte Carlo techniques and a star-triangle transformation method. In particular, the critical exponent of the surface magnetization is found to be close to that of the perfect model, Β₁ = 1/2. The crossover from surface to bulk critical properties is discussed.     

Critical Casimir effects in 2D Ising model with curved defect lines

This work is aimed at studying the influence of critical Casimir effects on energetic properties of curved defect lines in the frame of 2D Ising model. Two types of defect curves were investigated. We start with a simple task of globule formation from four-defect line. It was proved that an exothermic reaction of collapse occurs and the dependence of energy release on temperature was observed. Critical Casimir energy of extensive line of constant curvature was also examined. It was shown that its critical Casimir energy is proportional to curvature that leads to the tendency to radius decreasing under Casimir forces. The results obtained can be applied to proteins folding problem in polarized liquid.     

Monte Carlo study of dynamic phase transition in Ising metamagnet driven by oscillating magnetic field

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.     

Minor magnetization loops in two-dimensional dipolar Ising model

The two-dimensional dipolar Ising model is investigated for the relaxation and dynamics of minor magnetization loops. Monte Carlo simulations show that in a stripe phase an exponential decrease can be found for the magnetization maxima of the loops, M∼exp(−αN_l) where N_l is the number of loops. We discuss the limits of this behavior and its relation to the equilibrium phase diagram of the model.     

Monte Carlo study of the critical behavior and magnetic properties of La2/3Ca1/3MnO3 thin films

Critical exponents offer important information concerning the interaction mechanisms near the paramagnetic to ferromagnetic transition. In this work a Monte Carlo-Metropolis simulation of the critical behavior in La_2/3Ca_1/3MnO₃ thin films is addressed. Canonical ensemble averages for magnetization per site, magnetic susceptibility and specific heat of stoichiometric manganite within a three-dimensional classical Heisenberg model with nearest magnetic neighbor interactions are computed. The La_2/3Ca_1/3MnO₃ thin films were simulated addressing the thickness influence and thermal dependence. In the model, Mn magnetic ions are distributed on a simple cubic lattice according to the perovskite structure of this manganite. Ferromagnetic coupling for the bonds Mn³⁺-Mn³⁺(e_g-e_g′), Mn³⁺-Mn⁴⁺(e_g-d³) and Mn³⁺-Mn⁴⁺(e_g′-d³) were taken into account. On the basis of finite-size scaling theory, our best estimates of critical exponents, linked to the ferromagnetic to paramagnetic transition, for the correlation length, specific heat, magnetization and susceptibility are, respectively: v=0.56±0.01, α=0.16±0.03, β=0.34±0.04γ and γ=1.17±0.05. These theoretical results are consistent with the Rushbrooke equalitiy α+2β+γ=2.     

基于伊辛模型的Fe0.5Mn0.1Al0.4 合金磁化强度和磁熵变蒙特卡洛模拟

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