The BGK Model with External Confining Potential: Existence,Long-Time Behaviour and Time-Periodic Maxwellian Equilibria

We study global existence and long time behaviour for the inhomogeneous nonlinear BGK model for the Boltzmann equation with an external confining potential. For an initial datum f ₀≥0 with bounded mass, entropy and total energy we prove existence and strong convergence in L ¹ to a Maxwellian equilibrium state, by compactness arguments and multipliers techniques. Of particular interest is the case with an isotropic harmonic potential, in which Boltzmann himself found infinitely many time-periodic Maxwellian steady states. This behaviour is shared with the Boltzmann equation and other kinetic models. For all these systems we study the multistability of the time-periodic Maxwellians and provide necessary conditions on f ₀ to identify the equilibrium state, both in L ¹ and in Lyapunov sense. Under further assumptions on f, these conditions become also sufficient for the identification of the equilibrium in L ¹.     

Some Stability and Instability Criteria for Ideal Plane Flows

We investigate stability and instability of steady ideal plane flows for an arbitrary bounded domain. First, we obtain some general criteria for linear and nonlinear stability. Second, we find a sufficient condition for the existence of a growing mode to the linearized equation. Third, we construct a steady flow which is nonlinearly and linearly stable in the L² norm of vorticity but linearly unstable in the L² norm of velocity.     

Locating the minimum: Approach to equilibrium in a disordered, symmetric zero range process

We consider the dynamics of the disordered, one-dimensional, symmetric zero range process in which a particle from an occupied site k hops to its nearest neighbor with a quenched rate w(k). These rates are chosen randomly from the probability distribution f(w) ∼ (w − c) ⁿ , where c is the lower cutoff. For n>0, this model is known to exhibit a phase transition in the steady state from a low density phase with a finite number of particles at each site to a high density aggregate phase in which the site with the lowest hopping rate supports an infinite number of particles. In the latter case, it is interesting to ask how the system locates the site with globally minimum rate. We use an argument based on the local equilibrium, supported by Monte Carlo simulations, to describe the approach to the steady state. We find that at large enough time, regions with a smooth density profile are described by a diffusion equation with site-dependent rates, while the isolated points where the mass distribution is singular act as the boundaries of these regions. Our argument implies that the relaxation time scales with the system size L as L ^z with z = 2 + 1/(n + 1) for n>1 and suggests a different behavior for n<1.     

Canonical Analysis of Condensation in Factorised Steady States

We study the phenomenon of real space condensation in the steady state of a class of mass transport models where the steady state factorises. The grand canonical ensemble may be used to derive the criterion for the occurrence of a condensation transition but does not shed light on the nature of the condensate. Here, within the canonical ensemble, we analyse the condensation transition and the structure of the condensate, determining the precise shape and the size of the condensate in the condensed phase. We find two distinct condensate regimes: one where the condensate is gaussian distributed and the particle number fluctuations scale normally as L ^1/2 where L is the system size, and a second regime where the particle number fluctuations become anomalously large and the condensate peak is non-gaussian. Our results are asymptotically exact and can also be interpreted within the framework of sums of random variables. We further analyse two additional cases: one where the condensation transition is somewhat different from the usual second order phase transition and one where there is no true condensation transition but instead a pseudocondensate appears at superextensive densities. PACS numbers: 05.40.-a, 02.50.Ey, 64.60.-i.     

Finite-Size Scaling in the Energy-Entropy Plane for the 2D ± Ising Spin Glass

For L × L square lattices with L ≤ 20 the 2D Ising spin glass with +1 and −1 bonds is found to have a strong correlation between the energy and the entropy of its ground states. A fit to the data gives the result that each additional broken bond in the ground state of a particular sample of random bonds increases the ground state degeneracy by approximately a factor of 10/3. For x=0.5 (where x is the fraction of negative bonds), over this range of L, the characteristic entropy defined by the energy-entropy correlation scales with size as L ^1.78(2). Anomalous scaling is not found for the characteristic energy, which essentially scales as L ². When x=0.25, a crossover to L ² scaling of the entropy is seen near L=12. The results found here suggest a natural mechanism for the unusual behavior of the low temperature specific heat of this model, and illustrate the dangers of extrapolating from small L. PACS numbers: 75.10.Nr, 75.40.Mg, 75.50.Lk     

Stability of 3He2 4He M and 3He3 4He M L = 0,1 Clusters

We have studied the stability of mixed ³He/⁴He clusters in L = 0 and L = 1 states by the diffusion Monte Carlo method, employing the Tang-Toennies-Yiu (TTY) He-He potential. The clusters ³He⁴He_M ( ) and ³He₂⁴He_M (L = 0, S = 0) are stable for M > 1, while to bind two ³He in a triplet state the minimum number of ⁴He is four. Considering clusters with three ³He, ³He₃⁴He₄ is the smallest stable system in the L = 1 state, while ³He₃⁴He₈ is the smallest stable system in the L = 0 state.     

