Fractal nature of the phase space and energy landscape topology

1. Fractal nature of the phase space and energy landscape topology Gerardo G. Naumis Instituto de Física, UNAM. México D.F., Mexico. XVIII Meeting on Complex Fluids, San Luis Potosí, México.

2. Introduction Relaxation and flexibility in polymers, proteins, colloids and fluids. Energy landscape formalism • Topography of the phase space and energy landscape. • A modified Monte-Carlo method to test the topology and topography. • Applications to the most simple fluid • Conclusions

3. In many systems, we are interested in following the temporal evolution Failure in folding: Alzheimer's disease, cystic fibrosis, BSE (Mad Cow disease), an inherited form of emphysema, and even many cancers are believed to result from protein misfolding.

4. WHY IS PROTEIN FOLDING SO DIFFICULT TO UNDERSTAND? It's amazing that not only do proteins self-assemble -- fold -- but they do so amazingly quickly: some as fast as a millionth of a second. It takes about a day to simulate a nanosecond (1/1,000,000,000 of a second). Unfortunately, proteins fold on the tens of microsecond timescale (10,000 nanoseconds). Thus, it would take 10,000 CPU days to simulate folding -- i.e. it would take 30 CPU years! But is also important to understand the “mechanical flexibility”:

5. Dynamics and rheology in dense colloids, glasses, jamming in granular media, etc. Example: experimental results in coloids from the group of David Weitz High volume fraction supercooled fluid: volume fraction 0.56 Highlighted particles are slow by over a timestep of 3600 seconds. At this timestep, the largest slow cluster percolates. One timestep (18 seconds) later, the percolating supercooled fluid sample has broken up. Glass: volume fraction 0.60 The third sample is a glass at volume fraction 0.56. The highlighted particles are particles which are slow over an entire experiment, a timestep of 39,000 seconds. At this timestep, the largest slow cluster percolates. Over the experimentally accessible timescales, this percolating cluster never breaks up!

6. E Energy landscapes and rigidity There are many approaches to solve these problems, but in fact the Hamiltonian contains all this information… The mechanical state of the system is represented as a point in phase space: The allowed part of the phase space is determined by the “energy landscape” Basin Statistics of landscape

7. Decoy tree (protein:villin) 

9. Saddle points • Distribution of energy basins • Size of each basin FRACTALS! Some predictions were made about the range of the potential-roughness using catastophy theory. Short range: rough landspace Science, Vol 293, Issue 5537, 2067-2070 , 14 September 2001

10. FRACTALES: AUTOSIMILARIDAD

11. Rigidity Theory With N hindges, how many bars do I need to make the system rigid? Rigid Flexible Isostatic -1 1 0 4x (2 freedom degrees)-(# constraints)=# flexible movements f =(3N-constraints)/3NFraction of “floopy modes since: # flexible movements=# of normal modes of vibration with zero frequency f is a function of <r>. In the Maxwell approximation: f=2-5<r>/6.

12. Maxwell Counting (1860) 1/3 Rigid Flexible f f=0, <r>=2.4 2.0 2.2 2.4 <r>

13. Interpretation in terms of energy landscapes Floppy modes provide channels in the landscape: a lot of entropy!! Channels are not flat; there is a small curvature along the floopy coordinate, since floppy modes are not at non-zero frequency.

14. Energy landscapes and rigidity There are many approaches to solve these problems, but the Hamiltonian in fact contains all this information…

15. Adams-Gibbs equation: Tatsumisago et. al., Phys. Rev. Lett. 64, 1549 (1990).

16. Intermediate phase(P. Boolchand et. al., J. of Optoelectronics and Advanced Materials Vol. 3, 703 (2001)). <r> Self-organization: Thorpe, et. al., J. Non-Cryst. Solids 266, 859 (2000). Barré et. al., Phys. Rev. Lett. 94, 208701(2005) Boolchand et. al., J. of Non-Cryst. Solids 293, 348 (2001).

17. Rigidity of proteins and glasses Results of Mike Thorpe, Arizona State University To read more: G. Naumis, “Energy landscape and rigidity”, Phys. Rev. E71, 026114 (2005).

18. Topology of the phase space: the role of constraints… Hard-disks Restrictions: Boundary: Box, For N particles: Center of mass minimization observed in colloids!!! A. Van Blaaderen, Science 301, 471 (2003).

20. A simple example that explains the method…

22. Probabilidad de caer en un sitio de frontera: Sumando sobre toda la frontera: Definimos una coordinación y probabilidad promedio en los sitios de frontera:

23. The difference in dimensions between the phase space and the boundary is:

24. Application for the case of simple fluids: a hard disk system N=100 disks in a box, with hard core repulsion

27. CONCLUSIONS • The topology of the phase space and the topography of the energy landscape are important to understand several thermodynamical and relaxation phenomena. • This explains diverse features of simulations in associative fluids. • A method to obtain the fractal dimension using the Monte-Carlo rejection ratio was proposed. • The application of this method to a simple fluid shows the fractal nature of the phase space and that freezing occurs when the surface scales as the volume in phase space. • To read more: G.G. Naumis, Phys. Rev. E71, 056132 (2005).