The transverse spin-1 Ising model with random interactions

The phase diagrams of the transverse spin-1 Ising model with random interactions are investigated using a new technique in the effective field theory that employs a probability distribution within the framework of the single-site cluster theory based on the use of exact Ising spin identities. A model is adopted in which the nearest-neighbor exchange couplings are independent random variables distributed according to the law P(J_ij)=pδ(J_ij−J)+(1−p)δ(J_ij−αJ). General formulae, applicable to lattices with coordination number N, are given. Numerical results are presented for a simple cubic lattice. The possible reentrant phenomenon displayed by the system due to the competitive effects between exchange interactions occurs for the appropriate range of the parameter α.     

Upper and lower bounds to a correlation function for an Ising model with random interactions

We consider the general Ising model with random interactions Jp. We assume that the probability densities of the random interactions are statistically independent and that the averages of the absolute values of the random interactions, |Jp|, are finite. We then show that a correlation function for the regular Ising model with interactions |Jp| and the same quantity with opposite sign are an upper and a lower bound to the corresponding averaged correlation function of the random Ising model under consideration.     

Dynamics of the one-dimensional random transverse Ising model with next-nearest-neighbor interactions

The dynamics of the one-dimensional random transverse Ising model with both nearest-neighbor (NN) and next-nearest-neighbor (NNN) interactions is studied in the high-temperature limit by the method of recurrence relations. Both the time-dependent transverse correlation function and the corresponding spectral density are calculated for two typical disordered states. We find that for the case of bimodal disorder the dynamics of the system undergoes a crossover from a collective-mode behavior to a central-peak one and for the case of Gaussian disorder the dynamics is complex. For both cases, it is found that the central-peak behavior becomes more obvious and the collective-mode behavior becomes weaker as K_i increase, especially when K_i>J_i/2 (J_i and K_i are the exchange couplings of the NN and NNN interactions, respectively). However, the effects are small when the NNN interactions are weak (K_i<J_i/2).     

Renormalization group studies of phase diagrams of Ising and XY-model films of variable thickness

We compare the critical behavior of $\text{spin}\text{-}\text{1}\text{2}$ Ising model (ordinary phase transition in two dimensions) and classical XY-model (topological phase transition in two dimensions) films, with two flat surfaces and nearest neighbor couplings K_S between surface spins and K_B between all others, as a function of film thickness. We carry out a real space Migdal-style renormalization in two stages. In the bulk stage the film is first renormalized towards a double layer, with the renormalized parameters as inter- and intra-layer couplings. Then the double layer is renormalized with those couplings as initial parameters. From the RG-equations for the bulk stage we find a tricritical point, not only for the Ising model (in which case it is well known) but also for the XY-model. It signals the existence of distinct surface and bulk transitions for sufficiently large values of $\text{K}{}_{\mathrm{<mtext>S</mtext>}}\text{K}{}_{\mathrm{<mtext>B</mtext>}}$. For the Ising model the complete program can be carried out and the phase diagram for films of arbitrary thickness is constructed. For the double layer XY-model a sufficiently complex Migdal-style renormalization appears to be unfeasible, presumably due to the possible presence of strings between the layers. Therefore, in an alternative approach, two representations for the partition function of the double layer XY-model are given. The system can be described in terms of its topological excitations, i.e., vortices on each layer and strings, either closed or terminating in vortices, between them. The system is also written as three coupled Coulomb gases. Based on this representation a renormalization group is found, and used, together with a Griffiths inequality for the correlation functions, to obtain information on the phase diagram. If there is a single phase transition, there is one phase where the correlations exhibit power law behavior and another where they fall off exponentially. The transition temperature increases monotonously with interlayer coupling to twice its value in two dimensions, but the nature of the phase transition for any finite inter-layer coupling appears to be different from that at zero coupling. We suggest that this is associated with the behavior of the internal strings. These results should be relevant for the renormalization of the film with isotropic bulk couplings as well as for layered systems of finite thickness, with different inter- and intra-layer couplings.     

On the susceptibility of the one-dimensional Ising chain

The elements of the tensor of the susceptibility and the staggered susceptibility of a one-dimensional I Ising chain, with xx couplings between the spins are investigated in the presence of a magnetic field b parallel with or perpendicular to the x direction. Exact expressions are given for all susceptibilities, apart from the parallel susceptibilities in a transverse field which are evaluated by perturbation calculation. Special attention is paid to the conditions under which a susceptibility can have a minimum for b = 0. Furthermore, a system with weakly coupled Ising chains is considered on the basis of a model hamiltonian with separable interchain interactions.     

