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.

Causal and localizable quantum operations

We examine constraints on quantum operations imposed by relativistic causality. A bipartite superoperator is said to be localizable if it can be implemented by two parties (Alice and Bob) who share entanglement but do not communicate, it is causal if the superoperator does not convey information from Alice to Bob or from Bob to Alice. We characterize the general structure of causal complete-measurement superoperators, and exhibit examples that are causal but not localizable. We construct another class of causal bipartite superoperators that are not localizable by invoking bounds on the strength of correlations among the parts of a quantum system. A bipartite superoperator is said to be semilocalizable if it can be implemented with one-way quantum communication from Alice to Bob, and it is semicausal if it conveys no information from Bob to Alice. We show that all semicausal complete-measurement superoperators are semi localizable, and we establish a general criterion for semicausality. In the multipartite case, we observe that a measurement superoperator that projects onto the eigenspaces of a stabilizer code is localizable.

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.

Characterizing mixing and measurement in quantum mechanics

What fundamental constraints characterize the relationship between a mixture rho = Sigma (i)p(i)rho (i) of quantum states, the states rho (i) being mixed, and the probabilities p(i)? What fundamental constraints characterize the relationship between prior and posterior states in a quantum measurement? In this paper we show that then are many surprisingly strong constraints on these mixing and measurement processes that can be expressed simply in terms of the eigenvalues of the quantum states involved. These constraints capture in a succinct fashion what it means to say that a quantum measurement acquires information about the system being measured, and considerably simplify the proofs of many results about entanglement transformation.

Noise-free scattering of the quantized electromagnetic field from a dispersive linear dielectric

We study the scattering of the quantized electromagnetic field from a linear, dispersive dielectric using the scattering formalism for quantum fields. The medium is modeled as a collection of harmonic oscillators with a number of distinct resonance frequencies. This model corresponds to the Sellmeir expansion, which is widely used to describe experimental data for real dispersive media. The integral equation for the interpolating field in terms of the in field is solved and the solution used to find the out field. The relation between the ill and out creation and annihilation operators is found that allows one to calculate the S matrix for this system. In this model, we find that there are absorption bands, but the input-output relations are completely unitary. No additional quantum-noise terms are required.

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.

Calculation of product state distributions from resonance decay via Lanczos subspace filter diagonalization: Application to HO2

Resonance phenomena associated with the unimolecular dissociation of HO2 have been investigated quantum-mechanically by the Lanczos homogeneous filter diagonalization (LHFD) method. The calculated resonance energies, rates (widths), and product state distributions are compared to results from an autocorrelation function-based filter diagonalization (ACFFD) method. For calculating resonance wave functions via ACFFD, an analytical expression for the expansion coefficients of the modified Chebyshev polynomials is introduced. Both dissociation rates and product state distributions of O-2 show strong fluctuations, indicating the dissociation of HO2 is essentially irregular. (C) 2001 American Institute of Physics.

Quantum measurement of coherent tunnelling between quantum dots

We describe the conditional and unconditional dynamics of two coupled quantum dots when one dot is subjected to a measurement of its occupation number by coupling it to a third readout dot via the Coulomb interaction. The readout dot is coupled to source and drain leads under weak bias, and a tunnel current flows through a single bound state when energetically allowed. The occupation of the quantum dot near the readout dot shifts the bound state of the readout dot from a low conducting state to a high conducting state. The measurement is made by continuously monitoring the tunnel current through the readout dot. We show that there is a difference between the time scale for the measurement-induced decoherence between the localized states of the dots, and the time scale on which the system becomes localized due to the measurement.

Cavity QED analog of the harmonic-oscillator probability distribution function and quantum collapses

We establish a connection between the simple harmonic oscillator and a two-level atom interacting with resonant, quantized cavity and strong driving fields, which suggests an experiment to measure the harmonic-oscillator's probability distribution function. To achieve this, we calculate the Autler-Townes spectrum by coupling the system to a third level. We find that there are two different regions of the atomic dynamics depending on the ratio of the: Rabi frequency Omega (c) of the cavity field to that of the Rabi frequency Omega of the driving field. For Omega (c)

Effect of quantum interference on a three-level atom driven by two laser fields

We study the effect of quantum interference on the population distribution and absorptive properties of a V-type three-level atom driven by two lasers of unequal intensities and different angular frequencies. Three coupling configurations of the lasers to the atom are analysed: (a) both lasers coupled to the same atomic transition, (b) each laser coupled to different atomic transition and (c) each laser coupled to both atomic transitions. Dressed stales for the three coupling configurations are identified, and the population distribution and absorptive properties of the weaker field are interpreted in terms of transition dipole moments and transition frequencies among these dressed states. In particular, we find that in the first two cases there is no population inversion between the bare atomic states, but the population can be trapped in a superposition of the dressed states induced by quantum interference and the stronger held. We show that the trapping of the population, which results from the cancellation of transition dipole moments, does not prevent the weaker field to be coupled to the cancelled (dark) transitions. As a result, the weaker field can be strongly amplified on transparent transitions. In the case of each laser coupled to both atomic transitions the population can be trapped in a linear superposition of the excited bare atomic states leaving the ground state unpopulated in the steady state. Moreover, we find that the absorption rate of the weaker field depends on the detuning of the strong field from the atomic resonances and the splitting between the atomic excited states. When the strong held is resonant to one of the atomic transitions a quasi-trapping effect appears in one of the dressed states. In the quasi-trapping situation all the transition dipole moments are different from zero, which allows the weaker field to be amplified on the inverted transitions. When the strong field is tuned halfway between the atomic excited states, the population is completely trapped in one of the dressed states and no amplification is found for the weaker field.

Non-Markovian quantum optical processes and pseudomodes

Non-Markovian behaviour in atomic systems coupled to a structured reservoir of quantum EM field modes, such as in high Q cavities, is treated using a quasimode description, and the pseudo mode theory for single quantum reservoir excitations is obtained via Fano diagonalisation. The atomic transitions are coupled to a discrete set of (cavity) quasimodes, which are also coupled to a continuum set of (external) quasimodes with slowly varying coupling constants. Each pseudomode corresponds to a cavity quasimode, and the original reservoir structure is obtained in expressions for the equivalent atom-true mode coupling constants. Cases of multiple excitation of the reservoir are now treatable via Markovian master equations for the atom-discrete quasimode system.