Relaxation from a marginal state in optical bistability

Characteristic decay times for relaxation close to the marginal point of optical bistability are studied. A model-independent formula for the decay time is given which interpolates between Kramers time for activated decay and a deterministic relaxation time. This formula gives the decay time as a universal scaling function of the parameter which measures deviation from marginality. The standard deviation of the first-passage-time distribution is found to vary linearly with the decay time, close to marginality, with a slope independent of the noise intensity. Our results are substantiated by numerical simulations and their experimental relevance is pointed out.

Noise-induced scenario for inverted phase diagrams

We introduce a class of exactly solvable models exhibiting an ordering noise-induced phase transition in which order arises as a result of a balance between the relaxing deterministic dynamics and the randomizing character of the fluctuations. A finite-size scaling analysis of the phase transition reveals that it belongs to the universality class of the equilibrium Ising model. All these results are analyzed in the light of the nonequilibrium probability distribution of the system, which can be obtained analytically. Our results could constitute a possible scenario of inverted phase diagrams in the so-called lower critical solution temperature transitions.

Effects of external noise on the Swift-Hohenberg Equation

The Swift-Hohenberg equation is studied in the presence of a multiplicative noise. This stochastic equation could describe a situation in which a noise has been superimposed on the temperature gradient between the two plates of a Rayleigh-Bnard cell. A linear stability analysis and numerical simulations show that, in constrast to the additive-noise case, convective structures appear in a regime in which a deterministic analysis predicts a homogeneous solution.

External Fluctuations in Front Propagation

We study the effects of external noise in a one-dimensional model of front propagation. Noise is introduced through the fluctuations of a control parameter leading to a multiplicative stochastic partial differential equation. Analytical and numerical results for the front shape and velocity are presented. The linear-marginal-stability theory is found to increase its range of validity in the presence of external noise. As a consequence noise can stabilize fronts not allowed by the deterministic equation.

Exported oscillator competition: A concept to analyse complex rhytms

We study the interaction between two independent nonlinear oscillators competing through a neutral excitable element. The first oscillator, completely deterministic, acts as a normal pacemaker sending pulses to the neutral element which fires when it is excited by these pulses. The second oscillator, endowed with some randomness, though unable to make the excitable element to beat, leads to the occasional suppression of its firing. The missing beats or errors are registered and their statistics analyzed in terms of the noise intensity and the periods of both oscillators. This study is inspired in some complex rhythms such as a particular class of heart arrhythmia.

Consistência hidrológica de modelos digitais de elevação (MDE) para definição da rede de drenagem na sub-bacia do horto florestal Terra Dura, Eldorado do Sul, RS

O desenvolvimento da erosão hídrica ocorre em resposta ao modo como a água se move através e sobre uma determinada paisagem. O modelo digital de elevação (MDE) deve, portanto, ser o mais preciso possível, uma vez que constitui a base para a análise do relevo. Este trabalho teve como objetivo definir um modelo digital de elevação hidrologicamente consistente (MDEHC) e o método de direção de fluxo mais adequado para a definição da rede de drenagem na sub-bacia do horto florestal Terra Dura, município de Eldorado do Sul, RS. Foram testados os modelos gerados com os interpoladores Topogrid e redes triangulares irregulares (Triangulated Irregular Network -TIN) linear (TIN L) e TIN natural neighbor (TIN NN). A qualidade em relação às análises hidrológicas foi avaliada por meio da comparação das curvas de nível geradas pelos modelos testados com as curvas originais da sub-bacia (escala 1:10.000); da avaliação da quantidade de áreas planas; e da comparação da drenagem gerada pelos modelos a partir dos métodos de direção de fluxo Deterministic (D8) e Deterministic infinity (D∞ ou D infinito) com a drenagem original. Entre os modelos avaliados, o Topogrid apresentou maior consistência hidrológica, verificada na melhor continuidade das curvas de nível (menos arestas) e maior detalhamento da área de drenagem e divisores, acarretando menor quantidade de áreas planas e caminhos de fluxo mais detalhados, independentemente do método de direção de fluxo utilizado. Em relação à rede de drenagem, o método distribuído D∞ obteve melhor desempenho na descrição dos caminhos de fluxo, comparado ao método de direção única D8. O MDEHC Topogrid associado ao método D∞ proporcionou a identificação mais precisa dos caminhos preferenciais do fluxo que formam a rede de drenagem.

Excitability transitions and wave dynamics under spatiotemporal structured noise

