Trapping,percolation, and anomalous diffusion of particles in a two-dimensional random field

We analyze from first principles the advection of particles in a velocity field with HamiltonianH(x, y)=¯ V ₁ y–¯ V ₂ x+W ₁ (y)-W ₂ (x), whereW _i , i=1, 2, are random functions with stationary, independent increments. In the absence of molecular diffusion, the particle dynamics are very sensitive to the streamline topology, which depends on the mean-to-fluctuations ratio=max(|¯V₁¦/; ¦¯V₂|/), with =|W ₁ |^21/2=rms fluctuations. Remarkably, the model is exactly solvable for 1 and well suited for Monte Carlo simulations for all , providing a nice setting for studying seminumerically the influence of streamline topology on large-scale transport. First, we consider the statistics of streamlines for=0, deriving power laws for p_nc(L) and (L), which are, respectively, the escape probability and the length of escaping trajectories for a box of sizeL, L » 1. We also obtain a characterization of the statistical topography of the HamiltonianH. Second, we study the large-scale transport of advected particles with > 0. For 0 < < 1, a fraction of particles is trapped in closed field lines and another fraction undergoes unbounded motions; while for 1 all particles evolve in open streamlines. The fluctuations of the free particle positions about their mean is studied in terms of the normalized variablest ^– v/2[x(t)–x(t)] andt ^–v/2 [y(t)-(t)]. The large-scale motions are shown to be either Fickian (=1), or superdiffusive (=3/2) with a non-Gaussian coarse-grained probability, according to the direction of the mean velocity relative to the underlying lattice. These results are obtained analytically for 1 and extended to the regime 0<<1 by Monte Carlo simulations. Moreover, we show that the effective diffusivity blows up for resonant values of ) for which stagnation regions in the flow exist. We compare the results with existing predictions on the topology of streamlines based on percolation theory, as well as with mean-field calculations of effective diffusivities. The simulations are carried out with a CM 200 massively parallel computer with 8192 SIMD processors.     

On log-Sobolev inequalities for infinite lattice systems

For a system on an infinite lattice, we show that a Gibbs measure for a smooth local specification ={E } satisfying the Dobrushin uniqueness theorem also satisfies log-Sobolev inequality, provided it is satisfied for one-dimensional measures E _l .     

Radius,perimeter, and density profile for percolation clusters and lattice animals

Monte Carlo simulation and series expansion shows the radius of gyration of large clusters withs sites each to vary ass with0.56 in two and0.47 in three dimensions at the percolation threshold, and with(d=2)0.65 and(d=3)0.53 for random lattice animals (zero concentration). Clusters up tos=100 were used. The perimeter of random animals approaches 2.8s for larges on the simple cubic lattice. Monte Carlo simulation of the Eden process (growing animals) up tos=5,000 indicates a systematic variation of about ±0.05 for the effective exponent=(s) and thus suggests that the true asymptotic exponents may be compatible with the predictions of hyper-scaling.     

Measurement of the refractive index and optical absorption spectra of epitaxial bismuth substituted yttrium iron garnet films at uv to near-ir wavelengths

Refractive-index and optical-absorption spectra of Bi-substituted yttrium iron garnet films, epitaxially grown by liquid-phase epitaxy, have been measured in the spectral regime 0.26 m1.9 m by thin-film interference for 0.52 m and by ellipsometry for0.52 m. The Y_3–x–y Bi _x Pb _y Fe_5–z Pt _z O₁₂ films contain bismuth in the range Ox 1.42, lead in the range 0.01 y0.08 and platinum in the range 0.005<=z0.03. There is satisfactory coincidence between the results from ellipsometry and thin-film interference in the overlapping wavelength region. The materials investigated are the same as reported earlier from this laboratory in ter mof their magnetic and magnetooptic properties.     

Weyl's geometry and physics

