On the Hamiltonian interpolation of near-to-the identity symplectic mappings with application to symplectic integration algorithms

We reconsider the problem of the Hamiltonian interpolation of symplectic mappings. Following Moser's scheme, we prove that for any mapping , analytic and -close to the identity, there exists an analytic autonomous Hamiltonian system, H such that its time-one mapping _H differs from by a quantity exponentially small in 1/. This result is applied, in particular, to the problem of numerical integration of Hamiltonian systems by symplectic algorithms; it turns out that, when using an analytic symplectic algorithm of orders to integrate a Hamiltonian systemK, one actually follows exactly, namely within the computer roundoff error, the trajectories of the interpolating Hamiltonian H, or equivalently of the rescaled Hamiltonian K=^-1H, which differs fromK, but turns out to be ⁵ close to it. Special attention is devoted to numerical integration for scattering problems.     

Continuity of type-I intermittency from a measure-theoretical point of view

We study a simple dynamical system which displays a so-called type-I intermittency bifurcation. We determine the Bowen-Ruelle measure and prove that the expectation (g) of any continuous functiong and the Kolmogoroff-Sinai entropyh() are continuous functions of the bifurcation parameter. Therefore the transition is continuous from a measure-theoretical point of view. Those results could be generalized to any similar dynamical system.     

Unmixing of binary alloys by a vacancy mechanism of diffusion: a computer simulation

The initial stages of phase separation are studied for a model binary alloy (AB) with pairwise interactions _AA , _AB , _BB between nearest neighbors, assuming that there is no direct interchange of neighboring atoms possible, but only an indirect one mediated by vacancies (V) occurring in the system at a concentrationc _v and which are strictly conserved, as are the concentrationsc _A andc _B of the two species.A-atoms may jump to vacant sites with jump rate _A , B-atoms with jump rate _B (in the absence of interactions). Particular attention is paid to the question to what extent nonuniform distribution of vacancies affects the unmixing kinetics. Our study focuses on the special case _A = _B on a square lattice, considering three different choices of interactions with the same = _AB – ( _AA + _BB )/2: (i) _AB =, _AA = _BB = 0; (ii) _AA = 0, _AA = _BB ; = ; (iii) _AB = _BB = 0, _AA = –2. We obtain both the time evolution of the structure factorS(k,t) following a quench from infinite temperature to the considered temperature, and the timedependence of the mean cluster size and the various neighborhood probabilities of a vacancy. While in case (i) forc _V 0.16 the distribution of vacancies in the system stays nearly random, in case (ii) the vacancies cluster in theA-B interfacial region, and in case (iii) they get nearly completely expelled from theA-rich regions. While phase separation proceeds in case (i) only slightly faster than in case (ii), a significant slowing down of the relaxation is observed for case (iii), which shows up in a strong reduction of the effective exponents describing the growth.     

The slow passage through a steady bifurcation: Delay and memory effects

We consider the following problem as a model for the slow passage through a steady bifurcation: dy/dt = (t) y – y³ +, where is a slowly increasing function oft given by= _i + t ( ⁱ,<0). Both and are small parameters. This problem is motivated by laser experiments as well as theoretical studies of laser problems. In addition, this equation is a typical amplitude equation for imperfect steady bifurcations with cubic nonlinearities. When=0, we have found that=0 is not the point where the bifurcation transition is observed. This transition appears at a value = ^j > 0. We call ^j the delay of the bifurcation transition. We study this delay as a function of ⁱ, the initial position of, and, the imperfection parameter. To this end, we propose an asymptotic study of this equation as 0, small but fixed. Our main objective is to describe this delay in terms of the relative magnitude of and. Since time-dependent imperfections are always present in experiments, we analyze in the second part of the paper the effect of a small-amplitude but time-periodic imperfection given by (t) = cos(t).     

Geometric expansion of the boundary free energy of a dilute gas

We consider a dilute classical gas in a volume ^–1 which tends to ^d by dilation as 0. We prove that the pressurep(^–1) isC ^q in at =0 (thermodynamic limit), for anyq, provided the boundary isC ^q and provided the Ursell functionsu _n (x ₁, ...,x _n) admit moments of degreeq and have nice derivatives.     

Optimal multigrid algorithms for the massive Gaussian model and path integrals

Multigrid algorithms are presented which, in addition to eliminating the critical slowing down, can also eliminate the volume factor. The elimination of the volume factor removes the need to produce many independent fine-grid configurations for averaging out their statistical deviations, by averaging over the many samples produced on coarse grids during the multigrid cycle. Thermodynamic limits of observables can be calculated to relative accuracy _r in justO( _r ^-2 ) computer operations, where _r is the error relative to the standard deviation of the observable. In this paper, we describe in detail the calculation of the susceptibility in the one-dimensional massive Gaussian model, which is also a simple example of path integrals. Numerical experiments show that the susceptibility can be calculated to relative accuracy _r in about 8 _r ^-2 random number generations, independent of the mass size.     

Some remarkable degenerations of quantum groups

