Extending the use of canonical and microcanonical Monte Carlo algorithms to spin models

We describe the canonical and microcanonical Monte Carlo algorithms for different systems that can be described by spin models. Sites of the lattice, chosen at random, interchange their spin values, provided they are different. The canonical ensemble is generated by performing exchanges according to the Metropolis prescription whereas in the microcanonical ensemble, exchanges are performed as long as the total energy remains constant. A systematic finite size analysis of intensive quantities and a comparison with results obtained from distinct ensembles are performed and the quality of results reveal that the present approach may be an useful tool for the study of phase transitions, specially first-order transitions. (C) 2009 Elsevier B.V. All rights reserved.

Phase diagram and spectral properties of a new exactly integrable spin-1 quantum chain

The spectral properties and phase diagram of the exactly integrable spin-1 quantum chain introduced by Alcaraz and Bariev are presented. The model has a U(1) symmetry and its integrability is associated with an unknown R-matrix whose dependence on the spectral parameters is not of a different form. The associated Bethe ansatz equations that fix the eigenspectra are distinct from those associated with other known integrable spin models. The model has a free parameter t(p). We show that at the special point t(p) = 1, the model acquires an extra U(1) symmetry and reduces to the deformed SU(3) Perk-Schultz model at a special value of its anisotropy q = exp(i2 pi/3) and in the presence of an external magnetic field. Our analysis is carried out either by solving the associated Bethe ansatz equations or by direct diagonalization of the quantum Hamiltonian for small lattice sizes. The phase diagram is calculated by exploring the consequences of conformal invariance on the finite-size corrections of the Hamiltonian eigenspectrum. The model exhibits a critical phase ruled by the c = 1 conformal field theory separated from a massive phase by first-order phase transitions.

TEM, XRD and AFM study of poly(o-ethoxyaniline) films: new evidence for the formation of conducting islands

The existence of conducting islands in polyaniline films has long been proposed in the literature, which would be consistent with conducting mechanisms based on hopping. Obtaining direct evidence of conducting islands, however, is not straightforward. In this paper, conducting islands were visualized in poly(o-ethoxyaniline) (POEA) films prepared at low pH, using Transmission Electron Microscopy (TEM) and atomic force spectroscopy (AFS). The size of the islands varied between 67 and 470 angstrom for a pH=3.0, with a larger average being obtained with AFS, probably due to the finite size effect of the atomic force microscopy tip. In AFS, the conducting islands were denoted by regions with repulsive forces due to the double-layer forces. On the basis of X-ray diffraction (XRD) patterns for POEA in the powder form, we infer that the conducting islands are crystalline, and therefore a POEA film is believed to consist of conducting islands dispersed in an insulating, amorphous matrix. From conductivity measurements we inferred the charge transport to be governed by a typical quasi-one dimensional variable range hopping (VRH) mechanism.

Density-Profile Processes Describing Biological Signaling Networks: Almost Sure Convergence to Deterministic Trajectories

We introduce jump processes in R(k), called density-profile processes, to model biological signaling networks. Our modeling setup describes the macroscopic evolution of a finite-size spin-flip model with k types of spins with arbitrary number of internal states interacting through a non-reversible stochastic dynamics. We are mostly interested on the multi-dimensional empirical-magnetization vector in the thermodynamic limit, and prove that, within arbitrary finite time-intervals, its path converges almost surely to a deterministic trajectory determined by a first-order (non-linear) differential equation with explicit bounds on the distance between the stochastic and deterministic trajectories. As parameters of the spin-flip dynamics change, the associated dynamical system may go through bifurcations, associated to phase transitions in the statistical mechanical setting. We present a simple example of spin-flip stochastic model, associated to a synthetic biology model known as repressilator, which leads to a dynamical system with Hopf and pitchfork bifurcations. Depending on the parameter values, the magnetization random path can either converge to a unique stable fixed point, converge to one of a pair of stable fixed points, or asymptotically evolve close to a deterministic orbit in Rk. We also discuss a simple signaling pathway related to cancer research, called p53 module.

