Point-occurrence self-similarity in crackling-noise systems and in other complex systems

It has been recently found that a number of systems displaying crackling noise also show a remarkable behavior regarding the temporal occurrence of successive events versus their size: a scaling law for the probability distributions of waiting times as a function of a minimum size is fulfilled, signaling the existence on those systems of self-similarity in time-size. This property is also present in some non-crackling systems. Here, the uncommon character of the scaling law is illustrated with simple marked renewal processes, built by definition with no correlations. Whereas processes with a finite mean waiting time do not fulfill a scaling law in general and tend towards a Poisson process in the limit of very high sizes, processes without a finite mean tend to another class of distributions, characterized by double power-law waiting-time densities. This is somehow reminiscent of the generalized central limit theorem. A model with short-range correlations is not able to escape from the attraction of those limit distributions. A discussion on open problems in the modeling of these properties is provided.

On dome-shaped norms of reaction for size-to-age at maturity in fishes

The optimal size-to-age at maturity depends on growth and mortality rates, which vary with environment. Therefore, organisms in spatially or temporaly changing environments should develop adaptative phenotypic plasticity for this trait. Experimental work by Alm (1959) on several fish species shows a dome-shape norm of reaction for size-to-age at maturity: size at maturity is smaller in both fast-growing and slow-growing fishes, than it is in fish with a medium growth rate. Using computer simulations, we show that such a dome-shaped norm of reaction is optimal when assuming a finite life span and a negative relationship between production and survival rates. This latter assumption is supported by empirical data, as well as by physiological and emographic arguments.

Random mating with a finite number of matings.

Random mating is the null model central to population genetics. One assumption behind random mating is that individuals mate an infinite number of times. This is obviously unrealistic. Here we show that when each female mates a finite number of times, the effective size of the population is substantially decreased.

Scaling behavior and equilibrium lengths of knotted polymers

We use numerical simulations to investigate how the chain length and topology of freely fluctuating knotted polymer rings affect their various spatial characteristics such as the radius of the smallest sphere enclosing momentary configurations of simulated polymer chains. We describe how the average value of a characteristic changes with the chain size and how this change depends on the topology of the modeled polymers. Although the scaling profiles of a spatial characteristic for distinct knot types do not intersect (at least, in the range of our data), the profiles for nontrivial knots intersect the corresponding profile obtained for phantom polymers, i.e., those that are free to explore all available topological states. For each knot type, this point of intersection defines its equilibrium length with respect to the spatial characteristic. At this chain length, a polymer forming a given knot type will not tend to increase or decrease. on average, the value of the spatial characteristic when the polymer is released from its topological constraint. We show interrelations between equilibrium lengths defined with respect to spatial characteristics of different character and observe that they are related to the lengths of ideal geometric configurations of the corresponding knot types.

Some hypothesis tests for the covariance matrix when the dimension is large compared to the sample size

This paper analyzes whether standard covariance matrix tests work whendimensionality is large, and in particular larger than sample size. Inthe latter case, the singularity of the sample covariance matrix makeslikelihood ratio tests degenerate, but other tests based on quadraticforms of sample covariance matrix eigenvalues remain well-defined. Westudy the consistency property and limiting distribution of these testsas dimensionality and sample size go to infinity together, with theirratio converging to a finite non-zero limit. We find that the existingtest for sphericity is robust against high dimensionality, but not thetest for equality of the covariance matrix to a given matrix. For thelatter test, we develop a new correction to the existing test statisticthat makes it robust against high dimensionality.

Dpp signaling activity requires Pentagone to scale with tissue size in the growing Drosophila wing imaginal disc.

The wing of the fruit fly, Drosophila melanogaster, with its simple, two-dimensional structure, is a model organ well suited for a systems biology approach. The wing arises from an epithelial sac referred to as the wing imaginal disc, which undergoes a phase of massive growth and concomitant patterning during larval stages. The Decapentaplegic (Dpp) morphogen plays a central role in wing formation with its ability to co-coordinately regulate patterning and growth. Here, we asked whether the Dpp signaling activity scales, i.e. expands proportionally, with the growing wing imaginal disc. Using new methods for spatial and temporal quantification of Dpp activity and its scaling properties, we found that the Dpp response scales with the size of the growing tissue. Notably, scaling is not perfect at all positions in the field and the scaling of target gene domains is ensured specifically where they define vein positions. We also found that the target gene domains are not defined at constant concentration thresholds of the downstream Dpp activity gradients P-Mad and Brinker. Most interestingly, Pentagone, an important secreted feedback regulator of the pathway, plays a central role in scaling and acts as an expander of the Dpp gradient during disc growth.

Optimal growth strategies when mortality and production rates are size-dependent

Pontryagin's maximum principle from optimal control theory is used to find the optimal allocation of energy between growth and reproduction when lifespan may be finite and the trade-off between growth and reproduction is linear. Analyses of the optimal allocation problem to date have generally yielded bang-bang solutions, i.e. determinate growth: life-histories in which growth is followed by reproduction, with no intermediate phase of simultaneous reproduction and growth. Here we show that an intermediate strategy (indeterminate growth) can be selected for if the rates of production and mortality either both increase or both decrease with increasing body size, this arises as a singular solution to the problem. Our conclusion is that indeterminate growth is optimal in more cases than was previously realized. The relevance of our results to natural situations is discussed.