Non-Equilibrium Steady States of Finite¶Quantum Systems Coupled to Thermal Reservoirs

We study the non-equilibrium statistical mechanics of a 2-level quantum system, ?, coupled to two independent free Fermi reservoirs ?₁, ?₂, which are in thermal equilibrium at inverse temperatures β₁≠β₂. We prove that, at small coupling, the combined quantum system ?+?₁+?₂ has a unique non-equilibrium steady state (NESS) and that the approach to this NESS is exponentially fast. We show that the entropy production of the coupled system is strictly positive and relate this entropy production to the heat fluxes through the system. A part of our argument is general and deals with spectral theory of NESS. In the abstract setting of algebraic quantum statistical mechanics we introduce the new concept of the C-Liouvillean, L, and relate the NESS to zero resonance eigenfunctions of L ^*. In the specific model ?+?₁+?₂ we study the resonances of L ^* using the complex deformation technique developed previously by the authors in [JP1]. Dedicated to Jean Michel Combes on the occasion of his sixtieth birthday Received: 12 July 2001 / Accepted: 11 October 2001     

Airy Distribution Function: From the Area Under a Brownian Excursion to the Maximal Height of Fluctuating Interfaces

The Airy distribution function describes the probability distribution of the area under a Brownian excursion over a unit interval. Surprisingly, this function has appeared in a number of seemingly unrelated problems, mostly in computer science and graph theory. In this paper, we show that this distribution function also appears in a rather well studied physical system, namely the fluctuating interfaces. We present an exact solution for the distribution P(h_m,L) of the maximal height h_m (measured with respect to the average spatial height) in the steady state of a fluctuating interface in a one dimensional system of size L with both periodic and free boundary conditions. For the periodic case, we show that P(h_m,L)=L^−1/2f(h_m L^−1/2) for all L>0 where the function f(x) is the Airy distribution function. This result is valid for both the Edwards–Wilkinson (EW) and the Kardar–Parisi–Zhang interfaces. For the free boundary case, the same scaling holds P(h_m,L)=L^−1/2F(h_m L^−1/2), but the scaling function F(x) is different from that of the periodic case. We compute this scaling function explicitly for the EW interface and call it the F-Airy distribution function. Numerical simulations are in excellent agreement with our analytical results. Our results provide a rather rare exactly solvable case for the distribution of extremum of a set of strongly correlated random variables. Some of these results were announced in a recent Letter [S.N. Majumdar and A. Comtet, Phys. Rev. Lett. 92: 225501 (2004)].     

延时对色关联噪声诱导的逻辑生长过程的影响

应用小延时近似方法,研究了色关联噪声诱导的延时逻辑生长过程,得到了肿瘤细胞数的稳态概率分布P_st(x)的近似解析表达式,发现延时τ的变化可以使P_st(x)发生由多极值结构向单极值结构的转换,延时τ还可以使随机系统的平均值〈x〉、二阶矩〈x²〉、归一化涨落V_ar的极值位置和极值大小发生改变. 关键词： 逻辑生长过程 延时 关联色噪声 统计性质     

Asymptotic Stability of the Relativistic Boltzmann Equation for the Soft Potentials

In this paper it is shown that unique solutions to the relativistic Boltzmann equation exist for all time and decay with any polynomial rate towards their steady state relativistic Maxwellian provided that the initial data starts out sufficiently close in L^￥_l{L^\infty_\ell}. If the initial data are continuous then so is the corresponding solution. We work in the case of a spatially periodic box. Conditions on the collision kernel are generic in the sense of Dudyński and Ekiel-Jeżewska (Commun Math Phys 115(4):607–629, 1985); this resolves the open question of global existence for the soft potentials.     

Test of the stability of a three-phase contact line with negative line tension

The equilibrium structure of a three-phase contact line with negative line tension using a mean-field free-energy functional is calculated in the square-gradient approximation. The equilibrium density profiles are found by solving the Euler—Lagrange equations on a square grid of N ² points covering an area of L ². The fluctuations about the equilibrium structure are analysed via the spectrum of the free energy's second functional derivative. The equilibrium configuration is found to be stable with respect to fluctuations in the structure of the three-phase line and of the interfaces that meet at this line. In addition, the behaviour is investigated numerically of the lowest eigenvalue as the area of the grid is increased. The lowest eigenvalue is always positive and vanishes as 1/L ².     

Binding energies of negatively charged excitons in quantum disks

