title: Short-Time Effects on Eigenstate Structure in Sinai Billiards and Related Systems

abstract: There is much latitude between the requirements of Schnirelman's theorem regarding the ergodicity of individual high-energy eigenstates of classically chaotic systems on the one hand, and the extreme requirements of random matrix theory on the other. It seems likely that some eigenstate statistics and long-time transport behavior bear nonrandom imprints of the underlying classical dynamics while simultaneously obeying Schnirelman's theorem. Indeed this was shown earlier in the case of systems which approach classical ergodicity slowly, and is also realized in the scarring of eigenstates, even in the $\hbar \to 0$ limit, along unstable periodic orbits and their manifolds. Here we demonstrate the nonrandom character of eigenstates of Sinai-like systems. We show that mixing between channels in Sinai systems is dramatically deficient compared to random matrix theory predictions. The deficit {\it increases} as $|\log \hbar|$ for $\hbar\to 0$, and is due to the vicinity of the measure zero set of orbits which never collide with the Sinai obstruction. Coarse graining to macroscopic scales recovers the Schnirelman result. Three systems are investigated here: a Sinai-type billiard, a quantum map which possesses the essential properties of the Sinai billiard, and a unitary map corresponding to a quasirandom Hamiltonian. Various wavefunction and long-time transport statistics are defined, theoretically investigated, and compared to numerical data.


title: Chaotic enhancement of decay. The effect of classical phase space structure

abstract: We investigate the decay process from a time dependent potential well in the semiclassical regime. The classical dynamics is chaotic and the decay rate shows an irregular behavior as a function of the system parameters. By studying the weak-chaos regime we are able to connect the decay irregularities to the presence of nonlinear resonances in the classical phase space. A quantitative analytical prediction which accounts for the numerical results is obtained.


title: Truncations of random unitary matrices

abstract: We analyze properties of non-hermitian matrices of size M constructed as square submatrices of unitary (orthogonal) random matrices of size N>M, distributed according to the Haar measure. In this way we define ensembles of random matrices and study the statistical properties of the spectrum located inside the unit circle. In the limit of large matrices, this ensemble is characterized by the ratio M/N. For the truncated CUE we derive analytically the joint density of eigenvalues from which easily all correlation functions are obtained. For N-M fixed and N--> infinity the universal resonance-width distribution with N-M open channels is recovered.


title: Fractal escapes in Newtonian and relativistic multipole gravitational fields

abstract: We study the planar motion of test particles in gravitational fields produced by an external material halo, of the type found in many astrophysical systems, such as elliptical galaxies and globular clusters. Both the Newtonian and the general-relativistic dynamics are examined, and in the relativistic case the dynamics of both massive and massless particles are investigated. The halo field is given in general by a multipole expansion; we restrict ourselves to multipole fields of pure order, whose Newtonian potentials are homogeneous polynomials in cartesian coordinates. A pure (n)-pole field has (n) different escapes, one of which is chosen by the particle according to its initial conditions. We find that the escape has a fractal dependency on the initial conditions for (n>2) both in the Newtonian and the relativistic cases for massive test particles, but with important differences between them. The relativistic motion of massless particles, however, was found to be regular for all the fields we could study. The box-counting dimension was used in each case to quantify the sensitivity to initial conditions which arises from the fractality of the escape route.


title: Heat conduction in 2d nonlinear lattices

abstract: The divergence of the heat conductivity in the thermodynamic limit is investigated in 2d-lattice models of anharmonic solids with nearest-neighbour interaction from single-well potentials. Two different numerical approaches based on nonequilibrium and equilibrium simulations provide consistent indications in favour of a logarithmic divergence in "ergodic", i.e. highly chaotic, dynamical regimes. Analytical estimates obtained in the framework of linear-response theory confirm this finding, while tracing back the physical origin of this {\sl anomalous} transport to the slow diffusion of the energy of long-wavelength effective Fourier modes. Finally, numerical evidence of {\sl superanomalous} transport is given in the weakly chaotic regime, typically found below some energy density threshold.


title: Chaos and Fractals in Geodesic Motions Around a Non-Rotating Black-Hole with an External Halo

abstract: We investigate the occurrence chaos in the escape of test particles moving in the field of a Schwarzschild black hole surrounded by an external halo. The motion of both material particles and zero rest mass particles is considered. The chaos is characterized by the fractal dimension of boundary between the basins of the different escapes, which is a topologically invariant characterization. We find chaos in the motion of both material particles and null geodesics.


title: Action Correlations in Integrable Systems

abstract: In many problems of quantum chaos the calculation of sums of products of periodic orbit contributions is required. A general method of computation of these sums is proposed for generic integrable models where the summation over periodic orbits is reduced to the summation over integer vectors uniquely associated with periodic orbits. It is demonstrated that in multiple sums over such integer vectors there exist hidden saddle points which permit explicit evaluation of these sums. Saddle point manifolds consist of periodic orbits vectors which are almost mutually parallel. Different problems has been treated by this saddle point method, e.g. Berry's bootstrap relations, mean values of Green function products etc. In particular, it is obtained that suitably defined 2-point correlation form-factor for periodic orbit actions in generic integrable models is proportional to quantum density of states and has peaks at quantum eigenenergies.


title: Diffractive corrections in the trace formula for polygonal billiards

abstract: We derive contributions to the trace formula for the spectral density accounting for the role of diffractive orbits in two-dimensional polygonal billiards. In polygons, diffraction typically occurs at the boundary of a family of trajectories. In this case the first diffractive correction to the contribution of the family to the periodic orbit expansion is of order of the one of an isolated orbit, and gives the first $\sqrt{\hbar}$ correction to the leading semi-classical term. For treating these corrections Keller's geometrical theory of diffraction is inadequate and we develop an alternative approximation based on Kirchhoff's theory. Numerical checks show that our procedure allows to reduce the typical semi-classical error by about two orders of magnitude. The method permits to treat the related problem of flux-line diffraction with the same degree of accuracy.


title: Density-dependent diffusion in the periodic Lorentz gas

abstract: We study the deterministic diffusion coefficient of the two-dimensional periodic Lorentz gas as a function of the density of scatterers. Results obtained from computer simulations are compared to the analytical approximation of Machta and Zwanzig [Phys.Rev.Lett. 50, 1959 (1983)] showing that their argument is only correct in the limit of high densities. We discuss how the Machta-Zwanzig argument, which is based on treating diffusion as a Markovian hopping process on a lattice, can be corrected systematically by including microscopic correlations. We furthermore show that, on a fine scale, the diffusion coefficient is a non-trivial function of the density. We finally argue that, on a coarse scale and for lower densities, the diffusion coefficient exhibits a Boltzmann-like behavior, whereas for very high densities it crosses over to a regime which can be understood qualitatively by the Machta-Zwanzig approximation.


title: Relaxation and Diffusion for the Kicked Rotor

abstract: The dynamics of the kicked-rotor, that is a paradigm for a mixed system, where the motion in some parts of phase space is chaotic and in other parts is regular is studied statistically. The evolution (Frobenius-Perron) operator of phase space densities in the chaotic component is calculated in presence of noise, and the limit of vanishing noise is taken is taken in the end of calculation. The relaxation rates (related to the Ruelle resonances) to the invariant equilibrium density are calculated analytically within an approximation that improves with increasing stochasticity. The results are tested numerically. The global picture of relaxation to the equilibrium density in the chaotic component when the system is bounded and of diffusive behavior when it is unbounded is presented.


title: Relaxation to the Invariant Density for Kicked Rotor

abstract: The relaxation rates to the invariant density in the chaotic phase space component of the kicked rotor (standard map) are calculated analytically for a large stochasticity parameter, $K$. These rates are the logarithms of the poles of the matrix elements of the resolvent, $ \hat{R}(z)=(z-\hat{U})^{-1} $, of the classical evolution operator $\hat{U} $. The resolvent poles are located inside the unit circle. For hyperbolic systems this is a rigorous result, but little is known about mixed systems such as the kicked rotor. In this work, the leading relaxation rates of the kicked rotor are calculated in presence of noise, to the leading order in $1/\sqrt{K}$. Then the limit of vanishing noise is taken and the relaxation rates are found to be finite, corresponding to poles lying inside the unit circle. It is found that the slow relaxation rates, in essence, correspond to diffusion modes in the momentum direction. Faster relaxation modes intermix the motion in the momentum and the angle space. The slowest relaxation rate of distributions in the angle space is calculated analytically by studying the dynamics of inhomogeneities projected down to this space. The analytical results are verified by numerical simulations.


title: Different Traces of Quantum Systems Having the Same Classical Limit

abstract: Many quantum systems may have the same classical limit. We argue that in the classical limit their traces do not necessarily converge one to another. The trace formula allows to express quantum traces by means of classical quantities as sums over periodic orbits of the classical system. To explain the lack of convergence of the traces we need the quantum corrections to the classical actions of periodic orbits. The four versions of the quantum baker map on the sphere serve as an illustration of this problem.


title: Can Strange Nonchaotic Dynamics be induced through Stochastic Driving?

abstract: Upon addition of noise, chaotic motion in low-dimensional dynamical systems can sometimes be transformed into nonchaotic dynamics: namely, the largest Lyapunov exponent can be made nonpositive. We study this phenomenon in model systems with a view to understanding the circumstances when such behaviour is possible. This technique for inducing ``order'' through stochastic driving works by modifying the invariant measure on the attractor: by appropriately increasing measure on those portions of the attractor where the dynamics is contracting, the overall dynamics can be made nonchaotic, however {\it not} a strange nonchaotic attractor. Alternately, by decreasing measure on contracting regions, the largest Lyapunov exponent can be enhanced. A number of different chaos control and anticontrol techniques are known to function on this paradigm.


title: Structure of Quantum Chaotic Wavefunctions: Ergodicity, Localization, and Transport

