PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “statistical physics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6Linked to original sources

Mean field theory of spherical gravitating systems

Important gaps remain in our understanding of the thermodynamics and statistical physics of self-gravitating systems. Using mean field theory, here we investigate the equilibrium properties of several spherically symmetric model systems confined in a finite domain consisting of either point masses or rotating mass shells of different dimension. We establish a direct connection between the spherically symmetric equilibrium states of a self-gravitating point mass system and a shell model of dimension 3. We construct the equilibrium density functions by maximizing the entropy subject to the usual constraints of normalization and energy, but we also take into account the constraint on the sum of the squares of the individual angular momenta, which is also an integral of motion for these symmetric systems. Two statistical ensembles are introduced that incorporate the additional constraint. They are used to investigate the possible occurrence of a phase transition as the defining parameters for each ensemble are altered.

Journal Article↗

Statistically preserved structures in shell models of passive scalar advection.

It was conjectured recently that statiscally preserved structures underlie the statistical physics of turbulent transport processes. We analyze here in detail the time-dependent (noncompact) linear operator that governs the dynamics of correlation functions in the case of shell models of passive scalar advection. The problem is generic in the sense that the driving velocity field is neither Gaussian nor delta correlated in time. We show how to naturally discuss the dynamics in terms of an effective compact operator that displays "zero modes," which determine the anomalous scaling of the correlation functions. Since shell models have neither a Lagrangian structure nor "shape dynamics," this example differs significantly from standard passive scalar advection. Nevertheless, with the necessary modifications, the generality and efficacy of the concept of statistically preserved structures are further exemplified. In passing we point out a bonus of the present approach, in providing analytic predictions for the time-dependent correlation functions in decaying turbulent transport.

Journal Article↗

Bimodality for plant sizes and spatial pattern in cohorts: the role of competition and site conditions.

The well-known size bimodality in cohorts of plants is revisited with methods emphasizing generic modeling. A link between bimodality and interface growth in nonequilibrium statistical physics is emphasized: the development of bimodality is understood as a phase transition like interface roughening. A specific model is proposed, inspired by gap models. It is used to illustrate generic results, to understand the end point of mean field approximation and aggregation of variables, and to discuss the role of site and competition for the development of a social hierarchy within a forest stand, in a wider ecological context.

Cohort Studies↗

Statistical mechanics of combinatorial auctions.

Combinatorial auctions are formulated as frustrated lattice gases on sparse random graphs, allowing the determination of the optimal revenue by methods of statistical physics. Transitions between computationally easy and hard regimes are found and interpreted in terms of the geometric structure of the space of solutions. We introduce an iterative algorithm to solve intermediate and large instances, and discuss competing states of optimal revenue and maximal number of satisfied bidders. The algorithm can be generalized to the hard phase and to more sophisticated auction protocols.

Algorithms↗

Chaos and threshold for irreversibility in sheared suspensions.

Systems governed by time reversible equations of motion often give rise to irreversible behaviour. The transition from reversible to irreversible behaviour is fundamental to statistical physics, but has not been observed experimentally in many-body systems. The flow of a newtonian fluid at low Reynolds number can be reversible: for example, if the fluid between concentric cylinders is sheared by boundary motion that is subsequently reversed, then all fluid elements return to their starting positions. Similarly, slowly sheared suspensions of solid particles, which occur widely in nature and science, are governed by time reversible equations of motion. Here we report an experiment showing precisely how time reversibility fails for slowly sheared suspensions. We find that there is a concentration dependent threshold for the deformation or strain beyond which particles do not return to their starting configurations after one or more cycles. Instead, their displacements follow the statistics of an anisotropic random walk. By comparing the experimental results with numerical simulations, we demonstrate that the threshold strain is associated with a pronounced growth in the Lyapunov exponent (a measure of the strength of chaotic particle interactions). The comparison illuminates the connections between chaos, reversibility and predictability.

Journal Article↗

Principles of protein-protein recognition from structure to thermodynamics.

