Fluctuations of Power Injection in Randomly Driven Granular Gases

We investigate the large deviation function π_∞(w) for the fluctuations of the power W(t) = wt, integrated over a time t, injected by a homogeneous random driving into a granular gas, in the infinite time limit. Our analytical study starts from a generalized Liouville equation and exploits a Molecular Chaos-like assumption. We obtain an equation for the generating function of the cumulants μ(λ) which appears as a generalization of the inelastic Boltzmann equation and has a clear physical interpretation. Reasonable assumptions are used to obtain μ(λ) in a closed analytical form. A Legendre transform is sufficient to get the large deviation function π_∞(w). Our main result, apart from an estimate of all the cumulants of W(t) at large times t, is that π_∞ has no negative branch. This immediately results in the inapplicability of the Gallavotti-Cohen Fluctuation Relation (GCFR), that in previous studies had been suggested to be valid for injected power in driven granular gases. We also present numerical results, in order to discuss the finite time behavior of the fluctuations of W (t) . We discover that their probability density function converges extremely slowly to its asymptotic scaling form: the third cumulant saturates after a characteristic time τ larger than ∼50 mean free times and the higher order cumulants evolve even slower. The asymptotic value is in good agreement with our theory. Remarkably, a numerical check of the GCFR is feasible only at small times (at most τ/10), since negative events disappear at larger times. At such small times this check leads to the misleading conclusion that GCFR is satisfied for π_∞(w). We offer an explanation for this remarkable apparent verification. In the inelastic Maxwell model, where a better statistics can be achieved, we are able to numerically observe the “failure” of GCFR.     

Almost-Sure Central Limit Theorem for Directed Polymers and Random Corrections

We consider a general model of directed polymers on the lattice , weakly coupled to a random environment. We prove that the central limit theorem holds almost surely for the discrete time random walk X _T associated to the polymer. Moreover we show that the random corrections to the cumulants of X _T are finite, starting from some dimension depending on the index of the cumulants, and that there are corresponding random corrections of order , , in the asymptotic expansion of the expectations of smooth functions of X _T . Full proofs are carried out for the first two cumulants. We finally prove a kind of local theorem showing that the ratio of the probabilities of the events to the corresponding probabilities with no randomness, in the region of “moderate” deviations from the average drift bT, are, for almost all choices of the environment, uniformly close, as , to a functional of the environment “as seen from (T,y)$”. Received: 14 October 1996 / Accepted: 28 March 1997     

First order phase transitions in the canonical and the microcanonical ensemble

In the canonical ensemble any singularity of a thermodynamic function at a temperatureT _c is smeared over a temperature range of orderT _T /N. Therefore it is rather difficult to distinguish between a discontinuous and a continuous phase transition on the basis of numerical data obtained for finite systems in the canonical ensemble. It is demonstrated for four model systems that this problem cannot be circumvented by considering higher cumulants of the energy distribution or cumulant ratios. On the other hand, the distinction between first and a second order phase transition is rather direct if based on the microcanonical density of states which is readily obtainable in the dynamical ensemble.     

Theorem on the Distribution of Short-Time Particle Displacements with Physical Applications

The distribution of the initial short-time displacements of particles is considered for a class of classical systems under rather general conditions on the dynamics and with Gaussian initial velocity distributions, while the positions could have an arbitrary distribution. This class of systems contains canonical equilibrium of a Hamiltonian system as a special case. We prove that for this class of systems the nth order cumulants of the initial short-time displacements behave as the 2n-th power of time for all n > 2, rather than exhibiting an nth power scaling. This has direct applications to the initial short-time behavior of the Van Hove self-correlation function, to its non-equilibrium generalizations the Green's functions for mass transport, and to the non-Gaussian parameters used in supercooled liquids and glasses. PACS Number: 05.20.-y, 02.30.Mv, 66.10.-x, 78.70.Nx, 05.60.Cd     

Random Matrix Theory and ζ(1/2+it)

