Non Equilibrium Current Fluctuations in Stochastic Lattice Gases

We study current fluctuations in lattice gases in the macroscopic limit extending the dynamic approach for density fluctuations developed in previous articles. More precisely, we establish a large deviation principle for a space-time fluctuation j of the empirical current with a rate functional I(j). We then estimate the probability of a fluctuation of the average current over a large time interval; this probability can be obtained by solving a variational problem for the functional I. We discuss several possible scenarios, interpreted as dynamical phase transitions, for this variational problem. They actually occur in specific models. We finally discuss the time reversal properties of I and derive a fluctuation relationship akin to the Gallavotti-Cohen theorem for the entropy production.     

Poisson-Vlasov: stochastic representation and numerical codes

In this paper, we examine the effects of the gravitational field on the dynamical evolution of the cavity-field entropy and the creation of the Schr?dinger-cat state in the Jaynes-Cummings model. We consider a moving two-level atom interacting with a single mode quantized cavity-field in the presence of a classical homogeneous gravitational field. Based on an su(2) algebra, as the dynamical symmetry group of the model, we derive the reduced density operator of the cavity-field which includes the effects of the atomic motion and the gravitational field. Also, we obtain the exact solution and the approximate solution for the system-state vector, and examine the atomic dynamics. By considering the temporal evolution of the cavity-field entropy as well as the dynamics of the Q-function of the cavity-field we study the effects of the gravitational field on the generation of the Schr?dinger-cat states of the cavity-field by using the Q-function, field entropy and approximate solution for the system-state vector. The results show that the gravitational field destroys the generation of the Schr?dinger-cat state of the cavity-field.     

Two-dimensional Brownian vortices

We introduce a stochastic model of 2D Brownian vortices associated with the canonical ensemble. The point vortices evolve through their usual mutual advection but they experience in addition a random velocity and a systematic drift generated by the system as a whole. The statistical equilibrium state of this stochastic model is the Gibbs canonical distribution. We consider a single species system and a system made of two types of vortices with positive and negative circulations. At positive temperatures, like-sign vortices repel each other (“plasma” case) and at negative temperatures, like-sign vortices attract each other (“gravity” case). We derive the stochastic equation satisfied by the exact vorticity field and the Fokker-Planck equation satisfied by the N-body distribution function. We present the BBGKY-like hierarchy of equations satisfied by the reduced distribution functions and close the hierarchy by considering an expansion of the solutions in powers of 1/N, where N is the number of vortices, in a proper thermodynamic limit. For spatially inhomogeneous systems, we derive the kinetic equations satisfied by the smooth vorticity field in a mean field approximation valid for N→+∞. For spatially homogeneous systems, we study the two-body correlation function, in a Debye-Hückel approximation valid at the order O(1/N). The results of this paper can also apply to other systems of random walkers with long-range interactions such as self-gravitating Brownian particles and bacterial populations experiencing chemotaxis. Furthermore, for positive temperatures, our study provides a kinetic derivation, from microscopic stochastic processes, of the Debye-Hückel model of electrolytes.     

Structural response of stiffened plates in similitude under a turbulent boundary layer excitation

In structural dynamics, similitude laws usually deal with simple configurations as thin flat plates with point forces. Only recently, few papers have analyzed stiffened shells or stochastic pressure loads.This research activity extends the applicability of some similitude laws, developed for thin flat plates under a turbulent boundary layer load, to ribbed plates forced by the same wall pressure fluctuations.The work addresses the problem of designing a scaled experimental test-article and, successively, of re-modulating the measured data in order to get the structural response of an original (unscaled) configuration.Due to the complexity of the structural domain, the design of a scaled configuration leads to a distorted similitude. Then, a simple approach, to circumvent the distortion effects, is proposed.     

Stochastic method to control Mycobacterium tuberculosis epidemic

