Breaking conjugate pairing in thermostatted billiards by magnetic field

We demonstrate that in the thermostatted three-dimensional Lorentz gas the symmetry of the Lyapunov spectrum can be broken by adding to the system an external magnetic field not perpendicular to the electric field. For perpendicular field vectors, there is a Hamiltonian reformulation of the dynamics and the conjugate pairing rule still holds. This indicates that symmetric Lyapunov spectra has nothing to do with time reversal symmetry or reversibility; instead, it seems to be related to the existence of a Hamiltonian connection.Comment: 4 pages, 3 figure

Peeping at chaos: Nondestructive monitoring of chaotic systems by measuring long-time escape rates

One or more small holes provide non-destructive windows to observe corresponding closed systems, for example by measuring long time escape rates of particles as a function of hole sizes and positions. To leading order the escape rate of chaotic systems is proportional to the hole size and independent of position. Here we give exact formulas for the subsequent terms, as sums of correlation functions; these depend on hole size and position, hence yield information on the closed system dynamics. Conversely, the theory can be readily applied to experimental design, for example to control escape rates.Comment: Originally 4 pages and 2 eps figures incorporated into the text; v2 has more numerical results and discussion: now 6 pages, 4 figure

Open circular billiards and the Riemann hypothesis

A comparison of escape rates from one and from two holes in an experimental container (e.g. a laser trap) can be used to obtain information about the dynamics inside the container. If this dynamics is simple enough one can hope to obtain exact formulas. Here we obtain exact formulas for escape from a circular billiard with one and with two holes. The corresponding quantities are expressed as sums over zeroes of the Riemann zeta function. Thus we demonstrate a direct connection between recent experiments and a major unsolved problem in mathematics, the Riemann hypothesis.Comment: 5 pages, 4 embedded postscript figures; v2: more explicit on how the Reimann Hypothesis arises from a comparison of one and two hole escape rate

A Paradox of State-Dependent Diffusion and How to Resolve It

Consider a particle diffusing in a confined volume which is divided into two equal regions. In one region the diffusion coefficient is twice the value of the diffusion coefficient in the other region. Will the particle spend equal proportions of time in the two regions in the long term? Statistical mechanics would suggest yes, since the number of accessible states in each region is presumably the same. However, another line of reasoning suggests that the particle should spend less time in the region with faster diffusion, since it will exit that region more quickly. We demonstrate with a simple microscopic model system that both predictions are consistent with the information given. Thus, specifying the diffusion rate as a function of position is not enough to characterize the behaviour of a system, even assuming the absence of external forces. We propose an alternative framework for modelling diffusive dynamics in which both the diffusion rate and equilibrium probability density for the position of the particle are specified by the modeller. We introduce a numerical method for simulating dynamics in our framework that samples from the equilibrium probability density exactly and is suitable for discontinuous diffusion coefficients.Comment: 21 pages, 6 figures. Second round of revisions. This is the version that will appear in Proc Roy So

Microscopic chaos from Brownian motion?

A recent experiment on Brownian motion has been interpreted to exhibit direct evidence for microscopic chaos. In this note we demonstrate that virtually identical results can be obtained numerically using a manifestly microscopically nonchaotic system.Comment: 3 pages, 1 figure, Comment on P. Gaspard et al, Nature vol 394, 865 (1998); rewritten in a more popular styl

Anomalous heat conduction and anomalous diffusion in nonlinear lattices, single walled nanotubes, and billiard gas channels

We study anomalous heat conduction and anomalous diffusion in low dimensional systems ranging from nonlinear lattices, single walled carbon nanotubes, to billiard gas channels. We find that in all discussed systems, the anomalous heat conductivity can be connected with the anomalous diffusion, namely, if energy diffusion is $\sigma^2(t)\equiv =2Dt^{\alpha} (0<\alpha\le 2)$, then the thermal conductivity can be expressed in terms of the system size $L$ as $\kappa = cL^{\beta}$ with $\beta=2-2/\alpha$. This result predicts that a normal diffusion ($\alpha =1$) implies a normal heat conduction obeying the Fourier law ($\beta=0$), a superdiffusion ($\alpha>1$) implies an anomalous heat conduction with a divergent thermal conductivity ($\beta>0$), and more interestingly, a subdiffusion ($\alpha <1$) implies an anomalous heat conduction with a convergent thermal conductivity ($\beta<0$), consequently, the system is a thermal insulator in the thermodynamic limit. Existing numerical data support our theoretical prediction.Comment: 15 Revtex pages, 16 figures. Invited article for CHAOS focus issue commemorating the 50th anniversary of the Fermi-Pasta-Ulam (FPU) mode

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted

Gravity Waves, Chaos, and Spinning Compact Binaries

Spinning compact binaries are shown to be chaotic in the Post-Newtonian expansion of the two body system. Chaos by definition is the extreme sensitivity to initial conditions and a consequent inability to predict the outcome of the evolution. As a result, the spinning pair will have unpredictable gravitational waveforms during coalescence. This poses a challenge to future gravity wave observatories which rely on a match between the data and a theoretical template.Comment: Final version published in PR

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