We show that an irreducible representation of a quantized enveloping algebraU at a ^th root of 1 has maximal dimension (= ^N ) if the corresponding symplectic leaf has maximal dimension (=2N). The method of the proof consists of a construction of a sequence of degenerations ofU , the last one being aq-commutative algebraU ^(2N) . This allows us to reduce many problems concerningU to that concerningU ^(2N) .To Armand Borel on his 70^th birthdaySupported in part by the NSF grant DMS-9103792     

More inequalities for critical exponents

A variety of rigorous inequalities for critical exponents is proved. Most notable is the low-temperature Josephson inequalitydv +2 2–. Others are 1 1 +v, 1 1 , 1,d 1 + 1/ (for d),dv, ₃ + (for d), ₄ , and _{2m 2m+2} (form 2). The hypotheses vary; all inequalities are true for the spin-1/2 Ising model with nearest-neighbor ferromagnetic pair interactions.NSF Predoctoral Fellow (1976–1979). Research supported in part by NSF Grant PHY 78-23952.     

Energy and momentum in general relativity

In general relativity, conservation of energy and momentum is expressed by an equation of the form /x= 0, where –gT represents the total energy, momentum, and stress. This equation arises from the divergence formula dV _v = (/x ^v )d ⁴ d. Here we show that this formula fails to account properly for the system of basis vectors e(x). We obtain the (invariant) divergence formula e dV _v = e (/x ^v + )d ⁴ d. Conservation of energy and momentum is therefore expressed by the covariant equation (/x ^v ) + = 0. We go on to calculate the variation of the action under uniform displacements in space-time. This calculation yields the covariant equation of conservation, as well as the fully symmetric energy tensor . Finally, we discuss the transfer of energy and momentum, within the context of Einstein's theory of gravitation.     

Random shearing direction models for isotropic turbulent diffusion

Recently, a rigorous renormalization theory for various scalar statistics has been developed for special modes of random advection diffusion involving random shear layer velocity fields with long-range spatiotemporal correlations. New random shearing direction models for isotropic turbulent diffusion are introduced here. In these models the velocity field has the spatial second-order statistics of an arbitrary prescribed stationary incompressible isotropic random field including long-range spatial correlations with infrared divergence, but the temporal correlations have finite range. The explicit theory of renormalization for the mean and second-order statistics is developed here. With the spectral parameter, for –<<4 and measuring the strength of the infrared divergence of the spatial spectrum, the scalar mean statistics rigorously exhibit a phase transition from mean-field behavior for <2 to anomalous behavior for with 2<<4 as conjectured earlier by Avellaneda and the author. The universal inertial range renormalization for the second-order scalar statistics exhibits a phase transition from a covariance with a Gaussian functional form for with <2 to an explicit family with a non-Gaussian covariance for with 2<<4. These non-Gaussian distributions have tails that are broader than Gaussian as varies with 2<<4 and behave for large values like exp(–C _c |x|^4–), withC _c an explicit constant. Also, here the attractive general principle is formulated and proved that every steady, stationary, zero-mean, isotropic, incompressible Gaussian random velocity field is well approximated by a suitable superposition of random shear layers.     

Fine Structures of the Permittivity and Characteristic Losses of Electrons in Fluorite

Experimental-computational spectra of the permittivity and characteristic losses –Im^–1 for energies in the range 5–21 eV at a temperature of 4.2 K and theoretical spectra of and –Im^–1 of a fluorite crystal are resolved into elementary transition bands. The parameters of transition bands (energies of their maxima E _i, band halfwidths H _i and areas S _i, and oscillator forces f _i) are determined. A correlation of the spectral bands of and –Im^–1is established, and their specific features are elucidated.     

Band-gap structure of the spectrum of periodic Maxwell operators

We investigate the band-gap structure of some second-order differential operators associated with the propagation of waves in periodic two-component media. Particularly, the operator associated with the Maxwell equations with position-dependent dielectric constant (x),xR ³, is considered. The medium is assumed to consist of two components: the background, where (x) = _b , and the embedded component composed of periodically positioned disjoint cubes, where (x) = _a . We show that the spectrum of the relevant operator has gaps provided some reasonable conditions are imposed on the parameters of the medium. Particularly, we show that one can open up at least one gap in the spectrum at any preassigned point provided that the size of cubesL, the distancel=L betwen them, and the contrast = _b / _a are chosen in such a way thatL ^–2, and quantities ^-1-3/2 and ² are small enough. If these conditions are satisfied, the spectrum is located in a vicinity of widthw(^3/2)^-1 of the set {² L ^-2 k ²:kZ³}. This means, in particular, that any finite number of gaps between the elements of this discrete set can be opened simultaneously, and the corresponding bands of the spectrum can be made arbitrarily narrow. The method developed shows that if the embedded component consists of periodically positioned balls or other domains which cannot pack the space without overlapping, one should expect pseudogaps rather than real gaps in the spectrum.     

Power law growth for the resistance in the Fibonacci model

