BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

James P. Sethna

James P. Sethna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.mtrl-scicond-mat.softphysics.data-ancond-mat.dis-nnBiological Physicscond-mat.supr-conphysics.comp-phmath-phmath.MPphysics.acc-phcond-matcond-mat.mes-hallQuantitative MethodsBiomoleculesMolecular Networksphysics.flu-dynPopulations and EvolutionTissues and OrgansCell Behaviorcond-mat.str-elmath.APmath.NAmath.OCnlin.AOnlin.CDnlin.CGnlin.PSOther Quantitative Biologypatt-solphysics.soc-phSocial and Information NetworksSubcellular Processes

Catalog footprint

What is connected

54works
33topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2022arXiv

Avalanches and deformation in glasses and disordered systems

In this chapter, we discuss avalanches in glasses and disordered systems, and the macroscopic dynamical behavior that they mediate. We briefly review three classes of systems where avalanches are observed: depinning transition of disordered interfaces, yielding of amorphous materials, and the jamming transition. Without extensive formalism, we discuss results gleaned from theoretical approaches -- mean-field theory, scaling and exponent relations, the renormalization group, and a few results from replica theory. We focus both on the remarkably sophisticated physics of avalanches and on relatively new approaches to the macroscopic flow behavior exhibited past the depinning/yielding transition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Universal scaling for disordered viscoelastic matter II: Collapses, global behavior and spatio-temporal properties

Disordered viscoelastic materials are ubiquitous and exhibit fascinating invariant scaling properties. In a companion article, we have presented comprehensive new results for the critical behavior of the dynamic susceptibility of disordered elastic systems near the onset of rigidity. Here we provide additional details of the derivation of the singular scaling forms of the longitudinal response near both jamming and rigidity percolation. We then discuss global aspects associated with these forms, and make scaling collapse plots for both undamped and overdamped dynamics in both the rigid and floppy phases. We also derive critical exponents, invariant scaling combinations and analytical formulas for universal scaling functions of several quantities such as transverse and density responses, elastic moduli, viscosities, and correlation functions. Finally, we discuss tentative experimental protocols to measure these behaviors in colloidal suspensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of Magnetic Vortex Dissipation in Sn-Segregated Boundaries in Nb$_3$Sn Superconducting RF Cavities

We study mechanisms of vortex nucleation in Nb$_3$Sn Superconducting RF (SRF) cavities using a combination of experimental, theoretical, and computational methods. Scanning transmission electron microscopy (STEM) image and energy dispersive spectroscopy (EDS) of some Nb$_3$Sn cavities show Sn segregation at grain boundaries in Nb$_3$Sn with Sn concentration as high as $\sim$35 at.\% and widths $\sim$3 nm in chemical composition. Using ab initio calculations, we estimate the effect excess tin has on the local superconducting properties of the material. We model Sn segregation as a lowering of the local critical temperature. We then use time-dependent Ginzburg-Landau theory to understand the role of segregation on magnetic vortex nucleation. Our simulations indicate that the grain boundaries act as both nucleation sites for vortex penetration and pinning sites for vortices after nucleation. Depending on the magnitude of the applied field, vortices may remain pinned in the grain boundary or penetrate the grain itself. We estimate the superconducting losses due to vortices filling grain boundaries and compare with observed performance degradation with higher magnetic fields. We estimate that the quality factor may decrease by an order of magnitude ($10^{10}$ to $10^9$) at typical operating fields if 0.03\% of the grain boundaries actively nucleate vortices. We additionally estimate the volume that would need to be filled with vortices to match experimental observations of cavity heating.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reaction rates and the noisy saddle-node bifurcation: Renormalization group for barrier crossing

Barrier crossing calculations in chemical reaction-rate theory typically assume that the barrier is large compared to the temperature. When the barrier vanishes, however, there is a qualitative change in behavior. Instead of crossing a barrier, particles slide down a sloping potential. We formulate a renormalization group description of this noisy saddle-node transition. We derive the universal scaling behavior and corrections to scaling for the mean escape time in overdamped systems with arbitrary barrier height. We also develop an accurate approximation to the full distribution of barrier escape times by approximating the eigenvalues of the Fokker-Plank operator as equally spaced. This lets us derive a family of distributions that captures the barrier crossing times for arbitrary barrier height. Our critical theory draws links between barrier crossing in chemistry, the renormalization group, and bifurcation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Visualizing probabilistic models in Minkowski space with intensive symmetrized Kullback-Leibler embedding

We show that the predicted probability distributions for any $N$-parameter statistical model taking the form of an exponential family can be explicitly and analytically embedded isometrically in a $N{+}N$-dimensional Minkowski space. That is, the model predictions can be visualized as control parameters are varied, preserving the natural distance between probability distributions. All pairwise distances between model instances are given by the symmetrized Kullback-Leibler divergence. We give formulas for these intensive symmetrized Kullback Leibler (isKL) coordinate embeddings, and illustrate the resulting visualizations with the Bernoulli (coin toss) problem, the ideal gas, $n$ sided die, the nonlinear least squares fit, and the Gaussian fit. We highlight how isKL can be used to determine the minimum number of parameters needed to describe probabilistic data, and conclude by visualizing the prediction space of the two-dimensional Ising model, where we examine the manifold behavior near its critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Visualizing probabilistic models: Intensive Principal Component Analysis

