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.     

Holes in the probability density of strongly colored noise driven systems

A qualitative change in the topology of the joint probability densityP(,x), which occurs for strongly colored noise in multistable systems, has recently been observed first by analog simulation (F. Moss and F. Marchesoni,Phys. Lett. A 131:322 (1988)) and confirmed by matrix continued fraction methods (Th. Leiber and H. Riskin, unpublished), and by analytic theory (P. Hänggi, P. Jung, and F. Marchesoni,J. Stat. Phys., this issue). Systems studied were of the classx=–U(x)/x+(t,), whereU(x) is a multistable potential and (t, ) is a colored, Gaussian noise of intensityD, for which =0, and (t) (s)=(D/)exp(–t–s/). When the noise correlation time is smaller than some critical value ₀, which depends onD, the two-dimensional densityP(,x) has the usual topology [P. Jung and H. Risken,Z. Phys. B 61:367 (1985); F. Moss and P. V. E. McClintock,Z. Phys. B 61:381 (1985)]: a pair of local maxima ofP(,x), which correspond to a pair of adjacent local minima ofU(x), are connected by a single saddle point which lies on thex axis. When >₀, however,the single saddle disappears and is replaced by a pair of off-axis saddles. A depression, or hole, which is bounded by the saddles and the local maxima thus appears. The most probable trajectory connecting the two potential wells therefore does not pass through the origin for >₀, but instead must detour around the local barrier. This observation implies that successful mean-first-passage-time theories of strongly colored noise driven systems must necessarily be two dimensional (Hänggiet al.). We have observed these holes for several potentialsU(x): (1)a soft, bistable potential by analog simulation (Moss and Marchesoni); (2) a periodic potential [Th. Leiber, F. Marchesoni, and H. Risken,Phys. Rev. Lett. 59:1381 (1987)] by matrix continued fractions; (3) the usual hard, bistable potential,U(x)=–ax ²/2+bx ⁴/4, by analog simulations only; and (4) a random potential for which the forcingf(x)=–U(x)/x is an approximate Gaussian with nonzero correlation length, i.e., colored spatiotemporal noise, by analog simulation. There is a critical curve ₀(D) in the versusD plane which divides the two topological behaviors. For a fixed value ofD, this curve is shifted toward larger values of ₀ for progressively weaker barriers between the wells. Therefore, strong barriers favor the observation of this topological transformation at smaller values of . Recently, an analytic expression for the critical curve, valid asymptotically in the small-D limit, has been obtained (Hänggiet al.).This paper will appear in a forthcoming issue of theJournal of Statistical Physics.     

A two-level Kaluza-Klein theory

Let the Lie groups G and H act on the manifold P in such a way that P fibres as a principal G-bundle over P/G and as an H-bundle over P/H. We find that every pair (,) where is an H-invariant connection form in PP/G and is a G-invariant connection form in PP/H corresponds uniquely to a connection form in PP/(H×G) and a cross-section of a vector bundle with base P/(H×G).     

Carrier lifetime in single-crystal PbS films

The temperature dependence of the majority carrier lifetime was studied in single-crystal PbS films. In n- and p-type films in the range 300-200 ° K, increases exponentially with decreasing temperature, with an activation energy of e O. 17-0. 2 eV. This m dependence is assumed due to the trapping of minority carriers at deep levels. Below 160 ° K in the n-type films, does not depend on the temperature, while in the p-type films, decreases exponentially with an activation energy of 0. 11–0. 14 eV.Translated from Vysshikh Uchebnykh Zavedenii Fizika, No. 6, pp. 64–67, June, 1970.     

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.     

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.     

Effective Potential for \mathcal{P}\mathcal{T}-Symmetric Quantum Field Theories

