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

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

Multicanonical Spin Glass Simulations

We report a Monte Carlo simulation of the $2D$ Edwards-Anderson spin glass model within the recently introduced multicanonical ensemble. Replica on lattices of size $L^2$ up to $L=48$ are investigated. Once a true groundstate is found, we are able to give a lower bound on the number of statistically independent groundstates sampled. Temperature dependence of the energy, entropy and other quantities of interest are easily calculable. In particular we report the groundstate results. Computations involving the spin glass order parameter are more tedious. Our data indicate that the large $L$ increase of the ergodicity time is reduced to an approximately $V^3$ power law. Altogether the results suggest that the multicanonical ensemble improves the situation of simulations for spin glasses and other systems which have to cope with similar problems of conflicting constraints.Comment: 24 page

Measuring Workload Differences Between Short-term Memory and Long-term Memory Scenarios in a Simulated Flight Environment

Four highly experienced Air Force pilots each flew four simulated flight scenarios. Two scenarios required a great deal of aircraft maneuvering. The other two scenarios involved less maneuvering, but required remembering a number of items. All scenarios were designed to be equaly challenging. Pilot's Subjective Ratings for Activity-level, Complexity, Difficulty, Stress, and Workload were higher for the manuevering scenarios than the memory scenarios. At a moderate workload level, keeping the pilots active resulted in better aircraft control. When required to monitor and remember items, aircraft control tended to decrease. Pilots tended to weigh information about the spatial positioning and performance of their aircraft more heavily than other items

Structure of the Energy Landscape of Short Peptides

We have simulated, as a showcase, the pentapeptide Met-enkephalin (Tyr-Gly-Gly-Phe-Met) to visualize the energy landscape and investigate the conformational coverage by the multicanonical method. We have obtained a three-dimensional topographic picture of the whole energy landscape by plotting the histogram with respect to energy(temperature) and the order parameter, which gives the degree of resemblance of any created conformation with the global energy minimum (GEM).Comment: 17 pages, 4 figure

Normalized entropy density of the 3D 3-state Potts model

Using a multicanonical Metropolis algorithm we have performed Monte Carlo simulations of the 3D 3-state Potts model on $L^3$ lattices with L=20, 30, 40, 50. Covering a range of inverse temperatures from $\beta_{\min}=0$ to $\beta_{\max}=0.33$ we calculated the infinite volume limit of the entropy density $s(\beta)$ with its normalization obtained from $s(0)=\ln 3$. At the transition temperature the entropy and energy endpoints in the ordered and disordered phase are estimated employing a novel reweighting procedure. We also evaluate the transition temperature and the order-disorder interface tension. The latter estimate increases when capillary waves are taken into account.Comment: 5 pages, 4 figure

Recent Results of Multimagnetical Simulations of the Ising Model

To investigate order-order interfaces, we perform multimagnetical Monte Carlo simulations of the $2D$ and $3D$ Ising model. Stringent tests of the numerical methods are performed by reproducing with high precision exact $2D$ results. In the physically more interesting $3D$ case we estimate the amplitude $F^s_0$ of the critical interfacial tension.Comment: talk presented at the workshop "Dynamics of First Order Phase Transitions", Juelich June 1-3; FSU-SCRI-92C-87 preprint; 7 pages; sorry no figures; needs vanilla.st

Biased Metropolis-Heat-Bath Algorithm for Fundamental-Adjoint SU(2) Lattice Gauge Theory

For SU(2) lattice gauge theory with the fundamental-adjoint action an efficient heat-bath algorithm is not known so that one had to rely on Metropolis simulations supplemented by overrelaxation. Implementing a novel biased Metropolis-heat-bath algorithm for this model, we find improvement factors in the range 1.45 to 2.06 over conventionally optimized Metropolis simulations. If one optimizes further with respect to additional overrelaxation sweeps, the improvement factors are found in the range 1.3 to 1.8.Comment: 4 pages, 2 figures; minor changes and one reference added; accepted for publication in PR

Deep drawing simulation of Tailored Blanks

Tailored blanks are increasingly used in the automotive industry. A tailored blank consists of different metal parts, which are joined by a welding process. These metal parts usually have different material properties. Hence, the main advantage of using a tailored blank is to provide the right material properties at specific parts of the blank. The movement of the weld during forming is extremely important. Unwanted weld displacement can cause damage to both the product and the tool. This depends mainly on the original weld position and the process parameters. However experimental determination of the optimum weld position is quite expensive. Therefore a numerical tool has been developed for simulations of tailored blank forming. The Finite Element Code Dieka is used for the deep drawing simulations of some geometrically simple products. The results have been validated by comparing them with experimental data and show a satisfactory correlation