Search PubMed⌕ Search

Biomedical subjects

S Havlin

Publications and source records attributed to S Havlin.

At least 37 records · Page 2Linked to original sources

Dependence of conductance on percolation backbone mass

We study , the average conductance of the backbone, defined by two points separated by Euclidean distance r, of mass M(B) on two-dimensional percolation clusters at the percolation threshold. We find that with increasing M(B) and for fixed r, asymptotically decreases to a constant, in contrast with the behavior of homogeneous systems and nonrandom fractals (such as the Sierpinski gasket) in which conductance increases with increasing M(B). We explain this behavior by studying the distribution of shortest paths between the two points on clusters with a given M(B). We also study the dependence of conductance on M(B) above the percolation threshold and find that (i) slightly above p(c), the conductance first decreases and then increases with increasing M(B) and (ii) further above p(c), the conductance increases monotonically for all values of M(B), as is the case for homogeneous systems.

Journal Article↗

Multifractal behavior of linear polymers in disordered media.

The scaling behavior of linear polymers in disordered media modeled by self-avoiding random walks (SAWs) on the backbone of two- and three-dimensional percolation clusters at their critical concentrations p(c) is studied. All possible SAW configurations of N steps on a single backbone configuration are enumerated exactly. We find that the moments of order q of the total number of SAWs obtained by averaging over many backbone configurations display multifractal behavior; i.e., different moments are dominated by different subsets of the backbone. This leads to generalized coordination numbers mu(q) and enhancement exponents gamma(q), which depend on q. Our numerical results suggest that the relation mu(1)=p(c)mu between the first moment mu(1) and its regular lattice counterpart mu is valid.

Journal Article↗

Swelling-collapse transition of self-attracting walks

We study the structural properties of self-attracting walks in d dimensions using scaling arguments and Monte Carlo simulations. We find evidence of a transition analogous to the Theta transition of polymers. Above a critical attractive interaction u(c), the walk collapses and the exponents nu and k, characterizing the scaling with time t of the mean square end-to-end distance approximately t(2nu) and the average number of visited sites approximately t(k), are universal and given by nu=1/(d+1) and k=d/(d+1). Below u(c), the walk swells and the exponents are as with no interaction, i.e., nu=1/2 for all d, k=1/2 for d=1 and k=1 for d>/=2. At u(c), the exponents are found to be in a different universality class.

Journal Article↗

Multifractal properties of the random resistor network

We study the multifractal spectrum of the current in the two-dimensional random resistor network at the percolation threshold. We consider two ways of applying the voltage difference: (i) two parallel bars, and (ii) two points. Our numerical results suggest that in the infinite system limit, the probability distribution behaves for small i as P(i) approximately 1/i, where i is the current. As a consequence, the moments of i of order q</=q(c)=0 do not exist and all currents of value below the most probable one have the fractal dimension of the backbone. The backbone can thus be described in terms of only (i) blobs of fractal dimension d(B) and (ii) high current carrying bonds of fractal dimension going from 1/nu to d(B).

Journal Article↗

Flow between two sites on a percolation cluster

We study the flow of fluid in porous media in dimensions d=2 and 3. The medium is modeled by bond percolation on a lattice of L(d) sites, while the flow front is modeled by tracer particles driven by a pressure difference between two fixed sites ("wells") separated by Euclidean distance r. We investigate the distribution function of the shortest path connecting the two sites, and propose a scaling ansatz that accounts for the dependence of this distribution (i) on the size of the system L and (ii) on the bond occupancy probability p. We confirm by extensive simulations that the ansatz holds for d=2 and 3. Further, we study two dynamical quantities: (i) the minimal traveling time of a tracer particle between the wells when the total flux is constant and (ii) the minimal traveling time when the pressure difference is constant. A scaling ansatz for these dynamical quantities also includes the effect of finite system size L and off-critical bond occupation probability p. We find that the scaling form for the distribution functions for these dynamical quantities for d=2 and 3 is similar to that for the shortest path, but with different critical exponents. Our results include estimates for all parameters that characterize the scaling form for the shortest path and the minimal traveling time in two and three dimensions; these parameters are the fractal dimension, the power law exponent, and the constants and exponents that characterize the exponential cutoff functions.

Journal Article↗

Decomposition of heartbeat time series: scaling analysis of the sign sequence.

The cardiac interbeat (RR) increment time series can be decomposed into two sub-sequences: a magnitude series and a sign series. The authors show that the sign sequence, a simple binary representation of the original RR series, retains fundamental scaling properties of the original series, is robust with respect to outliers, and may provide useful information about neuroautonomic control mechanisms.

Adrenergic beta-Antagonists↗

Application of statistical physics to heartbeat diagnosis.

