On the connection between the macroscopical and microscopical evolution in an exactly soluble hopping model: I. Neutral particles

The hopping rate equation for neutral particles, on an arbitrary periodical lattice, can be solved exactly. It is shown that if one scales the time t and the distances x(t→λ²t, x→λx) then, in the λ→∞ limit, the particle density tends to the solution of the diffusion equation faster than λ^?3. The diffusion coefficient is the same as obtained from both Kubo and Miller-Abrahams theory via the Einstein relation.     

Cellular automata approach to site percolation on ℤ2. A numerical study

We present a cellular automata model as a new approach to Bernoulli site percolation on the square lattice. A new macroscopic quantity is defined and numerically computed at each level step of the automata dynamics. Its limit manifests a critical behavior at a value of the site occupancy probability quite close to those obtained for site percolation on ² with the best-known numerical methods.     

The quantum harmonic oscillator on the sphere and the hyperbolic plane

A nonlinear model of the quantum harmonic oscillator on two-dimensional space of constant curvature is exactly solved. This model depends on a parameter λ that is related with the curvature of the space. First, the relation with other approaches is discussed and then the classical system is quantized by analyzing the symmetries of the metric (Killing vectors), obtaining a λ-dependent invariant measure dμ_λ and expressing the Hamiltonian as a function of the Noether momenta. In the second part, the quantum superintegrability of the Hamiltonian and the multiple separability of the Schrödinger equation is studied. Two λ-dependent Sturm-Liouville problems, related with two different λ-deformations of the Hermite equation, are obtained. This leads to the study of two λ-dependent families of orthogonal polynomials both related with the Hermite polynomials. Finally the wave functions Ψ_m,n and the energies E_m,n of the bound states are exactly obtained in both the sphere S² and the hyperbolic plane H².     

Mound morphology of the 2+1 -dimensional Wolf-Villain model caused by the step-edge diffusion effect

The mound morphology of the 2+1-dimensional Wolf-Villain model is studied by numerical simulation. The diffusion rule of this model has an intrinsic mechanism, i.e., the step-edge diffusion, to create a local uphill particle current, which leads to the formation of the mound. In the simulation, a noise reduction technique is employed to enhance the local uphill particle current. Our results for the dynamic exponent 1/z and the roughness exponent α obtained from the surface width show a dependence on the strength of the step-edge diffusion. On the other hand, λ(t), which describes the separation of the mounds, grows as a function of time in a power-law form in the regime where the coalescence of mounds occurs, λ(t)∼tⁿ, with n≈0.23-0.25 for a wide range of the deposition conditions under the step-edge diffusion effect. For m=1, a noise reduction factor of unity, the behavior of λ(t) in the saturated regime is also simulated. We find that the evolution behavior of λ(t) in the whole process can be described by the standard Family-Vicsek scaling.     

From quantum cellular automata to quantum lattice gases

A natural architecture for nanoscale quantum computation is that of a quantum cellular automaton. Motivated by this observation, we begin an investigation of exactly unitary cellular automata. After proving that there can be no nontrivial, homogeneous, local, unitary, scalar cellular automaton in one dimension, we weaken the homogeneity condition and show that there are nontrivial, exactly unitary, partitioning cellular automata. We find a one-parameter family of evolution rules which are best interpreted as those for a one-particle quantum automaton. This model is naturally reformulated as a two component cellular automaton which we demonstrate to limit to the Dirac equation. We describe two generalizations of this automaton, the second, of which, to multiple interacting particles, is the correct definition of a quantum lattice gas.     

十三点格子Boltzmann模型仿真

格子气和格子Boltzmann方法的迅速发展提供了一类求解流体力学问题的新方法。格子Boltzmann方法在保留了格子气模型优点的同时,克服了它的不足之处。本文讨论了一种三迭加HPP十三点模型,通过选择适当的平衡分布及参数,并用Chapman－Enskog展开和多尺度技术导出了Navier－Stokes方程。在微机上模拟了空腔流的流动问题,并与传统方法的计算结果进行了比较,结果表明该模型能较好的模拟复杂流动现象,并具有较好的工程应用背景。     

Baryon–Baryon and Meson–Baryon Potentials with QCD-Based Short-Range Cores

