RUTGERS, THE STATE UNIVERSITY OF NEW JERSEY
HILL CENTER, ROOM 114
SUNDAY, MONDAY AND TUESDAY
DECEMBER 13-18, 2008


List of Short talk requests

Below is a list of short talk requests along with titles & abstracts. The short talk schedule will be posted the week of December 7th

________________________________


D. Adams, University of Michigan
Talk preference: Session A or B
Title: The barrier method: a new algorithm to measure rare transitions in non-equilibrium systems with applications to a model by Maier and Stein
Coauthors: L. M. Sander and R. M. Ziff
Abstract: In 1993, Maier and Stein introduced a simple non-equilibrium bistable stochastic model to study exit phenomena in the absence of detailed balance. A number of theoretical predictions were made for exit times and exit paths that have been difficult to verify because those predictions are only valid in the low-noise limit, which is not accessible to brute-force simulations. We present a new rare-event algorithm, which we call the barrier method. This algorithm allows us to verify some of the predictions in the low-noise limit by speeding up calculations by several orders of magnitude.

Y-Y. Anh, Session A
Talk preference: Session A
Title: Link communities reveal multi-scale complexity in networks
Coauthors: James P. Bagrow and Sune Lehmann
Abstract: Networks have become a key approach to understanding systems of interacting objects, unifying the study of diverse phenomena including biological organisms and human society. One crucial step when studying the structure and dynamics of networks is to identify communities; groups of related nodes that correspond to functional subunits such as protein complexes or social spheres. Communities in networks often overlap such that nodes simultaneously belong to several groups. Meanwhile, many networks are known to possess multi-scale, hierarchical organisation, where communities are recursively grouped into a hierarchical structure. However, the fact that many real networks have communities with pervasive overlap, where each and every node belongs to more than one group, has the consequence that a global hierarchy of nodes cannot capture the relationships between overlapping groups. Here we reinvent communities as groups of links rather than nodes and show that this unorthodox approach successfully reconciles the antagonistic organising principles of overlapping communities and hierarchy. In contrast to the existing literature, which has entirely focused on grouping nodes, link communities naturally incorporate overlap while revealing hierarchical organisation. We find biologically relevant link communities in protein-protein interaction and metabolic networks and show that a large social network contains hierarchically organised, community structures spanning inner-city to regional scales while maintaining pervasive overlap. Our results imply that link communities are fundamental building blocks that reveal overlap and multi-scale hierarchical organisation in networks to be two aspects of the same phenomenon.

M. Barbosa, Cornell University
Talk preference: Session B or A
Title: A linear relation between solvation free energy and the potential of mean force in a lattice model of fluid
Coauthors: Marco Barbosa* and Benjamin Widom, Cornell University
Abstract: The solvation of an apolar solute in a solvent medium represented by a simple lattice gas is investigated in the Bethe lattice. We follow a previous study on the hydrophobic effect [Widom et al., Phys. Chem. Chem. Phys. (2003)] which found an almost linear relation between the solvation free energy and the potential of mean force at contact. In the mean field limit, the model studied here obeys an exact linear relation for those two quantities.

A. Baule, Rockefeller University
Talk preference: Session A
Title: Path integral approach to random motion with nonlinear friction
Coauthors: E. G. D. Cohen, H. Touchette
Abstract: Using a path integral approach, we derive an analytical solution of a nonlinear and singular Langevin equation, which has been introduced previously by P.-G. de Gennes as a simple phenomenological model for the stick-slip motion of a solid object on a vibrating horizontal surface. We show that the optimal (or most probable) paths of this model can be divided into two classes of paths, which correspond physically to a sliding or slip motion, where the object moves with a non-zero velocity over the underlying surface, and a stick-slip motion, where the object is stuck to the surface for a finite time. These two kinds of basic motions underlie the behavior of many more complicated systems with solid/solid friction and appear naturally in de Gennes' model in the path-integral framework.

J. Bechhoefer, Simon Fraser University
Talk preference: Session B
Title: Defects in DNA Replication: A tale of two regimes?
Abstract: In higher organisms, a vast amount of DNA must be replicated in a short time. The cell thus initiates replication at many distinct "origins" that are dispersed throughout the genome. After initiation, "forks" (domain boundaries) spread out bi-directionally from the origin site until they eventually coalesce with another fork. Unfortunately, defects along the DNA (such as single-strand DNA lesions or double-strand breaks) can temporarily block replication forks. We propose a formalism to model the effects of defects such as fork blocks and find that there are two qualitatively different regimes: a low-defect-density regime, where defects perturb replication locally and a high-defect-density regime, where replication changes qualitatively. Experimental evidence suggests that in normal cells, the defect density is just below the crossover between the two regimes.

