Creep Motion in Two-Dimensional Fully Frustrated Coulomb Gas Model

The creep motion in a two-dimensional fully frustrated square lattice Coulomb gas model with disorders is studied by using the Monte Carlo technique. The dependence of charge current density J on electric field E is investigated at low temperature T and at low E. The results show that the creep obeys the Arrhenius law J - C（T） exp[-U（E）/T]. The prefactor C（T） increases with the temperature in a power law relation with an exponent about 3.0. The energy barrier U （ E） increases logarithmically with Ec,/ E as U （ E） - Uo ln（ Ec/ E） with Ec being the critical field at zero temperature.     

Short-Time Dynamical Behaviour of Depinning Transition in the Uniformly Frustrated Two-Dimensional XY Model

A short-time dynamic scaling approach is extended to study the depinning transition of the two-dimensional frustrated XY model driven by external currents. We investigate the scaling behaviour of depinning transition in the XY model with three different flux densities f = 1/2, 1/25, 1/30. The short-time scaling behaviour in the depinning transition of the two-dimensional XY model is clearly shown up. Besides the critical current, the exponent θ is obtained.     

On the universality class of the 3d Ising model with long-range-correlated disorder

We analyze a controversial topic about the universality class of the three-dimensional Ising model with long-range-correlated disorder. Whereas both theoretical and numerical studies agree on the validity of the extended Harris criterion [A. Weinrib, B.I. Halperin, Phys. Rev. B 27 (1983) 413] and indicate the existence of a new universality class, numerical values of the critical exponents found so far differ considerably. To resolve this discrepancy we perform extensive Monte Carlo simulations of a 3d Ising model with non-magnetic impurities being arranged in a form of lines along randomly chosen axes of a lattice. The Swendsen-Wang algorithm is used alongside with a histogram reweighting technique and finite-size scaling analysis to evaluate the values of critical exponents governing magnetic phase transition. Our estimates for these exponents differ from both previous numerical simulations and are in favor of a non-trivial dependency of the critical exponents on the peculiarities of long-range correlation decay.     

Majority-vote on directed Erd?s-Rényi random graphs

Through Monte Carlo Simulation, the well-known majority-vote model has been studied with noise on directed random graphs. In order to characterize completely the observed order-disorder phase transition, the critical noise parameter q_c, as well as the critical exponents β/ν, γ/ν and 1/ν have been calculated as a function of the connectivity z of the random graph.     

Compensation temperatures and magnetic susceptibility of a mixed ferro-ferrimagnetic ternary alloy

We study the compensation behaviour and magnetic susceptibility of a mixed ferro-ferrimagnetic ternary alloy of the type ABpC1−p on a cubic lattice by the use of Monte Carlo simulations. We consider that the A and X (X=B or C) ions are alternately connected and have different Ising spins SA=3/2, SB=2, and SC=5/2, respectively. In particular, we find that the magnetic susceptibility may exhibit an outstanding behaviour in a vicinity of the critical and compensation temperatures. The relation of the studied model to ternary Prussian blue analogs is also discussed.     

Monte Carlo analysis of critical properties of the two-dimensional randomly site-diluted Ising model via Wang-Landau algorithm

The influence of random site dilution on the critical properties of the two-dimensional Ising model on a square lattice was explored by Monte Carlo simulations with the Wang-Landau sampling. The lattice linear size was L=20-120 and the concentration of diluted sites q=0.1,0.2,0.3. Its pure version displays a second-order phase transition with a vanishing specific heat critical exponent α, thus, the Harris criterion is inconclusive, in that disorder is a relevant or irrelevant perturbation for the critical behaviour of the pure system. The main effort was focused on the specific heat and magnetic susceptibility. We have also looked at the probability distribution of susceptibility, pseudocritical temperatures and specific heat for assessing self-averaging. The study was carried out in appropriate restricted but dominant energy subspaces. By applying the finite-size scaling analysis, the correlation length exponent ν was found to be greater than one, whereas the ratio of the critical exponents (α/ν) is negative and (γ/ν) retains its pure Ising model value supporting weak universality.     

A comparative study of the dynamic critical behavior of the four-state Potts like models