We study the origin of TB (tuberculosis) epidemic and complex distributions of various populations of TB infection within the stochastic framework. The stochastic nature of this disease infection could be linked to the stochastic behaviour at genome level which is exhibited in SNP (single nucleotide polymorphism) distributions of experimentally identified hotspot driver genes. Our results show the emergence of random clusters, and well-defined discrete domains of the respective species populations in the model driven by demographic stochasticity and intrinsic complex species interaction. The multifractal analysis of the time series of the species populations indicate that TB epidemic could be mainly caused by contact communication and is directional. We propose that any TB epidemic may have high chance of approximately periodic recurrence and can be controlled by optimizing some of the parameters involved in the system modelling.     

Tracer diffusion in lattice gases

It has been proved that a tracer particle in a reversible lattice gas converges to Brownian motion. However, only in a few particular cases has a strictly positive self-diffusion coefficientD been established. Here we supply the missing piece and show thatD>0 in general. The exceptions are one-dimensional lattice gases with nearest neighbor jumps only, for whichD=0. The proof establishes a variational formula forD which could be used to obtain realistic bounds.     

The empirical behavior of sampling methods for stochastic programming

We investigate the quality of solutions obtained from sample-average approximations to two-stage stochastic linear programs with recourse. We use a recently developed software tool executing on a computational grid to solve many large instances of these problems, allowing us to obtain high-quality solutions and to verify optimality and near-optimality of the computed solutions in various ways. Research supported by the Mathematical, Information, and Computational Sciences Division subprogram of the Office of Advanced Scientific Computing Research, U.S. Department of Energy, under Contract W-31-109-Eng-38, and by the National Science Foundation under Grant 9726385. Research supported by the Mathematical, Information, and Computational Sciences Division subprogram of the Office of Advanced Scientific Computing Research, U.S. Department of Energy, under Contract W-31-109-Eng-38, and by the National Science Foundation under Grant DMS-0073770. Research supported by the Mathematical, Information, and Computational Sciences Division subprogram of the Office of Advanced Scientific Computing Research, U.S. Department of Energy, under Contract W-31-109-Eng-38, and by the National Science Foundation under Grants 9726385 and 0082065.     

Modelling the dynamics of stochastic local search on <Emphasis Type="Italic">k</Emphasis>-<Emphasis Type="SmallCaps">sat</Emphasis>

A new analytical tool is presented to provide a better understanding of the search space of k-sat. This tool, termed the local value distribution , describes the probability of finding assignments of any value q′ in the neighbourhood of assignments of value q. The local value distribution is then used to define a Markov model to model the dynamics of a corresponding stochastic local search algorithm for k-sat. The model is evaluated by comparing the predicted algorithm dynamics to experimental results. In most cases the fit of the model to the experimental results is very good, but limitations are also recognised.     

Infinite-Dimensional Black-Scholes Equation with Hereditary Structure

This paper considers the option pricing problem for contingent claims of the European type in a (B,S)-market in which the stock price and the asset in the riskless bank account both have hereditary structures. The Black-Scholes equation for the classical option pricing problem is generalized to an infinite-dimensional equation to include the effects of time delay in the evolution of the financial market as well as a very general payoff function. A computational algorithm for the solution is also obtained via a double sequence of polynomials of a certain bounded linear functional on a Banach space and the time variable.     

Improved Estimation and Uncertainty Quantification Using Monte Carlo-Based Optimization Algorithms

In this article, we present a novel method to obtain both improved estimates and reliable stopping rules for stochastic optimization algorithms such as the Monte Carlo EM (MCEM) algorithm. By characterizing a stationary point, θ*, of the algorithm as the solution to a fixed point equation, we provide a parameter estimation procedure by solving for the fixed point of the update mapping. We investigate various ways to model the update mapping, including the use of a local linear (regression) smoother. This simple approach allows increased stability in estimating the value of θ* as well as providing a natural quantification of the estimation uncertainty. These uncertainty measures can then also be used to construct convergence criteria that reflect the inherent randomness in the algorithm. We establish convergence properties of our modified estimator. In contrast to existing literature, our convergence results do not require the Monte Carlo sample size to go to infinity. Simulation studies are provided to illustrate the improved stability and reliability of our estimator.