R. Bundschuh, Ohio State University
Talk preference: Session A or B(Session C only if its early)
Title: Flexibility of short DNA
Coauthors: Robert Forties, Ralf Bundschuh*, and Michael Poirier Abstract: Protein-bound duplex DNA is often bent or kinked. Yet, quantification of intrinsic DNA bending that might lead to such protein interactions remains enigmatic. DNA cyclization experiments have indicated that DNA may form sharp bends more easily than predicted by the established worm-like chain (WLC) model. One proposed explanation suggests that local melting of a few base pairs introduces flexible hinges. We have expanded this model to incorporate sequence and temperature dependence of the local melting, and tested it for three sequences at temperatures from 23 degrees C to 42 degrees C. We find that small melted bubbles are significantly more flexible than double-stranded DNA and can alter DNA flexibility at physiological temperatures. However, these bubbles are not flexible enough to explain the recently observed very sharp bends in DNA.

J. England, Lewis-Sigler Institute, Princeton University
Talk preference: Session A or B
Title: An Exactly Solvable Model of Structure from Sequence: The Solution to a Gaussian Folding Problem
Abstract: We suggest a Hamiltonian for a harmonically-bonded polymer with heterogeneous solvophobicity along its length exploring different conformations in a collapsed globule state. By computing the partition function in the physically meaningful parameter regime, we show that the structural ensemble of states can be derived from the eigensystem of a one-dimensional Schroedinger equation for a particle in a potential determined by the heteropolymer's "amino acid sequence."

Y-J. Chen, Cornell University
Talk preference: Session A
Title: Merging theory with experiment: improving the accuracy of scaling theories.
Coauthors: Stefanos Papanikolaou, James P. Sethna (Cornell University), Gianfranco Durin (INRIM and ISI foundation, Torino, Italy), Stefano Zapperi(INFM-CNR, Modena and ISI foundation Torino, Italy)
Abstract: Motivated by the experimental problem of analyzing data of crackling noise collected through a limited field of view[1], we have developed a flexible software environment, SloppyScaling, which fits multi-variable scaling functions to both experimental and simulation data. We've used this to test our proposed two-variable scaling functions against simulations on interface depinning models, enabling experiments to make better predictions. Importance sampling algorithms[2] allow us to estimate exponents with honest error bars[3]and improved confidence. Furthermore, we've discovered its utility as a theorist's playground: it allows us to easily identify corrections to scaling, add them to our theory, and explore crossovers away from well-understood scaling behavior.
    [1]A. Magni, G. Durin, J. P. Sethna and S. Zapperi, "Visualization of avalanches in magnetic thin films: temporal processing", J. Stat. Mech. (2009) P01020
    [2] Ryan Gutenkunst. "Sloppiness, Modeling, and Evolution in Biochemical Networks." PhD thesis, Cornell University, 2007.
    [3] "Bayesian Ensemble Approach to Error Estimation of Interatomic Potentials", Søren L. Frederiksen, Karsten W. Jacobsen, Kevin S. Brown, and James P. Sethna, Phys. Rev. Letters 93, 165501 (2004).

B. Daniels, Cornell University
Talk preference: Session B or A
Title: Statistical mechanics of the DNA supercoiling transition Coauthors: James P. Sethna
Abstract: When overtwisted, DNA forms the same wound coils that are familiar from phone cords and water hoses. Since DNA lives at the nanometer scale, however, it is subject to significant thermal fluctuations, and the machinery of statistical mechanics becomes necessary to accurately describe its behavior. We are specifically interested in the transition that occurs as turns are added to straight fluctuating DNA, when it suddenly nucleates a coiled structure known as a plectoneme. Single molecule experiments have shown that this nucleation is thermally activated, with hopping near the transition between states with and without a plectoneme. Combining techniques from polymer physics and transition state theory to characterize motion over the free energy barrier at the transition, we aim to explain why the experimentally measured rate of hopping is so slow, happening on the human-sized timescale of about 1 Hz. Positions_Wanted: Seeking a postdoctoral position. I have experience in applying statistical physics ideas to biological problems, including DNA supercoiling and many-parameter complex network models. See http://www.physics.cornell.edu/~bdaniels