Unsupervised learning makes manifest the underlying structure of data without curated training and specific problem definitions. However, the inference of relationships between data points is frustrated by the `curse of dimensionality' in high-dimensions. Inspired by replica theory from statistical mechanics, we consider replicas of the system to tune the dimensionality and take the limit as the number of replicas goes to zero. The result is the intensive embedding, which is not only isometric (preserving local distances) but allows global structure to be more transparently visualized. We develop the Intensive Principal Component Analysis (InPCA) and demonstrate clear improvements in visualizations of the Ising model of magnetic spins, a neural network, and the dark energy cold dark matter (ΛCDM) model as applied to the Cosmic Microwave Background.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ginzburg-Landau theory of the superheating field anisotropy of layered superconductors

We investigate the effects of material anisotropy on the superheating field of layered superconductors. We provide an intuitive argument both for the existence of a superheating field, and its dependence on anisotropy, for $κ= λ/ ξ$ (the ratio of magnetic to superconducting healing lengths) both large and small. On the one hand, the combination of our estimates with published results using a two-gap model for MgB${}_2$ suggests high anisotropy of the superheating field near zero temperature. On the other hand, within Ginzburg-Landau theory for a single gap, we see that the superheating field shows significant anisotropy only when the crystal anisotropy is large and the Ginzburg-Landau parameter $κ$ is small. We then conclude that only small anisotropies in the superheating field are expected for typical unconventional superconductors near the critical temperature. Using a generalized form of Ginzburg Landau theory, we do a quantitative calculation for the anisotropic superheating field by mapping the problem to the isotropic case, and present a phase diagram in terms of anisotropy and $κ$, showing type I, type II, or mixed behavior (within Ginzburg-Landau theory), and regions where each asymptotic solution is expected. We estimate anisotropies for a number of different materials, and discuss the importance of these results for radio-frequency cavities for particle accelerators.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measuring nonlinear stresses generated by defects in 3D colloidal crystals

The mechanical, structural and functional properties of crystals are determined by their defects and the distribution of stresses surrounding these defects has broad implications for the understanding of transport phenomena. When the defect density rises to levels routinely found in real-world materials, transport is governed by local stresses that are predominantly nonlinear. Such stress fields however, cannot be measured using conventional bulk and local measurement techniques. Here, we report direct and spatially resolved experimental measurements of the nonlinear stresses surrounding colloidal crystalline defect cores, and show that the stresses at vacancy cores generate attractive interactions between them. We also directly visualize the softening of crystalline regions surrounding dislocation cores, and find that stress fluctuations in quiescent polycrystals are uniformly distributed rather than localized at grain boundaries, as is the case in strained atomic polycrystals. Nonlinear stress measurements have important implications for strain hardening, yield, and fatigue.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pinning Susceptibility: The effect of dilute, quenched disorder on jamming

We study the effect of dilute pinning on the jamming transition. Pinning reduces the average contact number needed to jam unpinned particles and shifts the jamming threshold to lower densities, leading to a pinning susceptibility, $χ_p$. Our main results are that this susceptibility obeys scaling form and diverges in the thermodynamic limit as $χ_p \propto |ϕ- ϕ_c^\infty|^{-γ_p}$ where $ϕ_c^\infty$ is the jamming threshold in the absence of pins. Finite-size scaling arguments yield these values with associated statistical (systematic) errors $γ_p = 1.018 \pm 0.026 (0.291) $ in $d=2$ and $γ_p =1.534 \pm 0.120 (0.822)$ in $d=3$. Logarithmic corrections raise the exponent in $d=2$ to close to the $d=3$ value, although the systematic errors are very large.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The weirdest martensite: Smectic liquid crystal microstructure and Weyl-Poincaré invariance

Smectic liquid crystals are remarkable, beautiful examples of materials microstructure, with ordered patterns of geometrically perfect ellipses and hyperbolas. The solution of the complex problem of filling three-dimensional space with domains of focal conics under constraining boundary conditions yields a set of strict rules, which are similar to the compatibility conditions in a martensitic crystal. Here we present the rules giving compatible conditions for the concentric circle domains found at two-dimensional smectic interfaces with planar boundary conditions. Using configurations generated by numerical simulations, we develop a clustering algorithm to decompose the planar boundaries into domains. The interfaces between different domains agree well with the smectic compatibility conditions. We also discuss generalizations of our approach to describe the full three-dimensional smectic domains, where the variant symmetry group is the Weyl-Poincaré group of Lorentz boosts, translations, rotations, and dilatations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

"Irregularization" of Systems of Conservation Laws

We explore new ways of regulating defect behavior in systems of conservation laws. Contrary to usual regularization schemes (such as a vanishing viscosity limit), which attempt to control defects by making them smoother, our schemes result in defects which are \textit{more singular}, and we thus refer to such schemes as "irregularizations". In particular, we seek to produce \textit{delta shock} defects which satisfy a condition of \textit{stationarity}. We are motivated to pursue such exotic defects by a physical example arising from dislocation dynamics in materials physics, which we describe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Crossover behavior in interface depinning

We study the crossover scaling behavior of the height-height correlation function in interface depinning in random media. We analyze experimental data from a fracture experiment and simulate an elastic line model with non-linear couplings and disorder. Both exhibit a crossover between two different universality classes. For the experiment, we fit a functional form to the universal crossover scaling function. For the model, we vary the system size and the strength of the non-linear term, and describe the crossover between the two universality classes with a multiparameter scaling function. Our method provides a general strategy to extract scaling properties in depinning systems exhibiting crossover phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Overshoot during phenotypic switching of cancer cell populations

The dynamics of tumor cell populations is hotly debated: do populations derive hierarchically from a subpopulation of cancer stem cells (CSCs), or are stochastic transitions that mutate differentiated cancer cells to CSCs important? Here we argue that regulation must also be important. We sort human melanoma cells using three distinct cancer stem cell (CSC) markers - CXCR6, CD271 and ABCG2 - and observe that the fraction of non-CSC-marked cells first overshoots to a higher level and then returns to the level of unsorted cells. This clearly indicates that the CSC population is homeostatically regulated. Combining experimental measurements with theoretical modeling and numerical simulations, we show that the population dynamics of cancer cells is associated with a complex miRNA network regulating the Wnt and PI3K pathways. Hence phenotypic switching is not stochastic, but is tightly regulated by the balance between positive and negative cells in the population. Reducing the fraction of CSCs below a threshold triggers massive phenotypic switching, suggesting that a therapeutic strategy based on CSC eradication is unlikely to succeed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Scaling theory for the jamming transition

We propose a scaling ansatz for the elastic energy of a system near the critical jamming transition in terms of three relevant fields: the compressive strain $Δϕ$ relative to the critical jammed state, the shear strain $ε$, and the inverse system size $1/N$. We also use $ΔZ$, the number of contacts relative to the minimum required at jamming, as an underlying control parameter. Our scaling theory predicts new exponents, exponent equalities and scaling collapses for energy, pressure and shear stress that we verify with numerical simulations of jammed packings of soft spheres. It also yields new insight into why the shear and bulk moduli exhibit different scalings; the difference arises because the shear stress vanishes as $1/\sqrt{N}$ while the pressure approaches a constant in the thermodynamic limit. The success of the scaling ansatz implies that the jamming transition exhibits an emergent scale invariance, and that it should be possible to develop a renormalization-group theory for jamming.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shielding superconductors with thin films

Determining the optimal arrangement of superconducting layers to withstand large amplitude AC magnetic fields is important for certain applications such as superconducting radiofrequency cavities. In this paper, we evaluate the shielding potential of the superconducting film/insulating film/superconductor (SIS') structure, a configuration that could provide benefits in screening large AC magnetic fields. After establishing that for high frequency magnetic fields, flux penetration must be avoided, the superheating field of the structure is calculated in the London limit both numerically and, for thin films, analytically. For intermediate film thicknesses and realistic material parameters we also solve numerically the Ginzburg-Landau equations. It is shown that a small enhancement of the superheating field is possible, on the order of a few percent, for the SIS' structure relative to a bulk superconductor of the film material, if the materials and thicknesses are chosen appropriately.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sloppiness and Emergent Theories in Physics, Biology, and Beyond

Large scale models of physical phenomena demand the development of new statistical and computational tools in order to be effective. Many such models are `sloppy', i.e., exhibit behavior controlled by a relatively small number of parameter combinations. We review an information theoretic framework for analyzing sloppy models. This formalism is based on the Fisher Information Matrix, which we interpret as a Riemannian metric on a parameterized space of models. Distance in this space is a measure of how distinguishable two models are based on their predictions. Sloppy model manifolds are bounded with a hierarchy of widths and extrinsic curvatures. We show how the manifold boundary approximation can extract the simple, hidden theory from complicated sloppy models. We attribute the success of simple effective models in physics as likewise emerging from complicated processes exhibiting a low effective dimensionality. We discuss the ramifications and consequences of sloppy models for biochemistry and science more generally. We suggest that the reason our complex world is understandable is due to the same fundamental reason: simple theories of macroscopic behavior are hidden inside complicated microscopic processes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Visualization, coarsening and flow dynamics of focal conic domains in simulated Smectic-A liquid crystals