Strict monotonicity for critical points in percolation and ferromagnetic models

When is the numerical value of the critical point changed by an enhancement of the process or of the interaction? Ferromagnetic spin models, independent percolation, and the contact process are known to be endowed with monotonicity properties in that certain enhancements are capable of shifting the corresponding phase transition in only an obvious direction, e. g., the addition of ferromagnetic couplings can only increase the transition temperature. The question explored here is whether enhancements do indeed change the value of the critical point. We present a generally applicable approach to this issue. For ferromagnetic Ising spin systems, with pair interactions of finite range ind?2 dimensions, it is shown that the critical temperatureT _c is strictly monotone increasing in each coupling, with the first-order derivatives bounded by positive functions which are continuous on the set of fullyd-dimensional interactions. For independent percolation, with 0<p _c<1, we prove that any “essential enhancement” of the process has an effect on the critical probability, a result with applications to the question of the existence of “entanglements” and to invasion percolation with trapping.     

Generalized transformation for decorated spin models

The paper discusses the transformation of decorated Ising models into an effective undecorated spin model, using the most general Hamiltonian for interacting Ising models including a long range and high order interactions. The inverse of a Vandermonde matrix with equidistant nodes [−s,s] is used to obtain an analytical expression of the transformation. This kind of transformation is very useful to obtain the partition function of decorated systems. The method presented by Fisher is also extended, in order to obtain the correlation functions of the decorated Ising models transforming into an effective undecorated Ising model. We apply this transformation to a particular mixed spin-(1/2, 1) and (1/2, 2) square lattice with only nearest site interaction. This model could be transformed into an effective uniform spin-S square lattice with nearest and next-nearest interaction, furthermore the effective Hamiltonian also includes combinations of three-body and four-body interactions; in particular we considered spin 1 and 2.     

Monte Carlo technique for very large ising models

Rebbi's multispin coding technique is improved and applied to the kinetic Ising model with size 600*600*600. We give the central part of our computer program (for a CDC Cyber 76), which will be helpful also in a simulation of smaller systems, and describe the other tricks necessary to go to large lattices. The magnetizationM atT=1.4*T _c is found to decay asymptotically as exp(-t/2.90) ift is measured in Monte Carlo steps per spin, and M(t = 0) = 1 initially.     

A simple model for metamagnetic percolation

An Ising model with both antiferromagnetic and randomly diluted, infinitely strong ferromagnetic interactions is studied on a Bethe lattice. The results are used to interpret recent experimental data on the phase diagram of Fe_xMg_1-xCl₂.     

Anyons in an exactly solved model and beyond

A spin-1/2 system on a honeycomb lattice is studied. The interactions between nearest neighbors are of XX, YY or ZZ type, depending on the direction of the link; different types of interactions may differ in strength. The model is solved exactly by a reduction to free fermions in a static Z₂ gauge field. A phase diagram in the parameter space is obtained. One of the phases has an energy gap and carries excitations that are Abelian anyons. The other phase is gapless, but acquires a gap in the presence of magnetic field. In the latter case excitations are non-Abelian anyons whose braiding rules coincide with those of conformal blocks for the Ising model. We also consider a general theory of free fermions with a gapped spectrum, which is characterized by a spectral Chern number ν. The Abelian and non-Abelian phases of the original model correspond to ν = 0 and ν = ±1, respectively. The anyonic properties of excitation depend on ν mod 16, whereas ν itself governs edge thermal transport. The paper also provides mathematical background on anyons as well as an elementary theory of Chern number for quasidiagonal matrices.     

Correlations in inhomogeneous Ising models

Using general methods developed in a previous treatment we study correlations in inhomogeneous Ising models on a square lattice. The nearest neighbour couplings can vary both in strength and sign such that the coupling distribution is translationally invariant in diagonal direction. We calculate correlations parallel to the layering in the diagonally layered model with periodv=2, the so-called “general square lattice” model (GS). If the model has a finite critical temperature,T _c>0, we have a spontaneous magnetization belowT _c vanishing atT _c with the Ising exponent β=1/8. AtT _c correlations decay algebraically with critical exponnet η=1/4 and exponentially forT>T _c. In the frustrated case we have oscillatory behaviour superposed on the exponential decay where the wavevector of the oscillations changes at some “disorder temperature”T _D(>T _c) from commensurate to temperature-dependent in commensurate periods. If the critical temperature vanishes,T _c=0 we always have exponential decay at finite temperatures, while atT=T _c=0 we encounter either long-range order or algebraic decay with critical index η=1/2, i.e.T=0 is thus a critical point.     

Semi-infinite Ising model