D. David-rus, Ecole Normale Superiere, Paris
Talk preference: Session B or A
Title: Understanding regulation of the states of DARPP-32 phosphorylations- a stochastic approach
Abstract: In this work, I study a stochastic process that describes DARPP-32 phosphorylation states. I solve the steady state of the model for a particular choice of states and transition rules that describes DARPP-32 phosphorylation states when the number of phosphorylated states is very large. Regulation of this states provides a mechanism for integrating information arriving at dopaminoceptive neurons, in multiple brain regions, via a variety of neurotransmitters. DARPP-32 has been established as a crucial mediator of the biochemical,electrophysiological, transcriptional, and behavioral effects of dopamine. Understanding the nature of dopaminergic neurotransmission is important for understanding Parkinson disease, Huntington's chorea, and virtually all antischizophrenic drugs that are influenced by dopamine receptors.
I thank you for suggestions in the analytical calculations to Prof. J.L.Lebowitz and Prof. Larry Shepp.

A. Davidson, Rutgers University
Talk preference: Session A or B
Title: Energy-dependent and pathway-specific transitions of RNA levles in budding yeast induced by glucose-galactose shift
Coauthors: Davidson*, A., Chin, P., Patel, D., Shah, R., So, K., and Ji, S.
Abstract: When glucose is replaced with galactose in the growth medium, budding yeast cells exhibit genome-wide changes in the intracellular levels of RNA molecules, which have been measured using DNA microarrays at 0, 5, 120, 360, 450 and 850 minutes after the nutritional shift (Garcia-Martinez et al., Mol. Cell 15, 303-313, 2004). The results can be displayed as a set of over 6000 RNA trajectories, each of which in turn can be represented as a point in a 6-dimensional RNA concentration space (6DRCS). The distance between any pair of the points in 6DRCS is thought to be inversely proportional to the phenotypic similarity between the paired RNA molecules. When all the possible phenotypic distances of a given metabolic pathway are plotted as a frequency vs. phenotypic distance histogram (FvsPDH), a relatively smooth distribution was obtained which was found to fit a Planck radiation law-type equation (see S. Ji and K. So, this Conference). The FvsPDH of the glycolytic pathway can be divided into two FvsPDHs, one belonging to the energy-poor early phase (0 to 120 min) and the other belonging to the energy-rich late phase (360-850 min). The former FvsPDH was found to contain more RNA pairs separated by short phenotypic distances than the latter, indicating that glycolytic pathway is more active in the late phase than in the early phase. Almost exactly opposite changes were found for the oxphos pathway. Thus, the energy- dependent and metabolic pathway-specific RNA level changes in whole cells measured with DNA microarrays can be conveniently analyzed in terms of the changes in the position and the shape of the FvsPDHs which may serve a role in whole-cell metabolism that is analogous to the role of atomic spectra which revolutionized physics in the early decades of the 20th century.

F. Family, Emory University
Talk preference:Session A
Title: A Statistical Physics Look at Macular Degeneration
Coauthors: Fereydoon Family*, Hans Grossniklaus, Miguel Arizmendi, Karina Mazitello, James Glazier
Abstract: Age-related macular degeneration (AMD) is the leading cause of blindness in the adult population. Choroidal neovascularization, which is the abnormal growth of blood vessels in the choroidal region, is the most common cause of AMD. CNV is produced with age by accumulation of residual material in the retinal pigment epithelium cells (RPE). With time, incompletely degraded membrane material build up in the RPE in the form of lipofuscin, cause abnormal growth of blood vessels that break through the Bruch's membrane, and raise the macula and eventually lead to blindness. The fact that a number of far from equilibrium dynamical processes are involved in the formation and growth of AMD makes this a rich field for application of many techniques of statistical mechanics. I will give some examples of the open problems and mention the results of a deposition and aggregation model of lipofuscin formation in the RPE cells, as well as both two and three-dimensional models of the formation of CNV, that we have recently developed.

B. Fernandez, CNRS & NYU
Talk preference: Session A or B
Title: Athermal dynamics of strongly coupled stochastic three-state oscillators
Coauthors: Bastien Fernandez* and Lev Tsimring
Abstract: We study the collective behavior of a globally coupled ensemble of N cyclic stochastic three-state systems with rates of transition from state i-1 to state i proportional to the number of systems already in state i. While the mean field theory predicts only decaying oscillations in this system, direct numerical simulations indicate that the mean field exhibits stochastic oscillations even in the large N limit. The order parameter characterizing the level of synchrony among oscillators, increases monotonously with the coupling strength. We derive the exact solution of the full master equation for the stationary probability distribution and find the analytical expression for the order parameter.

