Stiffly accurate Runge-Kutta methods for stiff stochastic differential equations

In this paper we discuss implicit methods based on stiffly accurate Runge-Kutta methods and splitting techniques for solving Stratonovich stochastic differential equations (SDEs). Two splitting techniques: the balanced splitting technique and the deterministic splitting technique, are used in this paper. We construct a two-stage implicit Runge-Kutta method with strong order 1.0 which is corrected twice and no update is needed. The stability properties and numerical results show that this approach is suitable for solving stiff SDEs. (C) 2001 Elsevier Science B.V. All rights reserved.

Stochastic gauges in quantum dynamics for many-body simulations

Quantum dynamics simulations can be improved using novel quasiprobability distributions based on non-orthogonal Hermitian kernel operators. This introduces arbitrary functions (gauges) into the stochastic equations. which can be used to tailor them for improved calculations. A possible application to full quantum dynamic simulations of BEC's is presented. (C) 2001 Elsevier Science B.V. All rights reserved.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

Transforming chaos to periodic oscillations

We demonstrate that the dynamics of an autonomous chaotic class C laser can be controlled to a periodic state via external modulation of the pump. In the absence of modulation, above the chaos threshold, the laser exhibits Lorenz-like chaotic pulsations. The average amplitude and frequency of these pulsations depend on the pump power. We find that there exist parameter windows where modulation of the pump power extinguishes the chaos in favor of simpler periodic behavior. Moreover we find a number of locking ratios between the pump and laser output follow the Farey sequence.

Xmds: eXtensible multi-dimensional simulator

A software package that efficiently solves a comprehensive range of problems based on coupled complex nonlinear stochastic ODEs and PDEs is outlined. Its input and output syntax is formulated as a subset of XML, thus making a step towards a standard for specifying numerical simulations.

Free field and parafermionic realizations of twisted su(3)(k)((2)) current algebra

Free field and twisted parafermionic representations of twisted su(3)(k)((2)) current algebra are obtained. The corresponding twisted Sugawara energy-momentum tensor is given in terms of three (beta, gamma) pairs and two scalar fields and also in terms of twisted parafermionic currents and one scalar field. Two screening currents of the first kind are presented in terms of the free fields.

Riddling and invariance for discontinuous maps preserving Lebesgue measure

In this paper we use the mixture of topological and measure-theoretic dynamical approaches to consider riddling of invariant sets for some discontinuous maps of compact regions of the plane that preserve two-dimensional Lebesgue measure. We consider maps that are piecewise continuous and with invertible except on a closed zero measure set. We show that riddling is an invariant property that can be used to characterize invariant sets, and prove results that give a non-trivial decomposion of what we call partially riddled invariant sets into smaller invariant sets. For a particular example, a piecewise isometry that arises in signal processing (the overflow oscillation map), we present evidence that the closure of the set of trajectories that accumulate on the discontinuity is fully riddled. This supports a conjecture that there are typically an infinite number of periodic orbits for this system.

Set-valued Markov chains and negative semitrajectories of discretized dynamical systems

Computer simulation of dynamical systems involves a phase space which is the finite set of machine arithmetic. Rounding state values of the continuous system to this grid yields a spatially discrete dynamical system, often with different dynamical behaviour. Discretization of an invertible smooth system gives a system with set-valued negative semitrajectories. As the grid is refined, asymptotic behaviour of the semitrajectories follows probabilistic laws which correspond to a set-valued Markov chain, whose transition probabilities can be explicitly calculated. The results are illustrated for two-dimensional dynamical systems obtained by discretization of fractional linear transformations of the unit disc in the complex plane.

Exact form factors for the Josephson tunneling current and relative particle number fluctuations in a model of two coupled Bose-Einstein condensates

Form factors are derived for a model describing the coherent Josephson tunneling between two coupled Bose-Einstein condensates. This is achieved by studying the exact solution of the model within the framework of the algebraic Bethe ansatz. In this approach the form factors are expressed through determinant representations which are functions of the roots of the Bethe ansatz equations.

Integrability and exact spectrum of a pairing model for nucleons

A pairing model for nucleons, introduced by Richardson in 1966, which describes proton-neutron pairing as well as proton-proton and neutron-neutron pairing, is re-examined in the context of the quantum inverse scattering method. Specifically, this shows that the model is integrable by enabling the explicit construction of the conserved operators. We determine the eigenvalues of these operators in terms of the Bethe ansatz, which in turn leads to an expression for the energy eigenvalues of the Hamiltonian.

Exact solution at integrable coupling of a model for the Josephson effect between small metallic grains

A model is introduced for two reduced BCS systems which are coupled through the transfer of Cooper pairs between the systems. The model may thus be used in the analysis of the Josephson effect arising from pair tunneling between two strongly coupled small metallic grains. At a particular coupling strength the model is integrable and explicit results are derived for the energy spectrum, conserved operators, integrals of motion, and wave function scalar products. It is also shown that form factors can be obtained for the calculation of correlation functions. Furthermore, a connection with perturbed conformal field theory is made.

Integrable coupling in a model for Josephson tunneling between non-identical BCS systems

We extend a recent construction for an integrable model describing Josephson tunneling between identical BCS systems to the case where the BCS systems have different single particle energy levels. The exact solution of this generalized model is obtained through the Bethe ansatz.

Integrable variant of the one-dimensional Hubbard model

A new integrable model which is a variant of the one-dimensional Hubbard model is proposed. The integrability of the model is verified by presenting the associated quantum R-matrix which satisfies the Yang-Baxter equation. We argue that the new model possesses the SO(4) algebra symmetry, which contains a representation of the eta-pairing SU(2) algebra and a spin SU(2) algebra. Additionally, the algebraic Bethe ansatz is studied by means of the quantum inverse scattering method. The spectrum of the Hamiltonian, eigenvectors, as well as the Bethe ansatz equations, are discussed. (C) 2002 American Institute of Physics.

Economic and environmental impacts of pollution control in a system of environment and economic interdependence

The importance of the rate of change of the pollution stock in determining the damage to the environment has been an issue of increasing concern in the literature. This paper uses a three-sector (economy, population and environment), non-linear, discrete time, calibrated model to examine pollution control. The model explicitly links economic growth to the health of the environment. The stock of natural resources is affected by the rate of pollution flows, through their impact on the regenerative capacity of the natural resource stock. This can shed useful insights into pollution control strategies, particularly in developing countries where environmental resources are crucial for production in many sectors of the economy. Simulation exercises suggested that, under plausible assumptions, it is possible to reverse undesirable transient dynamics through pollution control expenditure, but this is dependent upon the strategies used for control. The best strategy is to spend money fostering the development of production technologies that reduce pollution rather than spending money dealing with the effects of the pollution flow into the environment. (C) 2001 Elsevier Science Ltd. All rights reserved.

Controlling the complex Lorenz equations by modulation

We demonstrate that a system obeying the complex Lorenz equations in the deep chaotic regime can be controlled to periodic behavior by applying a modulation to the pump parameter. For arbitrary modulation frequency and amplitude there is no obvious simplification of the dynamics. However, we find that there are numerous windows where the chaotic system has been controlled to different periodic behaviors. The widths of these windows in parameter space are narrow, and the positions are related to the ratio of the modulation frequency of the pump to the average pulsation frequency of the output variable. These results are in good agreement with observations previously made in a far-infrared laser system.