Smectic liquid crystals vividly illustrate the subtle interplay of broken translational and orientational symmetries, by exhibiting defect structures forming geometrically perfect confocal ellipses and hyperbolas. Here, we develop and numerically implement an effective theory to study the dynamics of focal conic domains in smectic-A liquid crystals. We use the information about the smectic's structure and energy density provided by our simulations to develop several novel visualization tools for the focal conics. Our simulations accurately describe both simple and extensional shear, which we compare to experiments, and provide additional insight into the coarsening dynamics of focal conic domains.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fracture strength: Stress concentration, extreme value statistics and the fate of the Weibull distribution

The fracture strength distribution of materials is often described in terms of the Weibull law which can be derived by using extreme value statistics if elastic interactions are ignored. Here, we consider explicitly the interplay between elasticity and disorder and test the asymptotic validity of the Weibull distribution through numerical simulations of the two-dimensional random fuse model. Even when the local fracture strength follows the Weibull distribution, the global failure distribution is dictated by stress enhancement at the tip of the cracks and sometimes deviates from the Weibull law. Only in the case of a pre-existing power law distribution of crack widths do we find that the failure strength is Weibull distributed. Contrary to conventional assumptions, even in this case, the Weibull exponent can not be simply inferred from the exponent of the initial crack width distribution. Our results thus raise some concerns on the applicability of the Weibull distribution in most practical cases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mechanical Properties of Growing Melanocytic Nevi and the Progression to Melanoma