We study the characteristic polynomials Z(U, θ) of matrices U in the Circular Unitary Ensemble (CUE) of Random Matrix Theory. Exact expressions for any matrix size N are derived for the moments of |Z| and Z/Z ^*, and from these we obtain the asymptotics of the value distributions and cumulants of the real and imaginary parts of log Z as N→∞. In the limit, we show that these two distributions are independent and Gaussian. Costin and Lebowitz [15] previously found the Gaussian limit distribution for Im log Z using a different approach, and our result for the cumulants proves a conjecture made by them in this case. We also calculate the leading order N→∞ asymptotics of the moments of |Z| and Z/Z ^*. These CUE results are then compared with what is known about the Riemann zeta function ζ (s) on its critical line Re s= 1/2, assuming the Riemann hypothesis. Equating the mean density of the non-trivial zeros of the zeta function at a height T up the critical line with the mean density of the matrix eigenvalues gives a connection between N and T. Invoking this connection, our CUE results coincide with a theorem of Selberg for the value distribution of log ζ(1/2+iT) in the limit T→∞. They are also in close agreement with numerical data computed by Odlyzko [29] for large but finite T. This leads us to a conjecture for the moments of |ζ(1/2+it) |. Finally, we generalize our random matrix results to the Circular Orthogonal (COE) and Circular Symplectic (CSE) Ensembles. Received: 20 December 1999 / Accepted: 24 March 2000     

Monte Carlo investigation of phase transitions in the frustrated Heisenberg model on a triangular lattice

The critical properties and phase transitions of the three-dimensional frustrated antiferromagnetic Heisenberg model on a triangular lattice have been investigated using the Monte Carlo method with a replica algorithm. The critical temperature has been determined and the character of the phase transitions has been analyzed using the method of fourth-order Binder cumulants. A second-order phase transition has been found in the three-dimensional frustrated Heisenberg model on a triangular lattice. The static magnetic and chiral critical exponents of the heat capacity α, the susceptibility γ and γ _k , the magnetization β and β _k , the correlation length ν and ν _k , as well as the Fisher exponents η and η _k , have been calculated in terms of the finite-size scaling theory. It has been demonstrated that the three-dimensional frustrated antiferromagnetic Heisenberg model on a triangular lattice forms a new universality class of the critical behavior.     

Current Fluctuations in the One-Dimensional Symmetric Exclusion Process with Open Boundaries

We calculate the first four cumulants of the integrated current of the one-dimensional symmetric simple exclusion process of N sites with open boundary conditions. For large system size N, the generating function of the integrated current depends on the densities ρ _a and ρ _b of the two reservoirs and on the fugacity z, the parameter conjugated to the integrated current, through a single parameter. Based on our expressions for these first four cumulants, we make a conjecture which leads to a prediction for all the higher cumulants. In the case ρ _a =1 and ρ _b =0, our conjecture gives the same universal distribution as the one obtained by Lee, Levitov, and Yakovets for one-dimensional quantum conductors in the metallic regime.     

Quasiclassical Hamiltonians for large-spin systems

We extend and apply a previously developed method for a semiclassical treatment of a system with large spin S. A multisite Heisenberg Hamiltonian is transformed into an effective classical Hamilton function which can be treated by standard methods for classical systems. Quantum effects enter in form of multispin interactions in the Hamilton function. The latter is written in the form of an expansion in powers of J/(TS), where J is the coupling constant. Main ingredients of our method are spin coherent states and cumulants. Rules and diagrams are derived for computing cumulants of groups of operators entering the Hamiltonian. The theory is illustrated by calculating the quantum corrections to the free energy of a Heisenberg chain which were previously computed by a Wigner-Kirkwood expansion. Received 5 May 1999 and received in final form 24 September 1999     

Heat conduction in a random medium

The method of time-ordered cumulants is used to investigate the behavior of heat pulses in a one-dimensional medium in which the thermal conductivity is random. A partial differential equation is obtained for the average temperature profile; it is the heat equation modified by the addition of a fourth-order spatial derivative. A solution is obtained by asymptotic series. The first two spatial moments of the average temperature profile are evaluated and are shown to tend to those of a Gaussian whent is large. Finally, an equation is obtained for the covariance function.Alfred P. Sloan Fellow.     

Multipeak fitting analysis of Raman spectra on DLCH film

The hydrogenated diamond‐like carbon (DLCH) film with 1‐µm thickness is deposited by direct hydrocarbon gas ion beam method on silicon wafer and annealed at 400 °C. Detailed Raman spectra feature are fitted from nine sets of different peak fitting functions, including Gaussian, Lorentzian and Breit‐Wigner‐Fano (BWF) functions. These fitting results obtained from a two‐peak combination show some specific variances on the G peak position, FWHM_G and I_D/I_G ratio for as‐deposited and as‐annealed DLCH films. The most popular two‐peak fitting method with full Gaussian function tends to exhibit a higher ratio of the G peak position shift and higher I_D/I_G ratio than others fitting methods, the drastic difference among the most popular G (G) & G (D) and B (G) & L (D) schemes also have brought out in I_D/I_G ratio. However, for a more complex four‐peak Gaussian function fitting Raman spectra, the I_D/I_G ratio is close to that of a two‐peak fitting function with a mixture functions of BWF (G) and L (D). Furthermore, a series of systematic peak fitting procedures and comparisons of Raman spectra have been discussed in this study. Copyright © 2009 John Wiley & Sons, Ltd.     