It is proposed to remove the difficulty of nonitegrability of length in the Weyl geometry by modifying the law of parallel displacement and using standard vectors. The field equations are derived from a variational principle slightly different from that of Dirac and involving a parameter . For =0 one has the electromagnetic field. For <0 there is a vector meson field. This could be the electromagnetic field with finite-mass photons, or it could be a meson field providing the missing mass of the universe. In cosmological models the two natural gauges are the Einstein gauge and the cosmic gauge. With the latter the universe has a fixed size, but the sizes of small systems decrease with time and their masses and energies increase, thus producing the Hubble effect. The field of a particle in this gauge is investigated, and it leads to an interesting solution of the Einstein equations that raises a question about the Schwarzschild solution.     

Antiferromagnetic Potts Models on the Square Lattice: A High-Precision Monte Carlo Study

We study the antiferromagnetic q-state Potts model on the square lattice for q=3 and q=4, using the Wang–Swendsen–Kotecký (WSK) Monte Carlo algorithm and a powerful finite-size-scaling extrapolation method. For q=3 we obtain good control up to correlation length 5000; the data are consistent with ()=Ae ^{2 p} (1+a ₁ e ^– + ...) as , with p1. The staggered susceptibility behaves as _stagg ^5/3. For q=4 the model is disordered (2) even at zero temperature. In appendices we prove a correlation inequality for Potts antiferromagnets on a bipartite lattice, and we prove ergodicity of the WSK algorithm at zero temperature for Potts antiferromagnets on a bipartite lattice.     

On the maximization of the central gravitational red-shift of a schwarzschild sphere with a given surface gravitational red-shift

Following Bondi static, spherically symmetric equilibrium configurations with a core and an envelope have been considered. It has been shown that for any configurations with nonnegative pressure and density and with a surface red-shiftz _s 4.77 arbitrarily large central red-shiftsz _c are possible in the limiting case of arbitrarily large radius. The effects of imposition of further constraints in the form of a real speed of sound not exceeding the speed of light are also examined. It is seen that for a given limiting sound-to-light-speed ratio . (i) There exists a limiting surface red-shiftz _s() 1.71. (ii) A configuration withz _s >z _s() is not possible, (iii) A configuration withz _s=z _s() has a unique and finitez _c=z _c(). (iv) Forz _s<z _s() arbitrarily large central red-shifts can be obtained for configurations with arbitrarily large radii.     

Space-and time-resolved diffusion-limited binary reaction kinetics in capillaries: experimental observation of segregation,anomalous exponents,and depletion zone

An experimental investigation of one-dimensional, diffusion-limited A+BC chemical reactions is reported. The persistence of reactant segregation and the formation of a depletion zone is observed and expressed in terms of the universal time exponents: (motion of the boundary zone), (width of instantaneous product formation zone), (rate of instantaneous local formation of product), (rate of instantaneous global formation of product), etc. There is good agreement with the recently predicted and/or simulated values:=1/2,=1/6,=2/3,=1/2, in contrast to classical predictions (=0,=1/2,=0,=–1/2). Furthermore, classically the segregation would not be preserved and there would be no formation of a depletion zone and no motion (just dissipation) of the reaction zone. We also discuss the relations to electrode oxidation-reduction reactions, i.e., A+CC where C is a catalyst, electrode, or trap.     

Thes-f model with Coulomb repulsion: Thermodynamics of two atomic cluster. Spin 1/2

A cluster of two atoms described by thes-f model with Coulomb repulsion has been considered. The interaction between localized 4f electrons (S=1/2) is taken in the molecular field approximation. The thermodynamic quantities like magnetization, specific heat and correlation functions n , n , S ^z n , S ^z n , S ^z (n –n ), n n and S ⁺ a ⁺ a as functions of temperature are presented for different band fillingN=0, 0.5, 1, 1.5, 2. The dependence of Curie temperature onN is calculated. The phase diagram forN=1 (T=0K) shows the possibility of existence of two phases: paramagnetic and ferromagnetic.The Curie temperature and the specific heat as functions ofN exhibit similar trends as found in experiments on doped magnetic semiconductors.     

On the validity of an einstein relation in models of interface dynamics