We present an analytic and numerical study of the effects of external fluctuations in active media. Our analytical methodology transforms the initial stochastic partial differential equations into an effective set of deterministic reaction-diffusion equations. As a result we are able to explain and make quantitative predictions on the systematic and constructive effects of the noise, for example, target patterns created out of noise and traveling or spiral waves sustained by noise. Our study includes the case of realistic noises with temporal and spatial structures.

Front dynamics in the presence of spatio-temporal structured noises

Front dynamics modeled by a reaction-diffusion equation are studied under the influence of spatiotemporal structured noises. An effective deterministic model is analytical derived where the noise parameters, intensity, correlation time, and correlation length appear explicitly. The different effects of these parameters are discussed for the Ginzburg-Landau and Schlögl models. We obtain an analytical expression for the front velocity as a function of the noise parameters. Numerical simulation results are in a good agreement with the theoretical predictions.

Intensity correlation functions of dye lasers: Comparison of colored-gain-noise and colored-loss-noise models

The intensity correlation functions C(t) for the colored-gain-noise model of dye lasers are analyzed and compared with those for the loss-noise model. For correlation times ¿ larger than the deterministic relaxation time td, we show with the use of the adiabatic approximation that C(t) values coincide for both models. For small correlation times we use a method that provides explicit expressions of non-Markovian correlation functions, approximating simultaneously short- and long-time behaviors. Comparison with numerical simulations shows excellent results simultaneously for short- and long-time regimes. It is found that, when the correlation time of the noise increases, differences between the gain- and loss-noise models tend to disappear. The decay of C(t) for both models can be described by a time scale that approaches the deterministic relaxation time. However, in contrast with the loss-noise model, a secondary time scale remains for large times for the gain-noise model, which could allow one to distinguish between both models.

Estimating dormant and active hematopoietic stem cell kinetics through extensive modeling of bromodeoxyuridine label-retaining cell dynamics.

Bone marrow hematopoietic stem cells (HSCs) are responsible for both lifelong daily maintenance of all blood cells and for repair after cell loss. Until recently the cellular mechanisms by which HSCs accomplish these two very different tasks remained an open question. Biological evidence has now been found for the existence of two related mouse HSC populations. First, a dormant HSC (d-HSC) population which harbors the highest self-renewal potential of all blood cells but is only induced into active self-renewal in response to hematopoietic stress. And second, an active HSC (a-HSC) subset that by and large produces the progenitors and mature cells required for maintenance of day-to-day hematopoiesis. Here we present computational analyses further supporting the d-HSC concept through extensive modeling of experimental DNA label-retaining cell (LRC) data. Our conclusion that the presence of a slowly dividing subpopulation of HSCs is the most likely explanation (amongst the various possible causes including stochastic cellular variation) of the observed long term Bromodeoxyuridine (BrdU) retention, is confirmed by the deterministic and stochastic models presented here. Moreover, modeling both HSC BrdU uptake and dilution in three stages and careful treatment of the BrdU detection sensitivity permitted improved estimates of HSC turnover rates. This analysis predicts that d-HSCs cycle about once every 149-193 days and a-HSCs about once every 28-36 days. We further predict that, using LRC assays, a 75%-92.5% purification of d-HSCs can be achieved after 59-130 days of chase. Interestingly, the d-HSC proportion is now estimated to be around 30-45% of total HSCs - more than twice that of our previous estimate.

Giant acceleration of free diffusion by use of titled periodic potentials

The effective diffusion coefficient for the overdamped Brownian motion in a tilted periodic potential is calculated in closed analytical form. Universality classes and scaling properties for weak thermal noise are identified near the threshold tilt where deterministic running solutions set in. In this regime the diffusion may be greatly enhanced, as compared to free thermal diffusion with, for a realistic experimental setup, an enhancement of up to 14 orders of magnitude.

Difussion in a titled periodic potential: Enhancement, universality, and scaling

An exact analytical expression for the effective diffusion coefficient of an overdamped Brownian particle in a tilted periodic potential is derived for arbitrary potentials and arbitrary strengths of the thermal noise. Near the critical tilt (threshold of deterministic running solutions) a scaling behavior for weak thermal noise is revealed and various universality classes are identified. In comparison with the bare (potential-free) thermal diffusion, the effective diffusion coefficient in a critically tilted periodic potential may be, in principle, arbitrarily enhanced. For a realistic experimental setup, an enhancement by 14 orders of magnitude is predicted so that thermal diffusion should be observable on a macroscopic scale at room temperature.