abstract: We discuss recent developments in the study of quantum wavefunctions and transport in classically ergodic systems. Surprisingly, short-time classical dynamics leaves permanent imprints on long-time and stationary quantum behavior, which are absent from the long-time classical motion. These imprints can lead to quantum behavior on single-wavelength or single-channel scales which are very different from random matrix theory expectations. Robust and quantitative predictions are obtained using semiclassical methods. Applications to wavefunction intensity statistics and to resonances in open systems are discussed.


title: Quantum limit of the laser linewidth in chaotic cavities and statistics of residues of scattering matrix poles

abstract: The quantum-limited linewidth of a laser cavity is enhanced above the Schawlow-Townes value by the Petermann factor K, due to the non-orthogonality of the cavity modes. We derive the relation between the Petermann factor and the residues of poles of the scattering matrix and investigate the statistical properties of the Petermann factor for cavities in which the radiation is scattered chaotically. For a single scattering channel we determine the complete probability distribution of K and find that the average Petermann factor $<K>$ depends non-analytically on the area of the opening, and greatly exceeds the most probable value. For an arbitrary number N of scattering channels we calculate $<K>$ as a function of the decay rate $\Gamma$ of the lasing mode. We find for $N\gg 1$ that for typical values of $\Gamma$ the average Petermann factor $<K>\propto \sqrt{N}\gg 1$ is parametrically larger than unity.


title: The passive polymer problem

abstract: In this article, we introduce a generalization of the diffusive motion of point-particles in a turbulent convective flow with given correlations to a polymer or membrane. In analogy to the passive scalar problem we call this the passive polymer or membrane problem. We shall focus on the expansion about the marginal limit of velocity-velocity correlations which are uncorrelated in time and grow with the distance x as |x|^epsilon, and epsilon small. This relation gets modified for polymers and membranes (the marginal advecting flow has correlations which are shorter ranged.) The construction is done in three steps: First, we reconsider the treatment of the passive scalar problem using the most convenient treatment via field theory and renormalization group. We explicitly show why IR-divergences and thus the system-size appear in physical observables. In a second step, we reformulate the problem in terms of a Langevin equation. This is not only interesting in its own, but also gives an efficient algorithm to determine S^{2n} numerically. In a third and final step, we generalize the Langevin treatment of a particle to polymers and membranes, or more generally to an elastic object of inner dimension D with 0 <= D <= 2. These objects can intersect each other. We also analyze what happens when self-intersections are no longer allowed.


title: Two-Particle Dispersion in Model Velocity Fields

abstract: We consider two-particle dispersion in a velocity field, where the relative two-point velocity scales according to $v^{2}(r)\propto r^{\alpha}$ and the corresponding correlation time scales as $\tau (r)\propto r^{\beta}$, and fix $\alpha =2/3$, as typical for turbulent flows. We show that two generic types of dispersion behavior arize: For $\alpha /2+\beta < 1$ the correlations in relative velocities decouple and the diffusion approximation holds. In the opposite case, $\alpha /2+\beta >1$, the relative motion is strongly correlated. The case of Kolmogorov flows corresponds to a marginal, nongeneric situation.


title: Artificiality of multifractal phase transitions

abstract: A multifractal phase transition is associated to a nonanalyticity in the generalised dimensions. We show that its occurrence is an artifact of the asymptotic scaling behaviour of integral moments and that it is not observed in an analysis based on differential n-point correlation densities.


title: Spectral properties of Atoms in Fields: A Semiclassical Analysis

abstract: We develop a semiclassical theory for the spectral rigidity of non-hydrogenic Rydberg atoms in electric fields and evaluate the significant deviations from the well-known Poissonian behaviour in the hydrogenic case. The resulting formula is shown to be in excellent agreement with the exact quantal result. We also investigate diamagnetic atoms; we find that, in contrast to the classically integrable atoms, diffraction has a small effect on the spectral rigidity in the classically chaotic atom. We show our predictions can also be of use in the mixed phase space regime.


title: Semiclassical Gaussian matrix elements for chaotic quantum wells

abstract: We derive semiclassical expressions for spectra, weighted by matrix elements of a Gaussian observable, relevant to a range of molecular and mesoscopic systems. We apply the formalism to the particular example of the resonant tunneling diode (RTD) in tilted fields. The RTD is a new experimental realization of a mesoscopic system exhibiting a transition to chaos. It has generated much interest and several different semiclassical theories for the RTD have been proposed recently.Our formalism clarifies the relationship between the different approaches and to previous work on semiclassical theories of matrix elements. We introduce three possible levels of approximation in the application of the stationary phase approximation, depending on typical length scales of oscillations of the semiclassical Green's function, relative to the degree of localization of the observable. Different types of trajectories (periodic, normal, closed and saddle orbits) are shown to arise from such considerations. We propose a new type of trajectory (``minimal orbits'')and show they provide the best real approximation to the complex saddle points of the stationary phase approximation. We test the semiclassical formulae on quantum calculations and experimental data. We show that the new types of trajectories (saddle and minimal orbits) provide accurate results. We discuss a divergence of the contribution of saddle orbits, which suggests the existence of bifurcation-type phenomena affecting the complex and non-periodic saddle orbits.


title: Direct measurement of intersection angle of invariant manifolds for area preserving mappings

abstract: Intersection angles of stable and unstable manifolds for area preserving mappings are numerically calculated by extremely accurate computation. With the use of multiprecision library the values of angle as small as 10^{-400} are obtained. The singular dependence of the angle on the magnitude of hyperbolicity is confirmed. The power-law type prefactor with Stokes constant is also in good agreement with analytical estimation.


title: Polymer Stretching by Turbulence

abstract: The stretching of a polymer chain by a large scale chaotic flow is considered. The steady state which emerges as a balance of the turbulent stretching and anharmonic resistance of the chain is quantitatively described, i.e. the dependency on the flow parameters (Lyapunov exponent statistics) and the chain characteristics (the number of beads and the inter-bead elastic potential) is made explicit.


title: Forecasting the SST space-time variability of the Alboran Sea with genetic algorithms

abstract: We propose a nonlinear ocean forecasting technique based on a combination of genetic algorithms and empirical orthogonal function (EOF) analysis. The method is used to forecast the space-time variability of the sea surface temperature (SST) in the Alboran Sea. The genetic algorithm finds the equations that best describe the behaviour of the different temporal amplitude functions in the EOF decomposition and, therefore, enables global forecasting of the future time-variability.


title: Linear and nonlinear time series analysis of the black hole candidate Cygnus X-1

abstract: We analyze the variability in the X-ray lightcurves of the black hole candidate Cygnus X-1 by linear and nonlinear time series analysis methods. While a linear model describes the over-all second order properties of the observed data well, surrogate data analysis reveals a significant deviation from linearity. We discuss the relation between shot noise models usually applied to analyze these data and linear stochastic autoregressive models. We debate statistical and interpretational issues of surrogate data testing for the present context. Finally, we suggest a combination of tools from linear andnonlinear time series analysis methods as a procedure to test the predictions of astrophysical models on observed data.


title: On the Turbulent Dynamics of Polymer Solutions

abstract: We study properties of dilute polymer solutions which are known to depend strongly on polymer elongation. The probability density function (PDF) of polymer end-to-end extensions $R$ in turbulent flows is examined. We demonstrate that if the value of the Lyapunov exponent $\lambda$ is smaller than the inverse molecular relaxation time $1/\tau$ then the PDF has a strong peak at the equilibrium size $R_0$ and a power tail at $R\gg R_0$. This confirms and extends the results of \cite{Lumley72}. There is no essential influence of polymers on the flow in the regime $\lambda\tau<1$. At $\lambda>1/\tau$ the majority of molecules is stretched to the linear size $R_{\rm op}\gg R_0$. The value of $R_{\rm op}$ can be much smaller than the maximal length of the molecules because of back reaction of the polymers on the flow, which suppresses velocity gradients thus preventing the polymers from maximal possible stretching.

categories: chao-dyn cond-mat nlin.CD physics.flu-dyn


title: An approach to Mel'nikov theory in celestial mechanics

