Showing papers in "arXiv: Statistical Mechanics in 1998"

Inverse Cubic Law for the Probability Distribution of Stock Price Variations

Abstract: The probability distribution of stock price changes is studied by analyzing a database (the Trades and Quotes Database) documenting every trade for all stocks in three major US stock markets, for the two year period Jan 1994 -- Dec 1995. A sample of 40 million data points is extracted, which is substantially larger than studied hitherto. We find an asymptotic power-law behavior for the cumulative distribution with an exponent alpha approximately 3, well outside the Levy regime 0< alpha <2.

Discrete molecular dynamics studies of the folding of a protein-like model

Abstract: Background: Many attempts have been made to resolve in time the folding of model proteins in computer simulations. Different computational approaches have emerged. Some of these approaches suffer from the insensitivity to the geometrical properties of the proteins (lattice models), while others are computationally heavy (traditional MD). Results: We use a recently-proposed approach of Zhou and Karplus to study the folding of the protein model based on the discrete time molecular dynamics algorithm. We show that this algorithm resolves with respect to time the folding --- unfolding transition. In addition, we demonstrate the ability to study the coreof the model protein. Conclusion: The algorithm along with the model of inter-residue interactions can serve as a tool to study the thermodynamics and kinetics of protein models.

Spontaneous Stratification in Granular Mixtures

Abstract: Granular materials size segregate when exposed to external periodic perturbations such as vibrations. Moreover, mixtures of grains of different sizes spontaneously segregate in the absence of external perturbations: when a mixture is simply poured onto a pile, the large grains are more likely to be found near the base, while the small grains are more likely to be near the top. Here, we report a spontaneous phenomenon arising when we pour a mixture between two vertical plates: the mixture spontaneously stratifies into alternating layers of small and large grains whenever the large grains are rougher than the small grains. In contrast, we find only spontaneous segregation when the large grains are more rounded than the small grains. The stratification is related to the occurrence of avalanches; during each avalanche the grains comprising the avalanche spontaneously stratify into a pair of layers through a "kink" mechanism, with the small grains forming a sublayer underneath the layer of large grains.

First excitations of the spin 1/2 Heisenberg antiferromagnet on the kagom\'e lattice

Abstract: We study the exact low energy spectra of the spin 1/2 Heisenberg antiferromagnet on small samples of the kagom\'e lattice of up to N=36 sites. In agreement with the conclusions of previous authors, we find that these low energy spectra contradict the hypothesis of N\'eel type long range order. Certainly, the ground state of this system is a spin liquid, but its properties are rather unusual. The magnetic ($\Delta S=1$) excitations are separated from the ground state by a gap. However, this gap is filled with nonmagnetic ($\Delta S=0$) excitations. In the thermodynamic limit the spectrum of these nonmagnetic excitations will presumably develop into a gapless continuum adjacent to the ground state. Surprisingly, the eigenstates of samples with an odd number of sites, i.e. samples with an unsaturated spin, exhibit symmetries which could support long range chiral order. We do not know if these states will be true thermodynamic states or only metastable ones. In any case, the low energy properties of the spin 1/2 Heisenberg antiferromagnet on the kagom\'e lattice clearly distinguish this system from either a short range RVB spin liquid or a standard chiral spin liquid. Presumably they are facets of a generically new state of frustrated two-dimensional quantum antiferromagnets.

Self-Organization Phenomena in Pedestrian Crowds

Abstract: Pedestrian crowds can very realistically be simulated with a social force model which describes the different influences affecting individual pedestrian motion by a few simple force terms. The model is able to reproduce the emergence of several empirically observed collective patterns of motion. These self-organization phenomena can be utilized for new flow optimization methods which are indispensable for skilful town- and traffic-planning.

Identifying the protein folding nucleus using molecular dynamics

Abstract: Molecular dynamics simulations of folding in an off-lattice protein model reveal a nucleation scenario, in which a few well-defined contacts are formed with high probability in the transition state ensemble of conformations. Their appearance determines folding cooperativity and drives the model protein into its folded conformation.

Coherent Moving States in Highway Traffic (Originally: Moving Like a Solid Block)

