A rigorous implementation of the Jeans--Landau--Teller approximation

Rigorous bounds on the rate of energy exchanges between vibrational and translational degrees of freedom are established in simple classical models of diatomic molecules. The results are in agreement with an elementary approximation introduced by Landau and Teller. The method is perturbative theory ``beyond all orders'', with diagrammatic techniques (tree expansions) to organize and manipulate terms, and look for compensations, like in recent studies on KAM theorem homoclinic splitting.Comment: 23 pages, postscrip

Boundary effects in the stepwise structure of the Lyapunov spectra for quasi-one-dimensional systems

Boundary effects in the stepwise structure of the Lyapunov spectra and the corresponding wavelike structure of the Lyapunov vectors are discussed numerically in quasi-one-dimensional systems consisting of many hard-disks. Four kinds of boundary conditions constructed by combinations of periodic boundary conditions and hard-wall boundary conditions are considered, and lead to different stepwise structures of the Lyapunov spectra in each case. We show that a spatial wavelike structure with a time-oscillation appears in the spatial part of the Lyapunov vectors divided by momenta in some steps of the Lyapunov spectra, while a rather stationary wavelike structure appears in the purely spatial part of the Lyapunov vectors corresponding to the other steps. Using these two kinds of wavelike structure we categorize the sequence and the kinds of steps of the Lyapunov spectra in the four different boundary condition cases.Comment: 33 pages, 25 figures including 10 color figures. Manuscript including the figures of better quality is available from http://newt.phys.unsw.edu.au/~gary/step.pd

Steady-state conduction in self-similar billiards

The self-similar Lorentz billiard channel is a spatially extended deterministic dynamical system which consists of an infinite one-dimensional sequence of cells whose sizes increase monotonically according to their indices. This special geometry induces a nonequilibrium stationary state with particles flowing steadily from the small to the large scales. The corresponding invariant measure has fractal properties reflected by the phase-space contraction rate of the dynamics restricted to a single cell with appropriate boundary conditions. In the near-equilibrium limit, we find numerical agreement between this quantity and the entropy production rate as specified by thermodynamics

Spectral properties of quantum NN-body systems versus chaotic properties of their mean field approximations

We present numerical evidence that in a system of interacting bosons there exists a correspondence between the spectral properties of the exact quantum Hamiltonian and the dynamical chaos of the associated mean field evolution. This correspondence, analogous to the usual quantum-classical correspondence, is related to the formal parallel between the second quantization of the mean field, which generates the exact dynamics of the quantum $N$-body system, and the first quantization of classical canonical coordinates. The limit of infinite density and the thermodynamic limit are then briefly discussed.Comment: 15 pages RevTeX, 11 postscript figures included with psfig, uuencoded gz-compressed .tar fil

Tran-SAS v1.0: A numerical model to compute catchment-scale hydrologic transport using StorAge Selection functions

This paper presents the tran-SAS package, which includes a set of codes to model solute transport and water residence times through a hydrological system. The model is based on a catchment-scale approach that aims at reproducing the integrated response of the system at one of its outlets. The codes are implemented in MATLAB and are meant to be easy to edit, so that users with minimal programming knowledge can adapt them to the desired application. The problem of large-scale solute transport has both theoretical and practical implications. On the one side, the ability to represent the ensemble of water flow trajectories through a heterogeneous system helps unraveling streamflow generation processes and allows us to make inferences on plant–water interactions. On the other side, transport models are a practical tool that can be used to estimate the persistence of solutes in the environment. The core of the package is based on the implementation of an age master equation (ME), which is solved using general StorAge Selection (SAS) functions. The age ME is first converted into a set of ordinary differential equations, each addressing the transport of an individual precipitation input through the catchment, and then it is discretized using an explicit numerical scheme. Results show that the implementation is efficient and allows the model to run in short times. The numerical accuracy is critically evaluated and it is shown to be satisfactory in most cases of hydrologic interest. Additionally, a higher-order implementation is provided within the package to evaluate and, if necessary, to improve the numerical accuracy of the results. The codes can be used to model streamflow age and solute concentration, but a number of additional outputs can be obtained by editing the codes to further advance the ability to understand and model catchment transport processes

Persistent Chaos in High Dimensions

An extensive statistical survey of universal approximators shows that as the dimension of a typical dissipative dynamical system is increased, the number of positive Lyapunov exponents increases monotonically and the number of parameter windows with periodic behavior decreases. A subset of parameter space remains in which topological change induced by small parameter variation is very common. It turns out, however, that if the system's dimension is sufficiently high, this inevitable, and expected, topological change is never catastrophic, in the sense chaotic behavior is preserved. One concludes that deterministic chaos is persistent in high dimensions.Comment: 4 pages, 3 figures; Changes in response to referee comment