abstract: Using a completely analytic procedure - based on a suitable extension of a classical method - we discuss an approach to the Poincar\'e-Mel'nikov theory, which can be conveniently applied also to the case of non-hyperbolic critical points, and even if the critical point is located at the infinity. In this paper, we concentrate our attention on the latter case, and precisely on problems described by Kepler-like potentials in one or two degrees of freedom, in the presence of general time-dependent perturbations. We show that the appearance of chaos (possibly including Arnol'd diffusion) can be proved quite easily and in a direct way, without resorting to singular coordinate transformations, such as the McGehee or blowing-up transformations. Natural examples are provided by the classical Gyld\'en problem, originally proposed in celestial mechanics, but also of interest in different fields, and by the general 3-body problem in classical mechanics.


title: Asymptotics of high order noise corrections

abstract: We consider an evolution operator for a discrete Langevin equation with a strongly hyperbolic classical dynamics and noise with finite moments. Using a perturbative expansion of the evolution operator we calculate high order corrections to its trace in the case of a quartic map and Gaussian noise. The leading contributions come from the period one orbits of the map. The asymptotic behaviour is investigated and is found to be independent up to a multiplicative constant of the distribution of noise.


title: Lyapunov exponents, entropy production and decoherence

abstract: We establish that the entropy production rate of a classically chaotic Hamiltonian system coupled to the environment settles, after a transient, to a meta-stable value given by the sum of positive generalized Lyapunov exponents. A meta-stable steady state is generated in this process. This behavior also occurs in quantum systems close to the classical limit where it leads to the restoration of quantum-classical correspondence in chaotic systems coupled to the environment.


title: Classical statistical mechanics of a few-body interacting spin model

abstract: We study the emergence of Boltzmann's law for the "single particle energy distribution" in a closed system of interacting classical spins. It is shown that for a large number of particles Boltzmann's law may occur, even if the interaction is very strong. Specific attention is paid to classical analogs of the average shape of quantum eigenstates and "local density of states", which are very important in quantum chaology. Analytical predictions are then compared with numerical data.

categories: chao-dyn cond-mat.stat-mech nlin.CD physics.class-ph quant-ph


title: Dissipative tunneling in presence of classical chaos in a mixed quantum-classical system

abstract: We consider the tunneling of a wave packet through a potential barrier which is coupled to a nonintegrable classical system and study the interplay of classical chaos and dissipation in the tunneling dynamics. We show that chaos-assisted tunneling is further enhanced by dissipation, while tunneling is suppressed by dissipation when the classical subsystem is regular.


title: Billiard Sequences and the Property of Splittability of Integrable Hamilton Systems

abstract: The paper establishes the property of splittability of billiard boundary sequences in n dimensional cube into subsequences of fractional parts. This reveals a new property of integrable and weak perturbated Hamilton systems: under a simple assumption, the boundary motion of elliptic orbits on stable KAM tori, if considering in cartesian coordinates, can be splitted into a countable set of discrete rotations. The rate of the split process, expressed in terms of some exceptional sets density, in dependence of number-theoretical characteristics of the orbits frequencies, is also examined.


title: Almost Periodic Passive Tracer Dispersion

abstract: The authors investigate the impact of external sources on the pattern formation of concentration profiles of passive tracers in a two-dimensional shear flow. By using the pullback attractor technique for the associated nonautonomous dynamical system, it is shown that a unique time-almost periodic concentration profile exists for time-almost periodic external source.


title: Sensitive boundary condition dependence of noise-sustained structure

abstract: Sensitive boundary condition dependence (BCD) is reported in a convectively unstable system with noise, where the amplitude of generated oscillatory dynamics in the downstream depends sensitively on the boundary value. This BCD is explained in terms of the manner in which the co-moving Lyapunov exponent (characterizing the convective instability) decreases from upstream to downstream. It is shown that a fractal BCD appears if the dynamics that represent the spatial change of the fixed point includes transient chaotic dynamics. By considering as an example a one-way-coupled map lattice, this theory for BCD is demonstrated.

categories: chao-dyn adap-org cond-mat.dis-nn nlin.AO nlin.CD


title: First Experimental Evidence for Chaos-Assisted Tunneling in a Microwave Annular Billiard

abstract: We report on first experimental signatures for chaos-assisted tunneling in a two-dimensional annular billiard. Measurements of microwave spectra from a superconducting cavity with high frequency resolution are combined with electromagnetic field distributions experimentally determined from a normal conducting twin cavity with high spatial resolution to resolve eigenmodes with properly identified quantum numbers. Distributions of so-called quasi-doublet splittings serve as basic observables for the tunneling between whispering gallery type modes localized to congruent, but distinct tori which are coupled weakly to irregular eigenstates associated with the chaotic region in phase space.


title: Semiclassical singularities from bifurcating orbits

abstract: We study how the singular behaviour of classical systems at bifurcations is reflected by their quantum counterpart. The semiclassical contributions of individual periodic orbits to trace formulae of Gutzwiller type are known to diverge when orbits bifurcate, a situation characteristic for systems with a mixed phase space. This singular behaviour is reflected by the quantum system in the semiclassical limit, in which the individual contributions remain valid closer to a bifurcation, while the true collective amplitude at the bifurcation increases with some inverse power of Planck's constant (with an exponent depending on the type of bifurcation). We illustrate the interplay between the two competing limits (closeness to a bifurcation and smallness of $\hbar$) numerically for a generic dynamical system, the kicked top.


title: General properties of propagation in chaotic systems

abstract: We conjecture that in one-dimensional spatially extended systems the propagation velocity of correlations coincides with a zero of the convective Lyapunov spectrum. This conjecture is successfully tested in three different contexts: (i) a Hamiltonian system (a Fermi-Pasta-Ulam chain of oscillators); (ii) a general model for spatio-temporal chaos (the complex Ginzburg-Landau equation); (iii) experimental data taken from a CO_2 laser with delayed feedback. In the last case, the convective Lyapunov exponent is determined directly from the experimental data.


title: Measurement of Wave Chaotic Eigenfunctions in the Time-Reversal Symmetry-Breaking Crossover Regime

abstract: We present experimental results on eigenfunctions of a wave chaotic system in the continuous crossover regime between time-reversal symmetric and time-reversal symmetry-broken states. The statistical properties of the eigenfunctions of a two-dimensional microwave resonator are analyzed as a function of an experimentally determined time-reversal symmetry breaking parameter. We test four theories of one-point eigenfunction statistics and introduce a new theory relating the one-point and two-point statistical properties in the crossover regime. We also find a universal correlation between the one-point and two-point statistical parameters for the crossover eigenfunctions.


title: Unstable dimension variability and synchronization of chaotic systems

abstract: An aspect of the synchronization dynamics is investigated in this work. We argue analytically and confirm numerically that the chaotic dynamics on the synchronization manifold exhibits unstable dimension variability. Unstable dimension variability is a cause of severe modeling difficulty for physical phenomena, since trajectories obtained from the mathematical model may not be related to trajectories of the actual system. We present and example of unstable dimension variability occurring in a system of two coupled chaotic maps, considering the dynamics in the synchronization manifold and its corresponding transversal direction, where a tongue-like structure is formed. The unstable dimension variability is revealed in the statistical distribution of the finite-time transversal Lyapunov exponent, having both negative and positive values.


title: Fixed points, stability and intermittency in a shell model for advection of passive scalars

abstract: We investigate the fixed points of a shell model for the turbulent advection of passive scalars introduced Jensen, Paladin and Vulpiani. The passive scalar field is driven by the velocity field of the popular GOY shell model. The scaling behavior of the static solutions is found to differ significantly from Obukhov-Corrsin scaling theta_n ~ k_n^(-1/3) which is only recovered in the limit where the diffusivity vanishes, D -> 0. From the eigenvalue spectrum we show that any perturbation in the scalar will always damp out, i.e. the eigenvalues of the scalar are negative and are decoupled from the eigenvalues of the velocity. Furthermore we estimate Lyapunov exponents and the intermittency parameters using a definition proposed by Benzi et al. The full model is as chaotic as the GOY model, measured by the maximal Lyapunov exponent, but is more intermittent.


title: A New Feature in Some Quasi-discontinuous Systems

abstract: Many systems can display a very short, rapid changing stage (quasi-discontinuous region) inside a relatively very long and slowly changing process. A quantitative definition for the "quasi-discontinuity" in these systems has been introduced. We have shown by a simplified model that extra-large Feigenbaum constants can be found inside some period-doubling cascades due to the quasi-discontinuity. As an example, this phenomenon has also been observed in Rose-Hindmash model describing neuron activities.


title: Better Nonlinear Models from Noisy Data: Attractors with Maximum Likelihood

abstract: A new approach to nonlinear modelling is presented which, by incorporating the global behaviour of the model, lifts shortcomings of both least squares and total least squares parameter estimates. Although ubiquitous in practice, a least squares approach is fundamentally flawed in that it assumes independent, normally distributed (IND) forecast errors: nonlinear models will not yield IND errors even if the noise is IND. A new cost function is obtained via the maximum likelihood principle; superior results are illustrated both for small data sets and infinitely long data streams.


title: Phase Synchronization with Type-II Intermittency in Chaotic Oscillators

abstract: We study the phase synchronization (PS) with type-II intermittency showing $\pm 2 \pi$ irregular phase jumping behavior before the PS transition occurs in a system of two coupled hyperchaotic R\"{o}ssler oscillators. The behavior is understood as a stochastic hopping of an overdamped particle in a potential which has $2 \pi$-periodic minima. We characterize it as type-II intermittency with external noise through the return map analysis. In $\epsilon_{t} < \epsilon < \epsilon_{c}$ (where $\epsilon_{t}$ is the bifurcation point of type-II intermittency and $\epsilon_{c}$ is the PS transition point in coupling strength parameter space), the average length of the time interval between two successive jumps follows $<l> \sim \exp(\mid\epsilon_{t} - \epsilon\mid^{2})$, which agrees well with the scaling law obtained from the Fokker-Planck equation.


title: Fractional Fokker-Planck Equation for Nonlinear Stochastic Differential Equations Driven by Non-Gaussian Levy Stable Noises

abstract: The Fokker-Planck equation has been very useful for studying dynamic behavior of stochastic differential equations driven by Gaussian noises. However, there are both theoretical and empirical reasons to consider similar equations driven by strongly non-Gaussian noises. In particular, they yield strongly non-Gaussian anomalous diffusion which seems to be relevant in different domains of Physics. We therefore derive in this paper a Fractional Fokker-Planck equation for the probability distribution of particles whose motion is governed by a nonlinear Langevin-type equation, which is driven by a Levy-stable noise rather than a Gaussian. We obtain in fact a general result for a Markovian forcing. We also discuss the existence and uniqueness of the solution of the Fractional Fokker-Planck equation


title: Robust and Irreversible Development in Cell Society as a General Consequence of Intra-Inter Dynamics

abstract: A dynamical systems scenario for developmental cell biology is proposed, based on numerical studies of a system with interacting units with internal dynamics and reproduction. Diversification, formation of discrete and recursive types, and rules for differentiation are found as a natural consequence of such a system. "Stem cells" that either proliferate or differentiate to different types stochastically are found to appear when intra-cellular dynamics are chaotic. Robustness of the developmental process against microscopic and macroscopic perturbations is shown to be a natural consequence of such intra-inter dynamics, while irreversibility in developmental process is discussed in terms of the gain of stability, loss of diversity and chaotic instability.


title: Reconstruction of systems with delayed feedback: (II) Application

abstract: We apply a recently proposed method for the analysis of time series from systems with delayed feedback to experimental data generated by a CO_2 laser. The method is able to estimate the delay time with an error of the order of the sampling interval, while an approach based on the peaks of either the autocorrelation function, or the time delayed mutual information would yield systematically larger values. We reconstruct rather accurately the equations of motion and, in turn, estimate the Lyapunov spectrum even for rather high dimensional attractors. By comparing models constructed for different 'embedding dimensions' with the original data, we are able to find the minimal faitfhful model. For short delays, the results of our procedure have been cross-checked using a conventional Takens time-delay embedding. For large delays, the standard analysis is inapplicable since the dynamics becomes hyperchaotic. In such a regime we provide the first experimental evidence that the Lyapunov spectrum, rescaled according to the delay time, is independent of the delay time itself. This is in full analogy with the independence of the system size found in spatially extended systems.


title: An exit-time approach to $\epsilon$-entropy

abstract: An efficient approach to the calculation of the $\epsilon$-entropy is proposed. The method is based on the idea of looking at the information content of a string of data, by analyzing the signal only at the instants when the fluctuations are larger than a certain threshold $\epsilon$, i.e., by looking at the exit-time statistics. The practical and theoretical advantages of our method with respect to the usual one are shown by the examples of a deterministic map and a self-affine stochastic process.


title: Dynamic bifurcations: hysteresis, scaling laws and feedback control

abstract: We review some properties of dynamical systems with slowly varying parameters, when a parameter is moved through a bifurcation point of the static system. Bifurcations with a single zero eigenvalue may create hysteresis cycles, whose area scales in a nontrivial way with the adiabatic parameter. Hopf bifurcations lead to the delayed appearance of oscillations. Feedback control theory motivates the study of a bifurcation with double zero eigenvalue, in which this delay is suppressed.


title: Quantum Resonances and Regularity Islands in Quantum Maps

abstract: We study analytically as well as numerically the dynamics of a quantum map near a quantum resonance of an order q. The map is embedded into a continuous unitary transformation generated by a time-independent quasi-Hamiltonian. Such a Hamiltonian generates at the very point of the resonance a local gauge transformation described the unitary unimodular group SU(q). The resonant energy growth of is attributed to the zero Liouville eigenmodes of the generator in the adjoint representation of the group while the non-zero modes yield saturating with time contribution. In a vicinity of a given resonance, the quasi-Hamiltonian is then found in the form of power expansion with respect to the detuning from the resonance. The problem is related in this way to the motion along a circle in a (q^2-1)-component inhomogeneous "magnetic" field of a quantum particle with $q$ intrinsic degrees of freedom described by the SU(q) group. This motion is in parallel with the classical phase oscillations near a non-linear resonance. The most important role is played by the resonances with the orders much smaller than the typical localization length, q << l. Such resonances master for exponentially long though finite times the motion in some domains around them. Explicit analytical solution is possible for a few lowest and strongest resonances.


title: From Generalized Synchrony to Topological Decoherence: Emergent Sets in Coupled Chaotic Systems

abstract: We consider the evolution of the unstable periodic orbit structure of coupled chaotic systems. This involves the creation of a complicated set outside of the synchronization manifold (the emergent set). We quantitatively identify a critical transition point in its development (the decoherence transition). For asymmetric systems we also describe a migration of unstable periodic orbits that is of central importance in understanding these systems. Our framework provides an experimentally measurable transition, even in situations where previously described bifurcation structures are inapplicable.


title: Universality classes for self-similarity of noiseless multi-dimensional Burgers turbulence and interface growth

abstract: The present work is devoted to the evolution of random solutions of the unforced Burgers and KPZ equations in d-dimensions in the limit of vanishing viscosity. We consider a cellular model and as initial condition assign a value for the velocity potential chosen independently within each cell. We show that the asymptotic behavior of the turbulence at large times is determined by the tail of the initial potential probability distribution function. Three classes of initial distribution leading to self-similar evolution are identified: (a) distributions with a power-law tail, (b) compactly supported potential, (c) stretched exponential tails. In class (c) we find that the mean potential (mean height of the surface) increases logarithmically with time and the 'turbulence energy' E(t) (mean square gradient of the surface) decays as 1/t times a logarithmic correction. In classes (a) and (b) we find that the changes in the mean potential and energy have a power-law time dependence. In class (c) the roughness of the surface, measured by its mean--square gradient, may either decrease or increase with time. We discuss also the influence of finite viscosity and long range correlation on the late stage evolution of the Burgers turbulence


title: Algebraic Entropy of Birational Maps with Invariant Curves

abstract: The relevance of the algebraic entropy in the study of birational discrete time dynamical systems highlights the need to relate it to other characteristics of these systems. In this letter, two complementary proofs are given that the foliation of the space by invariant curves implies that the algebraic entropy is zero.

categories: chao-dyn math-ph math.AG math.DS math.MP nlin.CD


title: Quantum-classical correspondence and nonclassical states generation in dissipative quantum optical systems

abstract: We develop a semiclassical method for the determination of the nonlinear dynamics of dissipative quantum optical systems in the limit of large number of photons N, based on the 1/N-expansion and the quantum-classical correspondence. The method has been used to tackle two problems: to study the dynamics of nonclassical state generation in higher-order anharmonic dissipative oscillators and to establish the difference between the quantum and classical dynamics of the second-harmonic generation in a self-pulsing regime. In addressing the first problem, we have obtained an explicit time dependence of the squeezing and the Fano factor for an arbitrary degree of anharmonism in the short-time approximation. For the second problem, we have established analytically a characteristic time scale when the quantum dynamics differs insignificantly from the classical one.

categories: chao-dyn cond-mat.stat-mech nlin.CD quant-ph


title: Bifurcation analysis of the plane sheet pinch

abstract: A numerical bifurcation analysis of the electrically driven plane sheet pinch is presented. The electrical conductivity varies across the sheet such as to allow instability of the quiescent basic state at some critical Hartmann number. The most unstable perturbation is the two-dimensional tearing mode. Restricting the whole problem to two spatial dimensions, this mode is followed up to a time-asymptotic steady state, which proves to be sensitive to three-dimensional perturbations even close to the point where the primary instability sets in. A comprehensive three-dimensional stability analysis of the two-dimensional steady tearing-mode state is performed by varying parameters of the sheet pinch. The instability with respect to three-dimensional perturbations is suppressed by a sufficiently strong magnetic field in the invariant direction of the equilibrium. For a special choice of the system parameters, the unstably perturbed state is followed up in its nonlinear evolution and is found to approach a three-dimensional steady state.


title: The effect of forcing on the spatial structure and spectra of chaotically advected passive scalars

abstract: The stationary distribution of passive tracers chaotically advected by a two-dimensional large-scale flow is investigated. The tracer field is force by resetting the value of the tracer in certain localised regions. This problem is mathematically equivalent to advection in open flows and results in a fractal tracer structure. The spectral exponent of the tracer field is different from that for a passive tracer with the usual additive forcing (the so called Batchelor spectrum) and is related to the fractal dimension of the set of points that have never visited the forcing regions. We illustrate this behaviour by considering a time-periodic flow whose effect is equivalent to a simple two-dimensional area-preserving map. We also show that similar structure in the tracer field is found when the flow is aperiodic in time.


title: Generalized Stable Multivariate Distribution and Anisotropic Dilations

abstract: After having closely re-examined the notion of a L\'evy's stable vector, it is shown that the notion of a stable multivariate distribution is more general than previously defined. Indeed, a more intrinsic vector definition is obtained with the help of non isotropic dilations and a related notion of generalized scale. In this framework, the components of a stable vector may not only have distinct Levy's stability indices $\alpha$'s, but the latter may depend on its norm. Indeed, we demonstrate that the Levy's stability index of a vector rather correspond to a linear application than to a scalar, and we show that the former should satisfy a simple spectral property.


title: Differential light scattering: probing the sonoluminescence collapse

abstract: We have developed a light scattering technique based on differential measurement and polarization (differential light scattering, DLS) capable in principle of retrieving timing information with picosecond resolution without the need for fast electronics. DLS was applied to sonoluminescence, duplicating known results (sharp turnaround, self-similar collapse); the resolution was limited by intensity noise to about 0.5 ns. Preliminary evidence indicates a smooth turnaround on a time scale of a few hundred picoseconds, and suggests the existence of subnanosecond features within a few nanoseconds of the turnaround.

categories: chao-dyn nlin.CD physics.flu-dyn physics.optics


title: The Nose-hoover thermostated Lorentz gas

abstract: We apply the Nose-Hoover thermostat and three variations of it, which control different combinations of velocity moments, to the periodic Lorentz gas. Switching on an external electric field leads to nonequilibrium steady states for the four models with a constant average kinetic energy of the moving particle. We study the probability density, the conductivity and the attractor in nonequilibrium and compare the results to the Gaussian thermostated Lorentz gas and to the Lorentz gas as thermostated by deterministic scattering.


title: Characteristic Relations of Type-I Intermittency in the Presence of Noise

abstract: Near the point of tangent bifurcation, the scaling properties of the laminar length of type-I intermittency are investigated in the presence of noise. Based on analytic and numerical studies, we show that the scaling relation of the laminar length is dramatically deformed from $\frac{1}{\sqrt{\epsilon}}$ for $\epsilon >0$ to $\exp\{\frac{1}{D}|\epsilon|^{3/2}\}$ for $\epsilon<0$ as $\epsilon$ passes the bifurcation point $(\epsilon=0)$. The results explain why two coupled R\"ossler oscillators exhibit deformation of the scaling relation of the synchronous length in the nearly synchronous regime.


title: Extraction of information about periodic orbits from scattering functions

abstract: As a contribution to the inverse scattering problem for classical chaotic systems, we show that one can select sequences of intervals of continuity, each of which yields the information about period, eigenvalue and symmetry of one unstable periodic orbit.


title: Turbulent 4-wave Interaction of Two Type of Waves

abstract: We consider turbulent 4-wave interaction of two types of waves: acoustic waves (dispersion $\omega = k$) and electromagnetic-type waves (dispersion $\Omega^2 = m^2 + p^2$). For large wave vectors ($ k \gg m$), when the dispersion of EM-type waves becomes quasiacoustic, the stationary spectra are obtained following a standard Zakharov approach. In the small wave number region ($ k, p \ll m$) we derive nonlinear differential Kompaneets-type kinetic equation for arbitrary distributions of the interacting fields and find when this equation has Kolmogorov-type power law solutions.


title: The boundary integral method for magnetic billiards

abstract: We introduce a boundary integral method for two-dimensional quantum billiards subjected to a constant magnetic field. It allows to calculate spectra and wave functions, in particular at strong fields and semiclassical values of the magnetic length. The method is presented for interior and exterior problems with general boundary conditions. We explain why the magnetic analogues of the field-free single and double layer equations exhibit an infinity of spurious solutions and how these can be eliminated at the expense of dealing with (hyper-)singular operators. The high efficiency of the method is demonstrated by numerical calculations in the extreme semiclassical regime.


title: Ultrasound scattering and the study of vortex correlations in disordered flows

abstract: In an idealized way, some turbulent flows can be pictured by assemblies of many vortices characterized by a set of particle distribution functions. Ultrasound provide an useful, nonintrusive, tool to study the spatial structure of vorticity in flows. This is analogous to the use of elastic neutron scattering to determine liquid structure. We express the dispersion relation, as well as the scattering cross section, of sound waves propagating in a ``liquid'' of identical vortices as a function of vortex pair correlation functions. In two dimensions, formal analogies with ionic liquids are pointed out.


title: Critical behavior of vorticity in two-dimensional turbulence

abstract: We point out some similitudes between the statistics of high Reynolds number turbulence and critical phenomena. An analogy is developed for two-dimensional decaying flows, in particular by studying the scaling properties of the two-point vorticity correlation function within a simple phenomenological framework. The inverse of the Reynolds number is the analogue of the small parameter that separates the system from criticality. It is possible to introduce a set of three critical exponents; for the correlation length, the autocorrelation function and a so-called susceptibility, respectively. The exponents corresponding to the well-known enstrophy cascade theory of Kraichnan and Batchelor are, remarkably, the same as the Gaussian approximation exponents for spin models. The limitations of the analogy, in particular the lack of universal scaling functions, are also discussed.


title: Anomalous scaling in a shell model of helical turbulence

abstract: In a helical flow there is a subrange of the inertial range in which there is a cascade of both energy and helicity. In this range the scaling exponents associated with the cascade of helicity can be defined. These scaling exponents are calculated from a simulation of the GOY shell model. The scaling exponents for even moments are associated with the scaling of the symmetric part of the probability density functions while the odd moments are associated with the anti-symmetric part of the probability density functions.


title: Transition to Chaotic Phase Synchronization through Random Phase Jumps

abstract: Phase synchronization is shown to occur between opposite cells of a ring consisting of chaotic Lorenz oscillators coupled unidirectionally through driving. As the coupling strength is diminished, full phase synchronization cannot be achieved due to random generation of phase jumps. The brownian dynamics underlying this process is studied in terms of a stochastic diffusion model of a particle in a one-dimensional medium.


title: Clustering of inertial particles in turbulent flows

abstract: We consider inertial particles suspended in an incompressible turbulent flow. Due to inertia of particles, their velocity field acquires small compressible component. Its presence leads to a new qualitative effect --- possibility of clustering. We show that this effect is significant for heavy particles, leading to strong fluctuations of the concentration.


title: Thermodynamic Cross-Effects from Dynamical Systems

abstract: We give a thermodynamically consistent description of simultaneous heat and particle transport, as well as of the associated cross-effects, in the framework of a chaotic dynamical system, a generalized multibaker map. Besides the density, a second field with appropriate source terms is included in order to mimic, after coarse graining, a spatial temperature distribution and its time evolution. A new expression is derived for the irreversible entropy production in a steady state, as the average of the growth rate of the relative density, a unique combination of the two fields.


title: Experimental improvement of chaotic synchronization due to multiplicative time-correlated Gaussian noise

abstract: The effect of time-correlated zero-mean Gaussian noise on chaotic synchronization is analyzed experimentally in small-size arrays of Chua's circuits. Depending on the correlation time, an improvement of the synchronization is found for different values of the noise amplitude and coupling diffusion between circuits.


title: Differentiability Implies Continuity in Neuronal Dynamics

abstract: Recent work has identified nonlinear deterministic structure in neuronal dynamics using periodic orbit theory. Troublesome in this work were the significant periods of time where no periodic orbits were extracted - "dynamically dark" regions. Tests for periodic orbit structure typically require that the underlying dynamics are differentiable. Since continuity of a mathematical function is a necessary but insufficient condition for differentiability, regions of observed differentiability should be fully contained within regions of continuity. We here verify that this fundamental mathematical principle is reflected in observations from mammalian neuronal activity. First, we introduce a null Jacobian transformation to verify the observation of differentiable dynamics when periodic orbits are extracted. Second, we show that a less restrictive test for deterministic structure requiring only continuity demonstrates widespread nonlinear deterministic structure only partially appreciated with previous approaches.


title: Decimation and Harmonic Inversion of Periodic Orbit Signals

abstract: We present and compare three generically applicable signal processing methods for periodic orbit quantization via harmonic inversion of semiclassical recurrence functions. In a first step of each method, a band-limited decimated periodic orbit signal is obtained by analytical frequency windowing of the periodic orbit sum. In a second step, the frequencies and amplitudes of the decimated signal are determined by either Decimated Linear Predictor, Decimated Pade Approximant, or Decimated Signal Diagonalization. These techniques, which would have been numerically unstable without the windowing, provide numerically more accurate semiclassical spectra than does the filter-diagonalization method.


title: Chaotic Transport and Current Reversal in Deterministic Ratchets

abstract: We address the problem of the classical deterministic dynamics of a particle in a periodic asymmetric potential of the ratchet type. We take into account the inertial term in order to understand the role of the chaotic dynamics in the transport properties. By a comparison between the bifurcation diagram and the current, we identify the origin of the current reversal as a bifurcation from a chaotic to a periodic regime. Close to this bifurcation, we observed trajectories revealing intermittent chaos and anomalous deterministic diffusion.


title: Longitudinal Structure Functions in Decaying and Forced Turbulence

abstract: In order to reliably compute the longitudinal structure functions in decaying and forced turbulence, local isotropy is examined with the aid of the isotropic expression of the incompressible conditions for the second and third order structure functions. Furthermore, the Karman-Howarth-Kolmogorov relation is investigated to examine the effects of external forcing and temporally decreasing of the second order structure function. On the basis of these investigations, the scaling range and exponents $\zeta_n$ of the longitudinal structure functions are determined for decaying and forced turbulence with the aid of the extended-self-similarity (ESS) method. We find that $\zeta_n$'s are smaller, for $n \geq 4$, in decaying turbulence than in forced turbulence. The reasons for this discrepancy are discussed. Analysis of the local slopes of the structure functions is used to justify the ESS method.


title: A MultiBaker Map for Thermodynamic Cross-Effects in Dynamical Systems

abstract: A consistent description of simultaneous heat and particle transport, including cross effects, and the associated entropy balance is given in the framework of a deterministic dynamical system. This is achieved by a multibaker map where, besides the phase-space density of the multibaker, a second field with appropriate source terms is included in order to mimic a spatial temperature distribution and its time evolution. Conditions are given to ensure consistency in an appropriately defined continuum limit with the thermodynamic entropy balance. They leave as the only free parameter of the model the entropy flux let directly into a surroundings. If it vanishes in the bulk, the transport properties of the model are described by the thermodynamic transport equations. Another choice leads to a uniform temperature distribution. It represents transport problems treated by means of a thermostatting algorithm, similar to the one considered in non-equilibrium molecular dynamics.


title: Semiclassical properties of eigenfunctions and occupation number distribution for a model of two interacting particles

abstract: Quantum-classical correspondence for the shape of eigenfunctions, local spectral density of states and occupation number distribution is studied in a chaotic model of two coupled quartic oscillators. In particular, it is shown that both classical quantities and quantum spectra determine global properties of occupation numbers and inverse participation ratio.


title: Exact Renormalization Scheme for Quantum Anosov Maps

abstract: An exact renormalization scheme is introduced for quantum Anosov maps (QAMs) on a torus for general boundary conditions (BCs), whose number is always finite. Given a QAM $\hat{U}$ with $k$ BCs and Planck's constant $\hbar =2\pi /p$ ($p$ integer), its $n$th renormalization iterate $\hat{U}^{(n)}={\cal R}^{n}(\hat{U})$ is associated with $k$ BCs for all $n$ and with a Planck's constant $\hbar ^{(n)}=\hbar /k^{n}$. It is shown that the quasienergy eigenvalue problem for $\hat{U}^{(n)}$ for {\em all} $k$ BCs is equivalent to that for $\hat{U}^{(n+1)}$ at some {\em fixed} BCs, corresponding, for $n>0$, to either strict {\em periodicity} for $kp$ even or {\em antiperiodicity} for $kp$ odd. The quantum cat maps are, in general, fixed points of either ${\cal R}$ or ${\cal R}^{2}$. The Hannay-Berry results turn out then to be significant also for general BCs.


title: Large Global Coupled Maps with Multiple Attractors

abstract: A system of N unidimensional global coupled maps (GCM), which support multiattractors is studied. We analize the phase diagram and some special features of the transitions (volume ratios and characteristic exponents), by controlling the number of elements of the initial partition that are in each basin of attraction. It was found important difference with widely known coupled systems with a single attractor.


title: A simple denoising algorithm using wavelet transform

abstract: Submission withdrawn because the authors erroneously submitted a revised version as a new submission, see nlin.CD/0002028.


title: Fourier Path Integral Monte Carlo Method for the Calculation of the Microcanonical Density of States

abstract: Using a Hubbard-Stratonovich transformation coupled with Fourier path integral methods, expressions are derived for the numerical evaluation of the microcanonical density of states for quantum particles obeying Boltzmann statistics. A numerical algorithmis suggested to evaluate the quantum density of states and illustrated on a one-dimensional model system.


title: Monte Carlo Studies of the Orientational Order-Disorder Phase Transition in Solid Ammonium Chloride

abstract: Monte Carlo methods are used to study the phase transition in ammonium chloride from the orientationally ordered $\delta$ phase to the orientationally disordered $\gamma$ phase. An effective pair potential is used to model the interaction between ions. Thermodynamic properties are computed in the canonical and isothermal-isobaric ensembles. Each ammonium ion is treated as a rigidly rotating body and the lattice is fixed in the low-temperature CsCl geometry. A simple extension of the Metropolis Monte Carlo method is used to overcome quasiergodicity in the rotational sampling. In the constant-$NVT$ calculations the lattice is held rigid; in the constant-$NpT$ calculations the lattice parameter is allowed to fluctuate. In both ensembles the order parameter rapidly falls to zero in the range (200 - 250)K, suggesting that the model disorders at a temperature in fair agreement with the experimental disordering temperature (243K). Peaks in the heat capacity and thermal expansivity curves are also found in the same temperature range.


title: Stability of the Hydrogen and Hydrogen-like Molecules

abstract: We present a simple proof of the stability of the hydrogen molecule $(M^+M^+m^-m^-)$. It does not rely on the proton-to-electron mass ratio $M/m$ being very large, and actually holds for arbitrary values of $M/m$. Some asymmetric molecules of the type $(m_1^+m_2^+m_3^-m_4^-)$ are also stable. Possible applications to molecules containing antiparticles and to exotic hadrons in the quark model are briefly outlined. Revtex Version 3.0, 4 Figures available on request by Fax or Mail.


title: Comparison of electron density functional models

abstract: Presented here are calculations of the distortion of the density of an electron gas due to the electrostatic field of a proton. Several models based upon the local density approximation (LDA) of density functional theory [linear response theory, Kohn-Sham (KS), optimized Thomas-Fermi theory (OTF), and OTF plus perturbation corrections] are compared with one another. These models, in turn, are compared with available results of quantum Monte Carlo (QMC) calculations for the same system. Comparison of the KS results with the OTF results shows a very reasonable agreement that seems to be progressively improvable. This provides encouragement for the application of the OTF model to condensed phase systems. The QMC calculations of the density do not agree well with the density functional results. The reasons for this poor agreement are not clear. This particular system is expected to be an ideal one for application of the LDA and, thus, the poor agreement is of fundamental importance. The lack of detail presented in the available QMC results leads us to conclude that the QMC calculations should be attempted again. in print: Molec. Phys. 82, 245-261(1994)


title: Density Functional Calculations On First-Row Transition Metals

abstract: The excitation energies and ionization potentials of the atoms in the first transition series are notoriously difficult to compute accurately. Errors in calculated excitation energies can range from 1--4 eV at the Hartree-Fock level, and errors as high as 1.5eV are encountered for ionization energies. In the current work we present and discuss the results of a systematic study of the first transition series using a spin-restricted Kohn-Sham density-functional method with the gradient-corrected functionals of Becke and Lee, Yang and Parr. Ionization energies are observed to be in good agreement with experiment, with a mean absolute error of approximately 0.15eV; these results are comparable to the most accurate calculations to date, the Quadratic Configuration Interaction (QCISD(T)) calculations of Raghavachari and Trucks. Excitation energies are calculated with a mean error of approximately 0.5eV, compared with $\sim 1\mbox{eV}$ for the local density approximation and 0.1eV for QCISD(T). These gradient-corrected functionals appear to offer an attractive compromise between accuracy and computational effort.


title: The Aqueous Solvation of Water: A Comparison of Continuum Methods with Molecular Dynamics

abstract: The calculation of the solvation properties of a single water molecule in liquid water is carried out in two ways. In the first, the water molecule is placed in a cavity and the solvent is treated as a dielectric continuum. This model is analyzed by numerically solving the Poisson equation using the DelPhi program. The resulting solvation properties depend sensitively on the shape and size of the cavity. In the second method, the solvent and solute molecules are treated explicitly in molecular dynamics simulations using Ewald boundary conditions. We find a 2 kcal/mole difference in solvation free energies predicted by these two methods when standard cavity radii are used. In addition, dielectric continuum theory assumes that the solvent reacts solely by realigning its electric moments linearly with the strength of the solute's electric field; the results of the molecular simulation show important non-linear effects. Non-linear solvent effects are generally of two types: dielectric saturation, due to solvent-solute hydrogen bonds, and electrostriction, a decrease in the solute cavity due to an increased electrostatic interaction. We find very good agreement between the two methods if the radii defining the solute cavity used in the continuum theory is decreased with the solute charges,


title: Magic Numbers of Silicon Clusters

abstract: A structural model for intermediate sized silicon clusters is proposed that is able to generate unique structures without any dangling bonds. This structural model consists of bulk-like core of five atoms surrounded by fullerene-like surface. Reconstruction of the ideal fullerene geometry results in the formation of crown atoms surrounded by $\pi$-bonded dimer pairs. This model yields unique structures for \Si{33}, \Si{39}, and \Si{45} clusters without any dangling bonds and hence explains why these clusters are least reactive towards chemisorption of ammonia, methanol, ethylene, and water. This model is also consistent with the experimental finding that silicon clusters undergo a transition from prolate to spherical shapes at \Si{27}. Finally, reagent specific chemisorption reactivities observed experimentally is explained based on the electronic structures of the reagents.


title: Vibrational States of the Hydrogen Isotopes on Pd(111)

abstract: The ground and excited vibrational states for the three hydrogen isotopes on the Pd(111) surface have been calculated. Notable features of these states are the high degree of anharmonicity, which is most prominently seen in the weak isotopic dependence of the parallel vibrational transition, and the narrow bandwidths of these states, which imply that atomic hydrogen is localized on a particular surface site on time scales of 100 picoseconds or more. Experiments to resolve ambiguities concerning the present system are suggested.


title: Li$_5$ as a pseudorotating planar cluster

abstract: A pseudorotating state of Li$_5$ observed in a recent EPR experiment is investigated using the local density functional method. The calculated isotropic spin population indicates that a planar $C_{2v}$ structure is more consistent with the experimental result than a suggested trigonal bipyramid. A new pseudorotating state based on the planar structure is proposed.


title: Ion Pair Potentials-of-Mean-Force in Water

abstract: Recent molecular simulation and integral equation results alkali-halide ion pair potentials-of-mean-force in water are discussed. Dielectric model calculations are implemented to check that these models produce that characteristic structure of contact and solvent-separated minima for oppositely charged ions in water under physiological thermodynamic conditions. Comparison of the dielectric model results with the most current molecular level information indicates that the dielectric model does not, however, provide an accurate description of these potentials-of-mean-force. We note that linear dielectric models correspond to modelistic implementations of second-order thermodynamic perturbation theory for the excess chemical potential of a distinguished solute molecule. Therefore, the molecular theory corresponding to the dielectric models is second-order thermodynamic perturbation theory for that excess chemical potential. The second-order, or fluctuation, term raises a technical computational issue of treatment of long-ranged interactions similar to the one which arises in calculation of the dielectric constant of the solvent. It is contended that the most important step for further development of dielectric models would be a separate assessment of the first-order perturbative term (equivalently the {\it potential at zero charge} ) which vanishes in the dielectric models but is generally nonzero. Parameterization of radii and molecular volumes should then be based of the second-order perturbative term alone. Illustrative initial calculations are presented and discussed.


title: Tests of Dielectric Model Descriptions of Chemical Charge Displacements in Water

abstract: A dielectric model of electrostatic solvation is applied to describe potentials of mean force in water along reaction paths for: a) formation of a sodium chloride ion pair; b) the symmetric SN2 exchange of chloride in methylchloride; and c) nucleophilic attack of formaldehyde by hydroxide anion. For these cases simulation and XRISM results are available for comparison. The accuracy of model predictions varies from spectacular to mediocre. It is argued that: a) dielectric models are physical models, even though simplistic and empirical; b) their successes suggest that second-order perturbation theory is a physically sound description of free energies of electrostatic solvation; and c) the most serious deficiency of the dielectric models lies in the definition of cavity volumes. Second-order perturbation theory should therefore be used to refine the dielectric models. These dielectric models make no attempt to assess the role of packing effects but for solvation of classical electrostatic interactions the dielectric models sometimes perform as well as the more detailed XRISM theory.