Propriedades críticas do processo epidêmico difusivo com interação de Lévy

The diffusive epidemic process (PED) is a nonequilibrium stochastic model which, exhibits a phase trnasition to an absorbing state. In the model, healthy (A) and sick (B) individuals diffuse on a lattice with diffusion constants DA and DB, respectively. According to a Wilson renormalization calculation, the system presents a first-order phase transition, for the case DA > DB. Several researches performed simulation works for test this is conjecture, but it was not possible to observe this first-order phase transition. The explanation given was that we needed to perform simulation to higher dimensions. In this work had the motivation to investigate the critical behavior of a diffusive epidemic propagation with Lévy interaction(PEDL), in one-dimension. The Lévy distribution has the interaction of diffusion of all sizes taking the one-dimensional system for a higher-dimensional. We try to explain this is controversy that remains unresolved, for the case DA > DB. For this work, we use the Monte Carlo Method with resuscitation. This is method is to add a sick individual in the system when the order parameter (sick density) go to zero. We apply a finite size scalling for estimates the critical point and the exponent critical =, e z, for the case DA > DB

Universaidade em sistemas da mecâniaestatística de não equilíbrio com estados absorventes e percolação geográfica

Complex systems have stimulated much interest in the scientific community in the last twenty years. Examples this area are the Domany-Kinzel cellular automaton and Contact Process that are studied in the first chapter this tesis. We determine the critical behavior of these systems using the spontaneous-search method and short-time dynamics (STD). Ours results confirm that the DKCA e CP belong to universality class of Directed Percolation. In the second chapter, we study the particle difusion in two models of stochastic sandpiles. We characterize the difusion through diffusion constant D, definite through in the relation h(x)2i = 2Dt. The results of our simulations, using finite size scalling and STD, show that the diffusion constant can be used to study critical properties. Both models belong to universality class of Conserved Directed Percolation. We also study that the mean-square particle displacement in time, and characterize its dependence on the initial configuration and particle density. In the third chapter, we introduce a computacional model, called Geographic Percolation, to study watersheds, fractals with aplications in various areas of science. In this model, sites of a network are assigned values between 0 and 1 following a given probability distribution, we order this values, keeping always its localization, and search pk site that percolate network. Once we find this site, we remove it from the network, and search for the next that has the network to percole newly. We repeat these steps until the complete occupation of the network. We study the model in 2 and 3 dimension, and compare the bidimensional case with networks form at start real data (Alps e Himalayas)

Expansion effects on back-to-back correlations

The back-to-back correlations (BBC) of particle-antiparticle pairs, signalling in-medium mass modification, are studied in a finite size thermalized medium. The width of BBC function is explicitly evaluated in the case of a nonrelativistic spherically symmetric expanding fireball. The effect of the flow is to reduce the BBC signal as compared to the case of non flow. Nevertheless, a significant signal survives finite-time emission plus expansion effects.

Chaotic thermalization in Yang-Mills-Higgs theory on a spacial lattice

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

An introduction to numerical methods in low-dimensional quantum systems

This is an introductory course to the Lanczos Method and Density Matrix Renormalization Group Algorithms (DMRG), two among the leading numerical techniques applied in studies of low-dimensional quantum models. The idea of studying the models on clusters of a finite size in order to extract their physical properties is briefly discussed. The important role played by the model symmetries is also examined. Special emphasis is given to the DMRG.

Monte Carlo renormalization group with evolution in the space of parameters

Using data from a single simulation we obtain Monte Carlo renormalization-group information in a finite region of parameter space by adapting the Ferrenberg-Swendsen histogram method. Several quantities are calculated in the two-dimensional N 2 Ashkin-Teller and Ising models to show the feasibility of the method. We show renormalization-group Hamiltonian flows and critical-point location by matching of correlations by doing just two simulations at a single temperature in lattices of different sizes to partially eliminate finite-size effects.

