Source author record

Nigel Goldenfeld

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
24topics
4close collaborators

Actions

Connect this record

Log in to claim

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

42 published item(s)

preprint2025arXiv

Beyond chaos: fluctuations, anomalies and spontaneous stochasticity in fluid turbulence

In this perspective, we consider the development of statistical hydrodynamics, focusing on the way in which the intrinsic stochasticity of turbulent phenomena was identified and is being explored. A major purpose of our discussion is to bring out the role of anomalies in turbulent phenomena, in ways that are not usually done, and to emphasize how the description of turbulent phenomena requires delicate considerations of asymptotic limits. The scope of our narrative includes selected historical aspects that are not usually emphasized, primarily due to G.I. Taylor, as well as discussions of certain aspects of the laminar-turbulent transition, the behaviour of turbulent drag at intermediate Reynolds numbers, and the statistics of fully-developed turbulence that exhibit spontaneous stochasticity.

preprint2021arXiv

Stochastic social behavior coupled to COVID-19 dynamics leads to waves, plateaus and an endemic state

It is well recognized that population heterogeneity plays an important role in the spread of epidemics. While individual variations in social activity are often assumed to be persistent, i.e. constant in time, here we discuss the consequences of dynamic heterogeneity. By integrating the stochastic dynamics of social activity into traditional epidemiological models we demonstrate the emergence of a new long timescale governing the epidemic in broad agreement with empirical data. Our model captures multiple features of real-life epidemics such as COVID-19, including prolonged plateaus and multiple waves, which are transiently suppressed due to the dynamic nature of social activity. The existence of the long timescale due to the interplay between epidemic and social dynamics provides a unifying picture of how a fast-paced epidemic typically will transition to the endemic state.

preprint2021arXiv

Time-dependent heterogeneity leads to transient suppression of the COVID-19 epidemic, not herd immunity

Epidemics generally spread through a succession of waves that reflect factors on multiple timescales. On short timescales, super-spreading events lead to burstiness and overdispersion, while long-term persistent heterogeneity in susceptibility is expected to lead to a reduction in the infection peak and the herd immunity threshold (HIT). Here, we develop a general approach to encompass both timescales, including time variations in individual social activity, and demonstrate how to incorporate them phenomenologically into a wide class of epidemiological models through parameterization. We derive a non-linear dependence of the effective reproduction number Re on the susceptible population fraction S. We show that a state of transient collective immunity (TCI) emerges well below the HIT during early, high-paced stages of the epidemic. However, this is a fragile state that wanes over time due to changing levels of social activity, and so the infection peak is not an indication of herd immunity: subsequent waves can and will emerge due to behavioral changes in the population, driven (e.g.) by seasonal factors. Transient and long-term levels of heterogeneity are estimated by using empirical data from the COVID-19 epidemic as well as from real-life face-to-face contact networks. These results suggest that the hardest-hit areas, such as NYC, have achieved TCI following the first wave of the epidemic, but likely remain below the long-term HIT. Thus, in contrast to some previous claims, these regions can still experience subsequent waves.

preprint2020arXiv

A Field-Theoretic Model for Chemotaxis in Run and Tumble Particles

In this paper we develop a field-theoretic description for run and tumble chemotaxis, based on a density functional description of crystalline materials modified to capture orientational ordering. We show that this framework, with its in-built multi-particle interactions, soft-core repulsion and elasticity is ideal for describing continuum collective phases with particle resolution, but on diffusive timescales. We show that our model exhibits particle aggregation in an externally imposed constant attractant field, as is observed for phototactic or thermotactic agents. We also show that this model captures particle aggregation through self-chemotaxis, an important mechanism that aids quorum dependent cellular interactions.

preprint2020arXiv

Window of Opportunity for Mitigation to Prevent Overflow of ICU capacity in Chicago by COVID-19