Specific recognition is illustrated by X-ray structures of protease-inhibitor, antigen-antibody and other high affinity complexes including five electron transfer complexes. We attempt to give a physical definition to affinity and specificity on the basis of these data. In a protein-protein complex, specific recognition results from the assembly of complementary surfaces into well-packed interfaces that cover about 1500 A2 and contain about ten hydrogen bonds. These interfaces are larger than between molecules in protein crystals, and smaller than between subunits in oligomeric proteins. We relate the size and chemical nature of interfaces in complexes to the thermodynamical parameters that characterize affinity: the heat capacity and free enthalpy (Gibbs energy) of dissociation at equilibrium, the activation free enthalpy for the dissociation reaction. The same structural and thermodynamical parameters are inadequate for representing the specificity of recognition. We propose instead to describe specificity with the help of statistical physics, and we illustrate the application of the random energy model to antigen-antibody recognition by analyzing results of computer simulations by docking.

Chemical Phenomena↗

Protein threading by learning.

By using techniques borrowed from statistical physics and neural networks, we determine the parameters, associated with a scoring function, that are chosen optimally to ensure complete success in threading tests in a training set of proteins. These parameters provide a quantitative measure of the propensities of amino acids to be buried or exposed and to be in a given secondary structure and are a good starting point for solving both the threading and design problems.

Amino Acid Sequence↗

Transcription/replication collisions cause bacterial transcription units to be longer on the leading strand of replication.

MOTIVATION: The costs of head-on versus codirectional collisions between the replication complex and the much slower transcription complex on the circular bacterial chromosomes are much debated. Although it is established that the number of genes on the leading strand is higher than on the lagging strand of replication, the consequences of collisions on the length of transcription units are unknown. RESULTS: Here, we show that transcription units are generally longer on the leading strand, in rough proportion to the bias in number of units between the two strands. We propose a statistical physics model, based on the assumption that the proportion of interrupted transcripts for each unit is the major factor contributing to these biases. Its main prediction is that a large replication/transcription speed ratio implies a low leading/lagging bias for transcription unit length and number. This model is validated by an analysis of proven and predicted units in Escherichia coli and Bacillus subtilis. The results are consistent with an equal cost of head-on versus codirectional collisions.

Bacillus subtilis↗

A normal tissue dose response model of dynamic repair processes.

A model is presented for serial, critical element complication mechanisms for irradiated volumes from length scales of a few millimetres up to the entire organ. The central element of the model is the description of radiation complication as the failure of a dynamic repair process. The nature of the repair process is seen as reestablishing the structural organization of the tissue, rather than mere replenishment of lost cells. The interactions between the cells, such as migration, involved in the repair process are assumed to have finite ranges, which limits the repair capacity and is the defining property of a finite-sized reconstruction unit. Since the details of the repair processes are largely unknown, the development aims to make the most general assumptions about them. The model employs analogies and methods from thermodynamics and statistical physics. An explicit analytical form of the dose response of the reconstruction unit for total, partial and inhomogeneous irradiation is derived. The use of the model is demonstrated with data from animal spinal cord experiments and clinical data about heart, lung and rectum. The three-parameter model lends a new perspective to the equivalent uniform dose formalism and the established serial and parallel complication models. Its implications for dose optimization are discussed.

Animals↗

Stochastic reasoning, free energy, and information geometry.

Belief propagation (BP) is a universal method of stochastic reasoning. It gives exact inference for stochastic models with tree interactions and works surprisingly well even if the models have loopy interactions. Its performance has been analyzed separately in many fields, such as AI, statistical physics, information theory, and information geometry. This article gives a unified framework for understanding BP and related methods and summarizes the results obtained in many fields. In particular, BP and its variants, including tree reparameterization and concave-convex procedure, are reformulated with information-geometrical terms, and their relations to the free energy function are elucidated from an information-geometrical viewpoint. We then propose a family of new algorithms. The stabilities of the algorithms are analyzed, and methods to accelerate them are investigated.

Algorithms↗

Nonequilibrium statistical mechanical models for cytoskeletal assembly: towards understanding tensegrity in cells.

The cytoskeleton is not an equilibrium structure. To develop theoretical tools to investigate such nonequilibrium assemblies, we study a statistical physical model of motorized spherical particles. Though simple, it captures some of the key nonequilibrium features of the cytoskeletal networks. Variational solutions of the many-body master equation for a set of motorized particles accounts for their thermally induced Brownian motion as well as for the motorized kicking of the structural elements. These approximations yield stability limits for crystalline phases and for frozen amorphous structures. The methods allow one to compute the effects of nonequilibrium behavior and adhesion (effective cross-linking) on the mechanical stability of localized phases as a function of density, adhesion strength, and temperature. We find that nonequilibrium noise does not necessarily destabilize mechanically organized structures. The nonequilibrium forces strongly modulate the phase behavior and have comparable effect as the adhesion due to cross-linking. Modeling transitions such as these allows the mechanical properties of cytoskeleton to rapidly and adaptively change. The present model provides a statistical mechanical underpinning for a tensegrity picture of the cytoskeleton.