Abstract: Recent advances in multiagent simulations have made possible the study of realistic traffic patterns and allow to test theories based on driver behaviour. Such simulations also display various empirical features of traffic flows, and are used to design traffic controls that maximise the throughput of vehicles in heavily transited highways. In addition to its intrinsic economic value, vehicular traffic is of interest because it may throw light on some social phenomena where diverse individuals competitively try to maximise their own utilities under certain constraints. In this paper, we present simulation results that point to the existence of cooperative, coherent states arising from competitive interactions that lead to a new phenomenon in heterogeneous highway traffic. As the density of vehicles increases, their interactions cause a transition into a highly correlated state in which all vehicles practically move with the same speed, analogous to the motion of a solid block. This state is associated with a reduced lane changing rate and a safe, high and stable flow. It disappears as the vehicle density exceeds a critical value. The effect is observed in recent evaluations of Dutch traffic data.

Coherent versus Incoherent Dynamics during Bose-Einstein Condensation in Atomic Gases

Abstract: We review and extend the theory of the dynamics of Bose-Einstein condensation in weakly interacting atomic gases. We present in a unified way both the semiclassical theory as well as the full quantum theory. This is achieved by deriving a Fokker-Planck equation that incorporates both the coherent and incoherent effects of the interactions in a dilute Bose gas. In first instance we focus our attention on the nonequilibrium dynamics of a homogeneous Bose gas with a positive interatomic scattering length. After that we discuss how our results can be generalized to the inhomogeneous situation that exists in the present experiments with magnetically trapped alkali gases, and how we can deal with a negative interatomic scattering length in that case as well. We also show how to arrive at a discription of the collective modes of the gas that obeys the Kohn theorem at all temperatures. The theory is based on the many-body T-matrix approximation throughout, since this approximation has the correct physical behavior near the critical temperature and also treats the coherent and incoherent processes taking place in the gas on an equal footing.

A Mathematical Model for the Behavior of Individuals in a Social Field

Abstract: Related to an idea of Lewin, a mathematical model for behavioral changes under the influence of a social field is developed. The social field reflects public opinion, social norms and trends. It is not only given by external factors (the environment) but also by the interactions of individuals. Two important kinds of interaction processes are distinguished: Imitative and avoidance processes. Variations of individual behavior are taken into account by ``diffusion coefficients''.

Resonance between Noise and Delay

Abstract: We propose here a stochastic binary element whose transition rate depends on its state at a fixed interval in the past With this delayed stochastic transition this is one of the simplest dynamical models under the influence of ``noise'' and ``delay'' We demonstrate numerically and analytically that we can observe resonant phenomena between the oscillatory behavior due to noise and that due to delay

Finite size scaling in three-dimensional bootstrap percolation

Abstract: We consider the problem of bootstrap percolation on a three dimensional lattice and we study its finite size scaling behavior. Bootstrap percolation is an example of Cellular Automata defined on the $d$-dimensional lattice $\{1,2,...,L\}^d$ in which each site can be empty or occupied by a single particle; in the starting configuration each site is occupied with probability $p$, occupied sites remain occupied for ever, while empty sites are occupied by a particle if at least $\ell$ among their $2d$ nearest neighbor sites are occupied. When $d$ is fixed, the most interesting case is the one $\ell=d$: this is a sort of threshold, in the sense that the critical probability $p_c$ for the dynamics on the infinite lattice ${\Bbb Z}^d$ switches from zero to one when this limit is crossed. Finite size effects in the three-dimensional case are already known in the cases $\ell\le 2$: in this paper we discuss the case $\ell=3$ and we show that the finite size scaling function for this problem is of the form $f(L)={\mathrm{const}}/\ln\ln L$. We prove a conjecture proposed by A.C.D. van Enter.

Introduction to the Bethe ansatz I

Abstract: The Bethe ansatz for the one-dimensional s=1/2 Heisenberg ferromagnet is introduced at an elementary level. The presentation follows Bethe's original work very closely. A detailed description and a complete classification of all two-magnon scattering states and two-magnon bound states are given for finite and infinite chains. The paper is designed as a tutorial for beginning graduate students. It includes 10 problems for further study.

The approach to thermal equilibrium in quantized chaotic systems

Abstract: We consider many-body quantum systems that exhibit quantum chaos, in the sense that the observables of interest act on energy eigenstates like banded random matrices. We study the time-dependent expectation values of these observables, assuming that the system is in a definite (but arbitrary) pure quantum state. We induce a probability distribution for the expectation values by treating the zero of time as a uniformly distributed random variable. We show explicitly that if an observable has a nonequilibrium expectation value at some particular moment, then it is overwhelmingly likely to move towards equilibrium, both forwards and backwards in time. For deviations from equilibrium that are not much larger than a typical quantum or thermal fluctuation, we find that the time dependence of the move towards equilibrium is given by the Kubo correlation function, in agreement with Onsager's postulate. These results are independent of the details of the system's quantum state.