Melanocytic nevi are benign proliferations that sometimes turn into malignant melanoma in a way that is still unclear from the biochemical and genetic point of view. Diagnostic and prognostic tools are then mostly based on dermoscopic examination and morphological analysis of histological tissues. To investigate the role of mechanics and geometry in the morpholgical dynamics of melanocytic nevi, we study a computation model for cell proliferation in a layered non-linear elastic tissue. Numerical simulations suggest that the morphology of the nevus is correlated to the initial location of the proliferating cell starting the growth process and to the mechanical properties of the tissue. Our results also support that melanocytes are subject to compressive stresses that fluctuate widely in the nevus and depend on the growth stage. Numerical simulations of cells in the epidermis releasing matrix metalloproteinases display an accelerated invasion of the dermis by destroying the basal membrane. Moreover, we suggest experimentally that osmotic stress and collagen inhibit growth in primary melanoma cells while the effect is much weaker in metastatic cells. Knowing that morphological features of nevi might also reflect geometry and mechanics rather than malignancy could be relevant for diagnostic purposes

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Collective Motion of Moshers at Heavy Metal Concerts

Human collective behavior can vary from calm to panicked depending on social context. Using videos publicly available online, we study the highly energized collective motion of attendees at heavy metal concerts. We find these extreme social gatherings generate similarly extreme behaviors: a disordered gas-like state called a mosh pit and an ordered vortex-like state called a circle pit. Both phenomena are reproduced in flocking simulations demonstrating that human collective behavior is consistent with the predictions of simplified models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Parameter Space Compression Underlies Emergent Theories and Predictive Models

We report a similarity between the microscopic parameter dependance of emergent theories in physics and that of multiparameter models common in other areas of science. In both cases, predictions are possible despite large uncertainties in the microscopic parameters because these details are compressed into just a few governing parameters that are sufficient to describe relevant observables. We make this commonality explicit by examining parameter sensitivity in a hopping model of diffusion and a generalized Ising model of ferromagnetism. We trace the emergence of a smaller effective model to the development of a hierarchy of parameter importance quantified by the eigenvalues of the Fisher Information Matrix. Strikingly, the same hierarchy appears ubiquitously in models taken from diverse areas of science. We conclude that the emergence of effective continuum and universal theories in physics is due to the same parameter space hierarchy that underlies predictive modeling in other areas of science.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quasi-periodic events in crystal plasticity and the self-organized avalanche oscillator

When external stresses in a system - physical, social or virtual - are relieved through impulsive events, it is natural to focus on the attributes of these avalanches. However, during the quiescent periods in between, stresses may be relieved through competing processes, such as slowly flowing water between earthquakes or thermally activated dislocation flow between plastic bursts. Such unassuming, smooth responses can have dramatic effects on the avalanche properties. Here we report a thorough experimental investigation of slowly compressed Ni microcrystals, covering three orders of magnitude in nominal strain rate, that exhibits unconventional quasi-periodic avalanche bursts and higher critical exponents as the strain rate is decreased. Our analytic and computational study naturally extends dislocation avalanche modeling to incorporate dislocation relaxations and reveals the emergence of the self-organized avalanche oscillator, a novel critical state exhibiting oscillatory approaches toward a depinning critical point. We demonstrate that the predictions of our theory are faithfully exhibited in our experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling theory of continuum dislocation dynamics in three dimensions: Self-organized fractal pattern formation

We focus on mesoscopic dislocation patterning via a continuum dislocation dynamics theory (CDD) in three dimensions (3D). We study three distinct physically motivated dynamics which consistently lead to fractal formation in 3D with rather similar morphologies, and therefore we suggest that this is a general feature of the 3D collective behavior of geometrically necessary dislocation (GND) ensembles. The striking self-similar features are measured in terms of correlation functions of physical observables, such as the GND density, the plastic distortion, and the crystalline orientation. Remarkably, all these correlation functions exhibit spatial power-law behaviors, sharing a single underlying universal critical exponent for each type of dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Senescent Cells in Growing Tumors: Population Dynamics and Cancer Stem Cells

Tumors are defined by their intense proliferation, but sometimes cancer cells turn senescent and stop replicating. In the stochastic cancer model in which all cells are tumorigenic, senescence is seen as the result of random mutations, suggesting that it could represent a barrier to tumor growth. In the hierarchical cancer model a subset of the cells, the cancer stem cells, divide indefinitely while other cells eventually turn senescent. Here we formulate cancer growth in mathematical terms and obtain predictions for the evolution of senescence. We perform experiments in human melanoma cells which are compatible with the hierarchical model and show that senescence is a reversible process controlled by survivin. We conclude that enhancing senescence is unlikely to provide a useful therapeutic strategy to fight cancer, unless the cancer stem cells are specifically targeted

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theoretical Field Limits for Multi-Layer Superconductors