U. Harbola, University of California, San Diego
Talk preference: Session A or B
Title: Fluctuation Theorems and Electron Counting Statistics
Coauthors: Shaul Mukamel, Massimiliano Esposito
Abstract: Fluctuation Theorems (FTs) describe universal properties of non-equilibrium fluctuations. We present a unified approach to the FTs by introducing a two-point measurement process.Application to electron counting statistics are discussed.

S. Harvey, Georgia Institute of Technology
Talk preference: Session A or B
Title: The Entropic Penalty of Confining a Polymer into a Very Small Space
Coauthors: Mark R. Smyda
Abstract: Determination of the entropic penalty of confinement of a chain polymer into a very small space is an important unsolved problem in polymer statistical mechanics. We present a method for calculating ΔS for the confinement of an elastic polymer of persistence length P in the long-chain limit, when volume exclusion effects are ignored. We consider three geometries: (1) parallel planes separated by a distance d; (2) a circular tube of diameter d; and (3) a sphere of diameter d. We provide results over the range of d/P from 0.01 to 100.

C. Henley, Cornell University
Talk preference: Session B or C
Title: Possible mechanisms to determine macroscopic left-right asymmetry in animals and plants
Abstract: How can systematic L/R asymmetry of the body plan be brought up from molecular to macroscopic scales? Basic symmetry principles suggest that the usual ``biological'' mechanisms Diffusion and gene regulation are insufficient to implement the "right-hand rule"; physical mechanisms involving the cytoskeleton seem always to be involved. I will mention two kinds of dynamic arrays of fibers in the cell wall that might be involved, respectively, in the handedness of snails and of growing plants.

S. Ji, Rutgers University
Talk preference: Session A or B
Title: The universal law of thermal transistions applicable to blackbody radiation, single-molecule enzymology and whole-cell metabolism Coauthors: Sungchul Ji* and Kenneth So
Abstract

K. Korolev, Harvard University
Talk preference: Session A or B
Title: Genetic Waves under Strong Noise
Coauthors: Oskar Hallatschek
Abstract: The rate at which a new advantageous mutation takes over the population is an important quantity because it sets the time scale of evolutionary change. In a spatially extended habitat, a beneficial mutation creates a genetic wave expanding from the location where the mutation occurred. The spreading of beneficial mutations or infectious diseases can usually be tackled when randomness is small, but chance often plays a major role in nature. We studied how an advantageous mutation spreads in the limit of strong number fluctuations. Mathematically, this spreading is described by the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov equation, a classical model in nonequilibrium physics, which is also used in chemical kinetics, ecology, and nuclear physics. We developed a powerful technique to analyze systems in which the effects of chance dominate deterministic forces and calculated the expansion velocity in one and two spatial dimensions. We also analyzed the population structure in the hybrid zone, where both the mutants and the wild type are present. Instead of a stationary, smooth transition region predicted by the classical no-noise approximation, we found non-trivial dynamics of a few rugged kinks that diffuse, give birth by division, and annihilate upon encounter.

B. Machta, Cornell University
Talk preference: Session A or C
Title: Criticality in Biological Membranes
Coauthors: Stefanos Papanikolaou, Sarah Veatch, Jim Sethna
Abstract: All cells are surrounded by a lipid bilayer membrane. This two-dimensional liquid is composed of thousands of lipids and proteins and is home to a host of biological functions. Recent work in giant plasma membrane vesicles (GPMVs) isolated from living cells demonstrates that these GPMVs can be tuned with a single parameter (temperature) to liquid-liquid criticality in the 2D Ising universality class, not far from in vivo temperatures [1,2]. Criticality requires the fine-tuning of two parameters suggesting important biological function, and its presence resolves many of the paradoxes associated with putative lipid rafts. Here we look at the significance of this proximity to criticality, both for understanding surprising features in membrane experiments, as well as from a more theoretical perspective. Why would a cell want to have a nearly critical membrane?