We estimate the growth in demand for ICU beds in Chicago during the emerging COVID-19 epidemic, using state-of-the-art computer simulations calibrated for the SARS-CoV-2 virus. The questions we address are these: (1) Will the ICU capacity in Chicago be exceeded, and if so by how much? (2) Can strong mitigation strategies, such as lockdown or shelter in place order, prevent the overflow of capacity? (3) When should such strategies be implemented? Our answers are as follows: (1) The ICU capacity may be exceeded by a large amount, probably by a factor of ten. (2) Strong mitigation can avert this emergency situation potentially, but even that will not work if implemented too late. (3) If the strong mitigation precedes April 1st, then the growth of COVID-19 can be controlled and the ICU capacity could be adequate. The earlier the strong mitigation is implemented, the greater the probability that it will be successful. After around April 1 2020, any strong mitigation will not avert the emergency situation. In Italy, the lockdown occurred too late and the number of deaths is still doubling every 2.3 days. It is difficult to be sure about the precise dates for this window of opportunity, due to the inherent uncertainties in computer simulation. But there is high confidence in the main conclusion that it exists and will soon be closed. Our conclusion is that, being fully cognizant of the societal trade-offs, there is a rapidly closing window of opportunity to avert a worst-case scenario in Chicago, but only with strong mitigation/lockdown implemented in the next week at the latest. If this window is missed, the epidemic will get worse and then strong mitigation/lockdown will be required after all, but it will be too late.

preprint2016arXiv

Stochastic predator-prey dynamics of transposons in the human genome

Transposable elements, or transposons, are DNA sequences that can jump from site to site in the genome during the life cycle of a cell, usually encoding the very enzymes which perform their excision. However, some transposons are parasitic, relying on the enzymes produced by the regular transposons. In this case, we show that a stochastic model, which takes into account the small copy numbers of the transposons in a cell, predicts noise-induced predator-prey oscillations with a characteristic time scale that is much longer than the cell replication time, indicating that the state of the predator-prey oscillator is stored in the genome and transmitted to successive generations. Our work demonstrates the important role of number fluctuations in the expression of mobile genetic elements, and shows explicitly how ecological concepts can be applied to the dynamics and fluctuations of living genomes.

preprint2016arXiv

Turbulence as a problem in non-equilibrium statistical mechanics

The transitional and well-developed regimes of turbulent shear flows exhibit a variety of remarkable scaling laws that are only now beginning to be systematically studied and understood. In the first part of this article, we summarize recent progress in understanding the friction factor of turbulent flows in rough pipes and quasi-two-dimensional soap films, showing how the data obey a two-parameter scaling law known as roughness-induced criticality, and exhibit power-law scaling of friction factor with Reynolds number that depends on the precise form of the nature of the turbulent cascade. These results hint at a non-equilibrium fluctuation-dissipation relation that applies to turbulent flows. The second part of this article concerns the lifetime statistics in smooth pipes around the transition, showing how the remarkable super-exponential scaling with Reynolds number reflects deep connections between large deviation theory, extreme value statistics, directed percolation and the onset of coexistence in predator-prey ecosystems. Both these phenomena reflect the way in which turbulence can be fruitfully approached as a problem in non-equilibrium statistical mechanics.

preprint2015arXiv

A framework for analyzing ecological trait-based models in multi-dimensional niche spaces

We develop an theoretical approach for predicting biodiversity in multi-dimensional niche spaces, arising due to ecological drivers such as competitive exclusion. The novelty of our approach relies on the fact that ecological niches are described by sequences of strings, which allows us to describe multiple traits. We define the mathematical framework for analyzing pattern forming instabilities in these models, showing surprisingly that the analytic linear theory predicts the asymptotically long time population distributions of niches in the model. We propose a test for identifying ecological drivers in biodiversity distributions, based on representing ecosystem data by means of a certain transform introduced in the theory.

preprint2015arXiv

A noise-induced mechanism for biological homochirality of early life self-replicators