The SIS structure---a thin superconducting film on a bulk superconductor separated by a thin insulating film---was propsed as a method to protect alternative SRF materials from flux penetration by enhancing the first critical field $B_{c1}$. In this work, we show that in fact $B_{c1}$ = 0 for a SIS structure. We calculate the superheating field $B_{sh}$, and we show that it can be enhanced slightly using the SIS structure, but only for a small range of film thicknesses and only if the film and the bulk are different materials. We also show that using a multilayer instead of a single thick layer is detrimental, as this decreases $B_{sh}$ of the film. We calculate the dissipation due to vortex penetration above the $B_{sh}$ of the film, and find that it is unmanageable for SRF applications. However, we find that if a gradient in the phase of the order parameter is introduced, SIS structures may be able to shield large DC and low frequency fields. We argue that the SIS structure is not beneficial for SRF cavities, but due to recent experiments showing low-surface-resistance performance above $B_{c1}$ in cavities made of superconductors with small coherence lengths, we argue that enhancement of $B_{c1}$ is not necessary, and that bulk films of alternative materials show great promise.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Universal scaling function for the two-dimensional Ising model in an external field: A pragmatic approach

We report an effective functional form for the spin-spin correlation function of the 2D Ising model as a function of temperature and field. Although the Ising model has been well studied, no analytical result for the spin-spin correlation function exists for arbitrary magnetic fields and temperatures. We show the validity of our form by comparison with simulations using the Wolff algorithm, and obtain remarkable precision by including analytic corrections to scaling. Given recent interest in comparing biomembrane heterogeneity to Ising criticality, our spin-spin correlation function may be used as a predictive quantitative measure for FRET or NMR membrane experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Critical Casimir forces in cellular membranes

Recent experiments suggest that membranes of living cells are tuned close to a miscibility critical point in the 2D Ising universality class. We propose that one role for this proximity to criticality in live cells is to provide a conduit for relatively long-ranged critical Casimir forces. Using techniques from conformal field theory we calculate potentials of mean force between membrane bound inclusions mediated by their local interactions with the composition order parameter. We verify these calculations using Monte-Carlo where we also compare critical and off-critical results. Our findings suggest that membrane bound proteins experience weak yet long range forces mediated by critical composition fluctuations in the plasma membranes of living cells.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From damage percolation to crack nucleation through finite size criticality

We present a unified theory of fracture in disordered brittle media that reconciles apparently conflicting results reported in the literature. Our renormalization group based approach yields a phase diagram in which the percolation fixed point, expected for infinite disorder, is unstable for finite disorder and flows to a zero-disorder nucleation-type fixed point, thus showing that fracture has mixed first order and continuous character. In a region of intermediate disorder and finite system sizes, we predict a crossover with mean-field avalanche scaling. We discuss intriguing connections to other phenomena where critical scaling is only observed in finite size systems and disappears in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geodesic acceleration and the small-curvature approximation for nonlinear least squares

It has been shown numerically that the performance of the Levenberg-Marquardt algorithm can be improved by including a second order correction known as the geodesic acceleration. In this paper we give the method a more sound theoretical foundation by deriving the geodesic acceleration correction without using differential geometry and showing that the traditional convergence proofs can be adapted to incorporate geodesic acceleration. Unlike other methods which include second derivative information, the geodesic acceleration does not attempt to improve the Gauss-Newton approximate Hessian, but rather is an extension of the small-residual approximation to cubic order. In deriving geodesic acceleration, we note that the small-residual approximation is complemented by a small-curvature approximation. This latter approximation provides a much broader justification for the Gauss-Newton approximate Hessian and Levenberg-Marquardt algorithm. In particular, it is justifiable even if the best fit residuals are large, is dependent only on the model and not on the data being fit, and is applicable for the entire course of the algorithm and not just the region near the minimum.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Improvements to the Levenberg-Marquardt algorithm for nonlinear least-squares minimization

When minimizing a nonlinear least-squares function, the Levenberg-Marquardt algorithm can suffer from a slow convergence, particularly when it must navigate a narrow canyon en route to a best fit. On the other hand, when the least-squares function is very flat, the algorithm may easily become lost in parameter space. We introduce several improvements to the Levenberg-Marquardt algorithm in order to improve both its convergence speed and robustness to initial parameter guesses. We update the usual step to include a geodesic acceleration correction term, explore a systematic way of accepting uphill steps that may increase the residual sum of squares due to Umrigar and Nightingale, and employ the Broyden method to update the Jacobian matrix. We test these changes by comparing their performance on a number of test problems with standard implementations of the algorithm. We suggest that these two particular challenges, slow convergence and robustness to initial guesses, are complimentary problems. Schemes that improve convergence speed often make the algorithm less robust to the initial guess, and vice versa. We provide an open source implementation of our improvements that allow the user to adjust the algorithm parameters to suit particular needs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structural Susceptibility and Separation of Time Scales in the van der Pol Oscillator

We use an extension of the van der Pol oscillator as an example of a system with multiple time scales to study the susceptibility of its trajectory to polynomial perturbations in the dynamics. A striking feature of many nonlinear, multi-parameter models is an apparently inherent insensitivity to large magnitude variations in certain linear combinations of parameters. This phenomenon of "sloppiness" is quantified by calculating the eigenvalues of the Hessian matrix of the least-squares cost function which typically span many orders of magnitude. The van der Pol system is no exception: Perturbations in its dynamics show that most directions in parameter space weakly affect the limit cycle, whereas only a few directions are stiff. With this study we show that separating the time scales in the van der Pol system leads to a further separation of eigenvalues. Parameter combinations which perturb the slow manifold are stiffer and those which solely affect the transients in the dynamics are sloppier.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Avalanche Spatial Structure and Multivariable Scaling Functions; Sizes, Heights, Widths, and Views through Windows

