Unusual Yang-Lee edge singularity in the one-dimensional axial-next-to-nearest-neighbor Ising model

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Dynamic critical behavior of percolation observables in the 2d Ising model

We present preliminary results of our numerical study of the critical dynamics of percolation observables for the two-dimensional Ising model. We consider the (Monte-Carlo) short-time evolution of the system obtained with a local heat-bath method and with the global Swendsen-Wang algorithm. In both cases, we find qualitatively different dynamic behaviors for the magnetization and Omega, the order parameter of the percolation transition. This may have implications for the recent attempts to describe the dynamics of the QCD phase transition using cluster observables.

Dynamic behavior of cluster observables for the 2d Ising model

We present results of our numerical study of the critical dynamics of percolation observables for the two-dimensional Ising model. We consider the (Monte Carlo) short-time evolution of the system with small initial magnetization and heat-bath dynamics. We find qualitatively different dynamic behaviors for the magnetization M and for Ω, the so-called strength of the percolating cluster, which is the order parameter of the percolation transition. More precisely, we obtain a (leading) exponential form for Ω as a function of the Monte Carlo time t, to be compared with the power-law increase encountered for M at short times. Our results suggest that, although the descriptions in terms of magnetic or percolation order parameters may be equivalent in the equilibrium regime, greater care must be taken to interpret percolation observables at short times.

The Yang-Lee zeros of the 1D Blume-Capel model on connected and non-connected rings

We carry out a numerical and analytic analysis of the Yang-Lee zeros of the ID Blume-Capel model with periodic boundary conditions and its generalization on Feynman diagrams for which we include sums over all connected and nonconnected rings for a given number of spins. In both cases, for a specific range of the parameters, the zeros originally on the unit circle are shown to depart from it as we increase the temperature beyond some limit. The curve of zeros can bifurcate- and become two disjoint arcs as in the 2D case. We also show that in the thermodynamic limit the zeros of both Blume-Capel models on the static (connected ring) and on the dynamical (Feynman diagrams) lattice tend to overlap. In the special case of the 1D Ising model on Feynman diagrams we can prove for arbitrary number of spins that the Yang-Lee zeros must be on the unit circle. The proof is based on a property of the zeros of Legendre polynomials.

Propriedades magnéticas do modelo de Ising bidimensional em uma rede bipartida de spins 1/2 e 1: um estudo de ferrimagnetismo com vacâncias

Pós-graduação em Ciência e Tecnologia de Materiais - FC

Estudo de percolação de clusters de Monte Carlo para o modelo de Ising bidimensional

Pós-graduação em Física - IFT

Yang-Lee zeros of the two- and three-state Potts model defined on π3 Feynman diagrams

We present both analytical and numerical results on the position of partition function zeros on the complex magnetic field plane of the q=2 state (Ising) and the q=3 state Potts model defined on phi(3) Feynman diagrams (thin random graphs). Our analytic results are based on the ideas of destructive interference of coexisting phases and low temperature expansions. For the case of the Ising model, an argument based on a symmetry of the saddle point equations leads us to a nonperturbative proof that the Yang-Lee zeros are located on the unit circle, although no circle theorem is known in this case of random graphs. For the q=3 state Potts model, our perturbative results indicate that the Yang-Lee zeros lie outside the unit circle. Both analytic results are confirmed by finite lattice numerical calculations.

Short-time dynamics of percolation observables

We consider the critical short-time evolution of magnetic and droplet-percolation order parameters for the Ising model in two and three dimensions, through Monte Carlo simulations with the (local) heat-bath method. We find qualitatively different dynamic behaviors for the two types of order parameters. More precisely, we find that the percolation order parameter does not have a power-law behavior as encountered for the magnetization, but develops a scale (related to the relaxation time to equilibrium) in the Monte Carlo time. We argue that this difference is due to the difficulty in forming large clusters at the early stages of the evolution. Our results show that, although the descriptions in terms of magnetic and percolation order parameters may be equivalent in the equilibrium regime, greater care must be taken to interpret percolation observables at short times. In particular, this concerns the attempts to describe the dynamics of the deconfinement phase transition in QCD using cluster observables.

Percolation of Monte Carlo clusters

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A aproximação de campo médio de Bethe-Peierls

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Novo comportamento crítico da singularidade de Yang-Lee

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Using the Ferrenberg-Swendsen Monte Carlo method for nonuniversal exponents

We use the method of Monte Carlo evolution in the coupling constant space of Ferrenberg and Swendsen to evaluate the nonuniversal exponent η* associated to a linear defect in a 2d Ising model. © 1989.

The energy landscape for solvent dynamics in electron transfer reactions: A minimalist model

Energy fluctuations of a solute molecule embedded in a polar solvent are investigated to depict the energy landscape for solvation dynamics. The system is modeled by a charged molecule surrounded by two layers of solvent dipolar molecules with simple rotational dynamics. Individual solvent molecules are treated as simple dipoles that can point toward or away from the central charge (Ising spins). Single-spin-flip Monte Carlo kinetics simulations are carried out in a two-dimensional lattice for different central charges, radii of outer shell, and temperatures. By analyzing the density of states as a function of energy and temperatures, we have determined the existence of multiple freezing transitions. Each of them can be associated with the freezing of a different layer of the solvent. (C) 2002 American Institute of Physics.

Agent-based model to rural-urban migration analysis

In this paper, we analyze the rural-urban migration phenomenon as it is usually observed in economies which are in the early stages of industrialization. The analysis is conducted by means of a statistical mechanics approach which builds a computational agent-based model. Agents are placed on a lattice and the connections among them are described via an Ising-like model. Simulations on this computational model show some emergent properties that are common in developing economies, such as a transitional dynamics characterized by continuous growth of urban population, followed by the equalization of expected wages between rural and urban sectors (Harris-Todaro equilibrium condition), urban concentration and increasing of per capita income. (c) 2005 Elsevier B.V. All rights reserved.

A Monte Carlo study of the anisotropic N=3 Ashkin-Teller model

The phase diagram of an asymmetric N = 3 Ashkin-Teller model is obtained by a numerical analysis which combines Monte Carlo renormalization group and reweighting techniques. Present results reveal several differences with those obtained by mean-field calculations and a Hamiltonian approach. In particular, we found Ising critical exponents along a line where Goldschmidt has located the Kosterlitz-Thouless multicritical point. On the other hand, we did find nonuniversal exponents along another transition line. Symmetry breaking in this case is very similar to the N = 2 case, since the symmetries associated to only two of the Ising variables are broken. However, for large values of the coupling constant ratio XW = W/K, when the only broken symmetry is of a hidden variable, we detected first-order phase transitions giving evidence supporting the existence of a multicritical point, as suggested by Goldschmidt, but in a different region of the phase diagram. © 2002 Elsevier Science B.V. All rights reserved.