The observed single-handedness of biological amino acids and sugars has long been attributed to autocatalysis. However, the stability of homochiral states in deterministic autocatalytic systems relies on cross inhibition of the two chiral states, an unlikely scenario for early life self-replicators. Here, we present a theory for a stochastic individual-level model of autocatalysis due to early life self-replicators. Without chiral inhibition, the racemic state is the global attractor of the deterministic dynamics, but intrinsic multiplicative noise stabilizes the homochiral states, in both well-mixed and spatially-extended systems. We conclude that autocatalysis is a viable mechanism for homochirality, without imposing additional nonlinearities such as chiral inhibition.

preprint2015arXiv

Crystal plasticity as a mean field depinning transition: results from a phase field crystal model

Until now, most of our knowledge about the universality class of crystal plasticity has come from simulations using discrete dislocation dynamics. These are force-controlled, typically at zero temperature, and deal with the creation and annihilation of dislocations phenomenologically. In this work, we go beyond these limitations by using phase field crystal simulations in two dimensions at finite temperature to extract the avalanche statistics of a simulated crystal under constant shear velocity. In addition to the avalanche size and energy distributions we extract the avalanche duration distributions and power spectra. All exponents and scaling functions extracted here for the statics and dynamics of crystal plasticity, belong to the mean field elastic depinning universality class, confirming earlier findings based on discrete dislocation dynamics.

preprint2015arXiv

Ecological collapse and the emergence of traveling waves at the onset of shear turbulence

The transition to turbulence exhibits remarkable spatio-temporal behavior that continues to defy detailed understanding. Near the onset to turbulence in pipes, transient turbulent regions decay either directly or, at higher Reynolds numbers through splitting, with characteristic time-scales that exhibit a super-exponential dependence on Reynolds number. Here we report numerical simulations of transitional pipe flow, showing that a zonal flow emerges at large scales, activated by anisotropic turbulent fluctuations; in turn, the zonal flow suppresses the small-scale turbulence leading to stochastic predator-prey dynamics. We show that this "ecological" model of transitional turbulence reproduces the super-exponential lifetime statistics and phenomenology of pipe flow experiments. Our work demonstrates that a fluid on the edge of turbulence is mathematically analogous to an ecosystem on the edge of extinction, and provides an unbroken link between the equations of fluid dynamics and the directed percolation universality class.

preprint2015arXiv

Enhanced Microfluidic Mixing via a Tricritical Spiral Vortex Instability

Experimental measurements and numerical simulations are made on fluid flow through cross-slot devices with a range of aspect (depth:width) ratios, 0.4 < alpha < 3.87. For low Reynolds numbers Re, the flow is symmetric and a sharp boundary exists between fluid streams entering the cross-slot from opposite directions. Above an alpha-dependent critical value Re_c, the flow undergoes a symmetry-breaking bifurcation (though remains steady and laminar) and a spiral vortex structure develops about the central axis of the outflow channel. An order parameter characterizing the instability grows according to a sixth-order Landau potential, and shows a progression from second order to first order transitions as alpha increases. A tricritical point occurs for alpha ~ 0.55. The spiral vortex acts as a mixing region in the flow field and this phenomenon can be used to drive enhanced mixing in microfluidic devices.

preprint2014arXiv

Intermittent dislocation density fluctuations in crystal plasticity from a phase-field crystal model

Plastic deformation mediated by collective dislocation dynamics is investigated in the two-dimensional phase-field crystal model of sheared single crystals. We find that intermittent fluctuations in the dislocation population number accompany bursts in the plastic strain-rate fluctuations. Dislocation number fluctuations exhibit a power-law spectral density $1/f^2$ at high frequencies $f$. The probability distribution of number fluctuations becomes bimodal at low driving rates corresponding to a scenario where low density of defects alternate at irregular times with high population of defects. We propose a simple stochastic model of dislocation reaction kinetics that is able to capture these statistical properties of the dislocation density fluctuations as a function of shear rate.

preprint2013arXiv

Determination of the universality class of crystal plasticity