A unified potential model of baryon–baryon and meson–baryon interactions at low energies is proposed. In this model, the short-range cores which simulate recent lattice QCD calculations are introduced. In baryon–baryon sector, our potentials give a very good agreement with NN and YN scattering data and have very similar behavior to those from the lattice QCD calculations. The π N and KN interactions are also well described in our model.     

Macroscopic behaviour of a charged Boltzmann gas

We consider a classical charged gas (with self-consistent Coulomb interaction) described by a solvable linearized Boltzmann equation with thermalization on uniformly distributed scatterers. It is shown that if one scales the time t, the reciprocal space coordinate k and the Debye length l as λ²t, (1/λ)k, λl, respectively, in the λ → ∞ limit the charge density is equal to the solution of the corresponding diffusion-conduction (macroscopic) equation.     

Effects of mass media on opinion spreading in the Sznajd sociophysics model

In this work we consider the influence of mass media in the dynamics of the two-dimensional Sznajd model. This influence acts as an external field, and it is introduced in the model by means of a probability p of the agents to follow the media opinion. We performed Monte Carlo simulations on square lattices with different sizes, and our numerical results suggest a change on the critical behavior of the model, with the absence of the usual phase transition for p>∼0.18. Another effect of the probability p is to decrease the average relaxation times τ, that are log-normally distributed, as in the standard model. In addition, the τ values depend on the lattice size L in a power-law form, τ∼L^α, where the power-law exponent depends on the probability p.     

Similarity solutions of reaction–diffusion equation with space- and time-dependent diffusion and reaction terms

We consider solvability of the generalized reaction–diffusion equation with both space- and time-dependent diffusion and reaction terms by means of the similarity method. By introducing the similarity variable, the reaction–diffusion equation is reduced to an ordinary differential equation. Matching the resulting ordinary differential equation with known exactly solvable equations, one can obtain corresponding exactly solvable reaction–diffusion systems. Several representative examples of exactly solvable reaction–diffusion equations are presented.     

三维十五点格子Boltzmann模型仿真

格子气和格子Ｂｏｌｔｚｍａｎｎ方法的迅速发展提供了一类求解流体力学问题的新方法。格子Ｂｏｌｔｚｍａｎｎ方法在保留了格子气模型优点的同时,克服了它的不足之处。本文讨论了一种三维十五点格子Ｂｏｌｔｚｍａｎｎ模型,通过选择适当的平衡分布及参数,并用Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ展开和多尺度技术导出了Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程．在微机上模拟了工程中比较常见的管排绕流问题,并与实验观察到的结果进行了比较,结果表明该模型能较好的模拟复杂流动现象,并具有较好的工程应用背景。     

Diffusive Hydrodynamic Limits for Systems of Interacting Diffusions with Finite Range Random Interaction

We consider a system of interacting diffusive particles with finite range random interaction. The variables can be interpreted as charges at sites indexed by a periodic multidimensional lattice. The equilibrium states of the system are canonical Gibbs measures with finite range random interaction. Under the diffusive scaling of lattice spacing and time, we derive a deterministic nonlinear diffusion equation for the time evolution of the macroscopic charge density. This limit is almost sure with respect to the random environment. Received: 3 October 1996 / Accepted: 13 February 1997     

Pulses in a complex Ginzburg-Landau system: Persistence under coupling with slow diffusion

The Ginzburg-Landau (GL) equation is essential for understanding the dynamics of patterns in a wide variety of physical contexts. It governs the evolutions of small amplitude instabilities near criticality. If the instabilities are, however, driven by two coupled instability mechanisms, of which one corresponds with a neutrally stable mode, their evolution is described by a GL equation coupled to a diffusion equation.In this paper, we study the influence of an additional diffusion equation on the existence of pulse solutions in the complex GL equation. In light of recently developed insights into the effect of slow diffusion on the stability of pulses, we consider the case of slow diffusion, i.e., in which the additional diffusion equation acts on a long spatial scale.In previous work [A. Doelman, G. Hek, N. Valkhoff, Stabilization by slow diffusion in a real Ginzburg-Landau system, J. Nonlinear Sci. 14 (2004) 237-278; A. Doelman, G. Hek, N.J.M. Valkhoff, Algebraically decaying pulses in a Ginzburg-Landau system with a neutrally stable mode, Nonlinearity 20 (2007) 357-389], we restricted ourselves to a model with both real coefficients and, more importantly, a real amplitude A rather than the complex-valued A that is needed to completely describe the pattern formation near criticality. In this simpler setting, we proved that pulse solutions of the GL equation can both persist and be stabilized under coupling with a slow diffusion equation. In the current work, we no longer make these restrictions, so that the problem is higher-dimensional and intrinsically harder. By a combination of a geometrical approach and explicit perturbation analysis, we consider the persistence of the solitary pulse solution of the GL equation under coupling with the additional diffusion equation. In the two limiting situations of the nearly real GL equation and the near nonlinear Schrödinger equation, we show that the pulse solutions can indeed persist under this coupling.     

