Equilibrium distributions of discrete non-autonomous graphs

We introduce the notions of equilibrium distribution and time of convergence in discrete non-autonomous graphs. Under some conditions we give an estimate to the convergence time to the equilibrium distribution using the second largest eigenvalue of some matrices associated with the system.

Properties of Equilibrium Distributions of Order n

The present work is intended to discuss various properties and reliability aspects of higher order equilibrium distributions in continuous, discrete and multivariate cases, which contribute to the study on equilibrium distributions. At first, we have to study and consolidate the existing literature on equilibrium distributions. For this we need some basic concepts in reliability. These are being discussed in the 2nd chapter, In Chapter 3, some identities connecting the failure rate functions and moments of residual life of the univariate, non-negative continuous equilibrium distributions of higher order and that of the baseline distribution are derived. These identities are then used to characterize the generalized Pareto model, mixture of exponentials and gamma distribution. An approach using the characteristic functions is also discussed with illustrations. Moreover, characterizations of ageing classes using stochastic orders has been discussed. Part of the results of this chapter has been reported in Nair and Preeth (2009). Various properties of equilibrium distributions of non-negative discrete univariate random variables are discussed in Chapter 4. Then some characterizations of the geo- metric, Waring and negative hyper-geometric distributions are presented. Moreover, the ageing properties of the original distribution and nth order equilibrium distribu- tions are compared. Part of the results of this chapter have been reported in Nair, Sankaran and Preeth (2012). Chapter 5 is a continuation of Chapter 4. Here, several conditions, in terms of stochastic orders connecting the baseline and its equilibrium distributions are derived. These conditions can be used to rede_ne certain ageing notions. Then equilibrium distributions of two random variables are compared in terms of various stochastic orders that have implications in reliability applications. In Chapter 6, we make two approaches to de_ne multivariate equilibrium distribu- tions of order n. Then various properties including characterizations of higher order equilibrium distributions are presented. Part of the results of this chapter have been reported in Nair and Preeth (2008). The Thesis is concluded in Chapter 7. A discussion on further studies on equilib- rium distributions is also made in this chapter.