Although scaling phenomena have long been documented in crystalline plasticity, the universality class has been difficult to identify due to the rarity of avalanche events, which require large system sizes and long times in order to accurately measure scaling exponents and functions. Here we present comprehensive simulations of two-dimensional dislocation dynamics under shear, using finite-size scaling to extract scaling exponents and the avalanche profile scaling function from time-resolved measurements of slip-avalanches. Our results provide compelling evidence that both the static and dynamic universality classes are consistent with the mean-field interface depinning model.

preprint2013arXiv

Distribution of maximum velocities in avalanches near the depinning transition

We report exact predictions for universal scaling exponents and scaling functions associated with the distribution of the maximum collective avalanche propagation velocities $v_m$ in the mean field theory of the interface depinning transition. We derive the extreme value distribution $P(v_m|T)$ for the maximum velocities in avalanches of fixed duration $T$, and verify the results by numerical simulation near the critical point. We find that the tail of the distribution of maximum velocity for an arbitrary avalanche duration, $v_m$, scales as $P(v_m)\sim v_m^{-2}$ for large $v_m$. These results account for the observed power-law distribution of the maximum amplitudes in acoustic emission experiments of crystal plasticity, and are also broadly applicable to other systems in the mean-field interface depinning universality class, ranging from magnets to earthquakes.

preprint2013arXiv

Extrinsic noise driven phenotype switching in a self-regulating gene

Due to inherent noise in intracellular networks cellular decisions can be random, so genetically identical cells can display different phenotypic behavior even in identical environments. Most previous work in understanding the decision-making process has focused on the role of intrinsic noise in these systems. Yet, especially in the high copy-number regime, extrinsic noise has been shown to be much more significant. Here, using a prototypical example of a bistable self-regulating gene model, we develop a theoretical framework describing the combined effect of intrinsic and extrinsic noise on the dynamics of stochastic genetic switches. Employing our theory and Monte Carlo simulations, we show that extrinsic noise not only significantly alters the lifetimes of the phenotypic states, but can induce bistability in unexpected regions of parameter space, and may fundamentally change the escape mechanism. These results have implications for interpreting experimentally observed heterogeneity in cellular populations and for stochastic modeling of cellular decision processes.

preprint2013arXiv

Universal fluctuations and extreme statistics of avalanches near the depinning transtition

We derive exact predictions for universal scaling exponents and scaling functions associated with the statistics of maximum velocities vm during avalanches described by the mean field theory of the interface depinning transition. In particular, we find a robust power-law regime in the statistics of maximum events that can explain the observed distribution of the peak amplitudes in acoustic emission experiments of crystal plasticity. Our results are expected to be broadly applicable to a broad range of systems in the mean-field interface depinning universality class, ranging from magnets to earthquakes.

preprint2012arXiv

Anisotropic velocity statistics of topological defects under shear flow

We report numerical results on the velocity statistics of topological defects during the dynamics of phase ordering and non-relaxational evolution assisted by an external shear ow. We propose a numerically efficient tracking method for finding the position and velocity of defects, and apply it to vortices in a uniform field and dislocations in anisotropic stripe patterns. During relaxational dynamics, the distribution function of the velocity fuctuations is characterized by a dynamical scaling with a scaling function that has a robust algebraic tail with an inverse cube power law. This is characteristic to defects of codimension two, e.g. point defects in two dimensions and filaments in three dimensions, regardless of whether the motion is isotropic (as for vortices) or highly anisotropic (as for dislocations). However, the anisotropic dislocation motion leads to anisotropic statistical properties when the interaction between defects and their motion is in infuenced by the presence of an external shear flow transverse to the stripe orientation.

preprint2012arXiv

Avalanche statistics and intermittency in topological defect-mediated flows