title: Quantum Annealing: A New Method for Minimizing Multidimensional Functions

abstract: Quantum annealing is a new method for finding extrema of multidimensional functions. Based on an extension of classical, simulated annealing, this approach appears robust with respect to avoiding local minima. Further, unlike some of its predecessors, it does not require an approximation to a wavefunction. In this paper, we apply the technique to the problem of finding the lowest energy configurations of Lennard-Jones clusters of up to 19 particles (roughly 10$^5$ local minima). This early success suggests that this method may complement the widely implemented technique of simulated annealing.


title: Chemical Reactions of Silicon Clusters

abstract: Smalley and co-workers discovered that chemisorption reactivities of silicon clusters vary over three orders of magnitude as a function of cluster size. In particular, they found that \Si{33}, \Si{39}, and \Si{45} clusters are least reactive towards various reagents compared to their immediate neighbors in size. We explain these observations based on our stuffed fullerene model. This structural model consists of bulk-like core of five atoms surrounded by fullerene-like surface. Reconstruction of the ideal fullerene geometry gives rise to four-fold coordinated crown atoms and $\pi$-bonded dimer pairs. This model yields unique structures for \Si{33}, \Si{39}, and \Si{45} clusters without any dangling bonds and thus explains their lowest reactivity towards chemisorption of closed shell reagents. This model is also consistent with the experimental finding of Jarrold and Constant that silicon clusters undergo a transition from prolate to spherical shapes at \Si{27}. We justify our model based on an in depth analysis of the differences between carbon and silicon chemistry and bonding characteristics. Using our model, we further explain why dissociative chemisorption occurs on bulk surfaces while molecular chemisorption occurs on cluster surfaces. We also explain reagent specific chemisorption reactivities observed experimentally based on the electronic structures of the reagents. Finally, experiments on \Si{x}X$_y$ (X = B, Al, Ga, P, As, AlP, GaAs) are suggested as a means of verifying the proposed model. We predict that \Si{x}(AlP)$_y$ and \Si{x}(GaAs)$_y$ $(x = 25, 31, 37; y = 4)$ clusters will be highly inert and it may be possible to prepare macroscopic samples of these alloy clusters through high temperature reactions.


