A uniqueness result for scattering by infinite rough surfaces

Consider the Dirichlet boundary value problem for the Helmholtz equation in a non-locally perturbed half-plane with an unbounded, piecewise Lyapunov boundary. This problem models time-harmonic electromagnetic scattering in transverse magnetic polarization by one-dimensional rough, perfectly conducting surfaces. A radiation condition is introduced for the problem, which is a generalization of the usual one used in the study of diffraction by gratings when the solution is quasi-periodic, and allows a variety of incident fields including an incident plane wave to be included in the results obtained. We show in this paper that the boundary value problem for the scattered field has at most one solution. For the case when the whole boundary is Lyapunov and is a small perturbation of a flat boundary we also prove existence of solution and show a limiting absorption principle.

Averaging, slaving and balance dynamics in a simple atmospheric model

We report numerical results from a study of balance dynamics using a simple model of atmospheric motion that is designed to help address the question of why balance dynamics is so stable. The non-autonomous Hamiltonian model has a chaotic slow degree of freedom (representing vortical modes) coupled to one or two linear fast oscillators (representing inertia-gravity waves). The system is said to be balanced when the fast and slow degrees of freedom are separated. We find adiabatic invariants that drift slowly in time. This drift is consistent with a random-walk behaviour at a speed which qualitatively scales, even for modest time scale separations, as the upper bound given by Neishtadt’s and Nekhoroshev’s theorems. Moreover, a similar type of scaling is observed for solutions obtained using a singular perturbation (‘slaving’) technique in resonant cases where Nekhoroshev’s theorem does not apply. We present evidence that the smaller Lyapunov exponents of the system scale exponentially as well. The results suggest that the observed stability of nearly-slow motion is a consequence of the approximate adiabatic invariance of the fast motion.

A Wiener-Hopf type factorization for the exponential functional of Levy processes

For a Lévy process ξ=(ξt)t≥0 drifting to −∞, we define the so-called exponential functional as follows: Formula Under mild conditions on ξ, we show that the following factorization of exponential functionals: Formula holds, where × stands for the product of independent random variables, H− is the descending ladder height process of ξ and Y is a spectrally positive Lévy process with a negative mean constructed from its ascending ladder height process. As a by-product, we generate an integral or power series representation for the law of Iξ for a large class of Lévy processes with two-sided jumps and also derive some new distributional properties. The proof of our main result relies on a fine Markovian study of a class of generalized Ornstein–Uhlenbeck processes, which is itself of independent interest. We use and refine an alternative approach of studying the stationary measure of a Markov process which avoids some technicalities and difficulties that appear in the classical method of employing the generator of the dual Markov process.

A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids

In addition to the Hamiltonian functional itself, non-canonical Hamiltonian dynamical systems generally possess integral invariants known as ‘Casimir functionals’. In the case of the Euler equations for a perfect fluid, the Casimir functionals correspond to the vortex topology, whose invariance derives from the particle-relabelling symmetry of the underlying Lagrangian equations of motion. In a recent paper, Vallis, Carnevale & Young (1989) have presented algorithms for finding steady states of the Euler equations that represent extrema of energy subject to given vortex topology, and are therefore stable. The purpose of this note is to point out a very general method for modifying any Hamiltonian dynamical system into an algorithm that is analogous to those of Vallis etal. in that it will systematically increase or decrease the energy of the system while preserving all of the Casimir invariants. By incorporating momentum into the extremization procedure, the algorithm is able to find steadily-translating as well as steady stable states. The method is applied to a variety of perfect-fluid systems, including Euler flow as well as compressible and incompressible stratified flow.

Complex dynamics in simple systems with periodic parameter oscillations

We study systems with periodically oscillating parameters that can give way to complex periodic or nonperiodic orbits. Performing the long time limit, we can define ergodic averages such as Lyapunov exponents, where a negative maximal Lyapunov exponent corresponds to a stable periodic orbit. By this, extremely complicated periodic orbits composed of contracting and expanding phases appear in a natural way. Employing the technique of ϵ-uncertain points, we find that values of the control parameters supporting such periodic motion are densely embedded in a set of values for which the motion is chaotic. When a tiny amount of noise is coupled to the system, dynamics with positive and with negative nontrivial Lyapunov exponents are indistinguishable. We discuss two physical systems, an oscillatory flow inside a duct and a dripping faucet with variable water supply, where such a mechanism seems to be responsible for a complicated alternation of laminar and turbulent phases.

Chaos in foreign exchange markets: a sceptical view

This paper tests directly for deterministic chaos in a set of ten daily Sterling-denominated exchange rates by calculating the largest Lyapunov exponent. Although in an earlier paper, strong evidence of nonlinearity has been shown, chaotic tendencies are noticeably absent from all series considered using this state-of-the-art technique. Doubt is cast on many recent papers which claim to have tested for the presence of chaos in economic data sets, based on what are argued here to be inappropriate techniques.

Adsorption of organic molecules at the TiO2(110) surface: the effect of van der Waals interactions

Understanding the interaction of organic molecules with TiO2 surfaces is important for a wide range of technological applications. While density functional theory (DFT) calculations can provide valuable insight about these interactions, traditional DFT approaches with local exchange-correlation functionals suffer from a poor description of non-bonding van der Waals (vdW) interactions. We examine here the contribution of vdW forces to the interaction of small organic molecules (methane, methanol, formic acid and glycine) with the TiO2 (110) surface, based on DFT calculations with the optB88-vdW functional. The adsorption geometries and energies at different configurations were also obtained in the standard generalized gradient approximation (GGA-PBE) for comparison. We find that the optB88-vdW consistently gives shorter surface adsorbate-to-surface distances and slightly stronger interactions than PBE for the weak (physisorbed) modes of adsorption. In the case of strongly adsorbed (chemisorbed) molecules both functionals give similar results for the adsorption geometries, and also similar values of the relative energies between different chemisorption modes for each molecule. In particular both functionals predict that dissociative adsorption is more favourable than molecular adsorption for methanol, formic acid and glycine, in general agreement with experiment. The dissociation energies obtained from both functionals are also very similar, indicating that vdW interactions do not affect the thermodynamics of surface deprotonation. However, the optB88-vdW always predicts stronger adsorption than PBE. The comparison of the methanol adsorption energies with values obtained from a Redhead analysis of temperature programmed desorption data suggests that optB88-vdW significantly overestimates the adsorption strength, although we warn about the uncertainties involved in such comparisons.

Evolution PDEs and augmented eigenfunctions. half line.

The solution of an initial-boundary value problem for a linear evolution partial differential equation posed on the half-line can be represented in terms of an integral in the complex (spectral) plane. This representation is obtained by the unified transform introduced by Fokas in the 90's. On the other hand, it is known that many initial-boundary value problems can be solved via a classical transform pair, constructed via the spectral analysis of the associated spatial operator. For example, the Dirichlet problem for the heat equation can be solved by applying the Fourier sine transform pair. However, for many other initial-boundary value problems there is no suitable transform pair in the classical literature. Here we pose and answer two related questions: Given any well-posed initial-boundary value problem, does there exist a (non-classical) transform pair suitable for solving that problem? If so, can this transform pair be constructed via the spectral analysis of a differential operator? The answer to both of these questions is positive and given in terms of augmented eigenfunctions, a novel class of spectral functionals. These are eigenfunctions of a suitable differential operator in a certain generalised sense, they provide an effective spectral representation of the operator, and are associated with a transform pair suitable to solve the given initial-boundary value problem.