Yang-Lee zeros and the helix-coil transition in a continuum model of polyalanine

We calculate the Yang-Lee zeros for characteristic temperatures of the helix-coil transition in a continuum model of polyalanine. The distribution of these zeros differs from predictions of the Zimm-Bragg theory and supports recent claims that polyalanine exhibits a true phase transition. New estimates for critical exponents are presented and the relation of our results to the Lee-Yang theorem is discussed.Comment: 15 pages and 5 figure

Microcanonical thermostatistics analysis without histograms: cumulative distribution and Bayesian approaches

Microcanonical thermostatistics analysis has become an important tool to reveal essential aspects of phase transitions in complex systems. An efficient way to estimate the microcanonical inverse temperature $\beta(E)$ and the microcanonical entropy $S(E)$ is achieved with the statistical temperature weighted histogram analysis method (ST-WHAM). The strength of this method lies on its flexibility, as it can be used to analyse data produced by algorithms with generalised sampling weights. However, for any sampling weight, ST-WHAM requires the calculation of derivatives of energy histograms $H(E)$, which leads to non-trivial and tedious binning tasks for models with continuous energy spectrum such as those for biomolecular and colloidal systems. Here, we discuss two alternative methods that avoid the need for such energy binning to obtain continuous estimates for $H(E)$ in order to evaluate $\beta(E)$ by using ST-WHAM: (i) a series expansion to estimate probability densities from the empirical cumulative distribution function (CDF), and (ii) a Bayesian approach to model this CDF. Comparison with a simple linear regression method is also carried out. The performance of these approaches is evaluated considering coarse-grained protein models for folding and peptide aggregation.Comment: 9 pages, 11 figure

Structural transitions in biomolecules - a numerical comparison of two approaches for the study of phase transitions in small systems

We compare two recently proposed methods for the characterization of phase transitions in small systems. The usefulness of these techniques is evaluated for the case of structural transition in alanine-based peptides.Comment: Accepted for publication in Int. J. Mol. Sci., to appear in a special issue devoted to R.S. Berr

Helix Formation and Folding in an Artificial Peptide

We study the relation between $\alpha$-helix formation and folding for a simple artificial peptide, Ala$_{10}$-Gly$_5$-Ala$_{10}$. Our data rely on multicanonical Monte Carlo simulations where the interactions among all atoms are taken into account. The free-energy landscape of the peptide is evaluated for various temperatures. Our data indicate that folding of this peptide is a two-step process: in a first step two $\alpha$-helices are formed which afterwards re-arrange themselves into a U-like structure.Comment: 15 pages, with 9 eps figure

Spectral Density Study of the SU(3) Deconfining Phase Transition

We present spectral density reweighting techniques adapted to the analysis of a time series of data with a continuous range of allowed values. In a first application we analyze action and Polyakov line data from a Monte Carlo simulation on $L_t L^3 (L_t=2,4)$ lattices for the SU(3) deconfining phase transition. We calculate partition function zeros, as well as maxima of the specific heat and of the order parameter susceptibility. Details and warnings are given concerning i) autocorrelations in computer time and ii) a reliable extraction of partition function zeros. The finite size scaling analysis of these data leads to precise results for the critical couplings $\beta_c$, for the critical exponent $\nu$ and for the latent heat $\triangle s$. In both cases ($L_t=2$ and 4), the first order nature of the transition is substantiated