title: Decisions, Decisions: Noise and its Effects on Integral Monte Carlo Algorithms

abstract: In the present paper we examine the effects of noise on Monte Carlo algorithms, a problem raised previously by Kennedy and Kuti (Phys. Rev. Lett. {\bf 54}, 2473 (1985)). We show that the effects of introducing unbiased noise into the acceptance/rejection phase of the conventional Metropolis approach are surprisingly modest, and, to a significant degree, largely controllable. We present model condensed phase numerical applications to support these conclusions.


title: Magic Numbers for Classical Lennard-Jones Cluster Heat Capacities

abstract: Heat capacity curves as functions of temperature for classical atomic clusters bound by pairwise Lennard-Jones potentials were calculated for aggregate sizes from 4 to 24 using Monte Carlo methods. J-walking (or jump-walking) was used to overcome convergence difficulties due to quasi-ergodicity in the solid-liquid transition region. The heat capacity curves were found to differ markedly and nonmonotonically as functions of cluster size. Curves for N = 4, 5 and 8 consisted of a smooth, featureless, monotonic increase throughout the transition region, while curves for N = 7 and 15-17 showed a distinct shoulder in this region; the remaining clusters had distinguishable transition heat capacity peaks. The size and location of these peaks exhibited "magic number" behavior, with the most pronounced peaks occurring for magic number sizes of N = 13, 19 and 23. A comparison of the heat capacities with other cluster properties in the solid-liquid transition region that have been reported in the literature indicates partial support for the view that, for some clusters, the solid-liquid transition region is a coexistence region demarcated by relatively sharp, but separate, melting and freezing temperatures; some discrepancies, however, remain unresolved.