Transport of fluorobenzoate tracers in a vegetated hydrologic control volume: 2. Theoretical inferences and modeling

A theoretical analysis of transport in a controlled hydrologic volume, inclusive of two willow trees and forced by erratic water inputs, is carried out contrasting the experimental data described in a companion paper. The data refer to the hydrologic transport in a large lysimeter of different fluorobenzoic acids seen as tracers. Export of solute is modeled through a recently developed framework which accounts for nonstationary travel time distributions where we parameterize how output fluxes (namely, discharge and evapotranspiration) sample the available water ages in storage. The relevance of this work lies in the study of hydrologic drivers of the nonstationary character of residence and travel time distributions, whose definition and computation shape this theoretical transport study. Our results show that a large fraction of the different behaviors exhibited by the tracers may be charged to the variability of the hydrologic forcings experienced after the injection. Moreover, the results highlight the crucial, and often overlooked, role of evapotranspiration and plant uptake in determining the transport of water and solutes. This application also suggests that the ways evapotranspiration selects water with different ages in storage can be inferred through model calibration contrasting only tracer concentrations in the discharge. A view on upscaled transport volumes like hillslopes or catchments is maintained throughout the paper

Transition from regular to complex behaviour in a discrete deterministic asymmetric neural network model

We study the long time behaviour of the transient before the collapse on the periodic attractors of a discrete deterministic asymmetric neural networks model. The system has a finite number of possible states so it is not possible to use the term chaos in the usual sense of sensitive dependence on the initial condition. Nevertheless, at varying the asymmetry parameter, $k$, one observes a transition from ordered motion (i.e. short transients and short periods on the attractors) to a ``complex'' temporal behaviour. This transition takes place for the same value $k_{\rm c}$ at which one has a change for the mean transient length from a power law in the size of the system ($N$) to an exponential law in $N$. The ``complex'' behaviour during the transient shows strong analogies with the chaotic behaviour: decay of temporal correlations, positive Shannon entropy, non-constant Renyi entropies of different orders. Moreover the transition is very similar to that one for the intermittent transition in chaotic systems: scaling law for the Shannon entropy and strong fluctuations of the ``effective Shannon entropy'' along the transient, for $k > k_{\rm c}$.Comment: 18 pages + 6 figures, TeX dialect: Plain TeX + IOP macros (included

Localized behavior in the Lyapunov vectors for quasi-one-dimensional many-hard-disk systems

We introduce a definition of a "localization width" whose logarithm is given by the entropy of the distribution of particle component amplitudes in the Lyapunov vector. Different types of localization widths are observed, for example, a minimum localization width where the components of only two particles are dominant. We can distinguish a delocalization associated with a random distribution of particle contributions, a delocalization associated with a uniform distribution and a delocalization associated with a wave-like structure in the Lyapunov vector. Using the localization width we show that in quasi-one-dimensional systems of many hard disks there are two kinds of dependence of the localization width on the Lyapunov exponent index for the larger exponents: one is exponential, and the other is linear. Differences, due to these kinds of localizations also appear in the shapes of the localized peaks of the Lyapunov vectors, the Lyapunov spectra and the angle between the spatial and momentum parts of the Lyapunov vectors. We show that the Krylov relation for the largest Lyapunov exponent $\lambda\sim-\rho\ln\rho$ as a function of the density $\rho$ is satisfied (apart from a factor) in the same density region as the linear dependence of the localization widths is observed. It is also shown that there are asymmetries in the spatial and momentum parts of the Lyapunov vectors, as well as in their $x$ and $y$-components.Comment: 41 pages, 21 figures, Manuscript including the figures of better quality is available from http://www.phys.unsw.edu.au/~gary/Research.htm

Geometric dynamical observables in rare gas crystals

We present a detailed description of how a differential geometric approach to Hamiltonian dynamics can be used for determining the existence of a crossover between different dynamical regimes in a realistic system, a model of a rare gas solid. Such a geometric approach allows to locate the energy threshold between weakly and strongly chaotic regimes, and to estimate the largest Lyapunov exponent. We show how standard mehods of classical statistical mechanics, i.e. Monte Carlo simulations, can be used for our computational purposes. Finally we consider a Lennard Jones crystal modeling solid Xenon. The value of the energy threshold turns out to be in excellent agreement with the numerical estimate based on the crossover between slow and fast relaxation to equilibrium obtained in a previous work by molecular dynamics simulations.Comment: RevTeX, 19 pages, 6 PostScript figures, submitted to Phys. Rev.