A theoretical analysis on electronic structure of the (110) surface of TiO2-SnO2 mixed oxide

Mixed oxide compounds, such as TiO2-SnO2 system are widely used as gas sensors and should also provide varistor properties modifying the TiO2 surface. Therefore, a theoretical investigation has been carried out characterizing the effect of SnO2 on TiO2 addition on the electronic structure by means of ab initio SCF-LCAO calculations using all electrons. In order to take into account the finite size of the cluster, we have used the point charge model for the (TiO2)(15) cluster to study the effect on electronic structure of doping the TiO2 (110) Surface. The contracted basis set for titanium (4322/42/3), oxygen (33/3) and tin (43333/4333/43) atoms were used. The charge distributions, dipole moments, and density of states of doping TiO2 and vacancy formation are reported and analysed. (C) 2003 Elsevier B.V. All rights reserved.

Theoretical analysis on TiO(2)(110)/V surface

Theoretical analysis based on the Hartree-Fock method were performed in order to study the stoichiometric TiO(2) (110) surface and the vanadium substituted system. The Pople with polarization 3-21G* basis set level was used. The TiO(2) (110) surface was modeled using a (TiO(2))(15) cluster model. In order to take into account the finite size of the cluster, we have studied two different models: the point charge and the hydrogen saturated methodologies. The charge values used in the point charge calculations were optimized. The density of states, orbital self-consistend field (SCF) energies, and Mulliken charge values were analyzed. The method and model's dependence on the analyzed results are discussed. The theoretical results are compared with available experimental data. (C) 2001 John Wiley & Sons, Inc.

Vortex lattice and matching fields for a long superconducting wire

We investigate the flux penetration patterns and matching fields of a long cylindrical wire of circular cross section in the presence of an external magnetic field. For this study we write the London theory for a long cylinder both for the mixed and Meissner states, with boundary conditions appropriate for this geometry. Using the Monte Carlo simulated annealing method, the free energy of the mixed state is minimized with respect to the vortex position and we obtain the ground state of the vortex lattice for N=3 up to 18 vortices. The free energy of the Meissner and mixed states provides expressions for the matching fields. We find that, as in the case of samples of different geometry, the finite-size effect provokes a delay on the vortex penetration and a vortex accumulation in the center of the sample. The vortex patterns obtained are in good agreement with experimental results.

Update on the status of hadronic squeezed correlations at RHIC energies

In high energy heavy ion collisions a hot and dense medium is formed, where the hadronic masses may be shifted from their asymptotic values. If this mass modification occurs, squeezed back-to-back correlations (BBC) of particle-antiparticle pairs are predicted to appear, both in the femionic (fBBC) and in the bosonic (bBBC) sectors. Although they have unlimited intensity even for finite-size expanding systems, these hadronic squeezed correlations are very sensitive to their time emission distribution. Here we discuss results in case this time emission is parameterized by a Lévy-type distribution, showing that it reduces the signal even more dramatically than a Lorentzian distribution, which already reduces the intensity of the effect by orders of magnitude, as compared to the sudden emission. However, we show that the signal could still survive if the duration of the process is short, and if the effect is searched for lighter mesons, such as kaons. We compare some of our results to recent PHENIX preliminary data on squeezed correlations of K +K - pairs. © 2011 Pleiades Publishing, Ltd.

Influence of thermal gradient in vortex states of mesoscopic superconductors

In general, the studies of finite size effects in mesoscopic superconductors have been carried out in such a way that the temperature parameter is constant in the entire system. However, we could have situations where a real sample is near a heater source, as an example. In such situations, gradients of temperature are present. On the other hand, mesoscopic superconductors are interesting systems due to the fact that they present confinement effects which influence all the vortex dynamics. Thus, in this work we studied the influence of thermal gradients on the vortex dynamics in mesoscopic superconductors. For this purposes, we used the time dependent Ginzburg-Landau equations. The thermal gradients produce an asymmetric distribution of the currents around the system which, in turn, yield interesting vortex configurations and difficult the formation of giant vortices.