Topological properties of the one dimensional exponential random geometric graph

In this article we study the one-dimensional random geometric (random interval) graph when the location of the nodes are independent and exponentially distributed. We derive exact results and limit theorems for the connectivity and other properties associated with this random graph. We show that the asymptotic properties of a graph with a truncated exponential distribution can be obtained using the exponential random geometric graph. © 2007 Wiley Periodicals, Inc. Random Struct. Alg., 2008.

Structure and dynamics of two-dimensional sheared granular flows

The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analyzed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analyzed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analyzed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

Performance reliability of high-maintenance systems

The time–history of the performance of a system is treated as a stochastic corrective process, in which deterioration due to aging is counteracted at brief maintenance checks. Using a diffusion approximation for the deterioration, simple models are proposed for describing maintenance either by component replacement or by performance restoration. Equilibrium solutions of the models show that the performance has a probability distribution with exponential tails: the uncritical use of Gaussians can grossly underestimate the probability of poor performance. The proposed models are supported by recent observational evidence on aircraft track-keeping errors, which are shown to follow the modified exponential distribution derived here. The analysis also brings out the relation between the deterioration characteristics of the system and the intensity of the maintenance effort required to achieve a given performance reliability.

Observation of trap-assisted space charge limited conductivity in short channel MoS2 transistor

We present temperature dependent I-V measurements of short channel MoS2 field effect devices at high source-drain bias. We find that, although the I-V characteristics are ohmic at low bias, the conduction becomes space charge limited at high V-DS, and existence of an exponential distribution of trap states was observed. The temperature independent critical drain-source voltage (V-c) was also determined. The density of trap states was quantitatively calculated from V-c. The possible origin of exponential trap distribution in these devices is also discussed. (C) 2013 AIP Publishing LLC.

Belief propagation for minimum weight many-to-one matchings in the random complete graph

In a complete bipartite graph with vertex sets of cardinalities n and n', assign random weights from exponential distribution with mean 1, independently to each edge. We show that, as n -> infinity, with n' = n/alpha] for any fixed alpha > 1, the minimum weight of many-to-one matchings converges to a constant (depending on alpha). Many-to-one matching arises as an optimization step in an algorithm for genome sequencing and as a measure of distance between finite sets. We prove that a belief propagation (BP) algorithm converges asymptotically to the optimal solution. We use the objective method of Aldous to prove our results. We build on previous works on minimum weight matching and minimum weight edge cover problems to extend the objective method and to further the applicability of belief propagation to random combinatorial optimization problems.

Criticality of the exponential rate of decay for the largest nearest-neighbor link in random geometric graphs

Let n points be placed independently in d-dimensional space according to the density f(x) = A(d)e(-lambda parallel to x parallel to alpha), lambda, alpha > 0, x is an element of R-d, d >= 2. Let d(n) be the longest edge length of the nearest-neighbor graph on these points. We show that (lambda(-1) log n)(1-1/alpha) d(n) - b(n) converges weakly to the Gumbel distribution, where b(n) similar to ((d - 1)/lambda alpha) log log n. We also prove the following strong law for the normalized nearest-neighbor distance (d) over tilde (n) = (lambda(-1) log n)(1-1/alpha) d(n)/log log n: (d - 1)/alpha lambda <= lim inf(n ->infinity) (d) over tilde (n) <= lim sup(n ->infinity) (d) over tilde (n) <= d/alpha lambda almost surely. Thus, the exponential rate of decay alpha = 1 is critical, in the sense that, for alpha > 1, d(n) -> 0, whereas, for alpha <= 1, d(n) -> infinity almost surely as n -> infinity.

On a Ratio of Functions of Exponential Random Variables and Some Applications