Estimation of the correlation structure of crustal velocity heterogeneity from seismic reflection data

Numerous sources of evidence point to the fact that heterogeneity within the Earth's deep crystalline crust is complex and hence may be best described through stochastic rather than deterministic approaches. As seismic reflection imaging arguably offers the best means of sampling deep crustal rocks in situ, much interest has been expressed in using such data to characterize the stochastic nature of crustal heterogeneity. Previous work on this problem has shown that the spatial statistics of seismic reflection data are indeed related to those of the underlying heterogeneous seismic velocity distribution. As of yet, however, the nature of this relationship has remained elusive due to the fact that most of the work was either strictly empirical or based on incorrect methodological approaches. Here, we introduce a conceptual model, based on the assumption of weak scattering, that allows us to quantitatively link the second-order statistics of a 2-D seismic velocity distribution with those of the corresponding processed and depth-migrated seismic reflection image. We then perform a sensitivity study in order to investigate what information regarding the stochastic model parameters describing crustal velocity heterogeneity might potentially be recovered from the statistics of a seismic reflection image using this model. Finally, we present a Monte Carlo inversion strategy to estimate these parameters and we show examples of its application at two different source frequencies and using two different sets of prior information. Our results indicate that the inverse problem is inherently non-unique and that many different combinations of the vertical and lateral correlation lengths describing the velocity heterogeneity can yield seismic images with the same 2-D autocorrelation structure. The ratio of all of these possible combinations of vertical and lateral correlation lengths, however, remains roughly constant which indicates that, without additional prior information, the aspect ratio is the only parameter describing the stochastic seismic velocity structure that can be reliably recovered.

Development of a Computer Controlled Underbody Plow: TR412, Final Report, January 2006

Underbody plows can be very useful tools in winter maintenance, especially when compacted snow or hard ice must be removed from the roadway. By the application of significant down-force, and the use of an appropriate cutting edge angle, compacted snow and ice can be removed very effectively by such plows, with much greater efficiency than any other tool under those circumstances. However, the successful operation of an underbody plow requires considerable skill. If too little down pressure is applied to the plow, then it will not cut the ice or compacted snow. However, if too much force is applied, then either the cutting edge may gouge the road surface, causing significant damage often to both the road surface and the plow, or the plow may ride up on the cutting edge so that it is no longer controllable by the operator. Spinning of the truck in such situations is easily accomplished. Further, excessive down force will result in rapid wear of the cutting edge. Given this need for a high level of operator skill, the operation of an underbody plow is a candidate for automation. In order to successfully automate the operation of an underbody plow, a control system must be developed that follows a set of rules that represent appropriate operation of such a plow. These rules have been developed, based upon earlier work in which operational underbody plows were instrumented to determine the loading upon them (both vertical and horizontal) and the angle at which the blade was operating.These rules have been successfully coded into two different computer programs, both using the MatLab® software. In the first program, various load and angle inputs are analyzed to determine when, whether, and how they violate the rules of operation. This program is essentially deterministic in nature. In the second program, the Simulink® package in the MatLab® software system was used to implement these rules using fuzzy logic. Fuzzy logic essentially replaces a fixed and constant rule with one that varies in such a way as to improve operational control. The development of the fuzzy logic in this simulation was achieved simply by using appropriate routines in the computer software, rather than being developed directly. The results of the computer testing and simulation indicate that a fully automated, computer controlled underbody plow is indeed possible. The issue of whether the next steps toward full automation should be taken (and by whom) has also been considered, and the possibility of some sort of joint venture between a Department of Transportation and a vendor has been suggested.

Option valuation as an expectation in the complex domain: The Black-Scholes case

It is very well known that the first succesful valuation of a stock option was done by solving a deterministic partial differential equation (PDE) of the parabolic type with some complementary conditions specific for the option. In this approach, the randomness in the option value process is eliminated through a no-arbitrage argument. An alternative approach is to construct a replicating portfolio for the option. From this viewpoint the payoff function for the option is a random process which, under a new probabilistic measure, turns out to be of a special type, a martingale. Accordingly, the value of the replicating portfolio (equivalently, of the option) is calculated as an expectation, with respect to this new measure, of the discounted value of the payoff function. Since the expectation is, by definition, an integral, its calculation can be made simpler by resorting to powerful methods already available in the theory of analytic functions. In this paper we use precisely two of those techniques to find the well-known value of a European call