We present several recent studies based on statistical physics concepts that can be used as diagnostic tools for heart failure. We describe the scaling exponent characterizing the long-range correlations in heartbeat time series as well as the multifractal features recently discovered in heartbeat rhythm. It is found that both features, the long-range correlations and the multifractility, are weaker in cases of heart failure.

Adult↗

Sleep-wake differences in scaling behavior of the human heartbeat: analysis of terrestrial and long-term space flight data.

We compare scaling properties of the cardiac dynamics during sleep and wake periods for healthy individuals, cosmonauts during orbital flight, and subjects with severe heart disease. For all three groups, we find a greater degree of anticorrelation in the heartbeat fluctuations during sleep compared to wake periods. The sleep-wake difference in the scaling exponents for the three groups is comparable to the difference between healthy and diseased individuals. The observed scaling differences are not accounted for simply by different levels of activity, but appear related to intrinsic changes in the neuroautonomic control of the heartbeat.

Astronauts↗

Scaling in nature: from DNA through heartbeats to weather.

The purpose of this report is to describe some recent progress in applying scaling concepts to various systems in nature. We review several systems characterized by scaling laws such as DNA sequences, heartbeat rates and weather variations. We discuss the finding that the exponent alpha quantifying the scaling in DNA in smaller for coding than for noncoding sequences. We also discuss the application of fractal scaling analysis to the dynamics of heartbeat regulation, and report the recent finding that the scaling exponent alpha is smaller during sleep periods compared to wake periods. We also discuss the recent findings that suggest a universal scaling exponent characterizing the weather fluctuations.

Algorithms↗

Optimizing the success of random searches.

We address the general question of what is the best statistical strategy to adapt in order to search efficiently for randomly located objects ('target sites'). It is often assumed in foraging theory that the flight lengths of a forager have a characteristic scale: from this assumption gaussian, Rayleigh and other classical distributions with well-defined variances have arisen. However, such theories cannot explain the long-tailed power-law distributions of flight lengths or flight times that are observed experimentally. Here we study how the search efficiency depends on the probability distribution of flight lengths taken by a forager that can detect target sites only in its limited vicinity. We show that, when the target sites are sparse and can be visited any number of times, an inverse square power-law distribution of flight lengths, corresponding to Lévy flight motion, is an optimal strategy. We test the theory by analysing experimental foraging data on selected insect, mammal and bird species, and find that they are consistent with the predicted inverse square power-law distributions.

Animals↗

Statistical physics and physiology: monofractal and multifractal approaches.

Even under healthy, basal conditions, physiologic systems show erratic fluctuations resembling those found in dynamical systems driven away from a single equilibrium state. Do such "nonequilibrium" fluctuations simply reflect the fact that physiologic systems are being constantly perturbed by external and intrinsic noise? Or, do these fluctuations actually, contain useful, "hidden" information about the underlying nonequilibrium control mechanisms? We report some recent attempts to understand the dynamics of complex physiologic fluctuations by adapting and extending concepts and methods developed very recently in statistical physics. Specifically, we focus on interbeat interval variability as an important quantity to help elucidate possibly non-homeostatic physiologic variability because (i) the heart rate is under direct neuroautonomic control, (ii) interbeat interval variability is readily measured by noninvasive means, and (iii) analysis of these heart rate dynamics may provide important practical diagnostic and prognostic information not obtainable with current approaches. The analytic tools we discuss may be used on a wider range of physiologic signals. We first review recent progress using two analysis methods--detrended fluctuation analysis and wavelets--sufficient for quantifying monofractual structures. We then describe recent work that quantifies multifractal features of interbeat interval series, and the discovery that the multifractal structure of healthy subjects is different than that of diseased subjects.

Biophysical Phenomena↗

Multifractality in human heartbeat dynamics.

There is evidence that physiological signals under healthy conditions may have a fractal temporal structure. Here we investigate the possibility that time series generated by certain physiological control systems may be members of a special class of complex processes, termed multifractal, which require a large number of exponents to characterize their scaling properties. We report on evidence for multifractality in a biological dynamical system, the healthy human heartbeat, and show that the multifractal character and nonlinear properties of the healthy heart rate are encoded in the Fourier phases. We uncover a loss of multifractality for a life-threatening condition, congestive heart failure.

Adult↗

Dynamics of plaque formation in Alzheimer's disease.

Plaques that form in the brains of Alzheimer patients are made of deposits of the amyloid-beta peptide. We analyze the time evolution of amyloid-beta deposition in immunostained brain slices from transgenic mice. We find that amyloid-beta deposits appear in clusters whose characteristic size increases from 14 microm in 8-month-old mice to 22 microm in 12-month-old mice. We show that the clustering has implications for the biological growth of amyloid-beta by presenting a growth model that accounts for the experimentally observed structure of individual deposits and predicts the formation of clusters of deposits and their time evolution.