Neutral evolution of model proteins: diffusion in sequence space and overdispersion

Abstract: We simulate the evolution of model protein sequences subject to mutations. A mutation is considered neutral if it conserves 1) the structure of the ground state, 2) its thermodynamic stability and 3) its kinetic accessibility. All other mutations are considered lethal and are rejected. We adopt a lattice model, amenable to a reliable solution of the protein folding problem. We prove the existence of extended neutral networks in sequence space -- sequences can evolve until their similarity with the starting point is almost the same as for random sequences. Furthermore, we find that the rate of neutral mutations has a broad distribution in sequence space. Due to this fact, the substitution process is overdispersed (the ratio between variance and mean is larger than one). This result is in contrast with the simplest model of neutral evolution, which assumes a Poisson process for substitutions, and in qualitative agreement with biological data.

Field theory fo charged fluids and colloids

Abstract: A systematic field theory is presented for charged systems. The one-loop level corresponds to the classical Debye-Huckel (DH) theory, and exhibits the full hierarchy of multi-body correlations determined by pair-distribution functions given by the screened DH potential. Higher-loop corrections can lead to attractive pair interactions between colloids in asymmetric ionic environments. The free energy follows as a loop-wise expansion in half-integer powers of the density; the resulting two-phase demixing region shows pronounced deviations from DH theory for strongly charged colloids.

Interrelations between Stochastic Equations for Systems with Pair Interactions

Abstract: Several types of stochastic equations are important in thermodynamics, chemistry, evolutionary biology, population dynamics and quantitative social science. For systems with pair interactions four different types of equations are derived, starting from a master equation for the state space: First, general mean value and (co)variance equations. Second, Boltzmann-like equations. Third, a master equation for the configuration space allowing transition rates which depend on the occupation numbers of the states. Fourth, a Fokker-Planck equation and a ``Boltzmann-Fokker-Planck equation''. The interrelations of these equations and the conditions for their validity are worked out clearly. A procedure for a selfconsistent solution of the nonlinear equations is proposed. Generalizations to interactions between an arbitrary number of systems are discussed.

Multiscale behaviour of volatility autocorrelations in a financial market

Abstract: We perform a scaling analysis on NYSE daily returns. We show that volatility correlations are power-laws on a time range from one day to one year and, more important, that they exhibit a multiscale behaviour.

Energy Transport in the Integrable System in Contact with Various Types of Phonon Reservoirs

Abstract: We study how energy transport in an integrable system is affected by the spectral densities of heat reservoirs. The model investigated here is the quantum harmonic chain whose both ends are in contact with two heat reservoirs at different temperatures. The master equation for the reduced density matrix is derived on the assumption that the reservoirs are composed of an infinite number of independent harmonic oscillators. We evaluate temperature profile and energy flux in the stationary state for the master equation and discuss how they depend on the types of spectral densities. When we attach the reservoirs of the same type of spectral density, we find that the temperature profile is independent of the types. On the other hand, when the two reservoirs have different types of spectral densities, the energy profile near the ends of the chain depends on the types. When the coupling is finite, the temperature profile near the ends shows wide variation of behavior dependent on spectral densities and temperatures of reservoirs. This dependence is discussed with the Fokker-Planck equations obtained in the classical limit.

A Stochastic Liouvillian Algorithm to Simulate Dissipative Quantum Dynamics With Arbitrary Precision

Abstract: An exact and efficient new method to simulate dynamics in dissipative quantum systems is presented. A stochastic Liouville equation, deduced from Feynman and Vernon's path-integral expression of the reduced density matrix, is used to describe the exact dynamics at any dissipative strength and for arbitrarily low temperatures. The utility of the method is demonstrated by applications to a damped harmonic oscillator and a double-well system immersed in an Ohmic bath at low temperatures.

Spontaneous Breaking of Translational Invariance and Spatial Condensation in Stationary States on a Ring: I. The Neutral System

