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

Athermal Phase Separation of Self-Propelled Particles with no Alignment

We study numerically and analytically a model of self-propelled polar disks on a substrate in two dimensions. The particles interact via isotropic repulsive forces and are subject to rotational noise, but there is no aligning interaction. As a result, the system does not exhibit an ordered state. The isotropic fluid phase separates well below close packing and exhibits the large number fluctuations and clustering found ubiquitously in active systems. Our work shows that this behavior is a generic property of systems that are driven out of equilibrium locally, as for instance by self propulsion.Comment: 5 pages, 4 figure

Wind speed statistics for Goldstone, California, anemometer sites

An exploratory wind survey at an antenna complex was summarized statistically for application to future windmill designs. Data were collected at six locations from a total of 10 anemometers. Statistics include means, standard deviations, cubes, pattern factors, correlation coefficients, and exponents for power law profile of wind speed. Curves presented include: mean monthly wind speeds, moving averages, and diurnal variation patterns. It is concluded that three of the locations have sufficiently strong winds to justify consideration for windmill sites

Entropy-based analysis of the number partitioning problem

In this paper we apply the multicanonical method of statistical physics on the number-partitioning problem (NPP). This problem is a basic NP-hard problem from computer science, and can be formulated as a spin-glass problem. We compute the spectral degeneracy, which gives us information about the number of solutions for a given cost $E$ and cardinality $m$. We also study an extension of this problem for $Q$ partitions. We show that a fundamental difference on the spectral degeneracy of the generalized ($Q>2$) NPP exists, which could explain why it is so difficult to find good solutions for this case. The information obtained with the multicanonical method can be very useful on the construction of new algorithms.Comment: 6 pages, 4 figure

Run-and-tumble particles with hydrodynamics: sedimentation, trapping and upstream swimming

We simulate by lattice Boltzmann the nonequilibrium steady states of run-and-tumble particles (inspired by a minimal model of bacteria), interacting by far-field hydrodynamics, subject to confinement. Under gravity, hydrodynamic interactions barely perturb the steady state found without them, but for particles in a harmonic trap such a state is quite changed if the run length is larger than the confinement length: a self-assembled pump is formed. Particles likewise confined in a narrow channel show a generic upstream flux in Poiseuille flow: chiral swimming is not required

Statistical mechanics of random two-player games

Using methods from the statistical mechanics of disordered systems we analyze the properties of bimatrix games with random payoffs in the limit where the number of pure strategies of each player tends to infinity. We analytically calculate quantities such as the number of equilibrium points, the expected payoff, and the fraction of strategies played with non-zero probability as a function of the correlation between the payoff matrices of both players and compare the results with numerical simulations.Comment: 16 pages, 6 figures, for further information see http://itp.nat.uni-magdeburg.de/~jberg/games.htm

Modeling and parameter uncertainties for aircraft flight control system design

Values of plant dynamic uncertainties for some recent aircraft design and development programs are given. Histories of pertinent aerodynamic, inertial, and structural parameter variations are given for a period of time from program initiation to aircraft certification. These data can be used as typical of future vehicles so that control system design concepts are evaluated with due consideration to their sensitivity to uncertainties in plant dynamics

Interstellar extinction and the distribution of stellar populations in the direction of the ultra-deep Chandra Galactic field

We studied the stellar population in the central 6.6x6.6arcmin,region of the ultra-deep (1Msec) Chandra Galactic field - the "Chandra bulge field" (CBF) approximately 1.5 degrees away from the Galactic Center - using the Hubble Space Telescope ACS/WFC blue (F435W) and red (F625W) images. We mainly focus on the behavior of red clump giants - a distinct stellar population, which is known to have an essentially constant intrinsic luminosity and color. By studying the variation in the position of the red clump giants on a spatially resolved color-magnitude diagram, we confirm the anomalous total-to-selective extinction ratio, as reported in previous work for other Galactic bulge fields. We show that the interstellar extinction in this area is = 4 on average, but varies significantly between ~3-5 on angular scales as small as 1 arcminute. Using the distribution of red clump giants in an extinction-corrected color-magnitude diagram, we constrain the shape of a stellar-mass distribution model in the direction of this ultra-deep Chandra field, which will be used in a future analysis of the population of X-ray sources. We also show that the adopted model for the stellar density distribution predicts an infrared surface brightness in the direction of the "Chandra bulge field" in good agreement (i.e. within ~15%) with the actual measurements derived from the Spitzer/IRAC observations.Comment: 9 pages, 9 figures. Accepted for publication in A&