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

Generalized-ensemble Monte carlo method for systems with rough energy landscape

We present a novel Monte Carlo algorithm which enhances equilibrization of low-temperature simulations and allows sampling of configurations over a large range of energies. The method is based on a non-Boltzmann probability weight factor and is another version of the so-called generalized-ensemble techniques. The effectiveness of the new approach is demonstrated for the system of a small peptide, an example of the frustrated system with a rugged energy landscape.Comment: Latex; ps-files include

Multi-Overlap Simulations for Transitions between Reference Configurations

We introduce a new procedure to construct weight factors, which flatten the probability density of the overlap with respect to some pre-defined reference configuration. This allows one to overcome free energy barriers in the overlap variable. Subsequently, we generalize the approach to deal with the overlaps with respect to two reference configurations so that transitions between them are induced. We illustrate our approach by simulations of the brainpeptide Met-enkephalin with the ECEPP/2 energy function using the global-energy-minimum and the second lowest-energy states as reference configurations. The free energy is obtained as functions of the dihedral and the root-mean-square distances from these two configurations. The latter allows one to identify the transition state and to estimate its associated free energy barrier.Comment: 12 pages, (RevTeX), 14 figures, Phys. Rev. E, submitte

A comparison of extremal optimization with flat-histogram dynamics for finding spin-glass ground states

We compare the performance of extremal optimization (EO), flat-histogram and equal-hit algorithms for finding spin-glass ground states. The first-passage-times to a ground state are computed. At optimal parameter of tau=1.15, EO outperforms other methods for small system sizes, but equal-hit algorithm is competitive to EO, particularly for large systems. Flat-histogram and equal-hit algorithms offer additional advantage that they can be used for equilibrium thermodynamic calculations. We also propose a method to turn EO into a useful algorithm for equilibrium calculations. Keywords: extremal optimization. flat-histogram algorithm, equal-hit algorithm, spin-glass model, ground state.Comment: 10 LaTeX pages, 2 figure

Transition Matrix Monte Carlo Reweighting and Dynamics

We study an induced dynamics in the space of energy of single-spin-flip Monte Carlo algorithm. The method gives an efficient reweighting technique. This dynamics is shown to have relaxation times proportional to the specific heat. Thus, it is plausible for a logarithmic factor in the correlation time of the standard 2D Ising local dynamics.Comment: RevTeX, 5 pages, 3 figure

Spin glass overlap barriers in three and four dimensions

For the Edwards-Anderson Ising spin-glass model in three and four dimensions (3d and 4d) we have performed high statistics Monte Carlo calculations of those free-energy barriers $F^q_B$ which are visible in the probability density $P_J(q)$ of the Parisi overlap parameter $q$. The calculations rely on the recently introduced multi-overlap algorithm. In both dimensions, within the limits of lattice sizes investigated, these barriers are found to be non-self-averaging and the same is true for the autocorrelation times of our algorithm. Further, we present evidence that barriers hidden in $q$ dominate the canonical autocorrelation times.Comment: 20 pages, Latex, 12 Postscript figures, revised version to appear in Phys. Rev.

Potentials of Mean Force for Protein Structure Prediction Vindicated, Formalized and Generalized

Understanding protein structure is of crucial importance in science, medicine and biotechnology. For about two decades, knowledge based potentials based on pairwise distances -- so-called "potentials of mean force" (PMFs) -- have been center stage in the prediction and design of protein structure and the simulation of protein folding. However, the validity, scope and limitations of these potentials are still vigorously debated and disputed, and the optimal choice of the reference state -- a necessary component of these potentials -- is an unsolved problem. PMFs are loosely justified by analogy to the reversible work theorem in statistical physics, or by a statistical argument based on a likelihood function. Both justifications are insightful but leave many questions unanswered. Here, we show for the first time that PMFs can be seen as approximations to quantities that do have a rigorous probabilistic justification: they naturally arise when probability distributions over different features of proteins need to be combined. We call these quantities reference ratio distributions deriving from the application of the reference ratio method. This new view is not only of theoretical relevance, but leads to many insights that are of direct practical use: the reference state is uniquely defined and does not require external physical insights; the approach can be generalized beyond pairwise distances to arbitrary features of protein structure; and it becomes clear for which purposes the use of these quantities is justified. We illustrate these insights with two applications, involving the radius of gyration and hydrogen bonding. In the latter case, we also show how the reference ratio method can be iteratively applied to sculpt an energy funnel. Our results considerably increase the understanding and scope of energy functions derived from known biomolecular structures