Alzheimer Disease↗

Scaling features of noncoding DNA.

We review evidence supporting the idea that the DNA sequence in genes containing noncoding regions is correlated, and that the correlation is remarkably long range--indeed, base pairs thousands of base pairs distant are correlated. We do not find such a long-range correlation in the coding regions of the gene, and utilize this fact to build a Coding Sequence Finder Algorithm, which uses statistical ideas to locate the coding regions of an unknown DNA sequence. Finally, we describe briefly some recent work adapting to DNA the Zipf approach to analyzing linguistic texts, and the Shannon approach to quantifying the "redundancy" of a linguistic text in terms of a measurable entropy function, and reporting that noncoding regions in eukaryotes display a larger redundancy than coding regions. Specifically, we consider the possibility that this result is solely a consequence of nucleotide concentration differences as first noted by Bonhoeffer and his collaborators. We find that cytosine-guanine (CG) concentration does have a strong "background" effect on redundancy. However, we find that for the purine-pyrimidine binary mapping rule, which is not affected by the difference in CG concentration, the Shannon redundancy for the set of analyzed sequences is larger for noncoding regions compared to coding regions.

Algorithms↗

Clustering of identical oligomers in coding and noncoding DNA sequences.

We develop a quantitative method for analyzing repetitions of identical short oligomers in coding and noncoding DNA sequences. We analyze sequences presently available in the GenBank separately for primate, mammal, vertebrate, rodent, invertebrate and plant taxonomic partitions. We find that some oligomers "cluster" more than they would if randomly distributed, while other oligomers "repel" each other. To quantify this degree of clustering, we define clustering measures. We find that (i) clustering significantly differs in coding and noncoding DNA; (ii) in most cases, monomers, dimers and tetramers cluster in noncoding DNA but appear to repel each other in coding DNA. (iii) The degree of clustering for different sources (primates, invertebrates, and plants) is more conserved among these sources in the case of coding DNA than in the case of noncoding DNA. (iv) In contrast to other oligomers, we find that trimers always prefer to cluster. (v) Clustering of each particular oligomer is conserved within the same organism.

Animals↗

Dynamic feedback in an aggregation-disaggregation model.

We study an aggregation-disaggregation model which is relevant to biological processes such as the growth of senile plaques in Alzheimer disease. In this model, during the aggregation each deposited particle has a probability of producing a new particle in its vicinity, while during disaggregation the particles are anihilated randomly. The model is held in a dynamic equilibrium by a feedback mechanism which changes the disaggregation probability in proportion to the change in the total number of particles. We also include surface diffusion which influences the morphology of growing aggregates and colonies. A colony includes the descendents of a single particle. We investigate the statistical properties of the model in two dimensions. We find that unlike the colonies, individual aggregates are fractals with a fractal dimension of D(f)=1.92+/-0.06 in the absence of surface diffusion. We show that the surface diffusion changes the fractal dimension of aggregates: at a small aggregation-disaggregation rate, D(f) is independent of the strength of the surface diffusion, D(f)=1.73+/-0.03. At larger aggregation-disaggregation rates and different strengths of surface diffusion, aggregates with fractal dimensions between D(f)=1.73 and 1.92 form. The steady-state distribution of aggregate sizes is shown to be power law if the aggregation-disaggregation process dominates over the surface diffusion. In the limit of weak aggregation-disaggregation and strong surface diffusion the size distribution is log-normal.

Alzheimer Disease↗

Quasicrystals in a monodisperse system.

We investigate the formation of a two-dimensional quasicrystal in a monodisperse system, using molecular dynamics simulations of hard-sphere particles interacting via a two-dimensional square-well potential. We find that more than one stable crystalline phase can form for certain values of the square-well parameters. Quenching the liquid phase at a very low temperature, we obtain an amorphous phase. By heating this amorphous phase, we obtain a quasicrystalline structure with fivefold symmetry. From estimations of the Helmholtz potentials of the stable crystalline phases and of the quasicrystal, we conclude that the observed quasicrystal phase can be the stable phase in a specific range of temperatures.

Journal Article↗

Traveling time and traveling length in critical percolation clusters.

We study traveling time and traveling length for tracer dispersion in two-dimensional bond percolation, modeling flow by tracer particles driven by a pressure difference between two points separated by Euclidean distance r. We find that the minimal traveling time t(min) scales as t(min) approximately r(1.33), which is different from the scaling of the most probable traveling time, t* approximately r(1.64). We also calculate the length of the path corresponding to the minimal traveling time and find l(min) approximately r(1.13) and that the most probable traveling length scales as l* approximately r(1.21). We present the relevant distribution functions and scaling relations.

Journal Article↗