Topological defects dominate the deformation response of materials in processes ranging from quantum turbulence to crystal plasticity. We calculate the probability distribution function for the fluctuations in velocity $v$, using scaling arguments and a systematic cluster expansion method to account for density correlations. We find that the distribution has power-law tails with an exponent that takes the value -3 for $v\rightarrow \infty$, but a value -2 for intermediate values of $v$. We relate these regimes to the theory of avalanches, by directly computing the known avalanche scaling exponents.

preprint2012arXiv

Effect of weak rotation on large-scale circulation cessations in turbulent convection

We investigate the effect of weak rotation on the large-scale circulation (LSC) of turbulent Rayleigh-Bénard convection, using the theory for cessations in a low-dimensional stochastic model of the flow previously studied. We determine the cessation frequency of the LSC as a function of rotation, and calculate the statistics of the amplitude and azimuthal velocity fluctuations of the LSC as a function of the rotation rate for different Rayleigh numbers. Furthermore, we show that the tails of the reorientation probability distribution function remain unchanged for rotating systems, while the distribution of the LSC amplitude and correspondingly the cessation frequency are strongly affected by rotation. Our results are in close agreement with experimental observations.

preprint2011arXiv

Directed percolation describes lifetime and growth of turbulent puffs and slugs

We show that Directed Percolation (DP) simulations in a pipe geometry in 3+1 dimensions fully capture the observed complex phenomenology of the transition to turbulence. At low Reynolds numbers (Re), turbulent puffs form and spontaneously relaminarize. At high Re, turbulent slugs expand uniformly into the laminar regions. In a spatiotemporally intermittent state between these two regimes of Re, puffs split and turbulent regions exhibit laminar patches. DP also captures the main quantitative features of the transition, with a superexponentially diverging characteristic lifetime below the transition. Above the percolation threshold, active (turbulent) clusters expand into the inactive (laminar) phase with a well-defined velocity whose scaling with control parameter (Reynolds number or percolation probability) is consistent with experimental results. Our results provide strong evidence in favor of a conjecture of Pomeau.

preprint2011arXiv

Dislocations Jam At Any Density

Crystalline materials deform in an intermittent way via dislocation-slip avalanches. Below a critical stress, the dislocations are jammed within their glide plane due to long-range elastic interactions and the material exhibits plastic response, while above this critical stress the dislocations are mobile (the unjammed phase) and the material fails. We use dislocation dynamics and scaling arguments in two dimensions to show that the critical stress grows with the square root of the dislocation density. Consequently, dislocations jam at any density, in contrast to granular materials, which only jam below a critical density.

preprint2011arXiv

Fluctuation-driven Turing patterns

Models of diffusion driven pattern formation that rely on the Turing mechanism are utilized in many areas of science. However, many such models suffer from the defect of requiring fine tuning of parameters or an unrealistic separation of scales in the diffusivities of the constituents of the system in order to predict the formation of spatial patterns. In the context of a very generic model of ecological pattern formation, we show that the inclusion of intrinsic noise in Turing models leads to the formation of "quasi-patterns" that form in generic regions of parameter space and are experimentally distinguishable from standard Turing patterns. The existence of quasi-patterns removes the need for unphysical fine tuning or separation of scales in the application of Turing models to real systems.

preprint2011arXiv

Rare fluctuations and large-scale circulation cessations in turbulent convection

In turbulent Rayleigh-Benard convection, a large-scale circulation (LSC) develops in a nearly vertical plane, and is maintained by rising and falling plumes detaching from the unstable thermal boundary layers. Rare but large fluctuations in the LSC amplitude can lead to extinction of the LSC (a cessation event), followed by the re-emergence of another LSC with a different (random) azimuthal orientation. We extend previous models of the LSC dynamics to include momentum and thermal diffusion in the azimuthal plane, and calculate the tails of the probability distributions of both the amplitude and azimuthal angle. Our analytical results are in very good agreement with experimental data.

preprint2011arXiv

The dynamics of gene duplication and transposons in microbial genomes following a sudden environmental change