The Hamiltonian of a negatively charged exciton X⁻ (trion) in a quantum disk with parabolic confinement has been diagonalized to obtain the binding eigenenergy values of the L1 states as a function of the electron-to-hole effective mass ratio and the disk radius. It is found that a negatively charged exciton X⁻ in a quantum disk may have the second bound state with orbital angular momentum L=1 and the triplet state of the two bound electrons.     

Is realization of Majorana neutrino and neutrinoless double beta decay possible in the framework of standard weak interactions?

Usually it is supposed that Majorana neutrino produced in the superposition state χ _L = ν _L + (ν _L ) ^c and then follows the neutrinoless double beta decay. But since the standard weak interactions are chiral invariant then neutrino at production has definite helicity (ν _L and (ν _L ) ^c have opposite spirality). Then these neutrinos are separately produced and their superposition state cannot appear. Thus we see that for unsuitable helicity the neutrinoless double β decay is not possible even if it is supposed that neutrino is a Majorana particle (i.e. there is not a lepton number which is conserved). Also transition of Majorana neutrino ν _L into antineutrino (ν _L ) ^c at their oscillations is forbidden since helicity in vacuum holds. Transition Majora neutrino ν _L into (ν _R ) ^c (i.e., ν _L → (ν _R ) ^c ) at oscillations is unobserved since it is supposed that mass of (ν _R ) ^c is very big. If neutrino is a Dirac particle there can be transition of ν _L neutrino into (sterile) antineutrino $ \bar v_R $ \bar v_R (i.e., ν _L → $ \bar v_R $ \bar v_R ) at neutrino oscillations if there takes place double violation of lepton number. It is necessary also to remark that introducing of a Majorana neutrino implies violation of global and local gauge invariance in the standard weak interactions.     

Rotation of a self-bound many-body system

We discuss the rotational excitations of highly quantum many-body systems (for example, polyatomic molecules or microdroplets of helium). For a general system, the states F?, where and ? is a rotationally invariant ground or vibrational state, are shown to be eigenfunctions of L ² and L_z , with eigenvalues L(L+1)? ² and L? (for even L). These wavefunctions preserve the translational invariance and the permutation and inversion symmetries of the N-particle state ?. For harmonic pair interactions, the f = 1 wavefunctions are shown to be exact eigenstates of the N-body hamiltonian. For large N, the states F?(f=1) represent surface oscillations of the type first proposed by Bohr. An inequality for the rotational excitation energy is obtained variationally; it depends on two, three, and four-particle correlations. Other translationally invariant angular momentum eigenfunctions are also discussed.     

Statistical dynamics at critical bifurcations in Duffing-van der Pol oscillator

We study the characteristic features of certain statistical quantities near critical bifurcations such as onset of chaos, sudden widening and band-merging of chaotic attractor and intermittency in a periodically driven Duffing-van der Pol oscillator. At the onset of chaos the variance of local expansion rate is found to exhibit a self-similar pattern. For all chaotic attractors the variance Σ_n(q) of fluctuations of coarse-grained local expansion rates of nearby orbits has a single peak. However, multiple peaks are found just before and just after the critical bifurcations. On the other hand, Σ_n (q) associated with the coarse-grained state variable is zero far from the bifurcations. The height of the peak of Σ_n(q) is found to increase as the control parameter approached the bifurcation point. It is maximum at the bifurcation point. Power-law variation of maximal Lyapunov exponent and the mean value of the state variablex is observed near sudden widening and intermittency bifurcations while linear variation is seen near band-merging bifurcation. The standard deviation of local Lyapunov exponent λ(X,L) and the local mean valuex(L) of the coordinatex calculated after everyL time steps are found to approach zero in the limitL → ∞ asL ^-Β. Β is sensitive to the values of control parameters. Further weak and strong chaos are characterized using the probability distribution of ak-step difference quantity δx_k = x_i+k x _i.     

Asymmetric Diffusion and the Energy Gap Above¶the 111 Ground State of the Quantum XXZ Model

We consider the anisotropic three dimensional XXZ Heisenberg ferromagnet in a cylinder with axis along the 111 direction and boundary conditions that induce ground states describing an interface orthogonal to the cylinder axis. Let L be the linear size of the basis of the cylinder. Because of the breaking of the continuous symmetry around the axis, the Goldstone theorem implies that the spectral gap above such ground states must tend to zero as L→∞. In [3] it was proved that, by perturbing in a sub-cylinder with basis of linear size R≪L the interface ground state, it is possible to construct excited states whose energy gap shrinks as R ^-2. Here we prove that, uniformly in the height of the cylinder and in the location of the interface, the energy gap above the interface ground state is bounded from above and below by const.L ^-2. We prove the result by first mapping the problem into an asymmetric simple exclusion process on ℤ³ and then by adapting to the latter the recursive analysis to estimate from below the spectral gap of the associated Markov generator developed in [7]. Along the way we improve some bounds on the equivalence of ensembles already discussed in [3] and we establish an upper bound on the density of states close to the bottom of the spectrum. Received: 9 August 2001 / Accepted: 29 October 2001     