We introduce a systematic method for extracting multivariable universal scaling functions and critical exponents from data. We exemplify our insights by analyzing simulations of avalanches in an interface using simulations from a driven quenched Kardar-Parisi-Zhang (qKPZ) equation. We fully characterize the spatial structure of these avalanches- we report universal scaling functions for size, height and width distributions, and also local front heights. Furthermore, we resolve a problem that arises in many imaging experiments of crackling noise and avalanche dynamics, where the observed distributions are strongly distorted by a limited field of view. Through artificially windowed data, we show these distributions and their multivariable scaling functions may be written in terms of two control parameters, the window size and the characteristic length scale of the dynamics. For the entire system and the windowed distributions we develop accurate parameterizations for the universal scaling functions, including corrections to scaling and systematic error bars, facilitated by a novel software environment SloppyScaling.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dielectric breakdown and avalanches at non-equilibrium metal-insulator transitions

Motivated by recent experiments on the finite temperature Mott transition in VO2 films, we propose a classical coarse-grained dielectric breakdown model where each degree of freedom represents a nanograin which transitions from insulator to metal with increasing temperature and voltage at random thresholds due to quenched disorder. We describe the properties of the resulting non-equilibrium metal-insulator transition and explain the universal characteristics of the resistance jump distribution. We predict that by tuning voltage, another critical point is approached, which separates a phase of "bolt"-like avalanches from percolation-like ones.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fracture strength of disordered media: Universality, interactions and tail asymptotics