We investigate the short-time critical dynamics of the Baxter-Wu (BW) and n=3 Turban (3TU) models to estimate their global persistence exponent θ_g. We conclude that this new dynamical exponent can be useful in detecting differences between the critical behavior of these models which are very difficult to obtain in usual simulations. In addition, we estimate again the dynamical exponents of the four-state Potts (FSP) model in order to compare them with results previously obtained for the BW and 3TU models and to decide between two sets of estimates presented in the current literature. We also revisit the short-time dynamics of the 3TU model in order to check if, as already found for the FSP model, the anomalous dimension of the initial magnetization x₀ could be equal to zero.     

Dynamic compensation temperature in the kinetic spin-1 Ising model in an oscillating external magnetic field on alternate layers of a hexagonal lattice

The dynamic behavior of a two-sublattice spin-1 Ising model with a crystal-field interaction (D) in the presence of a time-varying magnetic field on a hexagonal lattice is studied by using the Glauber-type stochastic dynamics. The lattice is formed by alternate layers of spins σ=1 and S=1. For this spin arrangement, any spin at one lattice site has two nearest-neighbor spins on the same sublattice, and four on the other sublattice. The intersublattice interaction is antiferromagnetic. We employ the Glauber transition rates to construct the mean-field dynamical equations. Firstly, we study time variations of the average magnetizations in order to find the phases in the system, and the temperature dependence of the average magnetizations in a period, which is also called the dynamic magnetizations, to obtain the dynamic phase transition (DPT) points as well as to characterize the nature (continuous and discontinuous) of transitions. Then, the behavior of the total dynamic magnetization as a function of the temperature is investigated to find the types of the compensation behavior. Dynamic phase diagrams are calculated for both DPT points and dynamic compensation effect. Phase diagrams contain the paramagnetic (p) and antiferromagnetic (af) phases, the p+af and nm+p mixed phases, nm is the non-magnetic phase, and the compensation temperature or the L-type behavior that strongly depend on the interaction parameters. For D<2.835 and H₀>3.8275, H₀ is the magnetic field amplitude, the compensation effect does not appear in the system.     

Dynamic critical phenomena in two-dimensional fully frustrated Coulomb gas model with disorder

The dynamic critical phenomena near depinning transition in two-dimensional fully frustrated square lattice Coulomb gas model with disorders was studied using Monte Carlo technique. The ground state of the model system with disorder σ=0.3 is a disordered state. The dependence of charge current density J on electric field E was investigated at low temperatures. The nonlinear J-E behavior near critical depinning field can be described by a scaling function proposed for three-dimensional flux line system [M.B. Luo, X. Hu, Phys. Rev. Lett. 98 (2007) 267002]. We evaluated critical exponents and found an Arrhenius creep motion for field region Ec/2<E<Ec. The scaling law of the depinning transition is also obtained from the scaling function.     

Nonequilibrium phase transition in the kinetic Ising model on a two-layer square lattice under the presence of an oscillating field

The nonequilibrium or dynamic phase transitions are studied, within a mean-field approach, in the kinetic Ising model on a two-layer square lattice consisting of spin- 1/2 ions in the presence of a time varying (sinusoidal) magnetic field has been studied by using Glauber-type stochastic dynamics. The dynamic equations of motion are obtained in terms of the intralayer coupling constants J₁ and J₂ for the first and second layer, respectively, and interlayer coupling constant J₃ between these two layers. The nature (first- or second-order) of the transitions is characterized by investigating the behavior of the thermal variations of the dynamic order parameters. The dynamic phase transitions are obtained and the dynamic phase diagrams are constructed in the plane of the reduced temperature versus the amplitude of the magnetic field and found fourteen fundamental types of phase diagrams. Phase diagrams exhibit one, two or three dynamic tricritical points for various values of J₂/|J₁| and J₃/|J₁|. Besides the paramagnetic (p), ferromagnetic (f) and compensated (c) phases, there were the f+c,f+sf,c+sf,af+p,m+p,f+m and c+af, where the af, sf and m are the antiferromagnetic, surface ferromagnetic and mixed phases respectively. Coexistence phase regions also exist in the system.     

Kinetics of a mixed spin-1/2 and spin-3/2 Ising ferrimagnetic model