S. Maslov, Brookhaven National Laboratory
Talk preference: Session A or B
Title: "Home depot" model of evolution of prokaryotic metabolic networks and their regulation
Coauthors: Sandeep Krishna,Tin-Yau Pang, Kim Sneppen
Abstract: It has been reported [1] that in prokaryotes the number of transcription factors scales approximately quadratically with the total number of genes. As a consequence the fraction of transcriptional regulators among all genes in small bacterial genomes (< 500 genes) is less than 0.5%, while in large genomes (~10,000 genes) it reaches as high as 10%.
We recently proposed [2] a general explanation of this empirical scaling law and illustrated it using a simple model in which metabolic and regulatory networks co-evolve together. In this model prokaryotic organisms acquire new metabolic functions by the virtue of horizontal gene transfer of entire co-regulated metabolic pathways from a shared gene pool (the "universal metabolic network") followed by removal of redundant enzymes. This process can be compared to a homeowner buying a tool set from a hardware store (hence our "Home Depot" metaphor) and later returning duplicate or unnecessary items.
We view the full repertoire of metabolic enzymes (or more generally any non-regulatory proteins) encoded in the genome of an organism as its collection of tools. Adapting to a new environmental condition (e.g. learning to use a new nutrient source) involves acquiring new enzymes as well as reusing some of the enzymes/tools that are already encoded in the genome. As the toolbox of an organism grows larger, it can reuse its existing tools more often, and thus needs to acquire fewer new enzymes to master each new regulated task. From this analogy it follows that, in general, the number of regulators in an organism should scale faster than linearly with its total number of proteins.
Our model faithfully reproduces the empirically observed [1] quadratic scaling between these two numbers. Furthermore, the distribution of lengths of co-regulated pathways in our model approximately agrees with that in real-life metabolic network of E coli. Thus, the toolbox analogy provides a conceptual explanation for the empirically observed broad distribution of regulon sizes. I will describe several possible regulatory architectures ensuring proper coordination of activity of metabolic pathways with each other. It remains to be determined which of them (if any) are realized in real-life prokaryotes.
References:
    [1] E van Nimwegen, "Scaling laws in the functional content of genomes", Trends Genet 19:479-84 2003.
    [2] S Maslov, S Krishna, T Y Pang, K Sneppen, "Toolbox model of evolution of prokaryotic metabolic networks and their regulation", PNAS 106, 9743-9748 2009.

J. Menche, Max Planck Institute
Talk preference: Session A or B
Title: Activity patterns on scale-free networks
Coauthors: Angelo Valleriani, Reinhard Lipowsky
Abstract: We study the activity patterns of a generic two-state system on scale-free networks. The states of neighboring vertices interact according to a simple majority rule that is equivalent to Glauber dynamics at zero-temperature in Ising spin systems. On uncorrelated networks, this dynamics only exhibits two stable fixed points, where all vertices are in the same state. This situation changes when correlations between the degrees of adjacent vertices are introduced: With increasing correlations a growing number of additional attractors emerges. Most attractors are found in maximally correlated network configurations. We characterize the properties of these attractors in terms of the underlying network structure and give estimates for their total number. In networks with positive correlations the number of attractors grows with network size. This is not the case in networks with negative correlations, where a maximal number of attractors is reached at intermediate network sizes.

B. Miller, Texas Christian University
Talk preference: Session A or B
Title: A Minimal Model for the Study of Polychronous Groups
Coauthors: Bruce N. Miller * and Willard Maier
Abstract: The concept of a polychromous group in a neural network was introduced by Izhikevich as a model to facilitate learning. Here we present a minimal model of polychronous groups. The model is computationally efficient and allows the study of polychronous groups independent of specific neuron models.Computational experiments were performed with the model in one- and two-dimensional neural architectures to determine the dependence of the number of polychronous groups on various connectivity options. The possibility of using polychronous groups as computational elements will also be discussed.

S. Mishra, Syracuse University
Talk preference: Session A or B
Title: Pattern formation and traveling bands in dense layers of self-propelled rods
Coauthors: Aparna Baskaran, M. Cristina Marchetti
Abstract: A collection of interacting self-propelled (SP) hard rods in two dimensions can be used as a minimal model for a variety of active systems, including bacterial suspensions and vibrated granular layers. Hydrodynamics equations for SP rods were derived in A. Baskaran et al.[PRL 101, 268101 (2008)]. In the bulk limit the rods exhibit a mean-field transition from an isotropic to a polarized state at a density \rho_c . Linear stability analysis indicates that the uniform polarized (hence moving) state in unstable above a critical value v_0 ≃ (\rho_0 − \rho_c )^{-1} of the self-propulsion speed, with \rho_0 the mean density of rods. The numerical solution of the full nonlinear equations shows that above a critical value v_c (\rho_0 ), the system phase separates into traveling bands of ordered and disordered bands. The ordered bands are polarized along the direction of motion, transverse to the long direction of the bands. Similar patterns have been obtained from numerical studies of Vicsek-type models and in actin motility assays.

