Thermodynamics of the coil to frozen globule transition in heteropolymers
|
|
- Benjamin Fitzgerald
- 6 years ago
- Views:
Transcription
1 Thermodynamics of the coil to frozen globule transition in heteropolymers Vijay S. Pande Department of Physics, University of California at Berkeley, Berkeley, California 94720, and the Department of Physics and Center for Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts Alexander Yu. Grosberg a) and Toyoichi Tanaka Department of Physics and Center for Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts Received 18 December 1996; accepted 26 June 1997 Recent analytic theories and computer simulations of heteropolymers have centered on the freezing transition of globular heteropolymers. We present a simple analytic theory to describe the coil to globule collapse in heteropolymers and compare this to the computer simulation of the exhaustive enumeration of all 18-mer cubic lattice polymer conformations. We find that the collapse transition from coil to frozen globule can either be first or second order. The relevance to protein folding is also discussed American Institute of Physics. S I. INTRODUCTION In the early days of protein studies, the denaturation of globular proteins was thought of as a globule decollapse transition, and this was one of the reasons that brought much attention to the homopolymer coil to globule transition. By now, this later phenomenon is well understood in many respects. In particular, thermodynamics of this transition has been studied theoretically in great detail 1,2 and the transition has been also clearly observed experimentally despite the difficulties with phase segregation in poor solvent conditions. 3,4 Except for the simplest universal picture of the collapse of a flexible uncharged homopolymer, several other model systems have been examined, including annealed heteropolymers, 5 charged pyroelectrolytes and polyampholytes, 6 and polymeric mesogenes. 7 The corresponding phase diagrams are pretty rich and in some respects indicate certain similarities with properties of proteins. There are even assumptions albeit disputable ones in the literature 8 that the molten globule phase of proteins 9 can be modelled on the level of a homopolymer. Nevertheless, it is clear that the final native conformation of each protein is encoded in its sequence and thus it is well understood that protein folding requires a more specific heteropolymer aspect, which is believed to be freezing We think it is fair to say that freezing in general is understood much less clearly at present than collapse, as only the simplest models of freezing have been examined. Even more than that, these two aspects collapse and freezing are not well matched in the emerging general picture of heteropolymers. There were several works in the literature examining phase diagrams that include coil, liquid-like random globule similar to a homopolymer globule, and frozen globule phases of a heteropolymer chain. 13,14 After the submission of the present work, a new study was published 15 which contains several new important methodical innovations and pro- a On leave from the Institute of Chemical Physics, Russian Academy of Sciences, Moscow vides for new insights into the relationships of three possible phases of a heteropolymer chain. One of the important achievements of the work 15 is that it goes beyond the uncorrelated energy landscape or REM 16 model and thus continues the truck toward understanding the kinetics of folding. Despite the important achievements of the previous works, we would like with the present work to once again bring the attention to the problem of a phase diagram of a heteropolymer chain. In doing so, our main motivation is to take a closer look at the phase transition between frozen globule and coil state and to give arguments showing that this transition may be, and in many realistic cases is, of the first order. This fact seems to be missing in the previous studies, while we view it as important for the overall understanding of the possible scenarios of folding. We give simple analytical arguments to support our claim and we also demonstrate a solid numerical evidence for that. The later one is based on complete enumeration of all, not just maximally compact, conformations of an 18-mer on the cubic lattice. This enumeration has not been performed before, and it is clearly advantageous for our purposes compared to the previous model 14 of a two dimensional 16-mer. The possibility of the first order phase transition is fundamentally due to the obvious fact that the nature of the phases of a heteropolymer chain must depend on its sequence. It is well understood at present that the ensembles of possible sequences can be described much like Gibbs ensembles in the statistical mechanics, in terms of either a microcanonical ensemble, in which one specifies the energy E NS of the native state and thus knows how well this native state is optimized energetically, or in terms of a canonical ensemble, in which one controls E NS in terms of some effective temperature T d. There are several constructive models that show how to control native state energy using a canonical ensemble with temperature T d, including sequence selection 17 and Imprinting. 12 We shall show that phase transitions in the chain of a heteropolymer depend indeed on the type of sequence, that is, on T d. Our paper is organized as follows. We begin with an 5118 J. Chem. Phys. 107 (13), 1 October /97/107(13)/5118/7/$ American Institute of Physics
2 Pande, Grosberg, and Tanaka: Coil to globule transition 5119 analytic treatment and then compare those results with the computer simulation. We conclude with a discussion of the relevance of these results to proteins and other heteorpolymeric systems. II. ANALYTIC TREATMENT A. The model What are the relevant order parameters in our study? As to the coil to globule transition, in which the system goes from a phase of swollen, random walk type conformations perturbed only by excluded volume to one in which it is dense and collapsed, it can be described by chain size or density; for each conformation one defines (x) I (x r I ), where r I is the position vector of the monomer I in the conformation. As long as the thermodynamic limit is concerned, is zero in the coil phase and non-vanishing in the globule. A more delicate order parameter is necessary to consider the freezing transition, in which the system goes from a phase consisting of an exponential number of conformations O (e N ) to one in which only O (1) states dominate. As we expect the ground state conformation to be unique, a good order parameter to examine arbitrary conformation of the chain that undergoes a freezing transition is the overlap between and the ground state : Q I J r I r J r I r J, where (r) is a function that is concentrated on nearest neighboring space points on a lattice, for lattice spacing a, (a) 1 and (r a) 0], and thus (r I r J ) is nonvanishing only when monomers I and J in conformation are neighbors in contact. It is even more convenient to define, for a conformation, its fractional overlap to the ground state, q Q /Q max, 2 where Q max is the maximum number of contacts possible we assume that is maximally compact and thus Q Q max ). Obviously, q ranges from 0 to 1. At first, these two order parameters, and Q, may seem to be completely unrelated, but as dense chains have the possibility for forming more contacts, we expect that there should be some relationship between the density and the number of contacts of a given conformation which is the same as the overlap of a conformation with itself Q I J (r I r J )]. To see this, one has to introduce the coarse grained density (x) (1/z) I (x r I ), where z is the number of nearest neighbors on which is localized for each conformation, and define also spatial averages as dr/v. Noting that Nz Q max, one arrives at q Q 2 Q 2. 3 max Thus, the average density of the chain in some conformation is almost the same as fractional overlap of this conformation 1 with itself, and as density increases, the approximation that q improves. Note that q and q are proportional to the total number of bonds between monomers in the conformation and the number of bonds that the given conformation has in common with the ground state conformation. B. Free energy We will next calculate the free energy in terms of order parameters q and q q. In our calculations, we employ the most general Hamiltonian of short-range pairwise interactions, 12 H I J B si,s J r I r J, where s I 1,...,n is the species of monomer number I, n is the number of species, B ij is the matrix of species species interactions. We implicitly assume here that position vectors r I are such that the condition of chain connectivity is met, but do not assume maximally compact conformations. Formally, the free energy of interest can be written in a general form in terms of the Hamiltonian 4 : F seq,q T ln e H /T q q, 5 where the subscript seq indicates that the free energy is written for a given sequence, and summation is performed over all conformations. The dependence on the sequence enters both through the Hamiltonian and through the ground state conformation for the given sequence. It is technically cumbersome to treat the free energy 5 in the general form. What we shall do instead is we shall resort to an interpolation expression which can be derived using the following two arguments. For a homopolymer, when q-dependence is not in question, the free energy per monomer can be described in terms of a truncated virial series: 1 N F B C 2, 6 where B and C are, respectively, two- and three-body virial coefficients. In the simplest case, B in the equation 6 is related to B ij in the equation 4 by B b 3 z exp( B ij ) 1, where z and b are the coordination number and the mesh size of the underlying lattice. For simplicity, we shall take b 1 in what follows, such that volumes and densities are unitless. For a heteropolymer whose density is evenly distributed in space and cannot fluctuate when dependence is trivial, the following expression has been derived elsewhere: 18 1 B2 F q q B2 N T p T 1 q Tq s 3 2 ln q T q ln q 1 q ln 1 q. 7 This expression is detailed in the Appendix. We make now the assumption that the heterogeneity of monomeric interactions affects only the two-body interaction 4
3 5120 Pande, Grosberg, and Tanaka: Coil to globule transition terms, combine equations 6 and 7, and thus arrive at the following phenomenological interpolation which is central for the present work: 1 N F,q B B2 q B2 T p T 1 q C 2 Tq s 3 2 ln q T q ln q 1 q ln 1 q. 8 To understand the approximation behind the expression 8, we note that we have omitted the loop-related entropy part in the homopolymeric free-energy 6. This may seem strange, as we did take into account a seemingly analogous loop factor in the heteropolymeric part 7. We think this procedure is legitimate, because homopolymer loops are known to contribute only to the surface energy of the globule, i.e., to the terms of the order of N 2/3, while all terms in 8 are of order N. Before proceeding to describe the resulting phase diagram, we stop to comment on the differences between our theory and previous works. 13,14 Our model goes beyond previous models in that corrections to the Random Energy Model are included see the Appendix for details. The reader might be concerned that our formalism does not utilize complicated mathematical formalisms previously employed, such as replica field theory, and thus the current theory is an oversimplification of previous work and therefore less powerful. In fact, by breaking away from the replica formalism, we are able to introduce the relevant physical ingredients and advance beyond these previous works. C. Phase diagram An inspection of our main equation for the free energy 8 indicates that the q-dependence allows for a minima at either q 1 orq 0, i.e., with either complete dominance of the native state, or with the mixture of the conformations statistically unrelated to the native state. This result coincides with that of an uncorrelated energy landscape, or Random Energy Model REM. Thus, this theory can be considered to be a simple derivation of REM for heteropolymer freezing. Note also that the q 1 minimum exists if and only if the corresponding is not too small, meaning that this minimum can be realized only for the globular polymer; obviously, it corresponds to the state of the frozen globule. By contrast, the q 0 state can be realized both for the coil ( 0) and for the globule, which in this case is the random globule. Thus, we arrive at the three-region phase diagram which includes a coil phase and two globular phases random and native frozen globules ; we call these phases C, RG, and FG, respectively. We now comment on the nature of the boundaries between the three phases: Random to frozen globule When the solvent is sufficiently poor overall (B is sufficiently negative, the transition in q will not be accompanied by a transition in. This occurs because while there would be an entropic benefit to have a random globule with smaller, the energetic penalty due to a negative B 0 overrules this. Thus, we compare the free energy with max, but with q 0 and q 1. This yields the following condition for a transition from a random to a frozen globule: RG FG: s B2 T 1 T 1 9 T p, which agrees not surprisingly with previous works describing the transition from the random to the frozen globule state. 12 Random globule to coil As q 0 in both of these phases, the only heteropolymeric complication comes from an effective second virial coefficient B eff B B 2 /T. In terms of B eff, the transition to this order is exactly the same as the homopolymer globule to coil transition, with densities B eff /2C in the globular phase and 0 in the coil phase. To this approximation, as we neglect loop-related surface corrections to the homopolymer free energy, 1 the transition happens at C RG in principle : B B T We have taken into account that the coil free energy vanishes. However, this is a second order transition even if surface energy is taken into account, it remains a very weak transition 1, and for this reason it has a rather large width. In practical terms, it causes considerable complications for computer simulation, as it obscures the position of the transition point. Moreover, in reality the transition region borders the place where globule reaches its maximal possible dense packing density max. Thus, practically the transition is seen at B eff /2C max, that is C RG practically : B B2 T 2C max. 11 Frozen globule to coil As we increase the quality of the solvent, eventually the energetic penalty for decreasing will not be enough and changes as well as q upon denaturation. Thus, the two phases in question correspond to q 0, 0 C and q 1, max FG. The transition occurs at where the free energies are equal, that is, at B FG C: s max T B 2 C 2 max TT max p T. 12 Thus, if we are interested in denaturation, that is, in the destruction of the frozen globule state, we have to ask whether the FG is transformed into C or RG. These two regimes cross into one another where all the lines 9 and 12 and 11 meet, that is, at the corresponding triple point. To bring our results into a more obvious form that is more convenient for comparisons to simulations, we have to choose control parameters, as generally there may be too many of them. In particular, if one wants to speak in terms of a temperature change, then one must take into account the temperature dependence of both B and B 2. This makes equations rather cumbersome. We choose to get rid of temperatures T and T p by redefinitions B B/T,
4 Pande, Grosberg, and Tanaka: Coil to globule transition 5121 FIG. 1. Phase diagram as predicted theoretically. There are three phases: 1 Coil: expanded, random walk-like configurations and thus this phase has maximum entropy. 2 Random Globule: conformations are globular condensed, high density, but an exponentially large number of conformations are still relevant. 3 Frozen Globule: target conformation dominates equilibrium. Phase transitions in density between the random globule and the coil are second order, whereas the density transition between the frozen globule and the coil are first order. Also, the transition in the overlap with the ground state is first order and occurs at boundaries of the frozen globule phase. B p B/T p and constructing phase diagram in terms of B, B with C and B p considered as parameters. Such a diagram is displayed in Fig. 1. This diagram can be described as follows. For B 0 B c 0, there is a line of first order transitions in q at B B fp and a second transition from maximally compact states at B B m. These two lines meet at a tricritical point (B c 0, B fp ). For B 0 B c 0, there is a single line of first order transitions in both q and. As the coupling constant of and q is B p, better designed sequences will exhibit more pronounced jumps in in the frozen globule to coil transition. Summarizing the results of our analytic model, we find three phases coil, random globule, and frozen globule. This is not new and has been previously described. 13,14 However, we predict that well designed sequences will have a first order transition in density between the coil and frozen globule phases. This first order transition has been neither predicted and nor seen in computer simulations. In the next section, we will verify these analytic predictions by exactly calculating the partition function for a short heteropolymer chain computationally. III. COMPUTER SIMULATIONS A. The model We now compare the results of the theory to computer simulations. In the study of heteropolymers, one must be careful as a few conformations can dominate the partition function. Indeed, this phenomenon is freezing and is one of the most interesting aspects of heteropolymers and the central aspect of this paper. To avoid computational pitfalls, the exhaustive enumeration of all conformations is performed in order to exactly calculate the partition function. Exhaustive enumeration for all maximally compact chains can only be performed for relatively short lengths 19 and the 27-mer on the has become a common model system. Unfortunately, the enumeration of all not just maximally compact 27-mer conformations is out of reach of current computational power. One approach is to consider only two dimensional conformations for example, this approach was used in Ref. 14, in which all 16-mer conformations in two-dimensional 2D were enumerated. However, we have been able to enumerate all 5,577,317, mer conformations. One may be suspicious that the 18-mer is simply too small to give representative information about heteropolymers. To test this, we examined several heteropolymer properties of 18-mers and compared these results to that of 27- mers. Upon examining the density of states, and nature of the freezing transition, 18-mers were not qualitatively different from 27-mers. Also, one must choose what types of interactions between monomers one wishes to model. As the theory includes terms only to lowest order, it cannot discriminate between differences in interactions and considers only the variance of the interaction matrix. Thus, differences in interactions modeled computationally cannot be described analytically. To make the best connection with the analytic theory, we designed sequences with a variety of degrees of optimization for the Independent Interaction Model IIM. We found that the range of ground state energies from the best optimized sequences to that of random sequences ranged from 12 to 2. In particular, we have examined 6 values of ground state energies. For comparison, we also designed sequences with 2 and 5 letter Potts model and found no major qualitative differences. B. Results The result of our enumeration for a given sequence and set of interactions is the density of states P(E,Q,q ), i.e., how many conformations have energy E, total number of contacts Q, and percent of contacts q in common with. From the density of states, we can calculate the free energy as a function of B, Q, and q. In Fig. 2, we examine the free energy for a poorly designed sequence in a poor solvent. We see that there is first a first order freezing transition i.e., a jump in q ) and then a second order globule to coil transition as we lower B. In a good solvent, the transition in q and Q go together and are both first order Fig. 3. In general, we find strong qualitative agreement with the theoretical prediction. First, consider the phase diagram of an optimally designed sequence Fig. 4 : we see that all elements of the analytic phase diagram are present. It is not surprising that our lowest order approximation to the free energy in the high temperature expansion does not completely capture all quantitative aspects of the curves, but certainly the nature of the transitions and the general shape of
5 5122 Pande, Grosberg, and Tanaka: Coil to globule transition FIG. 2. Free energy for a poorly designed sequence in a poor solvent. We see a first order transition from frozen globule a to random globule b and then a second order transition to the coil state c dark areas have lower free energy. Frame d shows F(q ) for these three regimes. the curves agree well with theory. Unlike the 2D case, 14 we computationally find a coil to frozen globule transition with one conformation dominating equilibrium for B 0 0, even for random sequences. Note that the assumption of constant density does not seem to lead to disastrous effects in the prediction of 18-mer behavior. While it is not clear whether this approximation would hold as well for larger systems, proteins are not very FIG. 4. The phase diagram from computer simulations. Shown are the results from enumeration of all 18-mer conformations with IIM interactions and B p 0.5. All aspects of the theory are seen in the computer simulations, including order of the transitions, etc. Note that the B 0 scale here is for mean interactions and does not include the excluded volume present in the nature of conformations used. Thus, B 0 0 corresponds to a very good solvent. long either and, based upon the number of non-local contacts present, 18-mers potentially correspond to short proteins, i.e., residues. In Fig. 5, we plot the phase boundaries for sequences with different ground state energies. We see that increasing B p decreases the stability of the ground state, but does not FIG. 3. Free energy for a well designed sequence in a good solvent. We see a first order transition from a frozen globule a to a random globule c with the coexistence between the two shown in b dark areas have lower free energy. Frame d shows F(q q) for these three regimes. FIG. 5. In addition to the data of Fig. 4, we now include several other values for B p the symbols are numbers referring to the relative ground state energy. We see that lowering B p makes the frozen state less stable and thus moves the boundary between the frozen and the random globule. On the other hand, changing B p does not have any effect of the homopolymer-like random globule to coil transition a globule is here defined to be maximally compact.
6 Pande, Grosberg, and Tanaka: Coil to globule transition 5123 FIG. 6. In addition to the data of Fig. 5, we now include curves delineating the random globule to coil phase boundary for different definitions of a globule the number of contacts defining a globule is encircled; maximum overlap has 16 contacts. For all definitions of the phase boundary, the random globule to coil transition is independent of the nature of design. have any significant effect on the globule to coil transition. Finally, in Fig. 6, we define the boundary between the coil and the globule in different ways with different threshold Q); again, we find that the aspects of the coil to either random or frozen globule transition are not affected by the nature of design. FIG. 7. Full phase diagram of the system, including the dependencies on B 0, B, and ground state energy E. IV. CONCLUSIONS We have described the phase transitions between the three phases for a single heteropolymer chain: coil, globule, and frozen globule. We found that heteropolymeric aspects have two effects: 1 due to good contacts, a heteropolymer can be collapsed even though the mean interaction is repulsive as was shown previously 11 ; 2 the freezing transition is coupled to the collapse transition and, for sufficiently good solvents, can lead to a first order transition in density Fig. 7. Heteropolymers are interesting both for their relevance to biological systems as well as biologically inspired systems. As for the latter, heteropolymer gels have attracted great interest as many microscopic properties of single polymer chains, such as the coil to globule transition, have macroscopic analogs in gels, which are macroscopic macromolecules. For example, in good solvents, the freezing transition will be accompanied by a first order transition in density; this may be a simple means to identify the freezing transition in these materials. As for proteins, recently protein folding kinetics from the random globule to the frozen globule has been described using a similar non-rem theory and physical arguments which link the nature of the free energy as a function of q i.e., a single free energy barrier with the nature of kinetics. However, if the transition is from the coil to the frozen globule, one must consider the free energy barrier not in the one dimensional q phase space, but the two dimensional (q,q) space. This adds complications to the theory, but may be important in the understanding of many proteins whose denatured states are coils, not globules. ACKNOWLEDGMENTS The work was supported by the National Science Foundation DMR VSP acknowledges support from the Miller Institute for Basic Research. AYG acknowledges the support of a Kao Fellowship. APPENDIX: ON THE EQUATION (8) As we expect freezing to occur at temperatures higher than that of the phase segregation of monomers, we can use a high temperature expansion to describe freezing. 12 This yields F seq,q H,q 1 2T H2,q H,q 2 T ln M, A1 where (...),q means the average over all conformations with the given density and overlap to the ground state, and M is the total number of conformations. We have kept terms to O (1/T 2 ), which is exact for the Independent Interaction Model in which the elements of B ij are taken from a Gaussian distribution ; 11 higher order terms may be necessary to quantitatively model a particular set of interactions B ij. As one can only analytically perform a statistical analysis, we have to average over the sequences. In principle, one could directly average over the ensemble of sequences with
7 5124 Pande, Grosberg, and Tanaka: Coil to globule transition the given value of the ground state energy. This ensemble of sequences is in a way similar to the microcanonical ensemble in regular statistical mechanics. It is, however, technically simpler to use analog of the canonical ensemble where native state energy E is not fixed, but rather controlled through a temperature T p equivalently, 1/T p is a Lagrange multiplier to constrain E ), and thus it has the Gibbs probability distribution: P seq exp E NS seq /T p 1 1 T p H, A2 where in the last transformation we have resorted to high T p expansion which is justified for the same reason as high temperature expansion above. Thus, we characterize a given ensemble of sequences by the value of T p : for lower T p,we model sequences whose native states are better optimized energetically. This prescription for energetic optimization of the native state has many incarnations, including minimal frustration, 10 sequence selection, 17 and Imprinting. 12 When one averages free energy A1 over the probability distribution A2, the result appears very simple, because, due to the structure of the Hamiltonian, the energy depends directly on q, and thus it is very easy to implement the condition of the fixed q, and the result can be written in the form F,q P seq F seq,q E,q TS,q, A3 seq where the q-dependent part of the energy E(,q) for the maximally compact conformations is given by E,q max B2 q B2 1 q, A4 T p T where B 2 ij p i p j (B ij kl p k p l B kl ) is the variance of the elements of the interaction matrix, 12 p i is the fraction of monomers of type i, and the entropy S(q) of conformations with q contacts in common with is S q, ln M q s 3 2 ln q q ln q 1 q ln 1 q. A5 The first energy term comes from the mean of the density of states H and describes how the energy is pulled down due to selection; indeed, to this order, the ground state energy is E B 2 /T p, and thus this term says that on average, the energy is given by the energy of the native state times the fraction of native contacts. The second energy term comes from the width of the density of states, resulting from the H 2 H 2 terms; the q dependence of this width enters from the correlator H 2, 20 as two conformations with a given overlap q with also have q contacts in common, on average. Higher order terms can modify the q dependence of the width, depending on the nature of interactions. Note that the q dependence of the mean and width of the density of states is a correction to REM, describing the nature of energy correlations. One can understand the energy terms in a simple manner: each native contact yields an energetic bonus of E /Q max and each non-native contact yields B 2 /T, i.e., is essentially annealed at temperature T. The first two entropy terms detail how many conformations exist for a fixed set of q contacts; as each common contact pins two parts of the chain together, for q contacts, the number of conformations is the loop entropy 1 of O (Q) loops of length O (1/Q). 21,22 The third entropy term results from how many ways one can choose q contacts out of the total Q max using Sterling s formula. 23 The parameter s is roughly the entropy per contact; its relevance to freezing has been detailed extensively elsewhere. 18 Finally, note that the contribution from the loop factor term in d dimensions is (d/2)q ln q. This term cancels exactly with the q ln q term from the native contact mixing entropy at d 2 and predicts very different behavior for two dimensions A. Yu. Grosberg and A. R. Khokhlov, Statistical Physics of Macromolecules American Institute of Physics, New York, A. Yu. Grosberg and D. V. Kuznetsov, Macromolecules 25, B. Chu, J. Polym. Sci., Part B Polym. Phys. 31, ; B. Chu, Q. Ying, and A. Yu. Grosberg, Macromolecules 28, C. Wu and S. Q. Zhou, Macromolecules 28, A. Yu. Grosberg, Biofizika 29, ; T. Garel, L. Leibler, and H. Orland, J. Phys. II 4, Y. Kantor and M. Kardar, Phys. Rev. E 52, ; A. V. Dobrynin and M. Rubinstein, J. Phys. II 5, L. Noirez, P. Keller, and J. P. Cotton, Liq. Cryst. 18, C. Wu and S. Zhou, Phys. Rev. Lett. 77, O. B. Ptitsyn, Adv. Protein Chem. 47, J. D. Bryngelson and P. G. Wolynes, Proc. Natl. Acad. Sci., USA 84, E. Shakhnovich and A. Gutin, Biophys. Chem. 34, V. S. Pande, A. Yu. Grosberg, and T. Tanaka, Macromolecules 28, J. D. Bryngelson and P. G. Wolynes, Biopolymers 30, A. Dinner, A. Sali, M. Karplus, and E. Shakhnovich, J. Chem. Phys. 101, S. Plotkin, J. Wang, and P. G. Wolynes, J. Chem. Phys. 106, B. Derrida, Phys. Rev. Lett. 45, E. I. Shakhnovich and A. M. Gutin, Proc. Natl. Acad. Sci. USA 90, V. S. Pande, A. Yu. Grosberg, and T. Tanaka, Folding & Design 2, V. S. Pande, C. Joerg, A. Yu. Grosberg, and T. Tanaka, J. Phys. A 27, V. S. Pande, A. Yu. Grosberg, C. Joerg, and T. Tanaka, Phys. Rev. Lett. 77, A. M. Gutin and E. I. Shakhnovich, J. Chem. Phys. 100, Note that the validity of the loop factor relies on the assumption of Gaussian statistics of the chains. For coiled conformations, this holds to a reasonable approximation and for globular conformations, the chains behave as chains in a melt, again with Gaussian statistics. Thus, the entropy A5 can be considered as the interpolation of the entropy for these two regimes and thus is most questionable in the region in between swollen and collapsed conformations; as this region is not of principle interest, we expect A5 to be sufficient. Another approximation implicit in A5 is that loops lengths are treated in a mean field manner and are assumed to be of the same length; this is valid for short chains, such as the 18-mers used in the computer simulation section. Inclusion of the fluctuation of loop lengths is an interesting correction to our theory, but beyond the scope of this paper. 23 S. Plotkin, J. Wang, and P. G. Wolynes, Phys. Rev. E 53, We are indebted to E. Shakhnovich for this point.
Freezing of compact random heteropolymers with correlated sequence fluctuations
Freezing of compact random heteropolymers with correlated sequence fluctuations Arup K. Chakraborty Department of Chemical Engineering and Department of Chemistry, University of California, Berkeley, California
More informationarxiv:cond-mat/ v1 2 Feb 94
cond-mat/9402010 Properties and Origins of Protein Secondary Structure Nicholas D. Socci (1), William S. Bialek (2), and José Nelson Onuchic (1) (1) Department of Physics, University of California at San
More informationRandom heteropolymer adsorption on disordered multifunctional surfaces: Effect of specific intersegment interactions
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 15 15 OCTOBER 1998 Random heteropolymer adsorption on disordered multifunctional surfaces: Effect of specific intersegment interactions Simcha Srebnik Department
More informationSwelling and Collapse of Single Polymer Molecules and Gels.
Swelling and Collapse of Single Polymer Molecules and Gels. Coil-Globule Transition in Single Polymer Molecules. the coil-globule transition If polymer chains are not ideal, interactions of non-neighboring
More informationA simple theory and Monte Carlo simulations for recognition between random heteropolymers and disordered surfaces
A simple theory and Monte Carlo simulations for recognition between random heteropolymers and disordered surfaces Arup K. Chakraborty Department of Chemical Engineering and Department of Chemistry, University
More informationProteins polymer molecules, folded in complex structures. Konstantin Popov Department of Biochemistry and Biophysics
Proteins polymer molecules, folded in complex structures Konstantin Popov Department of Biochemistry and Biophysics Outline General aspects of polymer theory Size and persistent length of ideal linear
More informationStatistical Physics of The Symmetric Group. Mobolaji Williams Harvard Physics Oral Qualifying Exam Dec. 12, 2016
Statistical Physics of The Symmetric Group Mobolaji Williams Harvard Physics Oral Qualifying Exam Dec. 12, 2016 1 Theoretical Physics of Living Systems Physics Particle Physics Condensed Matter Astrophysics
More informationA first-order transition in the charge-induced conformational changes of polymers
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 22 8 JUNE 2002 A first-order transition in the charge-induced conformational changes of polymers Yi Mao, Alexander L. Burin, a) Mark A. Ratner, and Martin
More informationChap. 2. Polymers Introduction. - Polymers: synthetic materials <--> natural materials
Chap. 2. Polymers 2.1. Introduction - Polymers: synthetic materials natural materials no gas phase, not simple liquid (much more viscous), not perfectly crystalline, etc 2.3. Polymer Chain Conformation
More informationarxiv: v1 [cond-mat.soft] 22 Oct 2007
Conformational Transitions of Heteropolymers arxiv:0710.4095v1 [cond-mat.soft] 22 Oct 2007 Michael Bachmann and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11,
More informationEngineering of stable and fast-folding sequences of model proteins
Proc. Natl. Acad. Sci. USA Vol. 90, pp. 7195-7199, August 1993 Biophysics Engineering of stable and fast-folding sequences of model proteins E. I. SHAKHNOVICH AND.A. M. GUTIN Department of Chemistry, Harvard
More informationMany proteins spontaneously refold into native form in vitro with high fidelity and high speed.
Macromolecular Processes 20. Protein Folding Composed of 50 500 amino acids linked in 1D sequence by the polypeptide backbone The amino acid physical and chemical properties of the 20 amino acids dictate
More informationFormation of microdomains in a quenched disordered heteropolymer
Formation of microdomains in a quenched disordered heteropolymer E.I. Shakhnovich, A.M. Gutin To cite this version: E.I. Shakhnovich, A.M. Gutin. Formation of microdomains in a quenched disordered heteropolymer.
More informationStatistical Mechanics of Simple Models of Protein Folding and Design
3192 Biophysical Joumal Volume 73 December 1997 3192-3210 Statistical Mechanics of Simple Models of Protein Folding and Design Vijay S. Pande,* Alexander Yu. Grosberg,# and Toyoichi Tanaka# *Physics Department,
More informationProtein design: a perspective from simple tractable models Eugene I Shakhnovich
Review R45 Protein design: a perspective from simple tractable models Eugene I Shakhnovich Recent progress in computational approaches to protein design builds on advances in statistical mechanical protein
More information2.4 DNA structure. S(l) bl + c log l + d, with c 1.8k B. (2.72)
2.4 DNA structure DNA molecules come in a wide range of length scales, from roughly 50,000 monomers in a λ-phage, 6 0 9 for human, to 9 0 0 nucleotides in the lily. The latter would be around thirty meters
More informationPolymer Solution Thermodynamics:
Polymer Solution Thermodynamics: 3. Dilute Solutions with Volume Interactions Brownian particle Polymer coil Self-Avoiding Walk Models While the Gaussian coil model is useful for describing polymer solutions
More informationPhase Transitions. µ a (P c (T ), T ) µ b (P c (T ), T ), (3) µ a (P, T c (P )) µ b (P, T c (P )). (4)
Phase Transitions A homogeneous equilibrium state of matter is the most natural one, given the fact that the interparticle interactions are translationally invariant. Nevertheless there is no contradiction
More informationTwo-Dimensional Polymers with Random. Eilon Brenner. for the M.Sc. degree. at Tel-Aviv University. School of Physics and Astronomy
TEL AVIV UNIVERSITY RAYMOND AND BEVERLY SACKLER FACULTY OF EXACT SCIENCES SCHOOL OF PHYSICS & ASTRONOMY a aia`-lz zhiqxaipe` miwiiecn mircnl dhlewtd xlw`q ilxaae cpeniix y"r dinepexhq`e dwiqitl xtqd zia
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 86 28 MAY 21 NUMBER 22 Mathematical Analysis of Coupled Parallel Simulations Michael R. Shirts and Vijay S. Pande Department of Chemistry, Stanford University, Stanford,
More informationS(l) bl + c log l + d, with c 1.8k B. (2.71)
2.4 DNA structure DNA molecules come in a wide range of length scales, from roughly 50,000 monomers in a λ-phage, 6 0 9 for human, to 9 0 0 nucleotides in the lily. The latter would be around thirty meters
More informationPolymer solutions and melts
Course M6 Lecture 9//004 (JAE) Course M6 Lecture 9//004 Polymer solutions and melts Scattering methods Effects of excluded volume and solvents Dr James Elliott Online teaching material reminder Overheads
More informationPredicting free energy landscapes for complexes of double-stranded chain molecules
JOURNAL OF CHEMICAL PHYSICS VOLUME 114, NUMBER 9 1 MARCH 2001 Predicting free energy landscapes for complexes of double-stranded chain molecules Wenbing Zhang and Shi-Jie Chen a) Department of Physics
More informationConfinement of polymer chains and gels
Confinement of polymer chains and gels Nefeli Georgoulia - Student number: 70732831 1 Introduction Confinement of polymer chains is significant in industrial as well as biological applications. For this
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 5 May 1998
Linking Rates of Folding in Lattice Models of Proteins with Underlying Thermodynamic Characteristics arxiv:cond-mat/9805061v1 [cond-mat.soft] 5 May 1998 D.K.Klimov and D.Thirumalai Institute for Physical
More informationCoil to Globule Transition: This follows Giant Molecules by Alexander Yu. Grosberg and Alexei R. Khokhlov (1997).
Coil to Globule Transition: This follows Giant Molecules by Alexander Yu. Grosberg and Alexei R. Khokhlov (1997). The Flory Krigbaum expression for the free energy of a self-avoiding chain is given by,
More informationPhysics 562: Statistical Mechanics Spring 2002, James P. Sethna Prelim, due Wednesday, March 13 Latest revision: March 22, 2002, 10:9
Physics 562: Statistical Mechanics Spring 2002, James P. Sethna Prelim, due Wednesday, March 13 Latest revision: March 22, 2002, 10:9 Open Book Exam Work on your own for this exam. You may consult your
More informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
More informationPotts And XY, Together At Last
Potts And XY, Together At Last Daniel Kolodrubetz Massachusetts Institute of Technology, Center for Theoretical Physics (Dated: May 16, 212) We investigate the behavior of an XY model coupled multiplicatively
More informationSpontaneous Symmetry Breaking
Spontaneous Symmetry Breaking Second order phase transitions are generally associated with spontaneous symmetry breaking associated with an appropriate order parameter. Identifying symmetry of the order
More informationThe glass transition as a spin glass problem
The glass transition as a spin glass problem Mike Moore School of Physics and Astronomy, University of Manchester UBC 2007 Co-Authors: Joonhyun Yeo, Konkuk University Marco Tarzia, Saclay Mike Moore (Manchester)
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 19 Mar 2001
Modeling two-state cooperativity in protein folding Ke Fan, Jun Wang, and Wei Wang arxiv:cond-mat/0103385v1 [cond-mat.soft] 19 Mar 2001 National Laboratory of Solid State Microstructure and Department
More informationLocal Interactions Dominate Folding in a Simple Protein Model
J. Mol. Biol. (1996) 259, 988 994 Local Interactions Dominate Folding in a Simple Protein Model Ron Unger 1,2 * and John Moult 2 1 Department of Life Sciences Bar-Ilan University Ramat-Gan, 52900, Israel
More informationFirst steps toward the construction of a hyperphase diagram that covers different classes of short polymer chains
Cent. Eur. J. Phys. 12(6) 2014 421-426 DOI: 10.2478/s11534-014-0461-z Central European Journal of Physics First steps toward the construction of a hyperphase diagram that covers different classes of short
More informationFREQUENCY selected sequences Z. No.
COMPUTER SIMULATIONS OF PREBIOTIC EVOLUTION V.I. ABKEVICH, A.M. GUTIN, and E.I. SHAKHNOVICH Harvard University, Department of Chemistry 12 Oxford Street, Cambridge MA 038 This paper is a review of our
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 23 Mar 2007
The structure of the free energy surface of coarse-grained off-lattice protein models arxiv:cond-mat/0703606v1 [cond-mat.soft] 23 Mar 2007 Ethem Aktürk and Handan Arkin Hacettepe University, Department
More informationA Deeper Look into Phase Space: The Liouville and Boltzmann Equations
A Deeper Look into Phase Space: The Liouville and Boltzmann Equations Brian Petko Statistical Thermodynamics Colorado School of Mines December 7, 2009 INTRODUCTION Phase Space, in a sense, can be a useful
More informationComputer simulation of polypeptides in a confinement
J Mol Model (27) 13:327 333 DOI 1.17/s894-6-147-6 ORIGINAL PAPER Computer simulation of polypeptides in a confinement Andrzej Sikorski & Piotr Romiszowski Received: 3 November 25 / Accepted: 27 June 26
More informationThermodynamics of nuclei in thermal contact
Thermodynamics of nuclei in thermal contact Karl-Heinz Schmidt, Beatriz Jurado CENBG, CNRS/IN2P3, Chemin du Solarium B.P. 120, 33175 Gradignan, France Abstract: The behaviour of a di-nuclear system in
More informationStretching lattice models of protein folding
Proc. Natl. Acad. Sci. USA Vol. 96, pp. 2031 2035, March 1999 Biophysics Stretching lattice models of protein folding NICHOLAS D. SOCCI,JOSÉ NELSON ONUCHIC**, AND PETER G. WOLYNES Bell Laboratories, Lucent
More informationGuessing the upper bound free-energy difference between native-like structures. Jorge A. Vila
1 Guessing the upper bound free-energy difference between native-like structures Jorge A. Vila IMASL-CONICET, Universidad Nacional de San Luis, Ejército de Los Andes 950, 5700- San Luis, Argentina Use
More informationUniversal correlation between energy gap and foldability for the random energy model and lattice proteins
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 14 8 OCTOBER 1999 Universal correlation between energy gap and foldability for the random energy model and lattice proteins Nicolas E. G. Buchler Biophysics
More informationarxiv:chem-ph/ v1 11 Nov 1994
chem-ph/9411008 Funnels, Pathways and the Energy Landscape of Protein Folding: A Synthesis arxiv:chem-ph/9411008v1 11 Nov 1994 Joseph D. Bryngelson, Physical Sciences Laboratory, Division of Computer Research
More informationEvolution of functionality in lattice proteins
Evolution of functionality in lattice proteins Paul D. Williams,* David D. Pollock, and Richard A. Goldstein* *Department of Chemistry, University of Michigan, Ann Arbor, MI, USA Department of Biological
More informationShort Announcements. 1 st Quiz today: 15 minutes. Homework 3: Due next Wednesday.
Short Announcements 1 st Quiz today: 15 minutes Homework 3: Due next Wednesday. Next Lecture, on Visualizing Molecular Dynamics (VMD) by Klaus Schulten Today s Lecture: Protein Folding, Misfolding, Aggregation
More informationDistance Constraint Model; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 11
Distance Constraint Model; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 11 Taking the advice of Lord Kelvin, the Father of Thermodynamics, I describe the protein molecule and other
More informationPhase Transition in a Bond. Fluctuating Lattice Polymer
Phase Transition in a Bond arxiv:1204.2691v2 [cond-mat.soft] 13 Apr 2012 Fluctuating Lattice Polymer Hima Bindu Kolli and K.P.N.Murthy School of Physics, University of Hyderabad Hyderabad 500046 November
More informationchem-ph/ Feb 95
LU-TP 9- October 99 Sequence Dependence of Self-Interacting Random Chains Anders Irback and Holm Schwarze chem-ph/9 Feb 9 Department of Theoretical Physics, University of Lund Solvegatan A, S- Lund, Sweden
More informationDesigning refoldable model molecules
Designing refoldable model molecules I. Coluzza, H. G. Muller, and D. Frenkel FOM Institute for Atomic and Molecular Physics, Kruislaan, 407 1098 SJ Amsterdam, The Netherlands Received 16 April 2003; published
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationTriangular Lattice Foldings-a Transfer Matrix Study.
EUROPHYSICS LETTERS Europhys. Lett., 11 (2)) pp. 157-161 (1990) 15 January 1990 Triangular Lattice Foldings-a Transfer Matrix Study. Y. KANT~R(*) and M. V. JARIC(**) (*) School of Physics and Astronomy,
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationTHEORY OF PROTEIN FOLDING: The Energy Landscape Perspective
Annu. Rev. Phys. Chem. 1997. 48:545 600 Copyright c 1997 by Annual Reviews Inc. All rights reserved THEORY OF PROTEIN FOLDING: The Energy Landscape Perspective José Nelson Onuchic Department of Physics,
More informationStacking and Hydrogen Bonding. DNA Cooperativity at Melting.
Stacking and Hydrogen Bonding. DNA Cooperativity at Melting. Vladimir F. Morozov *, Artem V. Badasyan, Arsen V. Grigoryan, Mihran A. Sahakyan, Evgeni Sh. Mamasakhlisov. 1 Department of Molecular Physics,
More informationNuclear Science Research Group
UR-NCUR 05-12 THE UNIVERSITY OF ROCHESTER Nuclear Science Research Group ROCHESTER, NEW YORK 14267-0216, USA Compound Nucleus with Surface Entropy; A Unified Phenomenology of Particle Evaporation, Fission,
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 25 Apr 2001
IL NUOVO CIMENTO Vol. 20D, N. 12bis, pp. 2265-2288 1998 arxiv:cond-mat/0104495v1 [cond-mat.soft] 25 Apr 2001 Self consistent Treatment of Copolymers with Arbitrary Sequences Yu.A. Kuznetsov ( ), E.G. Timoshenko(
More informationOlle Inganäs: Polymers structure and dynamics. Polymer physics
Polymer physics Polymers are macromolecules formed by many identical monomers, connected through covalent bonds, to make a linear chain of mers a polymer. The length of the chain specifies the weight of
More informationThe (magnetic) Helmholtz free energy has proper variables T and B. In differential form. and the entropy and magnetisation are thus given by
4.5 Landau treatment of phase transitions 4.5.1 Landau free energy In order to develop a general theory of phase transitions it is necessary to extend the concept of the free energy. For definiteness we
More informationPhase transitions and critical phenomena
Phase transitions and critical phenomena Classification of phase transitions. Discontinous (st order) transitions Summary week -5 st derivatives of thermodynamic potentials jump discontinously, e.g. (
More informationA simple technique to estimate partition functions and equilibrium constants from Monte Carlo simulations
A simple technique to estimate partition functions and equilibrium constants from Monte Carlo simulations Michal Vieth Department of Chemistry, The Scripps Research Institute, 10666 N. Torrey Pines Road,
More informationPhysica A. A semi-flexible attracting segment model of two-dimensional polymer collapse
Physica A 389 (2010) 1619 1624 Contents lists available at ScienceDirect Physica A journal homepage: www.elsevier.com/locate/physa A semi-flexible attracting segment model of two-dimensional polymer collapse
More informationTHE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION
THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION AND CALIBRATION Calculation of turn and beta intrinsic propensities. A statistical analysis of a protein structure
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 16 Nov 2002
Dependence of folding rates on protein length Mai Suan Li 1, D. K. Klimov 2 and D. Thirumalai 2 1 Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland 2 Institute
More informationLong Range Moves for High Density Polymer Simulations
arxiv:cond-mat/9610116v1 [cond-mat.soft] 15 Oct 1996 Long Range Moves for High Density Polymer Simulations J.M.Deutsch University of California, Santa Cruz, U.S.A. Abstract Monte Carlo simulations of proteins
More informationAmorphous Polymers: Polymer Conformation Laboratory 1: Module 1
D E P A R T M E N T O F M A T E R I A L S S C I E N C E A N D E N G I N E E R I N G M A S S A C H U S E T T S I N S T I T U T E O F T E C H N O L O G Y 3.014 Materials Laboratory Fall 2008 Amorphous Polymers:
More informationCollapse of a polymer chain: A Born Green Yvon integral equation study
Collapse of a polymer chain: A Born Green Yvon integral equation study Mark P. Taylor and J. E. G. Lipson Department of Chemistry, Dartmouth College, Hanover, New Hampshire 03755 Received 13 November 1995;
More informationMaster equation approach to finding the rate-limiting steps in biopolymer folding
JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 7 15 FEBRUARY 2003 Master equation approach to finding the rate-limiting steps in biopolymer folding Wenbing Zhang and Shi-Jie Chen a) Department of Physics
More informationScaling Law for the Radius of Gyration of Proteins and Its Dependence on Hydrophobicity
Scaling Law for the Radius of Gyration of Proteins and Its Dependence on Hydrophobicity LIU HONG, JINZHI LEI Zhou Pei-Yuan Center for Applied Mathematics, Tsinghua University, Beijing, People s Republic
More informationPhase transition and spontaneous symmetry breaking
Phys60.nb 111 8 Phase transition and spontaneous symmetry breaking 8.1. Questions: Q1: Symmetry: if a the Hamiltonian of a system has certain symmetry, can the system have a lower symmetry? Q: Analyticity:
More informationCompetition between face-centered cubic and icosahedral cluster structures
Competition between face-centered cubic and icosahedral cluster structures R. S. Berry University of Chicago, Chicago, IL 60637, USA B. M. Smyrnov, and A. Yu. Strizhev High-Temperatures Institute, Russian
More informationTransitions of tethered polymer chains: A simulation study with the bond fluctuation lattice model
THE JOURNAL OF CHEMICAL PHYSICS 128, 064903 2008 Transitions of tethered polymer chains: A simulation study with the bond fluctuation lattice model Jutta Luettmer-Strathmann, 1,a Federica Rampf, 2 Wolfgang
More informationResearch Paper 577. Correspondence: Nikolay V Dokholyan Key words: Go model, molecular dynamics, protein folding
Research Paper 577 Discrete molecular dynamics studies of the folding of a protein-like model ikolay V Dokholyan, Sergey V Buldyrev, H Eugene Stanley and Eugene I Shakhnovich Background: Many attempts
More informationLecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions: Van der Waals Interactions
More informationGrowth oscillations. LElTER TO THE EDITOR. Zheming Cheng and Robert Savit
J. Phys. A: Math. Gen. 19 (1986) L973-L978. Printed in Great Britain LElTER TO THE EDITOR Growth oscillations Zheming Cheng and Robert Savit Department of Physics, The University of Michigan, Ann Arbor,
More informationThe Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by
More informationPathways for protein folding: is a new view needed?
Pathways for protein folding: is a new view needed? Vijay S Pande 1, Alexander Yu Grosberg 2, Toyoichi Tanaka 2, and Daniel S Rokhsar 1;3 Theoretical studies using simplified models for proteins have shed
More information8.334: Statistical Mechanics II Spring 2014 Test 2 Review Problems
8.334: Statistical Mechanics II Spring 014 Test Review Problems The test is closed book, but if you wish you may bring a one-sided sheet of formulas. The intent of this sheet is as a reminder of important
More informationV.E Mean Field Theory of Condensation
V.E Mean Field heory of Condensation In principle, all properties of the interacting system, including phase separation, are contained within the thermodynamic potentials that can be obtained by evaluating
More informationarxiv: v1 [cond-mat.stat-mech] 6 Mar 2008
CD2dBS-v2 Convergence dynamics of 2-dimensional isotropic and anisotropic Bak-Sneppen models Burhan Bakar and Ugur Tirnakli Department of Physics, Faculty of Science, Ege University, 35100 Izmir, Turkey
More informationWhy Complexity is Different
Why Complexity is Different Yaneer Bar-Yam (Dated: March 21, 2017) One of the hardest things to explain is why complex systems are actually different from simple systems. The problem is rooted in a set
More informationCRITICAL SLOWING DOWN AND DEFECT FORMATION M. PIETRONI. INFN - Sezione di Padova, via F. Marzolo 8, Padova, I-35131, ITALY
CRITICAL SLOWING DOWN AND DEFECT FORMATION M. PIETRONI INFN - Sezione di Padova, via F. Marzolo 8, Padova, I-35131, ITALY E-mail: pietroni@pd.infn.it The formation of topological defects in a second order
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationarxiv: v1 [cond-mat.soft] 4 Sep 2018
Generalized Flory Theory for Rotational Symmetry Breaking of Complex Macromolecules Josh Kelly 1, Alexander Y. Grosberg, and Robijn Bruinsma 1,3 1 Department of Physics and Astronomy, University of California,
More informationPhysics 212: Statistical mechanics II Lecture XI
Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is
More informationIntroduction to Computational Structural Biology
Introduction to Computational Structural Biology Part I 1. Introduction The disciplinary character of Computational Structural Biology The mathematical background required and the topics covered Bibliography
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationUNDERSTANDING BOLTZMANN S ANALYSIS VIA. Contents SOLVABLE MODELS
UNDERSTANDING BOLTZMANN S ANALYSIS VIA Contents SOLVABLE MODELS 1 Kac ring model 2 1.1 Microstates............................ 3 1.2 Macrostates............................ 6 1.3 Boltzmann s entropy.......................
More informationEnsemble equivalence for non-extensive thermostatistics
Physica A 305 (2002) 52 57 www.elsevier.com/locate/physa Ensemble equivalence for non-extensive thermostatistics Raul Toral a;, Rafael Salazar b a Instituto Mediterraneo de Estudios Avanzados (IMEDEA),
More informationOn the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 12 22 SEPTEMBER 1998 On the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models Nobuyuki Matubayasi Institute
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationFrontiers in Physics 27-29, Sept Self Avoiding Growth Walks and Protein Folding
Frontiers in Physics 27-29, Sept. 2012 Self Avoiding Growth Walks and Protein Folding K P N Murthy, K Manasa, and K V K Srinath School of Physics, School of Life Sciences University of Hyderabad September
More informationarxiv:cond-mat/ v4 [cond-mat.stat-mech] 19 Jun 2007
arxiv:cond-mat/060065v4 [cond-mat.stat-mech] 9 Jun 007 Restoration of Isotropy in the Ising Model on the Sierpiński Gasket Naoto Yajima Graduate School of Human and Environmental Studies, Kyoto University,
More information3.091 Introduction to Solid State Chemistry. Lecture Notes No. 9a BONDING AND SOLUTIONS
3.091 Introduction to Solid State Chemistry Lecture Notes No. 9a BONDING AND SOLUTIONS 1. INTRODUCTION Condensed phases, whether liquid or solid, may form solutions. Everyone is familiar with liquid solutions.
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
1 Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 2, 24 March 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationAdsorption from a one-dimensional lattice gas and the Brunauer Emmett Teller equation
Proc. Natl. Acad. Sci. USA Vol. 93, pp. 1438 1433, December 1996 Chemistry Adsorption from a one-dimensional lattice gas and the Brunauer Emmett Teller equation (Ising modelreference systemgibbs ecess
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationAn Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory
in Free Massive Scalar Field Theory NCSR Demokritos National Technical University of Athens based on arxiv:1711.02618 [hep-th] in collaboration with Dimitris Katsinis March 28 2018 Entanglement and Entanglement
More informationQuantitative Stability/Flexibility Relationships; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 12
Quantitative Stability/Flexibility Relationships; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 12 The figure shows that the DCM when applied to the helix-coil transition, and solved
More informationPolymer Physics MSE 458 / CHEM 482 Spring 2018
Polymer Physics MSE 458 / CHEM 482 Spring 2018 Instructor: Prof. A.L. Ferguson 204 MSEB (217) 300-2354 alf@illinois.edu Grader: Class: Location: 4101 MSEB Time: 2:00 3:20 pm Days: T, Th Sections: A3 (CRN-38260)
More information