Abstract: We consider a model in which positive and negative particles with equal densities diffuse in an asymmetric, CP invariant way on a ring. The positive particles hop clockwise, the negative counter-clockwise and oppositely-charged adjacent particles may swap positions. The model depends on two parameters. Analytic calculations using quadratic algebras, inhomogeneous solutions of the mean-field equations and Monte-Carlo simulations suggest that the model has three phases. A pure phase in which one has three pinned blocks of only positive, negative particles and vacancies and in which translational invariance is broken. A mixed phase in which the current has a linear dependence on one parameter but is independent of the other one and of the density of the charged particles. In this phase one has a bump and a fluid. The bump (condensate) contains positive and negative particles only, the fluid contains charged particles and vacancies uniformly distributed. The mixed phase is separated from the disordered phase by a second-order phase-transition which has many properties of the Bose-Einstein phase-transition observed in equilibrium. Various critical exponents are found.

Wetting and Capillary Condensation in Symmetric Polymer Blends: A comparison between Monte Carlo Simulations and Self-Consistent Field Calculations

Abstract: We present a quantitative comparison between extensive Monte Carlo simulations and self-consistent field calculations on the phase diagram and wetting behavior of a symmetric, binary (AB) polymer blend confined into a film. The flat walls attract one component via a short range interaction. The critical point of the confined blend is shifted to lower temperatures and higher concentrations of the component with the lower surface free energy. The binodals close the the critical point are flattened compared to the bulk and exhibit a convex curvature at intermediate temperatures -- a signature of the wetting transition in the semi-infinite system. Investigating the spectrum of capillary fluctuation of the interface bound to the wall, we find evidence for a position dependence of the interfacial tension. This goes along with a distortion of the interfacial profile from its bulk shape. Using an extended ensemble in which the monomer-wall interaction is a stochastic variable, we accurately measure the difference between the surface energies of the components, and determine the location of the wetting transition via the Young equation. The Flory-Huggins parameter at which the strong first order wetting transition occurs is independent of chain length and grows quadratically with the integrated wall-monomer interaction strength. We estimate the location of the prewetting line. The prewetting manifests itself in a triple point in the phase diagram of very thick films and causes spinodal dewetting of ultrathin layers slightly above the wetting transition. We investigate the early stage of dewetting via dynamic Monte Carlo simulations.

Critical Exponents of the N-vector model

Abstract: Recently the series for two RG functions (corresponding to the anomalous dimensions of the fields phi and phi^2) of the 3D phi^4 field theory have been extended to next order (seven loops) by Murray and Nickel. We examine here the influence of these additional terms on the estimates of critical exponents of the N-vector model, using some new ideas in the context of the Borel summation techniques. The estimates have slightly changed, but remain within errors of the previous evaluation. Exponents like eta (related to the field anomalous dimension), which were poorly determined in the previous evaluation of Le Guillou--Zinn-Justin, have seen their apparent errors significantly decrease. More importantly, perhaps, summation errors are better determined. The change in exponents affects the recently determined ratios of amplitudes and we report the corresponding new values. Finally, because an error has been discovered in the last order of the published epsilon=4-d expansions (order epsilon^5), we have also reanalyzed the determination of exponents from the epsilon-expansion. The conclusion is that the general agreement between epsilon-expansion and 3D series has improved with respect to Le Guillou--Zinn-Justin.

Development of Stresses in Cohesionless Poured Sand