Equilibrium distributions of topological states in circular DNA: Interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and ΔLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(ΔLk, K), is related to two conditional distributions: P(ΔLk|K), the distribution of ΔLk for a particular K, and P(K|ΔLk) and also to two simple distributions: P(ΔLk), the distribution of ΔLk irrespective of K, and P(K). We explored the relationships between these distributions. P(ΔLk, K), P(ΔLk), and P(K|ΔLk) were calculated from the simulated distributions of P(ΔLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|ΔLk), the distribution of knot types for a particular value of ΔLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of ΔLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |ΔLk| > 12, only one or two knot types dominate the P(K|ΔLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations

We describe and test a Markov chain model of microsatellite evolution that can explain the different distributions of microsatellite lengths across different organisms and repeat motifs. Two key features of this model are the dependence of mutation rates on microsatellite length and a mutation process that includes both strand slippage and point mutation events. We compute the stationary distribution of allele lengths under this model and use it to fit DNA data for di-, tri-, and tetranucleotide repeats in humans, mice, fruit flies, and yeast. The best fit results lead to slippage rate estimates that are highest in mice, followed by humans, then yeast, and then fruit flies. Within each organism, the estimates are highest in di-, then tri-, and then tetranucleotide repeats. Our estimates are consistent with experimentally determined mutation rates from other studies. The results suggest that the different length distributions among organisms and repeat motifs can be explained by a simple difference in slippage rates and that selective constraints on length need not be imposed.

Convergence time to equilibrium distributions of autonomous and periodic non-autonomous graphs

We present some estimates of the time of convergence to the equilibrium distribution in autonomous and periodic non-autonomous graphs, with ergodic stochastic adjacency matrices, using the eigenvalues of these matrices. On this way we generalize previous results from several authors, that only considered reversible matrices.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Characterizations of some continuous distributions using partial moments

Inthis paper,we define partial moments for a univariate continuous random variable. A recurrence relationship for the Pearson curve using the partial moments is established. The interrelationship between the partial moments and other reliability measures such as failure rate, mean residual life function are proved. We also prove some characterization theorems using the partial moments in the context of length biased models and equilibrium distributions

Analysis of thermodynamic non-ideality in terms of protein solvation

The effects of thermodynamic non-ideality on the forms of sedimentation equilibrium distributions for several isoelectric proteins have been analysed on the statistical-mechanical basis of excluded volume to obtain an estimate of the extent of protein solvation. Values of the effective solvation. parameter delta are reported for ellipsoidal as well as spherical models of the proteins, taken to be rigid, impenetrable macromolecular structures. The dependence of the effective solvated radius upon protein molecular mass exhibits reasonable agreement with the relationship calculated for a model in which the unsolvated protein molecule is surrounded by a 0.52-nm solvation shell. Although the observation that this shell thickness corresponds to a double layer of water molecules may be of questionable relevance to mechanistic interpretation of protein hydration, it augurs well for the assignment of magnitudes to the second virial coefficients of putative complexes in the quantitative characterization of protein-protein interactions under conditions where effects of thermodynamic non-ideality cannot justifiably be neglected. (C) 2001 Elsevier Science B.V. All rights reserved.

Ultracentrifugal studies of the effect of molecular crowding by trimethylamine N-oxide on the self-association of muscle glycogen phosphorylase b

The suitability of sedimentation equilibrium for characterizing the self-association of muscle glycogen phosphorylase b has been reappraised. Whereas sedimentation equilibrium distributions for phosphorylase b in 40 mM Hepes buffer (pH 6.8) supplemented with 1 mM AMP signify a lack of chemical equilibrium attainment, those in buffer supplemented additionally with potassium sulfate conform with the requirements of a dimerizing system in chemical as we:ll as sedimentation equilibrium. Because the rate of attainment of chemical equilibrium under the former conditions is sufficiently slow to allow resolution of the dimeric and tetrameric enzyme species by sedimentation velocity, this procedure has been used to examine the effects of thermodynamic nonideality arising from molecular crowding try trimethylamine N-oxide on the self-association behaviour of phosphorylase b. In those terms the marginally enhanced extent of phosphorylase b self-association observed in the presence of high concentrations of the cosolute is taken to imply that the effects of thermodynamic nonideality on the dimer-tetramer equilibrium are being countered by those displacing the T reversible arrow R isomerization equilibrium for dimer towards the smaller, nonassociating T state. Because the R state is the enzymically active form, an inhibitory effect is the predicted consequence of molecular crowding by high concentrations of unrelated solutes. Thermodynamic nonideality thus provides an alternative explanation for the inhibitory effects of high concentrations of glycerol, sucrose and ethylene glycol on phosphorylase b activity, phenomena that have been attributed to extremely weak interaction of these cryoprotectants with the T state of the enzyme.

On the existence of pure-strategy equilibria in large games

Over the years, several formalizations and existence results for games with a continuum of players have been given. These include those of Schmeidler (1973), Rashid (1983), Mas-Colell (1984), Khan and Sun (1999) and Podczeck (2007a). The level of generality of each of these existence results is typically regarded as a criterion to evaluate how appropriate is the corresponding formalization of large games. In contrast, we argue that such evaluation is pointless. In fact, we show that, in a precise sense, all the above existence results are equivalent. Thus, all of them are equally strong and therefore cannot rank the different formalizations of large games.

Comparison of information structures

We introduce two ways of comparing information structures, say ${\cal I}$ and${\cal J}$. First we say that ${\cal I}$ is richer than ${\cal J}$ when forevery compact game $G$, all correlated equilibrium distributions of $G$ inducedby ${\cal J}$ are also induced by ${\cal I}$. Second, we say that ${\cal J}$is faithfully reproducable from ${\cal I}$ when all the players can computefrom their information in ${\cal I}$ ``new information'' that they could havereceived from ${\cal J}$. We prove that ${\cal I}$ is richer than ${\cal J}$if and only if ${\cal J}$ is faithfully reproducable from ${\cal I}$.

The role of lower partial moments in stochastic modeling

Lower partial moments plays an important role in the analysis of risks and in income/poverty studies. In the present paper, we further investigate its importance in stochastic modeling and prove some characterization theorems arising out of it. We also identify its relationships with other important applied models such as weighted and equilibrium models. Finally, some applications of lower partial moments in poverty studies are also examined

DIFFUSION OF MAGNETIC FIELD AND REMOVAL OF MAGNETIC FLUX FROM CLOUDS VIA TURBULENT RECONNECTION

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence depends on whether magnetic fields can change their topology via reconnection in highly conducting media. Recent progress in understanding fast magnetic reconnection in the presence of turbulence reassures that the magnetic field behavior in computer simulations and turbulent astrophysical environments is similar, as far as magnetic reconnection is concerned. This makes it meaningful to perform MHD simulations of turbulent flows in order to understand the diffusion of magnetic field in astrophysical environments. Our studies of magnetic field diffusion in turbulent medium reveal interesting new phenomena. First of all, our three-dimensional MHD simulations initiated with anti-correlating magnetic field and gaseous density exhibit at later times a de-correlation of the magnetic field and density, which corresponds well to the observations of the interstellar media. While earlier studies stressed the role of either ambipolar diffusion or time-dependent turbulent fluctuations for de-correlating magnetic field and density, we get the effect of permanent de-correlation with one fluid code, i.e., without invoking ambipolar diffusion. In addition, in the presence of gravity and turbulence, our three-dimensional simulations show the decrease of the magnetic flux-to-mass ratio as the gaseous density at the center of the gravitational potential increases. We observe this effect both in the situations when we start with equilibrium distributions of gas and magnetic field and when we follow the evolution of collapsing dynamically unstable configurations. Thus, the process of turbulent magnetic field removal should be applicable both to quasi-static subcritical molecular clouds and cores and violently collapsing supercritical entities. The increase of the gravitational potential as well as the magnetization of the gas increases the segregation of the mass and magnetic flux in the saturated final state of the simulations, supporting the notion that the reconnection-enabled diffusivity relaxes the magnetic field + gas system in the gravitational field to its minimal energy state. This effect is expected to play an important role in star formation, from its initial stages of concentrating interstellar gas to the final stages of the accretion to the forming protostar. In addition, we benchmark our codes by studying the heat transfer in magnetized compressible fluids and confirm the high rates of turbulent advection of heat obtained in an earlier study.

On exact time averages of a massive Poisson particle

In this work we study, under the Stratonovich definition, the problem of the damped oscillatory massive particle subject to a heterogeneous Poisson noise characterized by a rate of events, lambda(t), and a magnitude, Phi, following an exponential distribution. We tackle the problem by performing exact time averages over the noise in a similar way to previous works analysing the problem of the Brownian particle. From this procedure we obtain the long-term equilibrium distributions of position and velocity as well as analytical asymptotic expressions for the injection and dissipation of energy terms. Considerations on the emergence of stochastic resonance in this type of system are also set forth.

Dynamics and folding of single two-stranded coiled-coil peptides studied by fluorescent energy transfer confocal microscopy

We report single-molecule measurements on the folding and unfolding conformational equilibrium distributions and dynamics of a disulfide crosslinked version of the two-stranded coiled coil from GCN4. The peptide has a fluorescent donor and acceptor at the N termini of its two chains and a Cys disulfide near its C terminus. Thus, folding brings the two N termini of the two chains close together, resulting in an enhancement of fluorescent resonant energy transfer. End-to-end distance distributions have thus been characterized under conditions where the peptide is nearly fully folded (0 M urea), unfolded (7.4 M urea), and in dynamic exchange between folded and unfolded states (3.0 M urea). The distributions have been compared for the peptide freely diffusing in solution and deposited onto aminopropyl silanized glass. As the urea concentration is increased, the mean end-to-end distance shifts to longer distances both in free solution and on the modified surface. The widths of these distributions indicate that the molecules are undergoing millisecond conformational fluctuations. Under all three conditions, these fluctuations gave nonexponential correlations on 1- to 100-ms time scale. A component of the correlation decay that was sensitive to the concentration of urea corresponded to that measured by bulk relaxation kinetics. The trajectories provided effective intramolecular diffusion coefficients as a function of the end-to-end distances for the folded and unfolded states. Single-molecule folding studies provide information concerning the distributions of conformational states in the folded, unfolded, and dynamically interconverting states.