We consider models of interface dynamics derived from Ising systems with Kac interactions and we prove the validity of the Einstein relation=, where is the proportionality coefficient in the motion by curvature, is the interface mobility, and is the surface tension.     

Coherent states of theq-canonical commutation relations

For theq-deformed canonical commutation relationsa(f)a (g)=(1-q)f,g 1+qa (g)a(f) forf, g in some Hilbert space we consider representations generated from a vector satisfyinga(f)=<f, >, where . We show that such a representation exists if and only if 1. Moreover, for <1 these representations are unitarily equivalent to the Fock representation (obtained for =0). On the other hand representations obtained for different unit vectors are disjoint. We show that the universal C^*-algebra for the relations has a largest proper, closed, two-sided ideal. The quotient by this ideal is a naturalq-analogue of the Cuntz algebra (obtained forq=0). We discuss the conjecture that, ford<, this analogue should, in fact, be equal to the Cuntz algebra itself. In the limiting casesq=±1 we determine all irreducible representations of the relations, and characterize those which can be obtained via coherent states.Supported in part by the NSF(USA), and NATO Available by anonymous FTPfrom nostrom.physik.Uni-Osnabrueck.DE     

Proof of confinement of static quarks in 3-dimensionalU(1) lattice gauge theory for all values of the coupling constant

We study the 3-dimensional pureU(1) lattice gauge theory with Villain action which is related to the 3-dimensional -ferromagnet by an exact duality transformation (and also to a Coulomb system). We show that its string tension is nonzero for all values of the coupling constantg ², and obeys a bound const·m _D^–1 for smallag ², with =4²/g ² and . A continuum limita0,m _D fixed, exists and represents a scalar free field theory of massm _D. The string tension m _D ^–2 in physical units tends to in this limit. Characteristic differences in the behaviour of the model for large and small coupling constantag ² are found. Renormalization group aspects are discussed.Work supported in part by Deutsche Forschungsgemeinschaft     

Analytic study of power spectra of the tent maps near band-splitting transitions

Successive band-splitting transitions occur in the one-dimensional map x_i+1=g(x_i),i=0, 1, 2,... withg(x)=x, (0 x 1/2) –x +, (1/2 <x 1) as the parameter is changed from 2 to 1. The transition point fromN (=2ⁿ) bands to 2Nbands is given by=(2)^1/N (n=0, 1,2,...). The time-correlation function _i=x_ix₀/(x₀)²,x_i x_i–x_i is studied in terms of the eigenvalues and eigenfunctions of the Frobenius-Perron operator of the map. It is shown that, near the transition point=2, _i–[(10–42)/17] _i,0-[(102-8)/51]_i,1 + [(7 + 42)/17](–1)ⁱe^–yi, where2(–2) is the damping constant and vanishes at=2, representing the critical slowing-down. This critical phenomenon is in strong contrast to the topologically invariant quantities, such as the Lyapunov exponent, which do not exhibit any anomaly at=2. The asymptotic expression for _i has been obtained by deriving an analytic form of _i for a sequence of which accumulates to 2 from the above. Near the transition point=(2)^1/N, the damping constant of _i fori N is given by _N=2(^N-2)/N. Numerical calculation is also carried out for arbitrary a and is shown to be consistent with the analytic results.     

The spectrum of a quasiperiodic Schrödinger operator

The spectrum (H) of the tight binding Fibonacci Hamiltonian (H _mn= _m,n+1+ _m+1,n + _m,n v(n),v(n)= ((n–1)), 1/ is the golden number) is shown to coincide with the dynamical spectrum, the set on which an infinite subsequence of traces of transfer matrices is bounded. The point spectrum is absent for any , and (H) is a Cantor set for 4. Combining this with Casdagli's earlier result, one finds that the spectrum is singular continuous for 16.On leave from the Central Research Institute for Physics, Budapest, Hungary     

Theory of percolation in fluids of long molecules