We present a study, within a mean-field approach, of the kinetics of a mixed ferrimagnetic model on a square lattice in which two interpenetrating square sublattices have spins that can take two values, , alternated with spins that can take the four values, . We use the Glauber-type stochastic dynamics to describe the time evolution of the system with a crystal-field interaction in the presence of a time-dependent oscillating external magnetic field. The nature (continuous and discontinuous) of transition is characterized by studying the thermal behaviors of average order parameters in a period. The dynamic phase transition points are obtained and the phase diagrams are presented in the reduced magnetic field amplitude (h) and reduced temperature (T) plane, and in the reduced temperature and interaction parameter planes, namely in the (h, T) and (d, T) planes, d is the reduced crystal-field interaction. The phase diagrams always exhibit a tricritical point in (h, T) plane, but do not exhibit in the (d, T) plane for low values of h. The dynamic multicritical point or dynamic critical end point exist in the (d, T) plane for low values of h. Moreover, phase diagrams contain paramagnetic (p), ferromagnetic (f), ferrimagnetic (i) phases, two coexistence or mixed phase regions, (f+p) and (i+p), that strongly depend on interaction parameters.     

Phase diagrams of a nonequilibrium mixed spin-1/2 and spin-2 Ising ferrimagnetic system under a time-dependent oscillating magnetic field

We present phase diagrams for a nonequilibrium mixed spin-1/2 and spin-2 Ising ferrimagnetic system on a square lattice in the presence of a time dependent oscillating external magnetic field. We employ the Glauber transition rates to construct the mean-field dynamical 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 dynamic phase transition points are obtained and the phase diagrams are presented in two different planes. Phase diagrams contain paramagnetic (p) and ferrimagnetic (i) phases, and one coexistence or mixed phase region, namely the i+p, that strongly depend on interaction parameters. The system exhibits the dynamic tricritical point and the reentrant behaviors.     

Multicritical dynamical phase diagrams of the kinetic Blume–Emery–Griffiths model with repulsive biquadratic coupling in an oscillating field

We study, within a mean-field approach, the stationary states of the kinetic Blume–Emery–Griffiths model with repulsive biquadratic coupling under the presence of a time-varying (sinusoidal) magnetic field. We employ the Glauber-type stochastic dynamics to construct set of dynamic equations of motion. The behavior of the time dependence of the order parameters and the behavior of the average order parameters in a period, which is also called the dynamic order parameters, as functions of the reduced temperature are investigated. The dynamic phase transition points are calculated and phase diagrams are presented in the reduced magnetic field amplitude and reduced temperature plane. The dynamical transition from one regime to the other can be of first- or second order depending on the region in the phase diagram. According to the values of the crystal field interaction or single-ion anisotropy constant and biquadratic exchange constant, we find 20 fundamental types of phase diagrams which exhibit many dynamic critical points, such as tricritical points, zero-temperature critical points, double critical end points, critical end point, triple point and multicritical point. Moreover, besides a disordered and ordered phases, seven coexistence phase regions exist in the system.     

Dynamic phase transition in the kinetic spin-2 Blume–Emery–Griffiths model in an oscillating field

We extend our recent paper [M. Keskin, O. Canko, M. Erta?, J. Exp. Theor. Phys. (Sov. Phys. JETP) 105 (2007) 1190.] to present a study, within a mean-field approach, the stationary states of the kinetic spin-2 Blume–Emery–Griffiths model in the presence of a time-dependent oscillating magnetic field by using the Glauber-type of stochastic dynamics. We found 20 fundamental types of dynamic phase diagrams where exhibit more complex and richer phase diagrams than our recent paper. Especially, the obtained dynamic phase diagrams show the dynamic triple, quadruple and dynamic double critical end points besides dynamic tricritical points that depending on interaction parameters. The phase diagrams also exhibit a disordered (d) and the ferromagnetic-2 (f₂) phases, and the f₂+d, f₂+fq, fq+d, f₂+f₁+fq and f₂+fq+d, where f₁ are fq the ferromagnetic-1 and ferroquadrupolar or simply quadrupolar phases respectively, coexistence phase regions that strongly depend on interaction parameters.     

The Glauber dynamics for a spin-1 metamagnetic Ising system with bilinear and biquadratic interactions

We present a study, within a mean-field approximation, of the dynamics of a spin-1 metamagnetic Ising system with bilinear and biquadratic interactions in the presence of a time-dependent oscillating external magnetic field. First, we employ the Glauber transition rates to construct the set of mean-field dynamic equations. Then, we study the time variation of the average order parameters to find the phases in the system. We also investigate the thermal behavior of dynamic order parameters to characterize the nature (first- or second-order) of the dynamic transitions. The dynamic phase transitions are obtained and the phase diagrams are constructed in two different the planes. The phase diagrams contain a disordered and ordered phases, and four different mixed phases that strongly depend on interaction parameters. Phase diagrams also display one or two dynamic tricritical points, a dynamic double critical end and dynamic quadruple points. A comparison is made with the results of the other metamagnetic Ising systems.     