A variety of genome transformations can occur as a microbial population adapts to a large environmental change. In particular, genomic surveys indicate that, following the transition to an obligate, host-dependent symbiont, the density of transposons first rises, then subsequently declines over evolutionary time. Here, we show that these observations can be accounted for by a class of generic stochastic models for the evolution of genomes in the presence of continuous selection and gene duplication. The models use a fitness function that allows for partial contributions from multiple gene copies, is an increasing but bounded function of copy number, and is optimal for one fully adapted gene copy. We use Monte Carlo simulation to show that the dynamics result in an initial rise in gene copy number followed by a subsequent fall due to adaptation to the new environmental parameters. These results are robust for reasonable gene duplication and mutation parameters when adapting to a novel target sequence. Our model provides a generic explanation for the dynamics of microbial transposon density following a large environmental changes such as host restriction.

preprint2011arXiv

The effect of Coulombic friction on spatial displacement statistics

The phenomenon of Coulombic friction enters the stochastic description of dry friction between two solids and the statistic characterization of vibrating granular media. Here we analyze the corresponding Fokker-Planck equation including both velocity and spatial components, exhibiting a formal connection to a quantum mechanical harmonic oscillator in the presence of a delta potential. Numerical solutions for the resulting spatial displacement statistics show a crossover from exponential to Gaussian displacement statistics. We identify a transient intermediate regime that exhibits multiscaling properties arising from the contribution of Coulombic friction. The possible role of these effects during observations in diffusion experiments is shortly discussed.

preprint2011arXiv

Theory of cooperation in a micro-organismal snow-drift game

We present a mean field model for the phase diagram of a community of micro-organisms, interacting through their metabolism so that they are, in effect, engaging in a cooperative social game. We show that as a function of the concentration of the nutrients glucose and histidine, the community undergoes a phase transition separating a state in which one strain is dominant to a state which is characterized by coexisting populations. Our results are in good agreement with recent experimental results, correctly predicting quantitative trends and the phase diagram.

preprint2010arXiv

Emergence of foams from the breakdown of the phase field crystal model

The phase field crystal (PFC) model captures the elastic and topological properties of crystals with a single scalar field at small undercooling. At large undercooling, new foam-like behavior emerges. We characterize this foam phase of the PFC equation and propose a modified PFC equation that may be used for the simulation of foam dynamics. This minimal model reproduces von Neumann's rule for two-dimensional dry foams, and Lifshitz-Slyozov coarsening for wet foams. We also measure the coordination number distribution and find that its second moment is larger than previously-reported experimental and theoretical studies of soap froths, a finding that we attribute to the wetness of the foam increasing with time.

preprint2010arXiv

Fixed points and limit cycles in the population dynamics of lysogenic viruses and their hosts

Starting with stochastic rate equations for the fundamental interactions between microbes and their viruses, we derive a mean field theory for the population dynamics of microbe-virus systems, including the effects of lysogeny. In the absence of lysogeny, our model is a generalization of that proposed phenomenologically by Weitz and Dushoff. In the presence of lysogeny, we analyze the possible states of the system, identifying a novel limit cycle, which we interpret physically. To test the robustness of our mean field calculations to demographic fluctuations, we have compared our results with stochastic simulations using the Gillespie algorithm. Finally, we estimate the range of parameters that delineate the various steady states of our model.

preprint2010arXiv

Plasticity and Dislocation Dynamics in a Phase Field Crystal Model

The critical dynamics of dislocation avalanches in plastic flow is examined using a phase field crystal (PFC) model. In the model, dislocations are naturally created, without any \textit{ad hoc} creation rules, by applying a shearing force to the perfectly periodic ground state. These dislocations diffuse, interact and annihilate with one another, forming avalanche events. By data collapsing the event energy probability density function for different shearing rates, a connection to interface depinning dynamics is confirmed. The relevant critical exponents agree with mean field theory predictions.

preprint2010arXiv

Statistical Mechanics of Horizontal Gene Transfer in Evolutionary Ecology