Animals↗

Estimation of distribution algorithms with Kikuchi approximations.

The question of finding feasible ways for estimating probability distributions is one of the main challenges for Estimation of Distribution Algorithms (EDAs). To estimate the distribution of the selected solutions, EDAs use factorizations constructed according to graphical models. The class of factorizations that can be obtained from these probability models is highly constrained. Expanding the class of factorizations that could be employed for probability approximation is a necessary step for the conception of more robust EDAs. In this paper we introduce a method for learning a more general class of probability factorizations. The method combines a reformulation of a probability approximation procedure known in statistical physics as the Kikuchi approximation of energy, with a novel approach for finding graph decompositions. We present the Markov Network Estimation of Distribution Algorithm (MN-EDA), an EDA that uses Kikuchi approximations to estimate the distribution, and Gibbs Sampling (GS) to generate new points. A systematic empirical evaluation of MN-EDA is done in comparison with different Bayesian network based EDAs. From our experiments we conclude that the algorithm can outperform other EDAs that use traditional methods of probability approximation in the optimization of functions with strong interactions among their variables.

Algorithms↗

On the analysis of large-scale genomic structures.

We apply methods from statistical physics (histograms, correlation functions, fractal dimensions, and singularity spectra) to characterize large-scale structure of the distribution of nucleotides along genomic sequences. We discuss the role of the extension of noncoding segments ("junk DNA") for the genomic organization, and the connection between the coding segment distribution and the high-eukaryotic chromatin condensation. The following sequences taken from GenBank were analyzed: complete genome of Xanthomonas campestri, complete genome of yeast, chromosome V of Caenorhabditis elegans, and human chromosome XVII around gene BRCA1. The results are compared with the random and periodic sequences and those generated by simple and generalized fractal Cantor sets.

Algorithms↗

Bubbles and denaturation in DNA.

The local opening of DNA is an intriguing phenomenon from a statistical-physics point of view, but is also essential for its biological function. For instance, the transcription and replication of our genetic code cannot take place without the unwinding of the DNA double helix. Although these biological processes are driven by proteins, there might well be a relation between these biological openings and the spontaneous bubble formation due to thermal fluctuations. Mesoscopic models, like the Peyrard-Bishop-Dauxois (PBD) model, have fairly accurately reproduced some experimental denaturation curves and the sharp phase transition in the thermodynamic limit. It is, hence, tempting to see whether these models could be used to predict the biological activity of DNA. In a previous study, we introduced a method that allows to obtain very accurate results on this subject, which showed that some previous claims in this direction, based on molecular-dynamics studies, were premature. This could either imply that the present PBD model should be improved or that biological activity can only be predicted in a more complex framework that involves interactions with proteins and super helical stresses. In this article, we give a detailed description of the statistical method introduced before. Moreover, for several DNA sequences, we give a thorough analysis of the bubble-statistics as a function of position and bubble size and the so-called l-denaturation curves that can be measured experimentally. These show that some important experimental observations are missing in the present model. We discuss how the present model could be improved.

Computer Simulation↗

Topological approach to neural complexity.

Considerable effort in modern statistical physics is devoted to the study of networked systems. One of the most important example of them is the brain, which creates and continuously develops complex networks of correlated dynamics. An important quantity which captures fundamental aspects of brain network organization is the neural complexity C(X) introduced by Tononi et al. [Proc. Natl. Acad. Sci. USA 91, 5033 (1994)]. This work addresses the dependence of this measure on the topological features of a network in the case of a Gaussian stationary process. Both analytical and numerical results show that the degree of complexity has a clear and simple meaning from a topological point of view. Moreover, the analytical result offers a straightforward and faster algorithm to compute the complexity of a graph than the standard one.

Animals↗

The interplay between viscoelastic and thermodynamic properties determines the birefringence of F-actin gels.