Consider L independent and identically distributed exponential random variables (r.vs) X-1, X-2 ,..., X-L and positive scalars b(1), b(2) ,..., b(L). In this letter, we present the probability density function (pdf), cumulative distribution function and the Laplace transform of the pdf of the composite r.v Z = (Sigma(L)(j=1) X-j)(2) / (Sigma(L)(j=1) b(j)X(j)). We show that the r.v Z appears in various communication systems such as i) maximal ratio combining of signals received over multiple channels with mismatched noise variances, ii)M-ary phase-shift keying with spatial diversity and imperfect channel estimation, and iii) coded multi-carrier code-division multiple access reception affected by an unknown narrow-band interference, and the statistics of the r.v Z derived here enable us to carry out the performance analysis of such systems in closed-form.

Solvation dynamics, energy distribution and trapping of a light solute ion

In the theoretical treatments of the dynamics of solvation of a newly created ion in a dipolar solvent, the self-motion of the solute is usually ignored. Recently, it has been shown that for a light ion the translational motion of the ion can significantly enhance its own rate of solvation. Therefore, solvation itself may not be the rate determining step in the equilibration. Instead, the rate determining step is the search of the low energy configuration which serves to localize the light ion. In this article a microscopic calculation of the probability distribution of the interaction energy of the nascent charge with the dipolar solvent molecules is presented in order to address this problem of solute trapping. It is found that to a good approximation, this distribution is Gaussian and the second moment of this distribution is exactly equal to the half of its own solvation energy. It is shown that this is in excellent agreement with the simulation results that are available for the model Brownian dipolar lattice and for liquid acetonitrile. If the distortion of the solvent by the ion is negligible then the same relation gives the energy distribution for the solvated ion, with the average centered at the final equilibrium solvation energy. These results are expected to be useful in understanding various chemical processes in dipolar liquids. Another interesting outcome of the present study is a simple dynamic argument that supports Onsager's ''inverse snow-ball'' conjecture of solvation of a light ion. A simple derivation of the semi-phenomenological relation between the solvation time correlation function and the single particle orientation, reported recently by Maroncelli et al. (J. Phys. Chem. 97 (1993) 13), is also presented.

Inverse surface melting in confined clusters: Ar-13 in zeolite L

Monte Carlo and molecular dynamics simulations on an Ar-13 cluster in zeolite L have been carried out at a series of temperatures to understand the rigid-nonrigid transition corresponding to the solid-liquid transition exhibited by the free Ar-13 cluster. The icosahedral geometry of the free cluster is no longer preferred when the cluster is confined in the zeolite. The root-mean-squared pair distance fluctuation, delta, exhibits a sharp, well-defined rigid-nonrigid transition at 17 K as compared to 27 K for the free cluster. Multiple peaks in the distribution of short-time averages of the guest-host interaction energy indicate coexistence of two phases.; It is shown that this transition is associated with the inner atoms becoming mobile at 17 K even while the outer layer atoms, which are in close proximity to the zeolitic wall, continue to be comparatively immobile. This may be contrasted with the melting of large free clusters of 40 or more atoms which exhibit surface melting. Guest-host interactions seem to play a predominant role in determining the properties of confined clusters. We demonstrate that the volume of the cluster increases rather sharply at 17 and 27 K respectively for the confined and the free cluster. Power spectra suggest that the motion of the inner atoms is generally parallel to the atoms which form the cage wall.

Water quality parameter estimation in steady-state distribution system

The maintenance of chlorine residual is needed at all the points in the distribution system supplied with chlorine as a disinfectant. The propagation and level of chlorine in a distribution system is affected by both bulk and pipe wall reactions. It is well known that the field determination of wall reaction parameter is difficult. The source strength of chlorine to maintain a specified chlorine residual at a target node is also an important parameter. The inverse model presented in the paper determines these water quality parameters, which are associated with different reaction kinetics, either in single or in groups of pipes. The weighted-least-squares method based on the Gauss-Newton minimization technique is used for the estimation of these parameters. The validation and application of the inverse model is illustrated with an example pipe distribution system under steady state. A generalized procedure to handle noisy and bad (abnormal) data is suggested, which can be used to estimate these parameters more accurately. The developed inverse model is useful for water supply agencies to calibrate their water distribution system and to improve their operational strategies to maintain water quality.