Abstract: The pressure distribution beneath a conical sandpile, created by pouring sand from a point source onto a rough rigid support, shows a pronounced minimum below the apex (`the dip'). Recent work of the authors has attempted to explain this phenomenon by invoking local rules for stress propagation that depend on the local geometry, and hence on the construction history, of the medium. We discuss the fundamental difference between such approaches, which lead to hyperbolic differential equations, and elastoplastic models, for which the equations are elliptic within any elastic zones present .... This displacement field appears to be either ill-defined, or defined relative to a reference state whose physical existence is in doubt. Insofar as their predictions depend on physical factors unknown and outside experimental control, such elastoplastic models predict that the observations should be intrinsically irreproducible .... Our hyperbolic models are based instead on a physical picture of the material, in which (a) the load is supported by a skeletal network of force chains ("stress paths") whose geometry depends on construction history; (b) this network is `fragile' or marginally stable, in a sense that we define. .... We point out that our hyperbolic models can nonetheless be reconciled with elastoplastic ideas by taking the limit of an extremely anisotropic yield condition.

Computer Simulations of Pedestrian Dynamics and Trail Formation

Abstract: A simulation model for the dynamic behaviour of pedestrian crowds is mathematically formulated in terms of a social force model, that means, pedestrians behave in a way as if they would be subject to an acceleration force and to repulsive forces describing the reaction to borders and other pedestrians. The computational simulations presented yield many realistic results that can be compared with video films of pedestrian crowds. Especially, they show the self-organization of collective behavioural patterns. By assuming that pedestrians tend to choose routes that are frequently taken the above model can be extended to an active walker model of trail formation. The topological structure of the evolving trail network will depend on the disadvantage of building new trails and the durability of existing trails. Computer simulations of trail formation indicate to be a valuable tool for designing systems of ways which satisfy the needs of pedestrians best. An example is given for a non-directed trail network.

Stochastic Energetics of Non-uniform Temperature Systems

Abstract: We propose an energetic interpretation ofstochastic processes described by Langevin equations with non-uniform temperature. In order to avoid It\^{o}-Stratonovich dilemma, we start with a Kramers equation, and derive a Fokker-Plank equation by the renormalization group method. We give a proper definition of heat for the system. Based on our formulations,we analyze two examples, the Thomson effect and a Brownian motor which realizes the Carnot efficiency.

Fast spin dynamics algorithms for classical spin systems

Abstract: We have proposed new algorithms for the numerical integration of the equations of motion for classical spin systems. In close analogy to symplectic integrators for Hamiltonian equations of motion used in Molecular Dynamics these algorithms are based on the Suzuki-Trotter decomposition of exponential operators and unlike more commonly used algorithms exactly conserve spin length and, in special cases, energy. Using higher order decompositions we investigate integration schemes of up to fourth order and compare them to a well established fourth order predictor-corrector method. We demonstrate that these methods can be used with much larger time steps than the predictor-corrector method and thus may lead to a substantial speedup of computer simulations of the dynamical behavior of magnetic materials.

Condensate Fluctuations in Trapped Bose Gases: Canonical vs. Microcanonical Ensemble

Abstract: We study the fluctuation of the number of particles in ideal Bose-Einstein condensates, both within the canonical and the microcanonical ensemble. Employing the Mellin-Barnes transformation, we derive simple expressions that link the canonical number of condensate particles, its fluctuation, and the difference between canonical and microcanonical fluctuations to the poles of a Zeta function that is determined by the excited single-particle levels of the trapping potential. For the particular examples of one- and three-dimensional harmonic traps we explore the microcanonical statistics in detail, with the help of the saddle-point method. Emphasizing the close connection between the partition theory of integer numbers and the statistical mechanics of ideal Bosons in one-dimensional harmonic traps, and utilizing thermodynamical arguments, we also derive an accurate formula for the fluctuation of the number of summands that occur when a large integer is partitioned.

Short-time Critical Dynamics of the 3-Dimensional Ising Model

Abstract: Comprehensive Monte Carlo simulations of the short-time dynamic behaviour are reported for the three-dimensional Ising model at criticality. Besides the exponent $\theta$ of the critical initial increase and the dynamic exponent $z$, the static critical exponents $ u$ and $\beta$ as well as the critical temperature are determined from the power-law scaling behaviour of observables at the beginning of the time evolution. States of very high temperature as well as of zero temperature are used as initial states for the simulations.

Lagrangian method for multiple correlations in passive scalar advection

Abstract: A Lagrangian method is introduced for calculating simultaneous n-point correlations of a passive scalar advected by a random velocity field, with random forcing and finite molecular diffusivity kappa. The method, which is here presented in detail, is particularly well suited for studying the kappa tending to 0 limit when the velocity field is not smooth. Efficient Monte Carlo simulations based on this method are applied to the Kraichnan model of passive scalar and lead to accurate determinations of the anomalous intermittency corrections in the fourth-order structure function as a function of the scaling exponent xi of the velocity field in two and three dimensions. Anomalous corrections are found to vanish in the limits xi tending to 0 and xi tending to 2, as predicted by perturbation theory.

The XX--model with boundaries. Part I: Diagonalization of the finite chain

Abstract: This is the first of three papers dealing with the XX finite quantum chain with arbitrary, not necessarily hermitian, boundary terms. This extends previous work where the periodic or diagonal boundary terms were considered. In order to find the spectrum and wave-functions an auxiliary quantum chain is examined which is quadratic in fermionic creation and annihilation operators and hence diagonalizable. The secular equation is in general complicated but several cases were found when it can be solved analytically. For these cases the ground-state energies are given. The appearance of boundary states is also discussed and in view to the applications considered in the next papers, the one and two-point functions are expressed in terms of Pfaffians.