title: Dynamical Fluctuating Charge Force Fields: Application to Liquid Water

abstract: A new molecular dynamics model in which the point charges on atomic sites are allowed to fluctuate in response to the environment is developed and applied to water. The idea for treating charges as variables is based on the concept of electronegativity equalization according to which: (a) The electronegativity of an atomic site is dependent on the atom's type and charge and is perturbed by the electrostatic potential it experiences from its neighbors and (b) Charge is transferred between atomic sites in such a way that electronegativities are equalized. The charges are treated as dynamical variables using an extended Lagrangian method in which the charges are given a fictitious mass, velocities and kinetic energy and then propagated according to Newtonian mechanics along with the atomic degrees of freedom. Models for water with fluctuating charges are developed using the geometries of two common fixed-charge water potentials: the simple point charge (SPC) and the 4-point transferable intermolecular potential (TIP4P). Both fluctuating charge models give accurate predictions for gas-phase and liquid state properties, including radial distribution functions, the dielectric constant, and the diffusion constant. The method does not introduce any new intermolecular interactions beyond those already present in the fixed charge models and increases the computer time by only a factor of 1.1, making this method tractable for large systems.


title: Many-body trial wave functions for atomic systems and ground states of small noble gas clusters

abstract: Clusters of sizes ranging from two to five are studied by variational quantum Monte Carlo techniques. The clusters consist of Ar, Ne and hypothetical lighter (``$1 \over 2$-Ne") atoms. A general form of trial function is developed for which the variational bias is considerably smaller than the statistical error of currently available diffusion Monte Carlo estimates. The trial functions are designed by a careful analysis of long- and short-range behavior as a function of inter-atomic distance; at intermediate distances, on the order of the average nearest neighbor distance, the trial functions are constructed to have considerable variational freedom. A systematic study of the relative importance of $n$-body contributions to the quality of the optimized trial wave function is made with $2\le n \le 5$. Algebraic invariants are employed to deal efficiently with the many-body interactions.


title: The Construction of Double-Ended Classical Trajectories

abstract: In the present paper we describe relaxation methods for constructing double-ended classical trajectories. We illustrate our approach with an application to a model anharmonic system, the Henon-Heiles problem. Trajectories for this model exhibit a number of interesting energy-time relationships that appear to be of general use in characterizing the dynamics.


title: Optical response of C60 and C70 fullerenes: Exciton and lattice fluctuation effects

abstract: Molecular exciton effects in the neutral and maximally doped C60 (C70) are considered using a tight binding model with long-range Coulomb interactions and bond disorder. By comparing calculated and observed optical spectra, we conclude that relevant Coulomb parameters for the doped cases are about half of those of the neutral systems. The broadening of absorption peaks is well simulated by the bond disorder model.

categories: chem-ph cond-mat hep-lat physics.chem-ph


title: Theoretical calculation of the water ion product K_W

abstract: A dielectric solvation model is applied to the prediction of the equilibrium ionization of liquid water over a wide range of density and temperature with the objective of calibrating that model for the study of ionization in water of organic acids, {\it e.g.\/}, proteins and nucleic acids. The model includes an approximate description of the polarizability of the dissociating water molecule. The calculated pK$_W$ are very sensitive to the value of the radii that parameterize the model. The radii required for the spherical molecular volumes of the water molecule in order to fit the experimental ion product are presented and discussed. These radii are larger than those commonly used. They decrease with increasing density as would be guessed but the rate of decrease is slight. They increase with increasing temperature, a variation opposite to what would be guessed if radii were strictly viewed as a distance of closest approach. The molecular theoretical principles that might provide an explanation of the thermodynamic state dependence of these radii are discussed.


title: Methods for molecular dynamics with nonadiabatic transitions

abstract: We show how the dynamically nonlocal formulation of classical nuclear motion in the presence of quantal electronic transitions presented many years ago by Pechukas can be localized in time using time dependent perturbation theory to give an impulsive force which acts when trajectories hop between electronic surfaces. The action of this impulsive force is completely equivalent to adjusting the nuclear velocities in the direction of the nonadiabatic coupling vector so as to conserve energy, a procedure which is widely used in surface hopping trajectory methods. This is the first time the precise connection between these two formulations of the nonadiabatic dynamics problem has been considered. We also demonstrate that the stationary phase approximation to the reduced propagator at the heart of Pechukas' theory is not unitary due to its neglect of nonstationary paths. As such mixed quantum-classical evolution schemes based on this approximation are not norm conserving and in general must fail to give the correct branching between different competing electronic states. Tully's phase coherent, fewest switches branching algorithm is guaranteed to conserve the norm. The branching between different alternatives predicted by this approach, however, may be inaccurate, due to use of the approximate local dynamics. We explore the relative merits of these different approximations using Tully's 1D two state example scattering problems for which numerically exact results are easily obtained.


title: C20: Fulleren, Bowl or Ring? New Results from Coupled-Cluster Calculations

abstract: Contrary to recent experimental evidence suggesting that the monocyclic ring is the most stable 20-atom carbon species, highly accurate calculations convincingly predict that the smallest fullerene, the dodecahedron C$_{20}$, has the lowest energy. A related corannulene-like bowl is nearly degenerate in energy to the fullerene. Thermodynamic considerations suggest that at formation temperatures of around 700 K the bowl should be the dominant species. The recent application of gradient corrections to LDA which supported the ring structure is qualitatively in error. (RK-94-02)


title: Locating transition states using double-ended classical trajectories

abstract: In this paper we present a method for locating transition states and higher-order saddles on potential energy surfaces using double-ended classical trajectories. We then apply this method to 7- and 8-atom Lennard-Jones clusters, finding one previously unreported transition state for the 7-atom cluster and two for the 8-atom cluster.


title: On the Mechanism of the Primary Charge Separation in Bacterial Photosynthesis

abstract: We present a detailed analysis of the mechanism of the primary charge separation process in bacterial photosynthesis using real-time path integrals. Direct computer simulations as well as an approximate analytical theory have been employed to map out the dynamics of the charge separation process in many regions of the parameter space relevant to bacterial photosynthesis. Two distinct parameter regions, one characteristic of sequential transfer and the other characteristic of superexchange, have been found to yield charge separation dynamics in agreement with experiments. Nonadiabatic theory provides accurate rate estimates for low-lying and very high-lying bacteriochlorophyll state energies, but it breaks down in between these two regimes.


title: Transferable Tight-Binding Potential for Hydrocarbons

abstract: A transferable tight-binding potential has been constructed for heteroatomic systems containing carbon and hydrogen. The electronic degree of freedom is treated explicitly in this potential using a small set of transferable parameters which has been fitted to small hydrocarbons and radicals. Transferability to other higher hydrocarbons were tested by comparison with ab initio calculations and experimental data. The potential can correctly reproduce changes in the electronic configuration as a function of the local bonding geometry around each carbon atom. This type of potential is well suited for computer simulations of covalently bonded systems in both gas-phase and condensed-phase systems.


title: Many-Body effects and resonances in universal quantum sticking of cold atoms to surfaces

abstract: The role of shape resonances and many-body effects on universal quantum sticking of ultra cold atoms onto solid surfaces is examined analytically and computationally using an exactly solvable representation of the Dyson equation. We derive the self-energy renormalization of the the transition amplitude between an ultra cold scattering atom and the bound states on the surface in order to elucidate the role of virtual phonon exchanges in the limiting behavior of the sticking probability. We demonstrate that, to first order in the interactions for finite ranged atom-surface potentials, virtual phonons can only rescale the strength of the atom-surface coupling and do not rescale the range of the coupling. Thus, universal sticking behaviour at ultra-low energies is to be expected for all finite ranged potentials. We demonstrate that the onset of the universal sticking behavior depends greatly on the position of the shape resonance of the renormalized potential and for sufficiently low energy shape resonances, deviations from the universal $s(E)\propto\sqrt{E}$ can occur near these energies. We believe that this accounts for many of the low energy sticking trends observed in the scattering of sub-millikelvin H atoms from superfluid $^4$He films.


title: Packing Transitions in Nanosized Li Clusters

abstract: Packing transitions in the lowest energy structures of Li clusters as a function of size have been identified via simulated annealing. For N>21, the large $p$ character of Li leads to unexpected ionic structures. At N~25, a packing pattern based on interpenetrating 13-atom icosahedra and similar to that of Na and K appears. This pattern persists until at N=55, where another transition to a structure based on a Mackay icosahedron occurs. For clusters of size 55 and 147, the optimized FCC structure representative of the bulk is still slightly higher in energy than the optimal MIC. (RK-94-03)


title: Dimensional Perturbation Theory on the Connection Machine

abstract: A recently developed linear algebraic method for the computation of perturbation expansion coefficients to large order is applied to the problem of a hydrogenic atom in a magnetic field. We take as the zeroth order approximation the $D \rightarrow \infty$ limit, where $D$ is the number of spatial dimensions. In this pseudoclassical limit, the wavefunction is localized at the minimum of an effective potential surface. A perturbation expansion, corresponding to harmonic oscillations about this minimum and higher order anharmonic correction terms, is then developed in inverse powers of $(D-1)$ about this limit, to 30th order. To demonstrate the implicit parallelism of this method, which is crucial if it is to be successfully applied to problems with many degrees of freedom, we describe and analyze a particular implementation on massively parallel Connection Machine systems (CM-2 and CM-5). After presenting performance results, we conclude with a discussion of the prospects for extending this method to larger systems.


title: Theoretical Study of Time-Resolved Fluorescence Anisotropy from Coupled Chromophore Pairs

abstract: Calculations of time-resolved fluorescence anisotropy from a pair of chromophores coupled by an excitation transfer interaction are presented. For the purpose of investigating the effects of nuclear motion on the energy transfer and anisotropy, an illustrative model is developed that provides each chromophore with a single intramolecular vibrational mode. Account is taken of non-instantaneous excitation and time- and frequency-resolved detection. Effects of excitation pulse duration, detection window duration and frequency resolution, and excitation transfer coupling strength on the time-resolved anisotropy are examined in detail. Effects of vibrational relaxation and dephasing are also examined using a simplified Redfield description of the effects of coupling to a thermal bath.


title: Comparative Study of Multicanonical and Simulated Annealing Algorithms in the Protein Folding Problem

abstract: We compare a few variants of the recently proposed multicanonical method with the well known simulated annealing for the effectiveness in search of the energy global minimum of a biomolecular system. For this we study in detail Met-enkephalin, one of the simplest peptides. We show that the new method not only outperforms simulated annealing in the search of the energy groundstate but also provides more statistical-mechanical information about the system.


title: Low-temperature dynamical simulation of spin-boson systems

abstract: The dynamics of spin-boson systems at very low temperatures has been studied using a real-time path-integral simulation technique which combines a stochastic Monte Carlo sampling over the quantum fluctuations with an exact treatment of the quasiclassical degrees of freedoms. To a large degree, this special technique circumvents the dynamical sign problem and allows the dynamics to be studied directly up to long real times in a numerically exact manner. This method has been applied to two important problems: (1) crossover from nonadiabatic to adiabatic behavior in electron transfer reactions, (2) the zero-temperature dynamics in the antiferromagnetic Kondo region 1/2<K<1 where K is Kondo's parameter.


title: Comment on "Monte Carlo Simulation of a First-Order Transition"

abstract: We comment on a recent article by Hao and Scheraga.


title: Funnels, Pathways and the Energy Landscape of Protein Folding: A Synthesis

abstract: The understanding, and even the description of protein folding is impeded by the complexity of the process. Much of this complexity can be described and understood by taking a statistical approach to the energetics of protein conformation, that is, to the energy landscape. The statistical energy landscape approach explains when and why unique behaviors, such as specific folding pathways, occur in some proteins and more generally explains the distinction between folding processes common to all sequences and those peculiar to individual sequences. This approach also gives new, quantitative insights into the interpretation of experiments and simulations of protein folding thermodynamics and kinetics. Specifically, the picture provides simple explanations for folding as a two-state first-order phase transition, for the origin of metastable collapsed unfolded states and for the curved Arrhenius plots observed in both laboratory experiments and discrete lattice simulations. The relation of these quantitative ideas to folding pathways, to uni-exponential {\em vs.} multi-exponential behavior in protein folding experiments and to the effect of mutations on folding is also discussed. The success of energy landscape ideas in protein structure prediction is also described. The use of the energy landscape approach for analyzing data is illustrated with a quantitative analysis of some recent simulations, and a qualitative analysis of experiments on the folding of three proteins. The work unifies several previously proposed ideas concerning the mechanism protein folding and delimits the regions of validity of these ideas under different thermodynamic conditions.


title: A Hybrid Decomposition Parallel Implementation of the Car-Parrinello Method

abstract: We have developed a flexible hybrid decomposition parallel implementation of the first-principles molecular dynamics algorithm of Car and Parrinello. The code allows the problem to be decomposed either spatially, over the electronic orbitals, or any combination of the two. Performance statistics for 32, 64, 128 and 512 Si atom runs on the Touchstone Delta and Intel Paragon parallel supercomputers and comparison with the performance of an optimized code running the smaller systems on the Cray Y-MP and C90 are presented.


title: Radiative association and inverse predissociation of oxygen atoms

abstract: The formation of $\mbox{O}_2$ by radiative association and by inverse predissociation of ground state oxygen atoms is studied using quantum-mechanical methods. Cross sections, emission spectra, and rate coefficients are presented and compared with prior experimental and theoretical results. At temperatures below 1000~K radiative association occurs by approach along the $1\,{}^3\Pi_u$ state of $\mbox{O}_2$ and above 1000~K inverse predissociation through the $\mbox{B}\,{}^3\Sigma_u^-$ state is the dominant mechanism. This conclusion is supported by a quantitative comparison between the calculations and data obtained from hot oxygen plasma spectroscopy.


title: Molecular Dynamics Simulations of a Pressure-induced Glass Transition

abstract: We simulate the compression of a two-component Lennard-Jones liquid at a variety of constant temperatures using a molecular dynamics algorithm in an isobaric-isothermal ensemble. The viscosity of the liquid increases with pressure, undergoing a broadened transition into a structurally arrested, amorphous state. This transition, like the more familiar one induced by cooling, is correlated with a significant increase in icosahedral ordering. In fact, the structure of the final state, as measured by an analysis of the bonding, is essentially the same in the glassy, frozen state whether produced by squeezing or by cooling under pressure. We have computed an effective hard-sphere packing fraction at the transition, defining the transition pressure or temperature by a cutoff in the diffusion constant, analogous to the traditional laboratory definition of the glass transition by an arbitrary, low cutoff in viscosity. The packing fraction at this transition point is not constant, but is consistently higher for runs compressed at higher temperature. We show that this is because the transition point defined by a constant cutoff in the diffusion constant is not the same as the point of structural arrest, at which further changes in pressure induce no further structural changes, but that the two alternate descriptions may be reconciled by using a thermally activated cutoff for the diffusion constant. This enables estimation of the characteristic activation energy for diffusion at the point of structural arrest.


title: Reversible Work Transiton State Theory: Application to Dissociative Adsorption of Hydrogen

abstract: A practical method for finding free energy barriers for transitions in high-dimensional classical and quantum systems is presented and used to calculate the dissociative sticking probability of H2 on a metal surface within transition state theory (TST). The reversible work involved in shifting the system confined to a hyperplane from the reactant region towards products is evaluated directly. Quantum mechanical degrees of freedom are included by using Feynman Path Integrals with the hyperplane constraint applied to the centroid of the cyclic paths. An optimal dividing surface for the rate estimated by TST is identified naturally in the course of the reversible work evaluation. The free energy barrier is determined relative to the reactant state directly, so an estimate of the transition rate can be obtained without requiring a solvable reference model for the transition state. The method has been applied to calculations of the sticking probability of a thermalized hydrogen gas on a Cu(110) surface. The two hydrogen atoms were included quantum mechanically, and over two hundred atoms in the Cu crystal where included classically. The activation energy for adsorption and desorption was determined and found to be significantly lowered by tunneling at low temperature. The calculated values agree quite well with experimental estimates. Dynamical corrections to the classical TST rate estimate were evaluated and found to be small.


title: Microwave detected, microwave-optical double resonance of NH3, NH2D, NHD2, and ND3. II. Predissociation dynamics of the A state

abstract: Using microwave detected, microwave-optical double resonance, we have measured the homogeneous linewidths of individual rovibrational transitions in the \~A state of NH3, NH2D, NHD2, and ND3. We have used this excited state spectroscopic data to characterize the height of the dissociation barrier and the mechanisms by which the molecule uses its excess vibrational and rotational energies to help overcome this barrier. To interpret the observed vibronic widths, a one dimensional, local mode potential has been developed along a N-H(D) bond. These calculations suggest the barrier height is roughly 2100cm-1, approximately 1000cm-1 below the ab initio prediction. The observed vibronic dependence of levels containing two or more quanta in nu2 is explained by a Fermi resonance between nu2 and the N-H(D) stretch. This interaction also explains the observed trends due to isotopic substitution. The rotational enhancement of the predissociation rates in the NH3 2(1) level is dominated by Coriolis coupling while for the same level in ND3, centrifugal effects dominate.