Recently, a class of -invariant scalar quantum field theories described by the non-Hermitian Lagrangian = () ² +g ² (i) was studied. It was found that there are two regions of . For <0 the -invariance of the Lagrangian is spontaneously broken, and as a consequence, all but the lowest-lying energy levels are complex. For 0 the -invariance of the Lagrangian is unbroken, and the entire energy spectrum is real and positive. The subtle transition at =0 is not well understood. In this paper we initiate an investigation of this transition by carrying out a detailed numerical study of the effective potential V _eff (_c) in zero-dimensional spacetime. Although this numerical work reveals some differences between the <0 and the >0 regimes, we cannot yet see convincing evidence of the transition at =0 in the structure of the effective potential for -symmetric quantum field theories.     

Macroscopic short-time electrodynamics for conductors

In the macroscopic electrodynamics (MED) of good conductors (metals) based on Ohm's law j=E, the momentum relaxation time =m/ne² of the electrons limits the application to electromagnetic (EM) processes with characteristic timest. An interesting physical difficulty occurs in MED since the EM field damping time _R=/ of metals is very small compared with the minimum macroscopic time scale, _R. Consequently, the damping and propagation of EM waves and pulses in good conductors cannot be correctly described within the frame of conventional MED. New hyperbolic EM wave equations with relaxation and memory are proposed, which no longer exhibit the _R deficiency. The latter is caused by Ohm's law, which breaks down for short-time processes, due to neglect of electron inertia. The advantages of the proposed and the disadvantages of the conventional EM wave equations for good conductors are discussed in applications.     

A stochastic particle system modeling the Carleman equation

Two species of Brownian particles on the unit circle are considered; both have diffusion coefficient >0 but different velocities (drift), 1 for one species and –1 for the other. During the evolution the particles randomly change their velocity: if two particles have the same velocity and are at distance ( being a positive parameter), they both may simultaneously flip their velocity according to a Poisson process of a given intensity. The analogue of the Boltzmann-Grad limit is studied when goes to zero and the total number of particles increases like ^–1. In such a limit propagation of chaos and convergence to a limiting kinetic equation are proven globally in time, under suitable assumptions on the initial state. If, furthermore, depends on and suitably vanishes when goes to zero, then the limiting kinetic equation (for the density of the two species of particles) is the Carleman equation.Dedicated to the memory of Paola Calderoni.     

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.     

Transport of light radiation in a spatially finite three-dimensional scattering medium

We propose a new method for solving radiation transport problems, which permits including in analytic form for the case of normal incidence the effect of spatial finiteness of the scattering medium. The formation of the light field accompanying changes in the optical parameters and optical dimensions of the medium is analyzed.this paper, we examine the simplest case of a geometry of a scattering medium in the form of a parallelipiped with optical length _x, height _y, and width _z. The analysis is performed for the case =1, _y = _z with the latter varying in the range 0.1 to . The results obtained show that the light field depends strongly on the optical dimensions of the medium. The limiting values of the optical dimensions (_y = _z), beginning with which the spatial finiteness of the medium can be neglected, are determined.Translated from Izvestiya Vysshikh Uchebnykh Zavednii, Fizika, No. 8, pp. 82–85, August, 1982.     

A matrix method for estimating the Liapunov exponent of one-dimensional systems

Let : [0, 1][0, 1] be a piecewise monotonie expanding map. Then admits an absolutely continuous invariant measure. A result of Kosyakin and Sandler shows that can be approximated by a sequence of absolutely continuous measures _n invariant under piecewise linear Markov maps _itn. Each _itn is constructed on the inverse images of the turning points of . The easily computable measures _n are used to estimate the Liapunov exponent of . The idea of using Markov maps for estimating the Liapunov exponent is applied to both expanding and nonexpanding maps.     

Characteristics analysis of metal-clad and absorptive dielectric waveguides by a simple and accurate perturbation method