Reexamination of the long-range Potts model: a multicanonical approach

We investigate the critical behavior of the one-dimensional q-state Potts model with long-range (LR) interaction $1/r^{d+\sigma}$, using a multicanonical algorithm. The recursion scheme initially proposed by Berg is improved so as to make it suitable for a large class of LR models with unequally spaced energy levels. The choice of an efficient predictor and a reliable convergence criterion is discussed. We obtain transition temperatures in the first-order regime which are in far better agreement with mean-field predictions than in previous Monte Carlo studies. By relying on the location of spinodal points and resorting to scaling arguments, we determine the threshold value $\sigma_c(q)$ separating the first- and second-order regimes to two-digit precision within the range $3 \leq q \leq 9$. We offer convincing numerical evidence supporting $\sigma_c(q)Comment: 18 pages, 18 figure

Parallel Excluded Volume Tempering for Polymer Melts

We have developed a technique to accelerate the acquisition of effectively uncorrelated configurations for off-lattice models of dense polymer melts which makes use of both parallel tempering and large scale Monte Carlo moves. The method is based upon simulating a set of systems in parallel, each of which has a slightly different repulsive core potential, such that a thermodynamic path from full excluded volume to an ideal gas of random walks is generated. While each system is run with standard stochastic dynamics, resulting in an NVT ensemble, we implement the parallel tempering through stochastic swaps between the configurations of adjacent potentials, and the large scale Monte Carlo moves through attempted pivot and translation moves which reach a realistic acceptance probability as the limit of the ideal gas of random walks is approached. Compared to pure stochastic dynamics, this results in an increased efficiency even for a system of chains as short as $N = 60$ monomers, however at this chain length the large scale Monte Carlo moves were ineffective. For even longer chains the speedup becomes substantial, as observed from preliminary data for $N = 200$

The Impact of Global Warming and Anoxia on Marine Benthic Community Dynamics: an Example from the Toarcian (Early Jurassic)

The Pliensbachian-Toarcian (Early Jurassic) fossil record is an archive of natural data of benthic community response to global warming and marine long-term hypoxia and anoxia. In the early Toarcian mean temperatures increased by the same order of magnitude as that predicted for the near future; laminated, organic-rich, black shales were deposited in many shallow water epicontinental basins; and a biotic crisis occurred in the marine realm, with the extinction of approximately 5% of families and 26% of genera. High-resolution quantitative abundance data of benthic invertebrates were collected from the Cleveland Basin (North Yorkshire, UK), and analysed with multivariate statistical methods to detect how the fauna responded to environmental changes during the early Toarcian. Twelve biofacies were identified. Their changes through time closely resemble the pattern of faunal degradation and recovery observed in modern habitats affected by anoxia. All four successional stages of community structure recorded in modern studies are recognised in the fossil data (i.e. Stage III: climax; II: transitional; I: pioneer; 0: highly disturbed). Two main faunal turnover events occurred: (i) at the onset of anoxia, with the extinction of most benthic species and the survival of a few adapted to thrive in low-oxygen conditions (Stages I to 0) and (ii) in the recovery, when newly evolved species colonized the re-oxygenated soft sediments and the path of recovery did not retrace of pattern of ecological degradation (Stages I to II). The ordination of samples coupled with sedimentological and palaeotemperature proxy data indicate that the onset of anoxia and the extinction horizon coincide with both a rise in temperature and sea level. Our study of how faunal associations co-vary with long and short term sea level and temperature changes has implications for predicting the long-term effects of “dead zones” in modern oceans