We study the asymptotic properties of fracture strength distributions of disordered elastic media by a combination of renormalization group, extreme value theory, and numerical simulation. We investigate the validity of the `weakest-link hypothesis' in the presence of realistic long-ranged interactions in the random fuse model. Numerical simulations indicate that the fracture strength is well described by the Duxbury-Leath-Beale (DLB) distribution which is shown to flow asymptotically to the Gumbel distribution. We explore the relation between the extreme value distributions and the DLB type asymptotic distributions, and show that the universal extreme value forms may not be appropriate to describe the non-universal low-strength tail.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nucleation at the DNA supercoiling transition

Twisting DNA under a constant applied force reveals a thermally activated transition into a state with a supercoiled structure known as a plectoneme. Using transition state theory, we predict the rate of this plectoneme nucleation to be of order 10^4 Hz. We reconcile this with experiments that have measured hopping rates of order 10 Hz by noting that the viscosity of the bead used to manipulate the DNA limits the measured rate. We find that the intrinsic bending caused by disorder in the base-pair sequence is important for understanding the free energy barrier that governs the transition. Both analytic and numerical methods are used in the calculations. We provide extensive details on the numerical methods for simulating the elastic rod model with and without disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A minimal model of plasma membrane heterogeneity requires coupling cortical actin to criticality

We present a minimal model of plasma membrane heterogeneity that combines criticality with connectivity to cortical cytoskeleton. Our model is motivated by recent observations of micron-sized critical fluctuations in plasma membrane vesicles that are detached from their cortical cytoskeleton. We incorporate criticality using a conserved order parameter Ising model coupled to a simple actin cytoskeleton interacting through point-like pinning sites. Using this minimal model, we recapitulate several experimental observations of plasma membrane raft heterogeneity. Small (r~20nm) and dynamic fluctuations at physiological temperatures arise from criticality. Including connectivity to cortical cytoskeleton disrupts large fluctuations, prevents macroscopic phase separation at low temperatures (T<=22°C), and provides a template for long lived fluctuations at physiological temperature (T=37°C). Cytoskeleton-stabilized fluctuations produce significant barriers to the diffusion of some membrane components in a manner that is weakly dependent on the number of pinning sites and strongly dependent on criticality. More generally, we demonstrate that critical fluctuations provide a physical mechanism to organize and spatially segregate membrane components by providing channels for interaction over large distances.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Adaptive mutation of biochemical reaction constants: Fisher's geometrical model without pleiotropy

The distribution of fitness effects of adaptive mutations remains poorly understood, both empirically and theoretically. We study this distribution using a version of Fisher's geometrical model without pleiotropy, such that each mutation affects only a single trait. We are motivated by the notion of an organism's chemotype, the set of biochemical reaction constants that govern its molecular constituents. From physical considerations, we expect the chemotype to be of high dimension and to exhibit very little pleiotropy. Our model generically predicts striking cusps in the distribution of the fitness effects of arising and fixed mutations. It further predicts that a single element of the chemotype should comprise all mutations at the high-fitness ends of these distributions. Using extreme value theory, we show that the two cusps with the highest fitnesses are typically well-separated, even when the chemotype possesses thousands of elements; this suggests a means to observe these cusps experimentally. More broadly, our work demonstrates that new insights into evolution can arise from the chemotype perspective, a perspective between the genotype and the phenotype.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bending crystals: Emergence of fractal dislocation structures

We provide a minimal continuum model for mesoscale plasticity, explaining the cellular dislocation structures observed in deformed crystals. Our dislocation density tensor evolves from random, smooth initial conditions to form self-similar structures strikingly similar to those seen experimentally - reproducing both the fractal morphologies and some features of the scaling of cell sizes and misorientations analyzed experimentally. Our model provides a framework for understanding emergent dislocation structures on the mesoscale, a bridge across a computationally demanding mesoscale gap in the multiscale modeling program, and a new example of self-similar structure formation in non-equilibrium systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Beyond power laws: Universality in the average avalanche shape

We report the measurement of multivariable scaling functions for the temporal average shape of Barkhausen noise avalanches, and show that they are consistent with the predictions of simple mean-field theories. We bypass the confounding factors of time-retarded interactions (eddy currents) by measuring thin permal- loy films, and bypass thresholding effects and amplifier distortions by applying Wiener deconvolution. We find experimental shapes that are approximately symmetric, and track the evolution of the scaling function. We solve a mean- field theory for the magnetization dynamics and calculate the form of the scaling function in the presence of a demagnetizing field and a finite field ramp-rate, yielding quantitative agreement with the experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Superheating field of superconductors within Ginzburg-Landau theory

We study the superheating field of a bulk superconductor within Ginzburg-Landau theory, which is valid near the critical temperature. We calculate, as functions of the Ginzburg-Landau parameter $κ$, the superheating field $\Hsh$ and the critical momentum $k_c$ characterizing the wavelength of the instability of the Meissner state to flux penetration. By mapping the two-dimensional linear stability theory into a one-dimensional eigenfunction problem for an ordinary differential equation, we solve the problem numerically. We demonstrate agreement between the numerics and analytics, and show convergence to the known results at both small and large $κ$. We discuss the implications of the results for superconducting RF cavities used in particle accelerators.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The geometry of nonlinear least squares with applications to sloppy models and optimization

Parameter estimation by nonlinear least squares minimization is a common problem with an elegant geometric interpretation: the possible parameter values of a model induce a manifold in the space of data predictions. The minimization problem is then to find the point on the manifold closest to the data. We show that the model manifolds of a large class of models, known as sloppy models, have many universal features; they are characterized by a geometric series of widths, extrinsic curvatures, and parameter-effects curvatures. A number of common difficulties in optimizing least squares problems are due to this common structure. First, algorithms tend to run into the boundaries of the model manifold, causing parameters to diverge or become unphysical. We introduce the model graph as an extension of the model manifold to remedy this problem. We argue that appropriate priors can remove the boundaries and improve convergence rates. We show that typical fits will have many evaporated parameters. Second, bare model parameters are usually ill-suited to describing model behavior; cost contours in parameter space tend to form hierarchies of plateaus and canyons. Geometrically, we understand this inconvenient parametrization as an extremely skewed coordinate basis and show that it induces a large parameter-effects curvature on the manifold. Using coordinates based on geodesic motion, these narrow canyons are transformed in many cases into a single quadratic, isotropic basin. We interpret the modified Gauss-Newton and Levenberg-Marquardt fitting algorithms as an Euler approximation to geodesic motion in these natural coordinates on the model manifold and the model graph respectively. By adding a geodesic acceleration adjustment to these algorithms, we alleviate the difficulties from parameter-effects curvature, improving both efficiency and success rates at finding good fits.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Discontinuities at the DNA supercoiling transition

While slowly turning the ends of a single molecule of DNA at constant applied force, a discontinuity was recently observed at the supercoiling transition, when a small plectoneme is suddenly formed. This can be understood as an abrupt transition into a state in which stretched and plectonemic DNA coexist. We argue that there should be discontinuities in both the extension and the torque at the transition, and provide experimental evidence for both. To predict the sizes of these discontinuities and how they change with the overall length of DNA, we organize a theory for the coexisting plectonemic state in terms of four length-independent parameters. We also test plectoneme theories, including our own elastic rod simulation, finding discrepancies with experiment that can be understood in terms of the four coexisting state parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Sloppiness, robustness, and evolvability in systems biology

The functioning of many biochemical networks is often robust -- remarkably stable under changes in external conditions and internal reaction parameters. Much recent work on robustness and evolvability has focused on the structure of neutral spaces, in which system behavior remains invariant to mutations. Recently we have shown that the collective behavior of multiparameter models is most often 'sloppy': insensitive to changes except along a few 'stiff' combinations of parameters, with an enormous sloppy neutral subspace. Robustness is often assumed to be an emergent evolved property, but the sloppiness natural to biochemical networks offers an alternative non-adaptive explanation. Conversely, ideas developed to study evolvability in robust systems can be usefully extended to characterize sloppy systems.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Variational method for estimating the rate of convergence of Markov Chain Monte Carlo algorithms

We demonstrate the use of a variational method to determine a quantitative lower bound on the rate of convergence of Markov Chain Monte Carlo (MCMC) algorithms as a function of the target density and proposal density. The bound relies on approximating the second largest eigenvalue in the spectrum of the MCMC operator using a variational principle and the approach is applicable to problems with continuous state spaces. We apply the method to one dimensional examples with Gaussian and quartic target densities, and we contrast the performance of the Random Walk Metropolis-Hastings (RWMH) algorithm with a ``smart'' variant that incorporates gradient information into the trial moves. We find that the variational method agrees quite closely with numerical simulations. We also see that the smart MCMC algorithm often fails to converge geometrically in the tails of the target density except in the simplest case we examine, and even then care must be taken to choose the appropriate scaling of the deterministic and random parts of the proposed moves. Again, this calls into question the utility of smart MCMC in more complex problems. Finally, we apply the same method to approximate the rate of convergence in multidimensional Gaussian problems with and without importance sampling. There we demonstrate the necessity of importance sampling for target densities which depend on variables with a wide range of scales.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Universally Sloppy Parameter Sensitivities in Systems Biology

Quantitative computational models play an increasingly important role in modern biology. Such models typically involve many free parameters, and assigning their values is often a substantial obstacle to model development. Directly measuring \emph{in vivo} biochemical parameters is difficult, and collectively fitting them to other data often yields large parameter uncertainties. Nevertheless, in earlier work we showed in a growth-factor-signaling model that collective fitting could yield well-constrained predictions, even when it left individual parameters very poorly constrained. We also showed that the model had a `sloppy' spectrum of parameter sensitivities, with eigenvalues roughly evenly distributed over many decades. Here we use a collection of models from the literature to test whether such sloppy spectra are common in systems biology. Strikingly, we find that every model we examine has a sloppy spectrum of sensitivities. We also test several consequences of this sloppiness for building predictive models. In particular, sloppiness suggests that collective fits to even large amounts of ideal time-series data will often leave many parameters poorly constrained. Tests over our model collection are consistent with this suggestion. This difficulty with collective fits may seem to argue for direct parameter measurements, but sloppiness also implies that such measurements must be formidably precise and complete to usefully constrain many model predictions. We confirm this implication in our signaling model. Our results suggest that sloppy sensitivity spectra are universal in systems biology models. The prevalence of sloppiness highlights the power of collective fits and suggests that modelers should focus on predictions rather than on parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Phonons in Molecular Quantum Dots: Density Functional Calculation of Franck-Condon Emission Rates in External Fields