We present a simple and accurate method for characteristic analysis of metal-clad dielectric waveguides and absorptive waveguides. The real partN of the complex modal indexN=N + iN is obtained by solving the corresponding real eigenvalue equation, and the imaginary partN is given by (n/), where= + i is the complex dielectric constant of the absorptive layer, and N/ is obtained by numerical differentiation. The method is straightforward, and the cumbersome solution of complex transcendental equations is completely eliminated. Results for simple structures are in good agreement with those obtained by exact analysis.     

On the theory of amplification of meso-ultrasound in anisotropic monopolar semiconductors

The electronic absorption coefficient (_e) of meso-ultrasound (_p ^–1 _en, _t, _ee; ql 1, where _p, _en, _t are the times of pulse relaxation, energy, and current carrier traps; _ee ^–1 is the frequency of intraelectron collisions; , q are the frequency and wave vector of sound; l is the carrier mean free path) in the presence of a permanent external field E₀ e ₀ ^–1 q₀T is calculated for anisotropic single-valley semiconductors with piezoelectrical and potential-deformation acoustoelectronic interaction. Considered arbitrary are: 1) the anisotropy of the tensors * and _p (and other crystal parameters); 2) the degree of Fermi degeneration of the carriers; 3) the dependence of _p on the carrier energy . The acoustoheat nonlinearity is neglected. The possibility of meso-ultrasound amplification by the transverse field (E₀ - q) is predicted. By changing the orientation of the vector q relative to the crystal axes, the transverse threshold field (E ^thr ) can be controlled smoothly and within broad limits, This permits the production of a mechanical sound-amplification regulator or a modulator of its intensity. The ratio between E ^thr and E ^thr yields a measure of the relative anisotropy of the mobility.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 2, pp. 87–92, February, 1977.     

Role of Activator and Vacancy Centers in the Decay of the Gamma Scintillators of CsI–Na Crystals

The influence of the concentration of an activator (C _NaI) and of plastic deformation on a change in the contribution of a slow component to the decay of the -scintillations of CsI–Na crystals was investigated, as well as the influence of C _NaI on a change in the shape of the luminescence excitation spectrum in the region of absorption of activator centers (AC) and of vacancy-related centers (VRC) and also on a change in the intercenter time of deexcitation of the centers indicated. It is shown that AC and VRC participate in the photoluminescence and -scintillations of CsI–Na crystals. In the -scintillations, AC are responsible for the component ₁ = 370 nsec, whereas the components ₁ = 460 nsec and ₂ 2 sec are associated with VRC. The reduction of ₁ from 770 to 560 nsec with an increase in C _NaI from 2·10^–3 to 3·10^–2 mole% and from 570 to 470 nsec after plastic deformation of the crystals ( = 5%) along the 111 axis is caused by a decrease in the number of VRC. The mechanisms underlying the -scintillations of the CaI–Na crystals containing AC and VRC and also the decrease in the number of VRC are discussed.     

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.     

Attenuation of a narrow light beam in a turbid medium

In model scattering media, in which the extinction coefficient and the probability of photon lifetime vary, the relation between the illumination near a parallel light beam and the optical depth is studied. It is shown that in media of low turbidity, Bouguer's Law holds at small penetration depths, because the scattered radiation is insignificant in comparison to the total energy balance. For beams of finite width, with constant , the illumination increases with increasing values. The importance of in relations describing the attenuation of a parallel light beam in turbid media is clarified.     

Correlation tensors of the blackbody radiation in rectangular cavities

Electromagnetic equilibrium fluctuations in finite cavities filled with a dissipative medium (dielectric function ()=+i) and bounded by walls of infinite conductivity are considered. Expanding the fields in terms of a complete and orthonormal set of functions and solving the Maxwell equations the response of the EM field to external forces (polarization and magnetization) is obtained. With the aid of the fluctuation dissipation theorem and the linear response functions the 2nd order correlation tensors of the EM field are derived.For rectangular cavities explicit considerations are made. In the case of transparent media (=0) the spectral energy density of the EM radiation is calculated.     

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