Transient electric field distribution within the working volume of an NEMP simulator

Nuclear electro-magnetic pulse (NEMP) simulators which are used in the simulation of transient electromagnetic fields due to a high altitude nuclear detonation are generally excited with a double exponential high voltage pulse. This results in a current distribution on the wires of the simulator and hence a transient electric field in the working volume of the simulator where the test object is kept. It is found that for the simulator under study, the current distribution is non-uniform and so is the field distribution along the width of the simulator in the working volume. To make the current distribution uniform, several methods have been suggested and the results of these methods are analyzed and suitable conclusions are arrived at from those results.

Molecular Dynamics of Thermotropic Liquid Crystals: Anomalous relaxation dynamics of calamitic and discotic liquid crystals

Recent optical kerr effect (OKE) studies have demonstrated that orientational relaxation of rod-like nematogens exhibits temporal power law decay at intermediate times not only near the isotropic–nematic (I–N) phase boundary but also in the nematic phase. Such behaviour has drawn an intriguing analogy with supercooled liquids. We have investigated both collective and single-particle orientational dynamics of a family of model system of thermotropic liquid crystals using extensive computer simulations. Several remarkable features of glassy dynamics are on display including non-exponential relaxation, dynamical heterogeneity, and non-Arrhenius temperature dependence of the orientational relaxation time. Over a temperature range near the I–N phase boundary, the system behaves remarkably like a fragile glass-forming liquid. Using proper scaling, we construct the usual relaxation time versus inverse temperature plot and explicitly demonstrate that one can successfully define a density dependent fragility of liquid crystals. The fragility of liquid crystals shows a temperature and density dependence which is remarkably similar to the fragility of glass forming supercooled liquids. Energy landscape analysis of inherent structures shows that the breakdown of the Arrhenius temperature dependence of relaxation rate occurs at a temperature that marks the onset of the growth of the depth of the potential energy minima explored by the system. A model liquid crystal, consisting of disk-like molecules, has also been investigated in molecular dynamics simulations for orientational relaxation along two isobars starting from the high temperature isotropic phase. The isobars have been so chosen that the phase sequence isotropic (I)–nematic (N)–columnar (C) appears upon cooling along one of them and the sequence isotropic (I)–columnar(C) along the other. While the orientational relaxation in the isotropic phase near the I–N phase transition shows a power law decay at short to intermediate times, such power law relaxation is not observed in the isotropic phase near the I–C phase boundary. The origin of the power law decay in the single-particle second-rank orientational time correlation function (OTCF) is traced to the growth of the orientational pair distribution functions near the I–N phase boundary. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational OTCF follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.

Statistics of the inverse-cascade regime in two-dimensional magnetohydrodynamic turbulence

We present a detailed direct numerical simulation of statistically steady, homogeneous, isotropic, two-dimensional magnetohydrodynamic turbulence. Our study concentrates on the inverse cascade of the magnetic vector potential. We examine the dependence of the statistical properties of such turbulence on dissipation and friction coefficients. We extend earlier work significantly by calculating fluid and magnetic spectra, probability distribution functions (PDFs) of the velocity, magnetic, vorticity, current, stream-function, and magnetic-vector-potential fields, and their increments. We quantify the deviations of these PDFs from Gaussian ones by computing their flatnesses and hyperflatnesses. We also present PDFs of the Okubo-Weiss parameter, which distinguishes between vortical and extensional flow regions, and its magnetic analog. We show that the hyperflatnesses of PDFs of the increments of the stream function and the magnetic vector potential exhibit significant scale dependence and we examine the implication of this for the multiscaling of structure functions. We compare our results with those of earlier studies.

Homogeneous isotropic superfluid turbulence in two dimensions: Inverse and forward cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the simplified, Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normalfluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli.