Second-order phase transitions in triangular ising models with two- and three-spin interactions

For Ising models with pair and three-spin interactions on the triangular lattice the nature of the phase diagram in the temperature-field plane is studied. Second-order transitions are located by the interface method of Müller-Hartmann and Zittartz.     

Inverse problem of reflection of electromagnetic waves from layered inhomogeneous plasma-like structures with refractive-index discontinuity

To reconstruct the distribution of free carriers in a layered inhomogeneous plasma-like half-space during the electromagnetic-wave reflection, we derive the Gelfrand-Levitan-Marchenko integral equation for structures with refractive-index discontinuity by the inverse problem method. Exact model solutions of the generalized equations are analyzed. The applicability limits and scope of these exact solutions for reconstruction problems are shown. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 50, No. 1, pp. 31–39, January 2007     

Exact results of a doubly decorated ising model with four-spin interactions

A doubly decorated Ising model with the crystal-field, two—and four-spin interactions is studied by applying the standard decoration-iteration transformation. Exact results for the critical boundaries, compensation temperatures and magnetization of the system are discussed in detail.     

Generalized motion by mean curvature as a macroscopic limit of stochastic ising models with long range interactions and Glauber dynamics

We study themacroscopic limit of an appropriately rescaledstochastic Ising model withlong range interactions evolving withGlauber dynamics as well as the correspondingmean field equation, which is nonlinear and nonlocal. In the limit we obtain an interface evolving with normal velocity k, wherek isthe mean curvature and thetransport coefficient is identified by aneffective Green-Kubo type formula. The above assertions are valid for all positive times, the motion of the interface being interpreted in theviscosity sense after the onset of the geometric singularities.Supported by ONRPartially supported by NSF, ARO, ONR and the Alfred P. Sloan Foundation     

Inverse vibration problem for un-damped 3-dimensional multi-story shear building models

Various researchers have contributed to the identification of the mass and stiffness matrices of two dimensional (2-D) shear building structural models for a given set of vibratory frequencies. The suggested methods are based on the specific characteristics of the Jacobi matrices, i.e., symmetric, tri-diagonal and semi-positive definite matrices. However, in case of three dimensional (3-D) structural models, those methods are no longer applicable, since their stiffness matrices are not tri-diagonal. In this paper the inverse problem for a special class of vibratory structural systems, i.e., 3-D shear building models, is investigated. A practical algorithm is proposed for solving the inverse eigenvalue problem for un-damped, 3-D shear buildings. The problem is addressed in two steps. First, using the target frequencies, a so-called normalized eigenvector matrix, which is a banded matrix containing the information related to the frequencies and mode shapes of the target structural system, is determined. In this regard, similar to the solution of inverse problem for 2-D shear building structural models in which an auxiliary structure is constructed by adding constraints (or springs) to the original system, three auxiliary structures are proposed to solve the problem for 3-D cases. In the second step, the normalized eigenvector matrix is utilized to obtain the normalized stiffness matrix; in turn, this matrix is decomposed into the stiffness and mass matrices of the system. Finally, a numerical example is presented to demonstrate the efficiency of the proposed algorithm in determining the mass and stiffness matrices of a 3-D structural model for a given set of target vibrational frequencies.     

Two-dimensional ising models with competing interaction—a Monte Carlo study

Two-dimensional Ising models on a square lattice with competing interactions along one axis or both axes are studied primarily by the Monte Carlo method. Several commensurate-incommensurate transitions are found. Dislocation-like configurations are identified near the sinusoidal — paramagnetic transition in accordance with the idea that the transition might be of Kosterlitz-Thouless, XY-like character.     

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.     

The inhomogeneous Suslov problem

We consider the Suslov problem of nonholonomic rigid body motion with inhomogeneous constraints. We show that if the direction along which the Suslov constraint is enforced is perpendicular to a principal axis of inertia of the body, then the reduced equations are integrable and, in the generic case, possess a smooth invariant measure. Interestingly, in this generic case, the first integral that permits integration is transcendental and the density of the invariant measure depends on the angular velocities. We also study the Painlevé property of the solutions.     