Many one-dimensional quasiperiodic systems based on the Fibonacci rule, such as the tight-binding HamiltonianH(n)=(n+1)+(n–1)+v(n) (n),n,l ²(),, wherev(n)=[(n+1)]–[n],[x] denoting the integer part ofx and the golden mean , give rise to the same recursion relation for the transfer matrices. It is proved that the wave functions and the norm of transfer matrices are polynomially bounded (critical regime) if and only if the energy is in the spectrum of the Hamiltonian. This solves a conjecture of Kohmoto and Sutherland on the power-law growth of the resistance in a one-dimensional quasicrystal.     

Energy level crossings in quantum mechanics

From the eigenvalue equationH \ _n () =E _n ()\ _n () withH H ₀ +V one can derive an autonomous system of first order differential equations for the eigenvaluesE _n () and the matrix elementsV _mn () where is the independent variable. To solve the dynamical system we need the initial valuesE _n ( = 0) and \ _n ( = 0). Thus one finds the motion of the energy levelsE _n (). We discuss the question of energy level crossing. Furthermore we describe the connection with the stationary state perturbation theory. The dependence of the survival probability as well as some thermodynamic quantities on is derived. This means we calculate the differential equations which these quantities obey. Finally we derive the equations of motion for the extended caseH =H ₀ +V ₁ + ² V ₂ and give an application to a supersymmetric Hamiltonian.     

Are there critical points on the boundaries of singular domains?

Generalizing a result of E. Ghys, we prove a general theorem that implies that if a rational functionf of the Riemann sphere of degree 2 leaves invariant a singular domainC (a disk or a ring) on which the rotation number off satisfies a diophantine condition, provided that on f is injective, then each boundary component ofC contains critical point off. The injectivity condition is always satisfied for singular disks associated to linearizable periodic elliptic points off(z)=z ⁿ +a, withn,n2 anda. We also show that the singular disks, associated to periodic elliptic points off(z)=e ^az that satisfy a diophantine condition, are unbounded in . In the end of the paper, we give a survey of the theory of iteration of entire functions of .     

Collisionless kinetic equations for massive particles with spins 1/2 and 0 in gravitational and electromagnetic fields

We define a covariant and gauge-invariant generalization of the Wigner functions of particles with spins 1/2 and 0. The collisionless kinetic equations are obtained for these particles in external gravitational and electromagnetic fields in the quasiclassical approximation; also obtained are the momentum representations of the energy-momentum tensor, current, and spin tensor, taking into account the effects of the spin's interaction with the gravitational field an electromagnetic field. The following notation is used: e and m are the charge and mass of the particles; is Planck's constant; (x) are the covariant-fixed Dirac matrices; ^,=(1/4)[, ]: a(^b)=(1/2) (a b +ab ); [A, B]=A·B – B·A; A,B=A·B+B·A; g(x)=det(g (x));R = –...; the speed of light c=1.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 9, pp. 47–53, September, 1990.The author wishes the thank Yu. G. Ignat'ev and members of the seminar in General Relativistic Statistics and Cosmology of the Kazan' Pedagogical Institute for useful discussions.     

Feynman integral and complex classical trajectories

Let _t be an analytic solution of the Schrödinger equation with the initial condition . Let _t be the solution of the Schrödinger equation with the initial condition =, where is an analytic function. When 0, then _t (x) _t (x)¹ ( _t (x)), where _t (x) trajectory starting from x. We relate this result to Feynman's sum over trajectories and complex stochastic differential equations.     

Замечание К теории по степенного возникновения подви жных слоев в плазме инертных газов

[1] t _B , t _B . , t _B , . .

---

A note on the theory of the successive production of moving striations in the plasma of inert gases

Approximate expressions are derived on the basis of Pekárek's theoretical paper [1] for the period of the maximum t_B of a wave packet produced by the passage of a wave of stratification before the aperture of a photomultiplier, and for its time width in the half-height t_B. The relaxation time of a wave of stratification, following from the theory [1], can thus be calculated by means of the experimentally measured velocity of motion of the maximum of a wave packet u and its width t_B. The calculation is supplemented by numerical data on the magnitude of errors committed by using approximate expressions.

     

Semi-infinite Potts model and percolation at surfaces

The effects of surfaces on percolation are investigated near the bulk percolation threshold ind=6– dimensions. Using field-theoretic methods, this is done within the framework of a semi-infinite continuousq-state Potts model withq1. Renormalization-group equations are obtained which imply that the usual scaling laws for surface and bulk exponents are valid to all orders in , and the surface exponents at the ordinary and special transition are computed to order . Our result for ₁ ^ord is in conformity with the one by Carton.     

Colored noise in activated rate processes

We consider bistable systems driven by stationary wideband Gaussian colored noise. We construct uniform asymptotic expansions of the stationary probability density function and of the activation rate, for small intensity and short correlation time of the noise. We find that for different values of the total power output / of the noise, different terms in the asymptotic expansions become dominant. For we recover previously derived results, while for =O() and new results are obtained.