The biological world, especially its majority microbial component, is strongly interacting and may be dominated by collective effects. In this review, we provide a brief introduction for statistical physicists of the way in which living cells communicate genetically through transferred genes, as well as the ways in which they can reorganize their genomes in response to environmental pressure. We discuss how genome evolution can be thought of as related to the physical phenomenon of annealing, and describe the sense in which genomes can be said to exhibit an analogue of information entropy. As a direct application of these ideas, we analyze the variation with ocean depth of transposons in marine microbial genomes, predicting trends that are consistent with recent observations using metagenomic surveys.

preprint2009arXiv

Emergence of heterogeneity and political organization in information exchange networks

We present a simple model of the emergence of the division of labor and the development of a system of resource subsidy from an agent-based model of directed resource production with variable degrees of trust between the agents. The model has three distinct phases, corresponding to different forms of societal organization: disconnected (independent agents), homogeneous cooperative (collective state), and inhomogeneous cooperative (collective state with a leader). Our results indicate that such levels of organization arise generically as a collective effect from interacting agent dynamics, and may have applications in a variety of systems including social insects and microbial communities.

preprint2009arXiv

Macroscopic effects of the spectral structure in turbulent flows

Two aspects of turbulent flows have been the subject of extensive, split research efforts: macroscopic properties, such as the frictional drag experienced by a flow past a wall, and the turbulent spectrum. The turbulent spectrum may be said to represent the fabric of a turbulent state; in practice it is a power law of exponent α(the "spectral exponent") that gives the revolving velocity of a turbulent fluctuation (or "eddy") of size s as a function of s. The link, if any, between macroscopic properties and the turbulent spectrum remains missing. Might it be found by contrasting the frictional drag in flows with differing types of spectra? Here we perform unprecedented measurements of the frictional drag in soap-film flows, where the spectral exponent α= 3 and compare the results with the frictional drag in pipe flows, where the spectral exponent α= 5/3. For moderate values of the Reynolds number Re (a measure of the strength of the turbulence), we find that in soap-film flows the frictional drag scales as Re^{-1/2}, whereas in pipe flows the frictional drag scales as Re^{-1/4} . Each of these scalings may be predicted from the attendant value of αby using a new theory, in which the frictional drag is explicitly linked to the turbulent spectrum. Our work indicates that in turbulence, as in continuous phase transitions, macroscopic properties are governed by the spectral structure of the fluctuations.

preprint2009arXiv

Optimality Properties of a Proposed Precursor to the Genetic Code

We calculate the optimality of a doublet precursor to the canonical genetic code with respect to mitigating the effects of point mutations and compare our results to corresponding ones for the canonical genetic code. We find that the proposed precursor has much less optimality than that of the canonical code. Our results render unlikely the notion that the doublet precursor was an intermediate state in the evolution of the canonical genetic code. These findings support the notion that code optimality reflects evolutionary dynamics, and that if such a doublet code originally had a biochemical significance, it arose before the emergence of translation.

preprint2009arXiv

Robust ecological pattern formation induced by demographic noise

We demonstrate that demographic noise can induce persistent spatial pattern formation and temporal oscillations in the Levin-Segel predator-prey model for plankton-herbivore population dynamics. Although the model exhibits a Turing instability in mean field theory, demographic noise greatly enlarges the region of parameter space where pattern formation occurs. To distinguish between patterns generated by fluctuations and those present at the mean field level in real ecosystems, we calculate the power spectrum in the noise-driven case and predict the presence of fat tails not present in the mean field case. These results may account for the prevalence of large-scale ecological patterns, beyond that expected from traditional non-stochastic approaches.

preprint2009arXiv

The friction factor of two-dimensional rough-boundary turbulent soap film flows