We report the calculation of various phonon overlaps and their corresponding phonon emission probabilities for the problem of an electron tunneling onto and off of the buckyball-dimer molecular quantum dot $C_{72}$, both with and without the influence of an external field. We show that the stretch mode of the two balls of the dumbbell couples most strongly to the electronic transition, and in turn that a field in the direction of the bond between the two $C_{36}$ balls is most effective at further increasing the phonon emission into the stretch mode. As the field is increased, phonon emission increases in probability with an accompanying decrease in probability of the dot remaining in the ground vibrational state. We also present a simple model to gauge the effect of molecular size on the phonon emission of molecules similar to our $C_{72}$ molecule, including the experimentally tested $C_{140}$. In our approach we do not assume that the hessians of the molecule are identical for different charge states. Our treatment is hence a generalization of the traditional phonon overlap calculations for coupled electron-photon transition in solids.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Scaling in Plasticity-Induced Cell-Boundary Microstructure: Fragmentation and Rotational Diffusion

We develop a simple computational model for cell boundary evolution in plastic deformation. We study the cell boundary size distribution and cell boundary misorientation distribution that experimentally have been found to have scaling forms that are largely material independent. The cell division acts as a source term in the misorientation distribution which significantly alters the scaling form, giving it a linear slope at small misorientation angles as observed in the experiments. We compare the results of our simulation to two closely related exactly solvable models which exhibit scaling behavior at late times: (i) fragmentation theory and (ii) a random walk in rotation space with a source term. We find that the scaling exponents in our simulation agree with those of the theories, and that the scaling collapses obey the same equations, but that the shape of the scaling functions depend upon the methods used to measure sizes and to weight averages and histograms.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Defect--Defect Correlation Functions, Generic Scale Invariance, and the Complex Ginzburg--Landau Equation

Motivated by the idea of developing a ``hydrodynamic'' description of spatiotemporal chaos, we have investigated the defect--defect correlation functions in the defect turbulence regime of the two--dimensional, anisotropic complex Ginzburg--Landau equation. We compare our results with the predictions of generic scale invariance. Using the topological nature of the defects, we prove that defect--defect correlations cannot decay as slowly as predicted by generic scale invariance. We also present results on the fluctuations of the amplitude field $A$.

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

Logarithmically slow domain growth in nonrandomly frustrated systems--- A summary of the major results

Here, we summarize the most important results of our study of logarithmically slow growth of domains following a quench in two models without randomness in their Hamiltonians. This is a slightly updated version of a paper to appear in the Proceedings of the 1st Annual Tohwa University International Symposium, Fukuoka, Japan (American Institute of Physics, 1992). It is meant to serve as a brief summary of cond-mat/9204015 for those who do not wish to read all the details contained therein (and don't want to hassle with 2 MBytes of tex/ps files).

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

Melting a copper cluster: Critical droplet theory

We simulate the melting of a 71 A diameter cluster of Cu. At low temperatures the crystal exhibits facets. With increasing temperatures the open facets pre-melt, the melted regions coalesce into a liquid envelope containing a crystalline nucleus, and the nucleus finally goes unstable to the supercooled liquid. Using critical droplet theory and experimental data for Cu, we explain the thermodynamics of the coexistence region. The width of the transition scales as (Number of particles) to the power (-1/4).

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

Tweed in Martensites: A Potential New Spin Glass

We've been studying the ``tweed'' precursors above the martensitic transition in shape--memory alloys. These characteristic cross--hatched modulations occur for hundreds of degrees above the first--order shape--changing transition. Our two--dimensional model for this transition, in the limit of infinite elastic anisotropy, can be mapped onto a spin--glass Hamiltonian in a random field. We suggest that the tweed precursors are a direct analogy of the spin--glass phase. The tweed is intermediate between the high--temperature cubic phase and the low--temperature martensitic phase in the same way as the spin--glass phase can be intermediate between ferromagnet and antiferromagnet.

0 citations0 reviewsNo code linkedOpen access