Multiresolution analysis in statistical mechanics. II. The wavelet transform as a basis for Monte Carlo simulations on lattices

In this paper, we extend our analysis of lattice systems using the wavelet transform to systems for which exact enumeration is impractical. For such systems, we illustrate a wavelet-accelerated Monte Carlo (WAMC) algorithm, which hierarchically coarse-grains a lattice model by computing the probability distribution for successively larger block spins. We demonstrate that although the method perturbs the system by changing its Hamiltonian and by allowing block spins to take on values not permitted for individual spins, the results obtained agree with the analytical results in the preceding paper, and ``converge'' to exact results obtained in the absence of coarse-graining. Additionally, we show that the decorrelation time for the WAMC is no worse than that of Metropolis Monte Carlo (MMC), and that scaling laws can be constructed from data performed in several short simulations to estimate the results that would be obtained from the original simulation. Although the algorithm is not asymptotically faster than traditional MMC, because of its hierarchical design, the new algorithm executes several orders of magnitude faster than a full simulation of the original problem. Consequently, the new method allows for rapid analysis of a phase diagram, allowing computational time to be focused on regions near phase transitions.Comment: 11 pages plus 7 figures in PNG format (downloadable separately

Multiresolution analysis in statistical mechanics. I. Using wavelets to calculate thermodynamic properties

The wavelet transform, a family of orthonormal bases, is introduced as a technique for performing multiresolution analysis in statistical mechanics. The wavelet transform is a hierarchical technique designed to separate data sets into sets representing local averages and local differences. Although one-to-one transformations of data sets are possible, the advantage of the wavelet transform is as an approximation scheme for the efficient calculation of thermodynamic and ensemble properties. Even under the most drastic of approximations, the resulting errors in the values obtained for average absolute magnetization, free energy, and heat capacity are on the order of 10%, with a corresponding computational efficiency gain of two orders of magnitude for a system such as a $4\times 4$ Ising lattice. In addition, the errors in the results tend toward zero in the neighborhood of fixed points, as determined by renormalization group theory.Comment: 13 pages plus 7 figures (PNG

Quantitative linkage of physiology and gene expression through empirical model construction: an investigation of diabetes

A methodology for the construction of predictive empirical models of physiological characteristics from microarray data is presented. The method, applied here to the study of the development of diabetes and insulin resistance, can be further expanded to other cases and to also include a variety of other data, such as protein expression, or metabolic flux data. The importance of several of the genes identified by the modeling methodology can be verified by comparison with results from prior literature. This implies potentially significant roles in diabetes for several of the uncharacterized genes discovered during the modeling procedure.Singapore-MIT Alliance (SMA

Nonlinear software sensor for monitoring genetic regulation processes with noise and modeling errors

Nonlinear control techniques by means of a software sensor that are commonly used in chemical engineering could be also applied to genetic regulation processes. We provide here a realistic formulation of this procedure by introducing an additive white Gaussian noise, which is usually found in experimental data. Besides, we include model errors, meaning that we assume we do not know the nonlinear regulation function of the process. In order to illustrate this procedure, we employ the Goodwin dynamics of the concentrations [B.C. Goodwin, Temporal Oscillations in Cells, (Academic Press, New York, 1963)] in the simple form recently applied to single gene systems and some operon cases [H. De Jong, J. Comp. Biol. 9, 67 (2002)], which involves the dynamics of the mRNA, given protein, and metabolite concentrations. Further, we present results for a three gene case in co-regulated sets of transcription units as they occur in prokaryotes. However, instead of considering their full dynamics, we use only the data of the metabolites and a designed software sensor. We also show, more generally, that it is possible to rebuild the complete set of nonmeasured concentrations despite the uncertainties in the regulation function or, even more, in the case of not knowing the mRNA dynamics. In addition, the rebuilding of concentrations is not affected by the perturbation due to the additive white Gaussian noise and also we managed to filter the noisy output of the biological systemComment: 21 pages, 7 figures; also selected in vjbio of August 2005; this version corrects a misorder in the last three references of the published versio

Phase Separation of Crystal Surfaces: A Lattice Gas Approach

We consider both equilibrium and kinetic aspects of the phase separation (``thermal faceting") of thermodynamically unstable crystal surfaces into a hill--valley structure. The model we study is an Ising lattice gas for a simple cubic crystal with nearest--neighbor attractive interactions and weak next--nearest--neighbor repulsive interactions. It is likely applicable to alkali halides with the sodium chloride structure. Emphasis is placed on the fact that the equilibrium crystal shape can be interpreted as a phase diagram and that the details of its structure tell us into which surface orientations an unstable surface will decompose. We find that, depending on the temperature and growth conditions, a number of interesting behaviors are expected. For a crystal in equilibrium with its vapor, these include a low temperature regime with logarithmically--slow separation into three symmetrically--equivalent facets, and a higher temperature regime where separation proceeds as a power law in time into an entire one--parameter family of surface orientations. For a crystal slightly out of equilibrium with its vapor (slow crystal growth or etching), power--law growth should be the rule at late enough times. However, in the low temperature regime, the rate of separation rapidly decreases as the chemical potential difference between crystal and vapor phases goes to zero.Comment: 16 pages (RevTex 3.0); 12 postscript figures available on request ([email protected]). Submitted to Physical Review E. SFU-JDSDJB-94-0

Exploring the gap between dynamic and constraint-based models of metabolism

Systems biology provides new approaches for metabolic engineering through the development of models and methods for simulation and optimization of microbial metabolism. Here we explore the relationship between two modeling frameworks in common use namely, dynamic models with kinetic rate laws and constraint-based flux models. We compare and analyze dynamic and constraint-based formulations of the same model of the central carbon metabolism of E. coli. Our results show that, if unconstrained, the space of steady states described by both formulations is the same. However, the imposition of parameter-range constraints can be mapped into kinetically feasible regions of the solution space for the dynamic formulation that is not readily transferable to the constraint-based formulation. Therefore, with partial kinetic parameter knowledge, dynamic models can be used to generate constraints that reduce the solution space below that identied by constraint-based models, eliminating infeasible solutions and increasing the accuracy of simulation and optimization methods.This research was supported by PhD Grants SFRH/BD/35215/2007 and SFRH/BD/25506/2005 from the Fundacao para a Ciencia e a Tecnologia (FCT) and the MIT-Portugal Program through the project "Bridging Systems and Synthetic Biology for the Development of Improved Microbial Cell Factories" (MIT-Pt/BS-BB/0082/2008)