Explicit construction of steady state of a model of chemical turbulence

Using only the microscopic dynamics, the nonequilibrium steady state of a one-dimensional cellular automaton (CA) model of chemical turbulence is explicitly constructed. A coding is found which decomposes the CA into three interacting shift systems, each of which has an independent steady-state distribution. It was previously shown that the steady state of this model is a Gibbs state. Hence the steady state can be represented in the formZ ^–1 e ^–F , whereF is the conditional energy of the system such that all conditional probabilities are continuous. It is shown that the conditional energy of this model has an approximate expression in terms of familiar models from equilibrium statistical mechanics.     

Effective nonlinear optical properties of shape distributed composite media

The effective linear and nonlinear optical properties of metal/dielectric composite media, in which ellipsoidal metal inclusions are distributed in shape, are investigated. The shape distribution function P(L _x, L _y) is assumed to be 2Δ^-2θ(L _x - 1/3 + Δ/3)θ(L _y - 1/3 + Δ/3)θ(2/3 + Δ/3 - L _x - L _y), where θ( ^{. . .} ) is the Heaviside function, Δ is the shape variance and L_i are the depolarization factors of the ellipsoidal inclusions along i-symmetric axes (i = x, y). Within the spectral representation, we adopt Maxwell-Garnett type approximation to study the effect of shape variance Δ on the effective nonlinear optical properties. Numerical results show that both the effective linear optical absorption α ∼ ωIm() and the modulus of the effective third-order optical nonlinearity enhancement |χ⁽³⁾ _e|/χ⁽³⁾ ₁ exhibit the nonmonotonic behavior with Δ. Moreover, with increasing Δ, the optical absorption and the nonlinearity enhancement bands become broad, accompanied with the decrease of their peaks. The adjustment of Δ from 0 to 1 allows us to examine the crossover behavior from no separation to large separation between optical absorption and nonlinearity enhancement peaks. As Δ → 0, i.e., the ellipsoidal shape deviates slightly from the spherical one, the dependence of |χ⁽³⁾ _e|/χ⁽³⁾ ₁ on Δ becomes strong first and then weak with increasing the imaginary part of inclusions' dielectric constant. In the dilute limit, the exact formula for the effective optical nonlinearity is derived, and the present approximation characterizes the exact results better than old mean field one does. Received 10 December 2002 Published online 4 June 2003 RID="a" ID="a"e-mail: lgaophys@pub.sz.jsinfo.net     

Stability of the Boltzmann equation with potential forces on torus

In this paper, we are concerned with the stability of solutions to the Cauchy problem of the Boltzmann equation with potential forces on torus. It is shown that the natural steady state with the symmetry of origin is asymptotically stable in the Sobolev space with exponential rate in time for any initially smooth, periodic, origin symmetric small perturbation, which preserves the same total mass, momentum and mechanical energy. For the non-symmetric steady state, it is also shown that it is stable in L¹-norm for any initial data with the finite total mass, mechanical energy and entropy.     

Energy landscapes, lowest gaps, and susceptibility of elastic manifolds at zero temperature

We study the effect of an external field on (1 + 1) and (2 + 1) dimensional elastic manifolds, at zero temperature and with random bond disorder. Due to the glassy energy landscape the configuration of a manifold changes often in abrupt, “first order”-type of large jumps when the field is applied. First the scaling behavior of the energy gap between the global energy minimum and the next lowest minimum of the manifold is considered, by employing exact ground state calculations and an extreme statistics argument. The scaling has a logarithmic prefactor originating from the number of the minima in the landscape, and reads ΔE ₁∼L ^θ[ln(L _z L ^{- ζ})]^-1/2, where ζ is the roughness exponent and θ is the energy fluctuation exponent of the manifold, L is the linear size of the manifold, and L_z is the system height. The gap scaling is extended to the case of a finite external field and yields for the susceptibility of the manifolds ∼L ^{2D + 1 - θ}[(1 - ζ)ln(L)]^1/2. We also present a mean field argument for the finite size scaling of the first jump field, h ₁∼L ^{d - θ}. The implications to wetting in random systems, to finite-temperature behavior and the relation to Kardar-Parisi-Zhang non-equilibrium surface growth are discussed. Received December 2000 and Received in final form April 2001