For the semi-infinite Ising model in two or more dimensions, we prove analyticity properties of the surface free energy and map out the phase diagram in the absence of an external magnetic field. We prove that this phase diagram contains critical lines where the parallel and/or the transverse correlation lengths diverge. The critical exponent,v _⊥, of the transverse correlation length is shown to be equal to the exponentv of the Ising model on an infinite lattice. In a second paper, these results will be used to analyze the wetting transition.     

Calculation of critical properties for the anisotropic two-layer Ising model on the Kagome lattice: Cellular automata approach

The critical point (K_c) of the two-layer Ising model on the Kagome lattice has been calculated with a high precision, using the probabilistic cellular automata with the Glauber algorithm. The critical point is calculated for different values of the inter- and intra-layer couplings (K₁≠K₂≠K₃≠K_z), where K₁, K₂ and K₃ are the nearest-neighbor interactions within each layer in the 1, 2 and 3 directions, respectively, and K_z is the intralayer coupling. A general ansatz equation for the critical point is given as a function of the inter- and intra-layer interactions, ξ=K₃/K₁,σ=K₂/K₁ and ω=K_z/K₁ for the one- and two-layer Ising models on the Kagome lattice.     

Constant coupling approximation for antiferromagnets with mixed ferro- and antiferromagnetic interactions

The constant-coupling approximation is extended to an antiferromagnetic spin-$\text{1}\text{2}$ system with two distinct anisotropic exchange interactions. The thermodynamic properties such as the transition temperature, magnetization, susceptibility and specific heat are discussed for three special cases: (i) ferro- and antiferromagnetic Ising interactions, (ii) isotropic ferro- and antiferromagnetic Heisenberg interactions; and (iii) isotropic ferromagnetic Heisenberg interactions and antiferromagnetic Ising interactions, allowing in each case for two different nearest-neighbour coupling constants. Numerical calculations have been performed for a layer structure with z = 6 intraplanar and z' = 6 interplanar nearest neighbours and the results are compared with those obtained in other approximations. Applying the theory to FeCl₂, the exchange constants are evaluated. It is shown that the calculated magnitudes of the interactions strongly depend upon the exchange-interaction model assumed.     

Lee–Yang Theorems and the Complexity of Computing Averages

We study the complexity of computing average quantities related to spin systems, such as the mean magnetization and susceptibility in the ferromagnetic Ising model, and the average dimer count (or average size of a matching) in the monomer-dimer model. By establishing connections between the complexity of computing these averages and the location of the complex zeros of the partition function, we show that these averages are #P-hard to compute, and hence, under standard assumptions, computationally intractable. In the case of the Ising model, our approach requires us to prove an extension of the famous Lee–Yang Theorem from the 1950s.     

Anisotropic compressible Ising model for quasi-one-dimensional hydrogen-bonded ferroelectrics

A compressible pseudo-spin Ising model hamiltonian is used to calculate the pressure-temperature phase diagram of quasi-one-dimensional hydrogen-bonded ferroelectric crystals such as CsD₂PO₄. We assume strong effective interactions, which are treated exactly, along chains, and weak volume dependent interactions, which are treated in the mean field approximation, between chains. In agreement with the experimental findings, the phase diagram exhibits a triple point and transition lines between ferroelectric, antiferroelectric, and paraelectric phases. However, the results suggest that the Ising model may be too simple to account for the detailed form of the phase diagram of CsD₂PO₄.     

Anomalies in Ising metamagnets

Ising metamagnets in a field with weak ferromagnetic intralayer interactions and highly coordinated antiferromagnetic interlayer couplings are studied using mean field theory and Monte Carlo simulations. In the antiferromagnetic phase anomalies in the magnetization and the specific heat are observed, reflecting the competing ordering tendencies of the external field and the interlayer couplings. Results are compared to recent experimental findings on FeBr₂.     

Inference of Kinetic Ising Model on Sparse Graphs

Based on dynamical cavity method, we propose an approach to the inference of kinetic Ising model, which asks to reconstruct couplings and external fields from given time-dependent output of original system. Our approach gives an exact result on tree graphs and a good approximation on sparse graphs, it can be seen as an extension of Belief Propagation inference of static Ising model to kinetic Ising model. While existing mean field methods to the kinetic Ising inference e.g., naïve mean-field, TAP equation and simply mean-field, use approximations which calculate magnetizations and correlations at time t from statistics of data at time t?1, dynamical cavity method can use statistics of data at times earlier than t?1 to capture more correlations at different time steps. Extensive numerical experiments show that our inference method is superior to existing mean-field approaches on diluted networks.