We use the reference interaction site model (RISM) integral equation theory to study the percolation behavior of fluids composed of long molecules. We examine the roles of hard core size and of length-to-width ratio on the percolation threshold. The critical density _c is a nonmonotonic function of these parameters exhibiting competition of different effects. Comparisons with Monte Carlo calculations of others are reasonably good. For critical exponents, the theory yields =2=2 for molecules of any noninfinite lengthL. WhenL is very large, the theory yields _cL ^–2. These predictions compare favorably with observations of the conductivity for random assemblies of conductive fibers. The threshold region where asymptotic scaling holds requires the correlation length (/ _c ) ^–v to be much larger thanL. Evidently, the range of densities in this region diminishes asL increases, requiring that density deviations from _c be no larger thanL ^–2. Otherwise, crossover behavior will be observed.     

Perturbations of flows on Banach spaces and operator algebras

For automorphism groups of operator algebras we show how properties of the difference _t – ' _t are reflected in relations between the generators , . Indeed for a von Neumann algebraM with separable predual we show that if _t – '_t 0.28 for smallt, then = 0(+)°^-1 where is an inner automorphism ofM and is a bounded derivation ofM. If the difference _t – ' _t =O(t) ast ; 0, then = + and if _t – ' _t 0.28 for allt then =. We prove analogous results for unitary groups on a Hilbert space andC ₀,C ₀ ^* groups on a Banach space.This paper subsumes an earlier work of the same title which appeared as a report from Z.I.F. der Universität BielefeldWith partial support of the U.S. National Science Foundation     

On the stability of spin waves at finite temperatures

The temperature dependent exchange stiffnessD, defined for long wavelength spin waves by=Dq ², is calculated for very weak itinerant ferromagnets. It is found that for a general single band of d-electrons,D reduces to a very simple form givingD(T)=D(0)/ ₀ in the limit ₀0. The stability condition for spin waves at finite temperatures is thenD(0) > 0 and/ ₀ > 0, where is the relative occupation ± spin sub-bands and ₀ its value at 0°K.The authors are grateful to E. P. Wohlfarth for helpful discussions.     

Correlation functionals of infinite volume quantum spin systems

The existence and analyticity of the correlation functionals of a quantum lattice in the infinite volume limit is proved. The result is valid at sufficiently high temperatures and for a large class of interactions. Our method estimates the kernelK for a set of Kirkwood-Salzburg equations. While a naive estimate would indicate that K =, we take into account cancellations between different contributions toK in order to show that for sufficiently high temperatures K <1, and this estimate is independent of the volume of the system.Supported in part by the Air Force Office of Scientific Research.     

Zeta functions and transfer operators for piecewise monotone transformations

Given a piecewise monotone transformationT of the interval and a piecewise continuous complex weight functiong of bounded variation, we prove that the Ruelle zeta function (z) of (T, g) extends meromorphically to {z<^-1} (where =lim g^°T^n-1...g^°Tg ^1/n ) and thatz is a pole of if and only ifz ^–1 is an eigenvalue of the corresponding transfer operator L. We do not assume that L leaves a reference measure invariant.Research partially supported by the Fonds National Suisse     

Critical behavior of the site percolation model on the square lattice in aL×M geometry

Relevant aspects of the critical behavior of the site percolation model in aL×M geometry (LM) are studied. It is shown that this geometry favors the growth of percolating clusters in theL-direction with respect to those growing in theM-direction, causing pronounced finite-size effects on the percolation probabilities. Scaling functions have an additional parameter, namelyM, which introduces a dependence of these functions on the aspect ratioL/M. At criticality, the probability of a site belonging to the percolation clusters (P _L,M ) behaves likeP _L,M L ^/v (L/M) with =5/36 andv=4/3, where is a suitable scaling function. Using scaling arguments it is conjectured and then tested by means of Monte Carlo simulations, the following asymptotic behavior (L/M)(L/M), (L,M, =1), for the leading term. Systematic deviations of the Monte Carlo data from the conjectured behavior are due to second order corrections to the leading term which can also be under-stood on the basis of scaling ideas. Finite-size dependent critical probabilities are also functions ofL/M as it follows from scaling arguments which are corroborated by the simulations.Financially supported by the Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) de la República Argentina