Monte Carlo test of the Goldstone mode singularity in 3D XY model

Monte Carlo simulations of magnetization and susceptibility in the 3D XY model are performed for system sizes up to L=384 (significantly exceeding the largest size L=160 considered in work published previously), and fields h ≥ 0.0003125 at two different coupling constants β=0.5, and β=0.55 in the ordered phase. We examine the prediction of the standard theory that the longitudinal susceptibility χ_∥ has a Goldstone mode singularity such that χ_∥ ∝h^-1/2 holds when h↦0. Most of our results, however, support another theoretical prediction that the singularity is of a more general form χ_∥ ∝h^ρ-1, where 1/2<ρ<1 is a universal exponent related to the ∼h^ρ variation of the magnetization.     

Percolation on patchwise heterogeneous surfaces under equilibrium conditions

The site percolation problem on square lattices whose sites are grouped in two types of energetically different patches is studied. Several lattices formed by collections of either randomly or orderly localized and no overlapped patches of different sizes are generated. The system is characterized by two parameters, namely, the size of each patch, l, and the energy difference between the two kind of sites, ΔE. Particles are adsorbed at equilibrium on the lattice. The critical coverage is determined by means of Monte Carlo simulations and finite-size scaling analysis. The percolative behavior of the system as a function of the parameters characterizing the heterogeneity of the energetic surface topography is presented and discussed.     

Driven lattice gas with nearest-neighbor exclusion: shear-like drive

We study the lattice gas with nearest-neighbor exclusion on the square lattice and Kawasaki (hopping) dynamics, under the influence of a nonuniform drive, via Monte Carlo simulation. The drive, which favors motion along the +x direction and inhibits motion in the opposite direction, varies linearly with y. (The boundaries along the drive direction are periodic, so that the system is not described by an equilibrium Gibbs distribution.) As in the uniformly driven case [R. Dickman, Phys. Rev. E 64, 16124 (2001)], the onset of sublattice ordering occurs at a lower density than in equilibrium, but here an unexpected feature appears: particles migrate out of the high-drive region. For intermediate system sizes (L ≃100), the accumulation of particles is sufficient for the low-drive region to become ordered at a global density of about 0.3. Above this density we observe a surprising reversal in the density profile, with particles accumulating to the high-drive region, due to jamming. For larger systems (L≥200) particles quickly jam in the high-drive region, as occurs under uniform drive, and the accumulation of particles in the low-field region is severely reduced.     

Phase Transition of a Distance-Dependent Ising Model on the Barabasi-Albert Network

We report our investigation on the behaviour of distance-dependent Ising models, which are located on the BA model network. The interaction strength between two nodes （the spins） is considered to obey an exponential decay dependence on the geometrical distance. The Monte Carlo simulation shows a phase transition from ferromagnetism to paramagnetism, and the critical temperature approaches a constant temperature as the interaction decaying exponent increases.     

Phase diagrams of a nonequilibrium mixed spin-3/2 and spin-2 Ising system in an oscillating magnetic field

The phase diagrams of the nonequilibrium mixed spin-3/2 and spin-2 Ising ferrimagnetic system on square lattice under a time-dependent external magnetic field are presented by using the Glauber-type stochastic dynamics. The model system consists of two interpenetrating sublattices of spins σ=3/2 and S=2, and we take only nearest-neighbor interactions between pairs of spins. The system is in contact with a heat bath at absolute temperature T_abs and the exchange of energy with the heat bath occurs via one-spin flip of the Glauber dynamics. First, we investigate the time variations of average order parameters to find the phases in the system and then the thermal behavior of the dynamic order parameters to obtain the dynamic phase transition (DPT) points as well as to characterize the nature (first- or second-order) phase transitions. The dynamic phase diagrams are presented in two different planes. Phase diagrams contain paramagnetic (p), ferrimagnetic (i₁, i₂, i₃) phases, and three coexistence or mixed phase regions, namely i₁+p, i₂+p and i₃+p mixed phases that strongly depend on interaction parameters.