I. Nemenman, Emory University
Talk preference: Session A or B
Title: Simplicity of Completion Time Distributions of Kinetic Proofreading-Like Biochemical Process
Abstract: Biochemical processes typically involve huge numbers of individual reversible steps, each with its own dynamical rate constants. For example, kinetic proofreading processes rely upon numerous sequential reactions in order to guarantee the precise construction of specific macromolecules. Here we show that, for a wide range of parameters, as the system size grows, the completion time for such processes attains a simple, almost universal distribution: it becomes either deterministic or exponentially distributed, with a very narrow transition between the two regimes. These findings suggest not only that one may not be able to understand individual elementary reactions from completion time observations, but also that such understanding may be unnecessary.

M. Palassini, University of Barcelona
Talk preference: Session A
Title: Delay and noise in negative-feedback genetic regulatory loops
Coauthors: Marta Dies
Abstract: Stochastic effects are important in gene regulation due to the small number of molecules involved. Another important, often neglected aspect of regulatory dynamics is the large separation of time scales between slow processes, such as transcription and translation, and fast processes, such as protein degradation. We consider a generic birth-and-death stochastic model of a negative- feedback loop, in which the separation of time scales is incorporated via an explicit delay in one of the loop arms. We show, both from exact simulation of the delayed Master Equation and from a Van Kampen volume expansion, that noise-sustained oscillations occurs quite generally in such a model. We propose that this mechanism might explain the experimentally observed temporal oscillations in the concentration levels of proteins p53 and Mdm2 in cells subjected to DNA damage, as well as similar oscillations observed in other genetic regulatory negative-feedback loops.

S. Papanikolaou, Cornell Univesity
Talk preference: Session A or B
Title: Beyond scaling: The average avalanche shape
Coauthors: Felipe Bohn, Rubem Luis Sommer, Gianfranco Durin, Stefano Zapperi and James P. Sethna
Abstract: Universality, scaling, and the renormalization group claim to predict all be- havior on long length and time scales asymptotically close to critical points. In practice, large simulations and heroic experiments have been needed to un- ambiguously test and measure the critical exponents and scaling functions. We announce here the measurement and prediction of universal corrections to scaling, applied to the temporal average shape of Barkhausen noise avalanches. We bypass the confounding factors of time-retarded interactions (eddy currents) by measuring thin permalloy films, and bypass thresholding effects and amplifier distortions by applying Wiener deconvolution. We show experimental shapes that are approximately symmetric, and measure the leading corrections to scal- ing. We solve a mean-field theory for the magnetization dynamics and calculate the relevant demagnetizing-field correction to scaling, showing qualitative agree- ment with the experiment. In this way, we move toward a quantitative theory useful at smaller time and length scales and farther from the critical point.

D. Sisan, NIST
Talk preference: Session B or A
Title: Event ordering in live cell imaging determined from temporal cross correlation asymmetry
Coauthors: Defne Yarar, Clare M. Waterman, and Jeffrey S. Urbach
Abstract: We use the temporal asymmetry of the cross-correlation function to determine the temporal ordering of spatially localized cellular events in live cell multi-channel fluorescence imaging. Temporal ordering, a term commonly used in cell biology, requires a nonequilibrium biochemical flux, which is quantifiable through the cross-correlation asymmetry. I'll briefly describe how the approach was applied to extract the temporal ordering of three proteins in the endocytic pathway: actin, sorting nexin 9, and clathrin.

V. Tkachenko, Ben-Gurion University of the Negev, Israel
Talk preference: Session A
Title: An inverse problem for 1d ordinary differential operator of order 4
Abstract

D. Rabson, University of South Florida
Talk preference: Session A
Title: Sometimes the Noise is the Signal
Coauthors: Chun-Min Lo, Douglas Lovelady
Abstract: Since 1984, electric cell-substrate impedance sensing (ECIS) has been used to monitor cell behavior in culture and has proven sensitive to morphological changes and cell mobility. Several authors have associated fluctuations in the measured impedance with cellular micromotion; however we are unaware of any previous work applying statistical techniques in order to distinguish two different cell types. We have demonstrated a method for distinguishing cancerous from non-cancerous cultures of human ovarian surface epithelial cells [1]; applying similar ideas, we have also determined the presence and concentration of the toxin cytochalisin B in cultures of 3T3 fibroblasts at levels lower than the detection thresholds of other techniques [2]. Measures of short-time and long-time correlation confirm that the noise from non-cancerous cultures has a higher degree of temporal order, order which we argue, based on a statistical-mechanical model, might arise from greater coordination of motion between healthy cells than between cancerous ones. The moral of the story: apply easy statistics to a field where people haven't previously!
    [1] D.C. Lovelady, T.C. Richmond, A.N. Maggi, C.-M. Lo, D.A. Rabson, Phys. Rev. E 76, 041908 (2007).
    [2] D.C. Lovelady, J. Friedman, S. Patel, D.A. Rabson, C.-M. Lo, Biosensors and Bioelectronics 24, 2250 (2009).

