Exchange Monte Carlo Method and Application to Spin Glass Simulations

We propose an efficient Monte Carlo algorithm for simulating a ``hardly-relaxing" system, in which many replicas with different temperatures are simultaneously simulated and a virtual process exchanging configurations of these replica is introduced. This exchange process is expected to let the system at low temperatures escape from a local minimum. By using this algorithm the three-dimensional $\pm J$ Ising spin glass model is studied. The ergodicity time in this method is found much smaller than that of the multi-canonical method. In particular the time correlation function almost follows an exponential decay whose relaxation time is comparable to the ergodicity time at low temperatures. It suggests that the system relaxes very rapidly through the exchange process even in the low temperature phase.Comment: 10 pages + uuencoded 5 Postscript figures, REVTe

The Approximating Hamiltonian Method for the Imperfect Boson Gas

The pressure for the Imperfect (Mean Field) Boson gas can be derived in several ways. The aim of the present note is to provide a new method based on the Approximating Hamiltonian argument which is extremely simple and very general.Comment: 7 page

Parallelization of Markov chain generation and its application to the multicanonical method

We develop a simple algorithm to parallelize generation processes of Markov chains. In this algorithm, multiple Markov chains are generated in parallel and jointed together to make a longer Markov chain. The joints between the constituent Markov chains are processed using the detailed balance. We apply the parallelization algorithm to multicanonical calculations of the two-dimensional Ising model and demonstrate accurate estimation of multicanonical weights.Comment: 15 pages, 5 figures, uses elsart.cl

Exactness of the Bogoliubov approximation in random external potentials

We investigate the validity of the Bogoliubov c-number approximation in the case of interacting Bose-gas in a \textit{homogeneous random} media. To take into account the possible occurence of type III generalized Bose-Einstein condensation (i.e. the occurrence of condensation in an infinitesimal band of low kinetic energy modes without macroscopic occupation of any of them) we generalize the c-number substitution procedure to this band of modes with low momentum. We show that, as in the case of the one-mode condensation for translation-invariant interacting systems, this procedure has no effect on the exact value of the pressure in the thermodynamic limit, assuming that the c-numbers are chosen according to a suitable variational principle. We then discuss the relation between these c-numbers and the (total) density of the condensate

Simulated Tempering: A New Monte Carlo Scheme

We propose a new global optimization method ({\em Simulated Tempering}) for simulating effectively a system with a rough free energy landscape (i.e. many coexisting states) at finite non-zero temperature. This method is related to simulated annealing, but here the temperature becomes a dynamic variable, and the system is always kept at equilibrium. We analyze the method on the Random Field Ising Model, and we find a dramatic improvement over conventional Metropolis and cluster methods. We analyze and discuss the conditions under which the method has optimal performances.Comment: 12 pages, very simple LaTeX file, figures are not included, sorr

Multicanonical Recursions

The problem of calculating multicanonical parameters recursively is discussed. I describe in detail a computational implementation which has worked reasonably well in practice.Comment: 23 pages, latex, 4 postscript figures included (uuencoded Z-compressed .tar file created by uufiles), figure file corrected

On the Wang-Landau Method for Off-Lattice Simulations in the "Uniform" Ensemble

We present a rigorous derivation for off-lattice implementations of the so-called "random-walk" algorithm recently introduced by Wang and Landau [PRL 86, 2050 (2001)]. Originally developed for discrete systems, the algorithm samples configurations according to their inverse density of states using Monte-Carlo moves; the estimate for the density of states is refined at each simulation step and is ultimately used to calculate thermodynamic properties. We present an implementation for atomic systems based on a rigorous separation of kinetic and configurational contributions to the density of states. By constructing a "uniform" ensemble for configurational degrees of freedom--in which all potential energies, volumes, and numbers of particles are equally probable--we establish a framework for the correct implementation of simulation acceptance criteria and calculation of thermodynamic averages in the continuum case. To demonstrate the generality of our approach, we perform sample calculations for the Lennard-Jones fluid using two implementation variants and in both cases find good agreement with established literature values for the vapor-liquid coexistence locus.Comment: 21 pages, 4 figure

Scientific, institutional and personal rivalries among Soviet geographers in the late Stalin era

Scientific, institutional and personal rivalries between three key centres of geographical research and scholarship (the Academy of Sciences Institute of Geography and the Faculties of Geography at Moscow and Leningrad State Universities) are surveyed for the period from 1945 to the early 1950s. It is argued that the debates and rivalries between members of the three institutions appear to have been motivated by a variety of scientific, ideological, institutional and personal factors, but that genuine scientific disagreements were at least as important as political and ideological factors in influencing the course of the debates and in determining their final outcome

Grundstate Properties of the 3D Ising Spin Glass

We study zero--temperature properties of the 3d Edwards--Anderson Ising spin glass on finite lattices up to size $12^3$. Using multicanonical sampling we generate large numbers of groundstate configurations in thermal equilibrium. Finite size scaling with a zero--temperature scaling exponent $y = 0.74 \pm 0.12$ describes the data well. Alternatively, a descriptions in terms of Parisi mean field behaviour is still possible. The two scenarios give significantly different predictions on lattices of size $\ge 12^3$.Comment: LATEX 9pages,figures upon request ,SCRI-9

An efficient, multiple range random walk algorithm to calculate the density of states

We present a new Monte Carlo algorithm that produces results of high accuracy with reduced simulational effort. Independent random walks are performed (concurrently or serially) in different, restricted ranges of energy, and the resultant density of states is modified continuously to produce locally flat histograms. This method permits us to directly access the free energy and entropy, is independent of temperature, and is efficient for the study of both 1st order and 2nd order phase transitions. It should also be useful for the study of complex systems with a rough energy landscape.Comment: 4 pages including 4 ps fig