We use momentum transfer arguments to predict the friction factor $f$ in two-dimensional turbulent soap-film flows with rough boundaries (an analogue of three-dimensional pipe flow) as a function of Reynolds number Re and roughness $r$, considering separately the inverse energy cascade and the forward enstrophy cascade. At intermediate Re, we predict a Blasius-like friction factor scaling of $f\propto\textrm{Re}^{-1/2}$ in flows dominated by the enstrophy cascade, distinct from the energy cascade scaling of $\textrm{Re}^{-1/4}$. For large Re, $f \sim r$ in the enstrophy-dominated case. We use conformal map techniques to perform direct numerical simulations that are in satisfactory agreement with theory, and exhibit data collapse scaling of roughness-induced criticality, previously shown to arise in the 3D pipe data of Nikuradse.

preprint2008arXiv

Lambda-prophage induction modeled as a cooperative failure mode of lytic repression

We analyze a system-level model for lytic repression of lambda-phage in E. coli using reliability theory, showing that the repressor circuit comprises 4 redundant components whose failure mode is prophage induction. Our model reflects the specific biochemical mechanisms involved in regulation, including long-range cooperative binding, and its detailed predictions for prophage induction in E. coli under ultra-violet radiation are in good agreement with experimental data.

preprint2005arXiv

Renormalization group approach to multiscale modelling in materials science

Dendritic growth, and the formation of material microstructure in general, necessarily involves a wide range of length scales from the atomic up to sample dimensions. The phase field approach of Langer, enhanced by optimal asymptotic methods and adaptive mesh refinement, copes with this range of scales, and provides an effective way to move phase boundaries. However, it fails to preserve memory of the underlying crystallographic anisotropy, and thus is ill-suited for problems involving defects or elasticity. The phase field crystal (PFC) equation-- a conserving analogue of the Hohenberg-Swift equation --is a phase field equation with periodic solutions that represent the atomic density. It can natively model elasticity, the formation of solid phases, and accurately reproduces the nonequilibrium dynamics of phase transitions in real materials. However, the PFC models matter at the atomic scale, rendering it unsuitable for coping with the range of length scales in problems of serious interest. Here, we show that a computationally-efficient multiscale approach to the PFC can be developed systematically by using the renormalization group or equivalent techniques to derive appropriate coarse-grained coupled phase and amplitude equations, which are suitable for solution by adaptive mesh refinement algorithms.

preprint1998arXiv

Dynamical Scaling Behavior of the Swift-Hohenberg Equation Following a Quench to the Modulated State

We study the kinetics of phase transitions in a Rayleigh-Benard system after onset of convection using 2D Swift-Hohenberg equation. An initially uniform state evolves to one whose ground state is spatially periodic. We confirmed previous results which showed that dynamical scaling occurs at medium quench ($ε= 0.25$) with scaling exponents 1/5 and 1/4 under zero noise and finite noise respectively. We find logarithmic scaling behavior for a deep quench ($ε= 0.75$) at zero noise. A simple method is devised to measure the proxy of domain wall length. We find that the energy and domain wall length exhibit scaling behavior with the same exponent. For $ε= 0.25$, the scaling exponents are 1/4 and 0.3 at zero and finite noise respectively.

preprint1996arXiv

Experimental Constraints on the Pairing State of the Cuprate Superconductors: an Emerging Consensus

We present a critical discussion of recent experimental probes of the pairing state of the high temperature superconductors, focusing primarily, but not exclusively, on \Yba, where the best data currently exist. Penetration depth measurements near \Tc\ give no indication of an extra transition, indicating that the pairing state is a one-dimensional representation of the crystal symmetry. Penetration depth measurements at low temperatures provide strong evidence for a change in sign of the gap function over the Fermi surface. Quantum mechanical phase interference experiments generally confirm this and in addition show that the nodal positions are consistent with a \dx2y2\ pairing state. This pairing state is consistent with photoemission measurements of the gap function, Raman scattering, the effect on \Tc\ of impurities, and many other data (reviewed by two of us previously) which indicate the presence of low lying excitations in the superconducting state. We also discuss evidence that apparently does not fit in with a \dx2y2\ pairing state, and we describe possible alternative scenarios.