S. Redner, Boston University
Talk preference: Session A
Title: Distribution of Species Body Masses
Coauthors: A. Clauset
Abstract: We present a model for the evolution of body masses of related species, in which mass M evolves by speciation-driven branching, multiplicative diffusion, and an extinction probability that increases weakly with mass, leading to a convection-diffusion equation for ln M. The resulting steady-state behavior agrees well with empirical data on recent terrestrial mammals, and the time-dependent behavior also agrees with data on extinct mammal species between 95-50 Myr ago.

T. Reichenbach, Rockefeller University
Talk preference: Session A or B
Title: A ratchet mechanism for low-frequency hearing in mammals
Coauthors: A. J. Hudspeth
Abstract: The sensitivity and frequency selectivity of hearing result from tuned amplification by an active process in the mechanoreceptive hair cells. The nature of the active process in the mammalian cochlea is intensely debated, for outer hair cells exhibit two forms of mechanical activity, active hair-bundle motility and membrane-based electromotility. Here we show theoretically that active hair-bundle motility and electromotility can together implement an efficient mechanism for amplification that functions like a ratchet: sound-evoked forces acting on the basilar membrane are transmitted to the hair bundles while electromotility decouples the active hair-bundle forces from the basilar membrane. Through a combination of analytical and computational techniques we demonstrate that the ratchet mechanism can naturally account for a variety of unexplained experimental observations from low-frequency hearing.

O. Sariyer, Koc University
Talk preference: Session A
Title: Charge-Ordered Phases of the $d = 3$ Spinless Falicov-Kimball Model: Renormalization-Group Theory
Coauthors: Ozan S. Sariyer*, Michael Hinczewski, and A. Nihat Berker
Abstract: The global phase diagram of the spinless Falicov-Kimball model in $d = 3$ spatial dimensions has been obtained by renormalization-group theory. This global phase diagram exhibits five distinct phases. Four of these phases are charge-ordered (CO) phases, in which the electron hopping strength diverges under repeated renormalizations. In the other ($\delta$) phase, the hopping strength vanishes under repeated renormalizations. The phase boundaries are second-order, except for an intermediate-temperature regime, where a first-order phase boundary between two CO phases occurs. The first-order phase boundary is delimited by bicritical points. The cross-sections of the global phase diagram with respect to the chemical potentials of the localized and mobile electrons, at all representative temperatures and hopping strengths, are calculated and exhibit four distinct topologies.

E.-M. Schoetz, Princeton University
Talk preference: Session A
Title: Dynamics of asexual reproduction in flatworms
Coauthors: Jared Talbot and Joern Dunkel
Abstract: Flatworms can reproduce by transverse fission (i.e. they simply split themselves in two) because they can regenerate the missing body parts. Naively, one would think that this kind of reproduction could be captured by models for cell growth in bacteria or other simple organisms. However, we find that there is much more to the story by monitoring >10 generations of individuals as well as the behavior of worm populations under different environmental conditions, such as temperature, feeding frequency and crowding.

A. Toom, UFPE, Brazil
Talk preference: Session A or B
Title: Non-Ergodicity and Growth Are Compatible for 1-D local Interaction
Coauthors: Alex D. Ramos (UFPE)
Abstract: We present results of Monte Carlo simulation and chaos approximation of a class of Markov processes. Each of their states can be written as a finite or infinite in both direction sequence of pluses and minuses As continuous time goes on,our sequence undergoes the following three types of local transformations: The first one, called flip, changes any minus into plus and any plus into minus with a rate beta. Another, called annihilation, eliminates two neighbor components with a rate alpha whenever they are in different states. The third one, called mitosis, doubles any component with a rate gamma. All of them occur at any place of the sequence independently. Our simulations and approximations suggest that with approprate positive alpha, beta and gamma this process has the following two properties. Growth: In the finite case, as the process goes on, the length of the sequence tends to infinity with a probability, which tends to 1 as the length of the initial sequence tends to infinity. Non-ergodicity: the infinite system is ergodic and the finite system keeps most of the time at two extremes, occasionally swinging from one to the other.

