Multiple separatrix crossing in multi-degree-of-freedom Hamiltonian flows

Summary We study separatrix crossing in near-integrablek-degree-of-freedom Hamiltonian flows, 2 <k < , whose unperturbed phase portraits contain separatrices inn degrees of freedom, 1 <n <k. Each of the unperturbed separatrices can be recast as a codimension-one separatrix in the 2k-dimensional phase space, and the collection of these separatrices takes on a variety of geometrical possibilities in the reduced representation of a Poincaré section on the energy surface. In general 0 l n of the separatrices will be available to the Poincaré section, and each separatrix may be completely isolated from all other separatrices or intersect transversely with one or more of the other available separatrices. For completely isolated separatrices, transitions across broken separatrices are described for each separatrix by the single-separatrix crossing theory of Wiggins, as modified by Beigie. For intersecting separatrices, a possible violation of a normal hyperbolicity condition complicates the analysis by preventing the use of a persistence and smoothness theory for compact normally hyperbolic invariant manifolds and their local stable and unstable manifolds. For certain classes of multi-degree-of-freedom flows, however, a local persistence and smoothness result is straightforward, and we study the global implications of such a local result. In particular, we find codimension-one partial barriers and turnstile boundaries associated with each partially destroyed separatrix. From the collection of partial barriers and turnstiles follows a rich phase space partitioning and transport formalism to describe the dynamics amongst the various degrees of freedom. A generalization of Wiggins' higher-dimensional Melnikov theory to codimension-one surfaces in the multi-separatrix case allows one to uncover invariant manifold geometry. In the context of this perturbative analysis and detailed numerical computations, we study invariant manifold geometry, phase space partitioning, and phase space transport, with particular attention payed to the role of a vanishing frequency in the limit approaching the intersection of the partially destroyed separatrices. The class of flows under consideration includes flows of basic physical relevance, such as those describing scattering phenomena. The analysis is illustrated in the context of a detailed study of a 3-degree-of-freedom scattering problem.     

Asymptotic behavior of basic contact process with rapid stirring

In this paper we consider the basic contact process with infection rate λ and stirring rateD. We study the asymptotic behavior of the critical value and survival probability asD→∞.     

New results for diffusion in Lorentz lattice gas cellular automata

New calculations to over ten million time steps have revealed a more complex diffusive behavior than previously reported of a point particle on a square and triangular lattice randomly occupied by mirror or rotator scatterers. For the square lattice fully occupied by mirrors where extended closed particle orbits occur, anomalous diffusion was still found. However, for a not fully occupied lattice the superdiffusion, first noticed by Owczarek and Prellberg for a particular concentration, obtains for all concentrations. For the square lattice occupied by rotators and the triangular lattice occupied by mirrors or rotators, an absence of diffusion (trapping) was found for all concentrations, except on critical lines, where anomalous diffusion (extended closed orbits) occurs and hyperscaling holds for all closed orbits withuniversal exponentsd _f =7/4 and =15/7. Only one point on these critical lines can be related to a corresponding percolation problem. The questions arise therefore whether the other critical points can be mapped onto a new percolation-like problem and of the dynamical significance of hyperscaling.     

Critical exponents,hyperscaling, and universal amplitude ratios for two- and three-dimensional self-avoiding walks

We make a high-precision Monte Carlo study of two- and three-dimensional self-avoiding walks (SAWs) of length up to 80,000 steps, using the pivot algorithm and the Karp-Luby algorithm. We study the critical exponentsv and 2 ₄ – as well as several universal amplitude ratios; in particular, we make an extremely sensitive test of the hyperscaling relationdv = 2 ₄ –. In two dimensions, we confirm the predicted exponentv=3/4 and the hyperscaling relation; we estimate the universal ratios <R _g ² >/<R _e ² >=0.14026±0.00007, <R _m ² >/<R _e ² >=0.43961±0.00034, and ^*=0.66296±0.00043 (68% confidence limits). In three dimensions, we estimatev=0.5877±0.0006 with a correctionto-scaling exponent ₁=0.56±0.03 (subjective 68% confidence limits). This value forv agrees excellently with the field-theoretic renormalization-group prediction, but there is some discrepancy for ₁. Earlier Monte Carlo estimates ofv, which were 0.592, are now seen to be biased by corrections to scaling. We estimate the universal ratios <R _g ² >/<R _e ² >=0.1599±0.0002 and ^*=0.2471±0.0003; since ^*>0, hyperscaling holds. The approach to ^* is from above, contrary to the prediction of the two-parameter renormalization-group theory. We critically reexamine this theory, and explain where the error lies. In an appendix, we prove rigorously (modulo some standard scaling assumptions) the hyperscaling relationdv = 2 ₄ – for two-dimensional SAWs.     