Sequential partitioning: an alternative to understanding size distributions of avalanches in first-order phase transitions

We study the problem of the partition of a system of initial size V into a sequence of fragments s1,s2,s3 . . . . By assuming a scaling hypothesis for the probability p(s;V) of obtaining a fragment of a given size, we deduce that the final distribution of fragment sizes exhibits power-law behavior. This minimal model is useful to understanding the distribution of avalanche sizes in first-order phase transitions at low temperatures.

Scaling calculation of isoscalar giant resonances in relativistic Thomas-Fermi theory

We derive analytical expressions for the excitation energy of the isoscalar giant monopole and quadrupole resonances in finite nuclei, by using the scaling method and the extended ThomasFermi approach to relativistic mean-field theory. We study the ability of several nonlinear σω parameter sets of common use in reproducing the experimental data. For monopole oscillations the calculations agree better with experiment when the nuclear matter incompressibility of the relativistic interaction lies in the range 220260 MeV. The breathing-mode energies of the scaling method compare satisfactorily with those obtained in relativistic RPA and time-dependent mean-field calculations. For quadrupole oscillations, all the analyzed nonlinear parameter sets reproduce the empirical trends reasonably well.

Scaling in relativistic Thomas-Fermi approach for nuclei

By using the scaling method we derive the virial theorem for the relativistic mean field model of nuclei treated in the ThomasFermi approach. The ThomasFermi solutions statisfy the stability condition against scaling. We apply the formalism to study the excitation energy of the breathing mode in finite nuclei with several relativistic parameter sets of common use.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Integrated random processes exhibiting long tails, finite moments, and power-law spectra

A dynamical model based on a continuous addition of colored shot noises is presented. The resulting process is colored and non-Gaussian. A general expression for the characteristic function of the process is obtained, which, after a scaling assumption, takes on a form that is the basis of the results derived in the rest of the paper. One of these is an expansion for the cumulants, which are all finite, subject to mild conditions on the functions defining the process. This is in contrast with the Lévy distribution¿which can be obtained from our model in certain limits¿which has no finite moments. The evaluation of the spectral density and the form of the probability density function in the tails of the distribution shows that the model exhibits a power-law spectrum and long tails in a natural way. A careful analysis of the characteristic function shows that it may be separated into a part representing a Lévy process together with another part representing the deviation of our model from the Lévy process. This

Random mixtures with orientational order, and the anisotropic resistivity tensor of high Tc superconductors

By generalizing effective-medium theory to the case of orientationally ordered but positionally disordered two component mixtures, it is shown that the anisotropic dielectric tensor of oxide superconductors can be extracted from microwave measurements on oriented crystallites of YBa2Cu3O7¿x embedded in epoxy. Surprisingly, this technique appears to be the only one which can access the resistivity perpendicular to the copper¿oxide planes in crystallites that are too small for depositing electrodes. This possibility arises in part because the real part of the dielectric constant of oxide superconductors has a large magnitude. The validity of the effective-medium approach for orientationally ordered mixtures is corroborated by simulations on two¿dimensional anisotropic random resistor networks. Analysis of the experimental data suggests that the zero-temperature limit of the finite frequency resistivity does not vanish along the c axis, a result which would simply the existence of states at the Fermi surface, even in the superconducting state

Is it the grain size or the characteristic pore size that controls the induced polarization relaxation time of clean sands and sandstones?

There is a wide range of evidence to suggest that permeability can be constrained through of induced polarization measurements. For clean sands and sandstones, current mechanistic models of induced polarization predict a relationship between the low-frequency time constant inferred from induced polarization measurements and the grain diameter. A number of observations do, however, disagree with this and indicate that the observed relaxation behavior is rather governed by the so-called dynamic pore radius L. To test this hypothesis, we have developed a set of new scaling relationships, which allow the relaxation time to be computed from the pore size and the permeability to be computed from both the Cole-Cole time constant and the formation factor. Moreover, these new scaling relationships can be also used to predict the dependence of the Cole-Cole time constant as a function of the water saturation under unsaturated conditions. Comparative tests of the proposed new relationships with regard to various published experimental results for saturated clean sands and sandstones as well as for partially saturated clean sandstones, do indeed confirm that the dynamic pore radius L is a much more reliable indicator of the observed relaxation behavior than grain-size-based models.

Solid phase crystallization under continuous heating: kinetic and microstructure scaling laws

The kinetics and microstructure of solid-phase crystallization under continuous heating conditions and random distribution of nuclei are analyzed. An Arrhenius temperature dependence is assumed for both nucleation and growth rates. Under these circumstances, the system has a scaling law such that the behavior of the scaled system is independent of the heating rate. Hence, the kinetics and microstructure obtained at different heating rates differ only in time and length scaling factors. Concerning the kinetics, it is shown that the extended volume evolves with time according to αex = [exp(κCt′)]m+1, where t′ is the dimensionless time. This scaled solution not only represents a significant simplification of the system description, it also provides new tools for its analysis. For instance, it has been possible to find an analytical dependence of the final average grain size on kinetic parameters. Concerning the microstructure, the existence of a length scaling factor has allowed the grain-size distribution to be numerically calculated as a function of the kinetic parameters