# Difference between revisions of "Minimal energy packings and collapse of sticky tangent hard-sphere polymers"

(→MEPs for polymers) |
(→Collapse of Sticky Tangent Hard-Sphere Polymers) |
||

Line 33: | Line 33: | ||

== Collapse of Sticky Tangent Hard-Sphere Polymers == | == Collapse of Sticky Tangent Hard-Sphere Polymers == | ||

+ | |||

+ | The authors use molecular dynamics (MD) simulations to look at polymer dynamics during temperature quenches, as the polymer is cooled at different rates. | ||

+ | |||

+ | The MD parameters are as follows: monomers interact via the potential energy <math>U_{harm}(r) = -\epsilon + \frac{k}{2} (\frac{r}{D} - 1)^2</math> for <math>r < r_c</math> and 0 otherwise. | ||

== Discussion == | == Discussion == |

## Revision as of 03:53, 21 September 2010

Minimal energy packings and collapse of sticky tangent hard-sphere polymers. Hoy R.S. and O'Hern C.S. (2010). Phys. Rev. Lett. 105, 068001.

Leon Furchtgott, APP 225. Fall 2010.

## Contents

## Summary/Introduction

This paper is related to work investigating The Free Energy Landscape of Hard Sphere Clusters (experimental) and Minimal Energy Clusters of Hard Spheres with Short Ranged Attractions (theoretical). In these papers, Arkus et al. were interested in enumerating minimal energy packings (MEPs) for 10 sticky hard spheres and predicting and measuring the probabilities of different configurations.

This paper looks at packings of sticky tangent hard-sphere polymers, which can model polymer collapse, protein folding and protein interactions. This is particularly interesting and different from the case of sticky hard spheres because of the addition of covalent-bond and chain uncrossability constraints.

This paper looks at 2 overarching questions. (1) How do the probabilities for obtaining polymer MEPs differ from those for sticky hard-sphere MEPs? (2) How do the properties of single compact polymers depend on collapse dynamics?

(1) Polymer constraints reduce the number of MEP configurations. This is because of *dividing surfaces* which split the packings into disjoint regions.

(2) For slow temperature quenches, single chains undergo a sharp transition to crystallites with a close-packed core. As the quench rate increases, the core size decreases and the exterior becomes more disordered.

## MEPs for polymers

The authors use Arkus et al's geometric and graph-theory techniques to generate all of the macrostates and microstates for a given <math>N</math> number of spheres and <math>N_c</math> number of contacts. They start by enumerating adjacency matrices satisfying hard-sphere and minimal rigidity conditions, as well as the proper polymeric covalent-bond constraints.

Given a valid adjacency matrix, they solve the system of quadratic equations <math>|r_i - r_j|^2 = d^2_{ij}</math> for the sphere positions <math>r_i</math>.

They thus obtain a set of macrostates, each characterized by an adjacency matrix that is nonisomorphic, and a unique set of interparticle distances -- each macrostate is distinct and cannot be obtained through rotation or reflection from another macrostate.

Each macrostate corresponds to several microstates, which, for polymers, correspond to multiple possible paths through a given macrostate.

The authors then examine the fraction <math>f</math> of adjacency matrices with <math>N_c</math> contacts that satisfies hard-sphere constraints and minimal rigidity, and they compare this fraction for monomers <math>f_m</math>, linear polymers <math>f_p</math> and polymer rings <math>f_r</math>. They find that hard-sphere constraints are more difficult to satisfy for polymers than for monomers: <math>f_r < f_p < f_m</math>.

This is because of the occurrence of *dividing surfaces* in polymer packings. A dividing surface is a minimal subset of a connected cluster of contacting monomers that geometrically splits it into two. If monomers J and K are split by S, then any path that starts in J cannot cross S unless it passes through all of the monomers in S. This greatly reduces the number of possible arrangements. (See figure 1 in the paper).

Entropy suppression via blocking is strongest for structures of both high and low symmetry, therefore polymer MEPs favor intermediate structural symmetry. This is unlike monomer MEPs, which favor low structural symmetry for entropic reasons.

## Collapse of Sticky Tangent Hard-Sphere Polymers

The authors use molecular dynamics (MD) simulations to look at polymer dynamics during temperature quenches, as the polymer is cooled at different rates.

The MD parameters are as follows: monomers interact via the potential energy <math>U_{harm}(r) = -\epsilon + \frac{k}{2} (\frac{r}{D} - 1)^2</math> for <math>r < r_c</math> and 0 otherwise.