F-actin gels of increasing concentrations (25-300 microM) display in vitro a progressive onset of birefringence due to orientational ordering of actin filaments. At F-actin concentrations <100 microM, this birefringence can be erased and restored at will by sonication and gentle flow, respectively. Hence, the orientational ordering does not result from a thermodynamic transition to a nematic phase but instead is due to mechanical stresses stored in the gels. In contrast, at F-actin concentrations > or =100 microM, gels display spontaneous birefringence recovery, at rest, which is the sign of true nematic ordering, in good agreement with statistical physics models of the isotropic/nematic transition. Well-aligned samples of F-actin gels could be produced and their small-angle x-ray scattering patterns are quite anisotropic. These patterns show no sign of filament positional short-range order and could be modeled by averaging the form factor with the Maier-Saupe nematic distribution function. The derived nematic order parameter S of the gels ranged from S = 0.7 at 300 microM to S = 0.4 at 25 microM. Both birefringence and small-angle x-ray scattering data indicate that, even in absence of cross-linking proteins, spontaneous cooperative alignment of actin filaments may arise in motile regions of living cells where F-actin concentrations can reach values of a few 100 microM.

Actins↗

Long-range correlations in rabbit brain neural activity.

We have analyzed the presence of persistence properties in rabbit brain electrical signals by means of non-equilibrium statistical physics tools. To measure long-memory properties of these experimental signals, we have first determined whether the data are fractional Gaussian noise (fGn) or fractional Brownian motion (fBm) by calculating the slope of the power spectral density plot of the series. The results show that the series correspond to fBm. Then, the data were studied by means of the bridge detrended scaled windowed variance analysis, detecting long-term correlation. Three different types of experimental signals have been studied: neural basal activity without stimulation, the response induced by a single flash light stimulus and the average of the activity evoked by 200 flash light stimulations. Analysis of the series revealed the existence of persistent behavior in all cases. Moreover, the results also exhibited an increasing correlation in the level of long-term memory from recordings without stimulation, to one sweep recording or 200 sweeps averaged recordings. Thus, brain neural electrical activity is affected not only by its most recent states, but also by previous states much more distant in the past.

Action Potentials↗

A theory of cortical responses.

This article concerns the nature of evoked brain responses and the principles underlying their generation. We start with the premise that the sensory brain has evolved to represent or infer the causes of changes in its sensory inputs. The problem of inference is well formulated in statistical terms. The statistical fundaments of inference may therefore afford important constraints on neuronal implementation. By formulating the original ideas of Helmholtz on perception, in terms of modern-day statistical theories, one arrives at a model of perceptual inference and learning that can explain a remarkable range of neurobiological facts.It turns out that the problems of inferring the causes of sensory input (perceptual inference) and learning the relationship between input and cause (perceptual learning) can be resolved using exactly the same principle. Specifically, both inference and learning rest on minimizing the brain's free energy, as defined in statistical physics. Furthermore, inference and learning can proceed in a biologically plausible fashion. Cortical responses can be seen as the brain's attempt to minimize the free energy induced by a stimulus and thereby encode the most likely cause of that stimulus. Similarly, learning emerges from changes in synaptic efficacy that minimize the free energy, averaged over all stimuli encountered. The underlying scheme rests on empirical Bayes and hierarchical models of how sensory input is caused. The use of hierarchical models enables the brain to construct prior expectations in a dynamic and context-sensitive fashion. This scheme provides a principled way to understand many aspects of cortical organization and responses. The aim of this article is to encompass many apparently unrelated anatomical, physiological and psychophysical attributes of the brain within a single theoretical perspective. In terms of cortical architectures, the theoretical treatment predicts that sensory cortex should be arranged hierarchically, that connections should be reciprocal and that forward and backward connections should show a functional asymmetry (forward connections are driving, whereas backward connections are both driving and modulatory). In terms of synaptic physiology, it predicts associative plasticity and, for dynamic models, spike-timing-dependent plasticity. In terms of electrophysiology, it accounts for classical and extra classical receptive field effects and long-latency or endogenous components of evoked cortical responses. It predicts the attenuation of responses encoding prediction error with perceptual learning and explains many phenomena such as repetition suppression, mismatch negativity (MMN) and the P300 in electroencephalography. In psychophysical terms, it accounts for the behavioural correlates of these physiological phenomena, for example, priming and global precedence. The final focus of this article is on perceptual learning as measured with the MMN and the implications for empirical studies of coupling among cortical areas using evoked sensory responses.

Biophysical Phenomena↗