Diffusion in fluids with large mean free paths: Non-classical behavior between Knudsen and Fickian limits

Diffusion in fluids is analyzed at non-classical conditions, intermediate between the Knudsen and Fickian limits. The fluid is considered in the framework of the Einstein’s diffusion evolution equation involving expansions of the density distribution in powers of displacement and time. The standard truncation of these expansions results in the classical model of diffusion; however, higher-order terms lead to a departure from classical behavior. This has not been studied or discussed adequately in the literature previously.Here, we present an exact solution of the Einstein’s diffusion evolution equation without truncation of the density expansions. This solution illustrates limitations in the classical truncations and demonstrates non-classical effects due to large mean free paths, λ. In particular, this new solution shows that, at large λ, there are significant quantitative deviations from classical diffusion profiles. In addition, this solution demonstrates a dramatic change in the diffusion mechanism from the state where the molecular motions are predominantly ballistic to one of molecular chaos. This has implications for fundamentals of fluids between the Knudsen and Fickian limits, and for a variety of fields where evolution of a system includes random, multi-scale displacement of particles, such as nanotechnology, vacuum techniques, turbulence, and astrophysics.     

Simulation of euclidean quantum field theories by a random walk process

A new Monte Carlo method for euclidean lattice field theory is introduced by writing the Boltzmann distribution e^?s as a solution of a diffusion type equation and constructing the associated random walk process. It is practically tested for a quantum mechanical model and a non-compact version of lattice QCD. It is explained where the main interest in this algorithm lies: the diffusion process coming from an action that can be generalized to include non-conservative forces. This possibility is exploited in our QCD version to implement gauge fixing without Faddeev-Popov ghosts.     

Renormalization-group derivation of Navier-Stokes equation

The Navier-Stokes equation is proved from first principles (rotational symmetry and conservation of momentum, mass, and energy) using renormalization-group ideas. That is, we consider a system described by one (classical) conserved vector field and two conserved scalar fields, and demonstrate that on a large scale it obeys the Navier-Stokes equation. No assumptions about the physical meanings of the fields are required; in particular, no results from thermodynamics are used. The result comes about because the Euler equation is an exact fixed point of an appropriate scale-coarsening transformation, and the coefficients of the eigenvectors of the transformation with the largest (most relevant) eigenvalues include (in dimensiond>2) the thermal conductivity and the bulk and shear viscosities, leading to Navier-Stokes behavior on a large scale. Ford<2, the largest eigenvalue corresponds to a convection term, and the Navier-Stokes equation is incorrect. Our method differs from previous renormalization approaches in using time-coarsening as well as space-coarsening transformations. This allows renormalization trajectories to be determined exactly, and allows the determination of the macroscopic behavior of specific microscopic models. The Navier-Stokes equation we obtain is almost, but not exactly, the same as the conventional one; distinguishing between them experimentally would require measurement of the very small asymmetry of the Brillouin line in a simple fluid.     

Dependence of the cesium absorption D1 linewidth on the thickness of atomic vapor column

A new nanocell has been elaborated, where the thickness of the atomic vapor column varies smoothly in the range L = 350–5100 nm. The cell allows studying the behavior of the resonance absorption at the D₁ line of cesium atoms by varying the thickness from L = λ / 2 to L = 5 λ with the step λ / 2 (λ being the resonant wavelength of the laser, 894 nm) and the laser intensity. It is shown that at low laser intensities a narrowing of the resonance absorption spectrum is observed for L = (2n + 1)λ/2 (with an integer n) up to L = (7/2)λ, whereas for L = nλ the spectrum broadens. The developed theoretical model well describes the experiment.