Stability of nonlinear difference inclusions

The stability of difference inclusions x(k+1) is an element of F(x(k)) is studied, where F(x) = {F(x, gimel) : is an element of Lambda} and the selections F(., gimel) : E -->E assume values in a Banach space E, partially ordered by a cone K. It is assumed that the operators F(.,gimel) are heterotone or pseudoconcave. The main results concern asymptotically stable absorbing sets, and include the case of a single equilibrium point. The results are applied to a number of practical problems.

A new wavelet-based method for the solution of the population balance equation

A new wavelet-based method for solving population balance equations with simultaneous nucleation, growth and agglomeration is proposed, which uses wavelets to express the functions. The technique is very general, powerful and overcomes the crucial problems of numerical diffusion and stability that often characterize previous techniques in this area. It is also applicable to an arbitrary grid to control resolution and computational efficiency. The proposed technique has been tested for pure agglomeration, simultaneous nucleation and growth, and simultaneous growth and agglomeration. In all cases, the predicted and analytical particle size distributions are in excellent agreement. The presence of moving sharp fronts can be addressed without the prior investigation of the characteristics of the processes. (C) 2001 Published by Elsevier Science Ltd.

Quantum criticality

We investigate the theory of quantum fluctuations in non-equilibrium systems having large critical fluctuations. This allows us to treat the limits imposed by nonlinearities to quantum squeezing and noise reduction, and also to envisage future tests of quantum theory in regions of macroscopic quantum fluctuations. A long-term objective of this research is to identify suitable physical systems in which macroscopic 'Schrodinger cat'-like behaviour may be observed. We investigate two systems in particular of much current experimental interest, namely the degenerate parametric oscillator near threshold, and the evaporatively cooled (BEC). We compare the results obtained in the positive-P representation, as a fully quantum mechanical calculation, with the truncated Wigner phase space equation, also known as semi-classical theory. We show when these results agree and differ in calculations taken beyond the linearized approximation. In the region where the largest quantum fluctuations and Schrodinger cat-like behaviour might be expected, we find that the quantum predictions correspond very closely to the semi-classical theory. Nature abhors observing a Schrodinger car.

Continuous quantum measurement of two coupled quantum dots using a point contact: A quantum trajectory approach

We obtain the finite-temperature unconditional master equation of the density matrix for two coupled quantum dots (CQD's) when one dot is subjected to a measurement of its electron occupation number using a point contact (PC). To determine how the CQD system state depends on the actual current through the PC device, we use the so-called quantum trajectory method to derive the zero-temperature conditional master equation. We first treat the electron tunneling through the PC barrier as a classical stochastic point process (a quantum-jump model). Then we show explicitly that our results can be extended to the quantum-diffusive limit when the average electron tunneling rate is very large compared to the extra change of the tunneling rate due to the presence of the electron in the dot closer to the PC. We find that in both quantum-jump and quantum-diffusive cases, the conditional dynamics of the CQD system can be described by the stochastic Schrodinger equations for its conditioned state vector if and only if the information carried away from the CQD system by the PC reservoirs can be recovered by the perfect detection of the measurements.

Light guiding light: Nonlinear refraction in rubidium vapor

Recently there has been experimental and theoretical interest in cross-dispersion effects in rubidium vapor, which allows one beam of light to be guided by another. We present theoretical results which account for the complications created by the D line hyperfine structure of rubidium as well as the presence of the two major isotopes of rubidium. This allows the complex frequency dependence of the effects observed in our experiments to be understood and lays the foundation for future studies of nonlinear propagation.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

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.

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.

The Dubinin-Radushkevich equation and the underlying microscopic adsorption description

The Dubinin-Radushkevich (DR) equation is widely used for description of adsorption in microporous materials, especially those of a carbonaceous origin. The equation has a semi-empirical origin and is based on the assumptions of a change in the potential energy between the gas and adsorbed phases and a characteristic energy of a given solid. This equation yields a macroscopic behaviour of adsorption loading for a given pressure. In this paper, we apply a theory developed in our group to investigate the underlying mechanism of adsorption as an alternative to the macroscopic description using the DR equation. Using this approach, we are able to establish a detailed picture of the adsorption in the whole range of the micropore system. This is different from the DR equation, which provides an overall description of the process. (C) 2001 Elsevier Science Ltd. All rights reserved.

Thue's theorem and the diophantine equation x2 - Dy2 = ±N

A constructive version of a theorem of Thue is used to provide representations of certain integers as x(2) - Dy-2, where D = 2, 3, 5, 6, 7.

Critical quantum fluctuations in the degenerate parametric oscillator

We develop a systematic theory of critical quantum fluctuations in the driven parametric oscillator. Our analytic results agree well with stochastic numerical simulations. We also compare the results obtained in the positive-P representation, as a fully quantum-mechanical calculation, with the truncated Wigner phase-space equation, also known as the semiclassical theory. We show when these results agree and differ in calculations taken beyond the linearized approximation. We find that the optimal broadband noise reduction occurs just above threshold. In this region where there are large quantum fluctuations in the conjugate variance and macroscopic quantum superposition states might be expected, we find that the quantum predictions correspond very closely to the semiclassical theory.

On positive solutions of nonlinear elliptic equations involving concave and critical nonlinearities

We study the existence of nonnegative solutions of elliptic equations involving concave and critical Sobolev nonlinearities. Applying various variational principles we obtain the existence of at least two nonnegative solutions.