Generalized beam propagation factor of hard-edged circular apertured diffracted Bessel–Gaussian beams

Based on the truncated second-order moments method on the cylindrical coordinate systems and the incomplete gamma function, an analytical expression of the generalized beam propagation factor (M_G² factor) of hard-edged circular apertured diffracted Bessel–Gaussian beams is derived and illustrated numerically. It is shown that the M_G² factor of hard-edged diffracted BGBs mainly depends on the truncation parameter δ and the beam parameters m and η. The results can be reduced to that for the non-truncated Bessel–Gaussian beams case and that for the truncated fundamental Gaussian beams case under certain conditions, respectively. The power fraction is also discussed analytically and numerically.     

Synthesis and Mössbauer Study of Maghemite Nanowire Arrays

An analytical formula for the distance dependence of the electric field gradient produced by a Gaussian charge density distribution n(r) is derived. This charge density is displaced by z ₀ along the z-axis. The system has cylindrical symmetry; hence it suffices to calculate V _zz(0). It turns out that V _zz(0) is always smaller than the value with the total charge shrunk into a point. For distances larger than about four times the Gaussian width σ the expression approaches the point charge value. For z₀ → 0, i.e. a spherically symmetric charge distribution around the origin, V _zz(0) vanishes quadratically, as required by symmetry. A slab-wise calculation in cylindrical coordinates is presented which shows the contribution to V _zz(0) for infinitesimally thin slabs as a function of distance from the origin. This analytical formula allows for a fast computation of electric field gradients from a given charge density distribution for Gaussian expansions of Slater-type orbitals. An example for a hydrogen atom will be given.     

On the Distribution of Surface Extrema in Several One- and Two-dimensional Random Landscapes

We study here a standard next-nearest-neighbor (NNN) model of ballistic growth on one-and two-dimensional substrates focusing our analysis on the probability distribution function P(M,L) of the number M of maximal points (i.e., local “peaks”) of growing surfaces. Our analysis is based on two central results: (i) the proof (presented here) of the fact that uniform one-dimensional ballistic growth process in the steady state can be mapped onto “rise-and-descent” sequences in the ensemble of random permutation matrices; and (ii) the fact, established in Ref. [G. Oshanin and R. Voituriez, J. Phys. A: Math. Gen. 37:6221 (2004)], that different characteristics of “rise-and-descent” patterns in random permutations can be interpreted in terms of a certain continuous-space Hammersley-type process. For one-dimensional system we compute P(M,L) exactly and also present explicit results for the correlation function characterizing the enveloping surface. For surfaces grown on 2d substrates, we pursue similar approach considering the ensemble of permutation matrices with long-ranged correlations. Determining exactly the first three cumulants of the corresponding distribution function, we define it in the scaling limit using an expansion in the Edgeworth series, and show that it converges to a Gaussian function as L → ∞.     

Spectra of Hamiltonians with generalized single-site dynamical disorder

Starting from the deformed commutation relationsa _q (t) a ^q (s)–q a ^q (s) a _q (t)=(t–s)1, –1q1 with a covariance (t–s) and a parameterq varying between –1 and 1, a stochastic process is constructed which continuously deforms the classical Gaussian and classical compound Poisson process. The moments of these distinguished stochastic processes are identified with the Hilbert space vacuum expectation values of products of with fixed parametersq, and . Thereby we can interpolate between dichotomic, random matrix and classical Gaussian and compound Poisson processes. The spectra of Hamiltonians with single-site dynamical disorder are calculated for an exponential covariance (coloured noise) by means of the time convolution generalized master equation formalism (TC-GME) and the partial cumulants technique. The final result for the spectral function is given as aq-dependent infinite continued fraction. In the case of the random matrix processes the infinite continued fraction can be summed up yielding a self-consistent equation for the one-particle Green function.     

Formal structure of kinetic theory