M. Transtrum, Cornell University
Talk preference: Session A or B
Title: Differential Geometric Approach to Fitting Data
Coauthors: Benjamin Machta, James Sethna
Abstract: Fitting nonlinear functions to experimental data can be a very difficult task. Standard algorithms are usually unreliable to automatically find fits, often requiring many adjustments by hand in order to converge. By considering the manifold of model predictions in data space, we find that the ideal path an algorithm should follow is a geodesic. The standard Levenberg-Marquardt algorithm approximately follows a geodesic path, but will often fail to converge because the geodesic intersects the boundaries of the manifold. When this happens, the algorithm "evaporates" parameters, pushing them to unphysical values, without finding a good fit. Typical multi-parameter models have boundaries with a hierarchy of widths, forming a long, narrow hyper-ribbon which is difficult for standard algorithms to navigate. We explain the origin of these boundaries in terms of the analyticity of the fitting function and suggest that they should guide the design of algorithms to improve the fitting process.

R. Vandiver, Bryn Mawr College
Talk preference: Session A or B
Title: On the mechanical stability of growing arteries
Coauthors: Alain Goriely
Abstract: Arteries are modeled, within the framework of nonlinear elasticity, as incompressible two-layer cylindrical structures that are residually stressed through differential growth. These structures are loaded by an axial force, internal pressure and have nonlinear, anisotropic, hyperelastic response to stresses. Parameters for this model are directly related to experimental observations. The mechanical stability of growing arteries and the role of residual stresses is investigated. It is shown that residual stress lowers the critical internal pressure leading to buckling and that a reduction of axial loading may lead to a buckling instability which may eventually lead to arterial tortuosity.

B. Vollmayr-Lee, Bucknell University
Talk preference: None given
Title: Coarsening Dynamics in a Chaotic Flow
Coauthors: *Benjamin Vollmayr-Lee, Bucknell and Daniel Beller, Brandeis
Abstract: We study the phase separation dynamics of a binary fluid system advected by a chaotic flow. The chaotic flow stretches and rips apart domains, and thus competes with the thermodynamic coarsening to yield a nonequilibrium steady state. We calculate local measures such as the finite-time Lyapunov exponent field and the local free energy density to characterize this steady state.

R. Weinkamer, Max Planck Institute of Colloids and Interfaces, Department of Biomaterials
Talk preference: Session A
Title: Control of bone remodeling
Coauthors: M. Rusconi, A. Valleriani, J.W.C. Dunlop, J. Kurths
Abstract: Trabecular bone is the network-like bone built from struts called trabeculae found inside vertebrae and close to joints. This spongy bone is continuously renewed during life by means of resorption and deposition of bone packets from and onto its surface. We developed a stochastic model, allowing us to extract information about the control of this remodeling process based on experimentally measured frequency distributions of the thickness of the trabeculae.

J. Xing,Virginia Tech.
Talk preference: Session A or B
Title: Mapping between stochastic dissipative and Hamiltonian systems
Abstract: Biological systems are away from equilibrium. Here we prove that a system described by stochastic differential equations can be mapped to a Hamiltonian system, at least for the case the stationary distribution exists. We used the result to obtain a generalized fluctuation-dissipation relation, and to derive the Zwanzig-Mori projection formula for general non-Hamiltonian systems.

Y. Zhang, Fundan University/University of Maryland
Talk preference: Session A
Title: A tug of war model for organelle transport can display three stable steady states
Abstract: Motion of organelles and vesicles moved by motor proteins in cells can be modeled by a tug-of-war model developed by M¨¹ller, Klumpp and Lipowsky [1]. By detailed theoretical analysis, we find that this model can, depending on the single-motor parameters and external force, exhibit one, two or three stable steady states [2]. The steady state motion of the cargo is determined by the initial numbers of the motors bound to the track. The three states correspond, respectively, to the cargo moving to the right, to the left or remaining stationary. Thus our study indicates that the possible motions of the cargo are determined by the intrinsic parameters of the tug-of-war model, while the final motion is determined by the initial conditions. Monte Carlo simulations confirm that our results are accurate when there are a large number of motors. When the number of motors is small, the cargo motion may change from one steady state to another. It is planned to study the transition times between the different steady states.
    [1] M. J. I. Muller, S. Klumpp and R. Lipowsky, PNAS, 105, 4609 (2008).
    [2] Y. Zhang, Phys. Rev. E. 79, 061918 (2009).

R. Zia, Virginia Tech
Talk preference: Session A or B
Title: Convection Cells driven by Spontaneous Symmetry Breaking
Coauthors: M.F.J. Pleimling and B. Schmittmann
Abstract: TBA
________________________________