Modulated phase of an ising system with competing interactions on a Cayley tree

A system with competing nearest-neighbour and next-nearest-neighbour interactions is considered on a Cayley tree. The phase diagram contains a modulated phase, as found for similar models on periodic lattices, but the multicritical Lifshitz point is at zero temperature. The variation of the wavevector with temperature in the modulated phase is studied in detail, it shows narrow commensurate steps between incommensurate regions (incomplete devil's staircase). The behaviour of the coherence length near the different transitions is also analyzed.     

Efficient algorithm for random-bond ising models in 2D

We present an efficient algorithm for calculating the properties of Ising models in two dimensions, directly in the spin basis, without the need for mapping to fermion or dimer models. The algorithm computes the partition function and correlation functions at a single temperature on any planar network of N Ising spins in O(N;{3/2}) time or less. The method can handle continuous or discrete bond disorder and is especially efficient in the case of bond or site dilution, where it executes in O(NlnN) time near the percolation threshold. We demonstrate its feasibility on the ferromagnetic Ising model and the +/-J random-bond Ising model and discuss the regime of applicability in cases of full frustration such as the Ising antiferromagnet on a triangular lattice.     

Evidence for spinodal singularities in high-dimensional nearest-neighbor ising models

Conventional theories of nucleation predict that the metastable state has an average lifetime which monotonically decreases as the system is quenched further from the condensation point. However, theories based on the coarsegrained Ginzburg-Landau free energy functional seem to indicate that for systems above six dimensions there is a sharp spinodal dividing the metastable and unstable regimes where the lifetime of the metastable state diverges. Monte Carlo simulations are used to investigate this discrepency. Both nucleation rates and bulk susceptibility measurements seem to support the prediction of the Ginzburg-Landau theories.     

Two-dimensional periodic frustrated ising models in a transverse field

We investigate the interplay of classical degeneracy and quantum dynamics in a range of periodic frustrated transverse field Ising systems at zero temperature. We find that such dynamics can lead to unusual ordered phases and phase transitions or to a quantum spin liquid (cooperative paramagnetic) phase as in the triangular and kagome lattice antiferromagnets, respectively. For the latter, we further predict passage to a bond-ordered phase followed by a critical phase as the field is tilted. These systems also provide exact realizations of quantum dimer models introduced in studies of high temperature superconductivity.     

Properties of metastable ising models evolving under the Swendsen-Wang dynamics

The behavior of the metastable nearest neighbor Ising model governed by Swendsen-Wang dynamics (SW) is investigated ind=2. The results are compared to those obtained in standard Metropolis dynamics. Both the SW and Metropolis systems are observed to decay from the metastable state via the formation of nucleating droplets. Nucleation rates are measured and found to agree with those predicted by classical nucleation theory. The growth rates of the droplets are observed to differ between the two dynamics. In addition, the dynamic critical exponentz is measured in a mean-field (Curie-Weiss) metastable Ising model at the spinodal. It is found that for SW dynamics,z=2. Since this is the same value as that obtained in the Metropolis case, this result shows that SW does not change the dynamical universality class at the spinodal.     

Inverse problem in acoustoplasma

Using an experimentally obtained shape of distribution of the optical emission along the discharge tube an attempt is made to determine working and independent parameters of the acoustic oscillations and acoustoplasma by means of solving the inverse problem.     

Critical surface of the ising model with first-neighbor,second-neighbor,and four-spin interactions

Exact values are obtained for the slopesK ₁ ^c (0, 0)/K ₂,K ₁ ^c (0, 0)/K ₄ of the critical surface of paramagnetic-ferromagnetic transitionsK ₁ ^c (K ₂,K ₄) for the two-dimensional Ising model on a square lattice with first-neighbor, second-neighbor, and four-spin couplingsK ₁,K ₂, andK ₄, respectively. The results are obtained using universality arguments to relate the slopes to known spin-spin correlation functions forK ₂=K ₄=0. The equivalence of different expressions for the slopes in terms of correlation functions yields sum rules for the divergent part of certain sums over the second-neighbor and four-spin energy-energy correlation functions.     

Exactly solvable one-dimensional inhomogeneous models

We present a simple way of constructing one-dimensional inhomogeneous models (random or quasiperiodic) which can be solved exactly. We treat the example of an Ising chain in a varying magnetic field, but our procedure can easily be extended to other one-dimensional inhomogeneous models. For all the models we can construct, the free energy and its derivatives with respect to temperature can be computed exactly at one particular temperature.     

Ordering and phase transitions in ising systems with competing short range and dipolar interactions

As a simple model of order-disorder ferroelectrics or dipolar magnets we consider a simple cubic Ising-system with nearest neighbor exchangeJ and dipolar interaction of strengthµ ²/a ³. ForJa ³/µ ²Ja ³/µ ²<0.1270 the ground state consists of ferromagnetic rows (in spin direction) arranged antiferromagnetically in the plane perpendicular to it. AtJa ³/µ ²=0.1279 the structure changes to a layered antiferromagnetic structure with a twocomponent order parameter, while forJa ³/µ ²>0.16429 the ferromagnetic phase becomes stable (with domain arrangements depending on the shape of the sample). For all critical values ofJa ³/µ ² where the bulk energies of two phases become equal also the interface energy between these phases is found to be zero. The ordering at nonzero temperature is studied by means of mean-field approximations (MFA) and Monte Carlo (MC) calculations. It turns out that forJa ³/µ ² of order unity the MFA overestimates ordering temperatures by about a factor of two, and predicts multicritical points (between the disordered and two ordered phases) at nonzero temperature, including two biaxial Lifshitz points which the MC work suggests to occur atT=0. In contrast to MFA the layered antiferromagnetic structure is found to be stable only at extremely lowT, because a metastable spin-glass phase (with random arrangement of ferromagnetic rows in the spin direction) has only slightly higher energy. The MFA also yields two regimes of helical phases which are “locked in” to the antiferromagnetic phases at uniaxial Lifshitz points occurring at the Brillouin zone boundary. In the MC-work various methods of treating the long-range interaction are investigated. While all kinds of truncations as well as compensating field methods are rather unsatisfactory in our case, Ewald summation techniques yield satisfactory results. Nevertheless strong fluctuations as well as strong finite size effects prevent us from making accurate exponent estimates, but arguments are given that there is no regime of broad visibility of Landaulike critical behavior. Finally the extension of our results to other lattices as well as experimental applications are briefly discussed.     

Correlations in inhomogeneous Ising models

We study correlations in inhomogeneous Ising models on a square lattice. The nearest neighbour couplings are allowed to be of arbitrary strength and sign such that the coupling distribution is translationally invariant either in horizontal or in diagonal direction, i.e. the models have a layered structure. By using transfer matrix techniques the spin-spin correlations are calculated parallel to the layering and are expressed as Toeplitz determinants. After working out the general methods we discuss two special examples in detail: the fully frustrated square lattice (FFS) and the chessboard model, both having no phase transition. At zero temperature correlations in the chessboard model decay exponentially, while in the FFS model one has algebraic decay with a critical index =1/2, i.e.T=0 is a critical point. At finite temperature we find exponential decay in both models with a correlation length determined by the excitation gap in the fermion spectrum. Due to frustration correlations may develop on oscillatory structure and spins separated by an odd diagonal distance are totally uncorrelated at all temperatures.Work performed within the research program of the Sonderforschungsbereich 125 Aachen-Jülich-Köla     

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 horizontal direction. We calculate correlations parallel to the layering in the horizontally layered model with periodv=2. If the model has a finite critical temperature,T _c>0, the order parameter in the frustrated case may become discontinuous forT0. Correlations atT=T _c decay algebraically with critical exponent =1/4 and exponentially forT>T _c. 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.Work performed within the research program of the Sonder forschungsbereich 125 Aachen-Jülich-Köln