We study the Liouville equation in the domain of small deviations from absolute equilibrium. The solution is expressed in terms of amplitudes ofn-body additive functions which are orthogonal with respect to the Gibbs weight factor. In the memory operator approach the memory operators are formally exact continued fractions. We show that with the isolation in the Liouville operator of a one-body additive operatorL _o, any memory operator can be written alternatively as an exact infinite series, each term of which can be calculated exactly. This yields improvements of the dressed particle approximation. We discuss the choice ofL _o, which is in general time dependent. The theory is developed both for smooth potentials and for hard spheres, where we use pseudo-Liouville operators. The theory can be formulated in an equivalent manner by introducing modified cumulant distributions, which are closely related to the amplitudes. The modified cumulants have the same spatial asymptotic properties as ordinary cumulants, but have superior short-time and small-distance behavior.Work supported by the National Science Foundation.     

Self‐focusing of a cosh‐Gaussian laser beam in magnetized plasma under relativistic‐ponderomotive regime

The combined effect of relativistic and ponderomotive nonlinearities on the self‐focusing of an intense cosh‐Gaussian laser beam (CGLB) in magnetized plasma have been investigated. Higher‐order paraxial‐ray approximation has been used to set up the self‐focusing equations, where higher‐order terms in the expansion of the dielectric function and the eikonal are taken into account. The effects of various lasers and plasma parameters viz. laser intensity (a₀), decentred parameter (b), and magnetic field (ω_c) on the self‐focusing of CGLB have been explored. The results are compared with the Gaussian profile of laser beams and relativistic nonlinearity. Self‐focusing can be enhanced by optimizing and selecting the appropriate laser‐plasma parameters. It is observed that the focusing of CGLB is fast in a nonparaxial region in comparison with that of a Gaussian laser beam and in a paraxial region in magnetized plasma. In addition, strong self‐focusing of CGLB is observed at higher values of a₀, b, and ω_c. Numerical results show that CGLB can produce ultrahigh laser irradiance over distances much greater than the Rayleigh length, which can be used for various applications.     

Invariant mass dependence of particle correlations in π+ p andK + p interactions at 250 GeV/c

We study two- and three-particle correlations as a function of invariant mass. Using data on π⁺ p andK ⁺ p collisions at 250 GeV/c, we compare correlation functions and normalised factorial cumulants for various charge combinations. Strong positive correlations are observed only at small invariant masses. The normalised cumulants for “exotic” [(??), (++)] and “nonexotic” pairs (+?) and triplets decrease in power-like fashion with increasing invariant mass. The mass dependence is not incompatible with the power-law behaviour as expected in a Dual Mueller-Regge framework. Comparison with FRITIOF reveals strong disagreements, which are due to too large production rates of resonances, such as ρ⁰ and η, and the absence of a Bose-Einstein low-mass enhancement in JETSET.

     

Factorial cumulant moments and correlation integrals in pp collisions at 400 GeV/c

The pseudorapidity distribution of charged particles produced in pp collisions at 400 GeV/c was measured by using LEBC films offered by CERN NA27 collaboration. The scaled factorial cumulant moments have been calculated. The results show that the second order cumulants have positive values, while the cumulants of higher order are consistent with zero except for the situation ofn _ch ≥ 4 events, where the third-order cumulants have positive values beyond the statistical uncertainties. It means that the observed increase of the higherorder factorial momentsF _q is almost due to the dynamical two-particle correlation in pp collisions at 400 GeV/c. The cumulant moments also have been calculated by star correlation integrals. It significantly reduced the statistical errors, especially for higher order cumulants. From Monte Carlo events with the same single particle spectrum and no correlations, we observed that for broad mixed-multiplicity distributions, a significant part ofK ₂ is coming from the single-particle fluctuation due to the fluctuating multiplicity.     

Beam wander of J 0- and I 0-Bessel Gaussian beams propagating in turbulent atmosphere

Root mean square (rms) beam wander of J ₀-Bessel Gaussian and I ₀-Bessel Gaussian beams, normalized by the rms beam wander of the fundamental Gaussian beam, is evaluated in atmospheric turbulence. Our formulation is based on the first and the second statistical moments obtained from the Rytov series. It is found that after propagating in atmospheric turbulence, the collimated J ₀-Bessel Gaussian and the I ₀-Bessel Gaussian beams have smaller rms beam wander than that of the Gaussian beam, regardless of the choice of Bessel width parameter. However, the extent of such an advantage depends on the chosen width parameter, Gaussian source size, propagation distance and the wavelength. Focusing at finite distances of the considered beams causes the rms beam wander to decrease sharply at the propagation distances equal to the focusing parameter.