Superexponential stability of KAM tori

We study the dynamics in the neighborhood of an invariant torus of a nearly integrable system. We provide an upper bound to the diffusion speed, which turns out to be of superexponentially small size exp[-exp(1/)], being the distance from the invariant torus. We also discuss the connection of this result with the existence of many invariant tori close to the considered one.     

Phase separation of symmetrical polymer mixtures in thin-film geometry

Monte Carlo simulations of the bond fluctuation model of symmetrical polymer blends confined between two neutral repulsive walls are presented for chain lengthN _A=N _B=32 and a wide range of film thicknessD (fromD=8 toD=48 in units of the lattice spacing). The critical temperaturesT _c (D) of unmixing are located by finite-size scaling methods, and it is shown that , wherev ₃0.63 is the correlation length exponent of the three-dimensional Ising model universality class. Contrary to this result, it is argued that the critical behavior of the films is ruled by two-dimensional exponents, e.g., the coexistence curve (difference in volume fraction of A-rich and A-poor phases) scales as , where ₂ is the critical exponent of the two-dimensional Ising universality class ( ₂=1/8). Since for largeD this asymptotic critical behavior is confined to an extremely narrow vicinity ofT _c (D), one observes in practice effective exponents which gradually cross over from ₂ to ₃ with increasing film thickness. This anomalous flattening of the coexistence curve should be observable experimentally.     

Multifractal analysis of Brownian zero set

We present a method for the derivation of the generating function and computation of critical exponents for several cluster models (staircase, bar-graph, and directed column-convex polygons, as well as partially directed self-avoiding walks), starting with nonlinear functional equations for the generating function. By linearizing these equations, we first give a derivation of the generating functions. The nonlinear equations are further used to compute the thermodynamic critical exponents via a formal perturbation ansatz. Alternatively, taking the continuum limit leads to nonlinear differential equations, from which one can extract the scaling function. We find that all the above models are in the same universality class with exponents _u =-1/2, _i =-1/3, and =2/3. All models have as their scaling function the logarithmic derivative of the Airy function.     

Lattice-BGK approach to simulating granular flows

Many continuum theories for granular flow produce an equation of motion for the fluctuating kinetic energy density (granular temperature) that accounts for the energy lost in inelastic collisions. Apart from the presence of an extra dissipative term, this equation is very similar in form to the usual temperature equation in hydrodynamics. It is shown how a lattice-kinetic model based on the Bhatnagar-Gross-Krook (BGK) equation that was previously derived for a miscible two-component fluid may be modified to model the continuum equations for granular flow. This is done by noting that the variable corresponding to the concentration of one species follows an equation that is essentially analogous to the granular temperature equation. A simulation of an unforced granular fluid using the modified model reproduces the phenomenon of clustering instability, namely the spontaneous agglomeration of particles into dense clusters, which occurs generically in all granular flows. The success of the continuum theory in capturing the gross features of this basic phenomenon is discussed. Some shear flow simulations are also presented.     

An experiment to test an explanation of a possible violation of quantum theory

An experiment to test a possible explanation of the Schmidt backwards causation results is suggested. The experiment might distinguish between many- and one- world interpretations of quantum theory.     

Necessary conditions for characterization of laws via mixed sums

SupposeX andY are independent and identically distributed, and independent ofU which satisfies 0U1. Recent work has centered on finding the lawsL(X) for whichX U(X+Y) where denotes equality in law. We show that this equation corresponds to a certain projective invariance property under random rotations. Implicitly or explicitly, it has been assumed that the characteristic function ofX has an expansion property near the origin. We show that solutions may be admitted in the absence of this condition when –logU has a lattice law. A continuous version of the basic problem replaces sums with a Lévy process. Instead we consider self-similar processes, showing that a solution exists only whenU is constant, and then all processes of a given order are admitted.This research was in part supported by NSERC grant A-8466.