Theory of pseudomodes in quantum optical processes

This paper deals with non-Markovian behavior in atomic systems coupled to a structured reservoir of quantum electromagnetic field modes, with particular relevance to atoms interacting with the field in high-Q cavities or photonic band-gap materials. In cases such as the former, we show that the pseudomode theory for single-quantum reservoir excitations can be obtained by applying the Fano diagonalization method to a system in which the atomic transitions are coupled to a discrete set of (cavity) quasimodes, which in turn are coupled to a continuum set of (external) quasimodes with slowly varying coupling constants and continuum mode density. Each pseudomode can be identified with a discrete quasimode, which gives structure to the actual reservoir of true modes via the expressions for the equivalent atom-true mode coupling constants. The quasimode theory enables cases of multiple excitation of the reservoir to now be treated via Markovian master equations for the atom-discrete quasimode system. Applications of the theory to one, two, and many discrete quasimodes are made. For a simple photonic band-gap model, where the reservoir structure is associated with the true mode density rather than the coupling constants, the single quantum excitation case appears to be equivalent to a case with two discrete quasimodes.

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.

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.

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.

Addressing shortcomings in the Household Responsibility System - Empirical analysis of the Two-Farmland System in Shandong Province

Despite apparent overwhelming benefits, implementation of the Household Responsibility System (HRS) in China contained a number of flaws. The Two-Farmland System (TFS), which originated in Pingdu City in Shandong Province, sought to address the twin problems of land fragmentation and economies of size. A stochastic frontier production function analysis that isolates the impacts of land allocation reforms suggests that the TFS increased efficiency by around 7%. This article highlights the need for empirical analysis to assess objectively the merits or otherwise of particular reforms. (C) 2002 Elsevier Science Inc. All rights reserved.

Feedback-stabilization of an arbitrary pure state of a two-level atom

Unit-efficiency homodyne detection of the resonance fluorescence of a two-level atom collapses the quantum state of the atom to a stochastically moving point on the Bloch sphere. Recently, Hofmann, Mahler, and Hess [Phys. Rev. A 57, 4877 (1998)] showed that by making part of the coherent driving proportional to the homodyne photocurrent one can stabilize the state to any point on the bottom-half of the sphere. Here we reanalyze their proposal using the technique of stochastic master equations, allowing their results to be generalized in two ways. First, we show that any point on the upper- or lower-half, but not the equator, of the sphere may be stabilized. Second, we consider nonunit-efficiency detection, and quantify the effectiveness of the feedback by calculating the maximal purity obtainable in any particular direction in Bloch space.

Estimating economic relationships subject to firm- and time-varying equality and inequality constraints

Applied econometricians often fail to impose economic regularity constraints in the exact form economic theory prescribes. We show how the Singular Value Decomposition (SVD) Theorem and Markov Chain Monte Carlo (MCMC) methods can be used to rigorously impose time- and firm-varying equality and inequality constraints. To illustrate the technique we estimate a system of translog input demand functions subject to all the constraints implied by economic theory, including observation-varying symmetry and concavity constraints. Results are presented in the form of characteristics of the estimated posterior distributions of functions of the parameters. Copyright (C) 2001 John Wiley Sons, Ltd.

Path integrals for continuous-time Markov chains

This note presents a method of evaluating the distribution of a path integral for Markov chains on a countable state space.

Two-link approximation schemes for loss networks with linear structure and trunk reservation

Loss networks have long been used to model various types of telecommunication network, including circuit-switched networks. Such networks often use admission controls, such as trunk reservation, to optimize revenue or stabilize the behaviour of the network. Unfortunately, an exact analysis of such networks is not usually possible, and reduced-load approximations such as the Erlang Fixed Point (EFP) approximation have been widely used. The performance of these approximations is typically very good for networks without controls, under several regimes. There is evidence, however, that in networks with controls, these approximations will in general perform less well. We propose an extension to the EFP approximation that gives marked improvement for a simple ring-shaped network with trunk reservation. It is based on the idea of considering pairs of links together, thus making greater allowance for dependencies between neighbouring links than does the EFP approximation, which only considers links in isolation.

A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Efficient simulation of a Tandem Jackson Network

The two-node tandem Jackson network serves as a convenient reference model for the analysis and testing of different methodologies and techniques in rare event simulation. In this paper we consider a new approach to efficiently estimate the probability that the content of the second buffer exceeds some high level L before it becomes empty, starting from a given state. The approach is based on a Markov additive process representation of the buffer processes, leading to an exponential change of measure to be used in an importance sampling procedure. Unlike changes of measures proposed and studied in recent literature, the one derived here is a function of the content of the first buffer. We prove that when the first buffer is finite, this method yields asymptotically efficient simulation for any set of arrival and service rates. In fact, the relative error is bounded independent of the level L; a new result which is not established for any other known method. When the first buffer is infinite, we propose a natural extension of the exponential change of measure for the finite buffer case. In this case, the relative error is shown to be bounded (independent of L) only when the second server is the bottleneck; a result which is known to hold for some other methods derived through large deviations analysis. When the first server is the bottleneck, experimental results using our method seem to suggest that the relative error is bounded linearly in L.

Quasimode theory of quantum optical processes in photonic band gap materials

This paper deals with atomic systems coupled to a structured reservoir of quantum EM field modes, with particular relevance to atoms interacting with the field in photonic band gap materials. The case of high Q cavities has been treated elsewhere using Fano diagonalization based on a quasimode approach, showing that the cavity quasimodes are responsible for pseudomodes introduced to treat non-Markovian behaviour. The paper considers a simple model of a photonic band gap case, where the spatially dependent permittivity consists of a constant term plus a small spatially periodic term that leads to a narrow band gap in the spectrum of mode frequencies. Most treatments of photonic band gap materials are based on the true modes, obtained numerically by solving the Helmholtz equation for the actual spatially periodic permittivity. Here the field modes are first treated in terms of a simpler quasimode approach, in which the quasimodes are plane waves associated with the constant permittivity term. Couplings between the quasimodes occur owing to the small periodic term in the permittivity, with selection rules for the coupled modes being related to the reciprocal lattice vectors. This produces a field Hamiltonian in quasimode form. A matrix diagonalization method may be applied to relate true mode annihilation operators to those for quasimodes. The atomic transitions are coupled to all the quasimodes, and the true mode atom-EM field coupling constants (one-photon Rabi frequencies) are related to those for the quasimodes and also expressions are obtained for the true mode density. The results for the one-photon Rabi frequencies differ from those assumed in other work. Expressions for atomic decay rates are obtained using the Fermi Golden rule, although these are valid only well away from the band gaps.

Atom lasers, coherent states, and coherence. I. Physically realizable ensembles of pure states

A laser, be it an optical laser or an atom laser, is an open quantum system that produces a coherent beam of bosons (photons or atoms, respectively). Far above threshold, the stationary state rho(ss) of the laser mode is a mixture of coherent-field states with random phase, or, equivalently, a Poissonian mixture of number states. This paper answers the question: can descriptions such as these, of rho(ss) as a stationary ensemble of pure states, be physically realized? Here physical realization is as defined previously by us [H. M. Wiseman and J. A. Vaccaro, Phys. Lett. A 250, 241 (1998)]: an ensemble of pure states for a particular system can be physically realized if, without changing the dynamics of the system, an experimenter can (in principle) know at any time that the system is in one of the pure-state members of the ensemble. Such knowledge can be obtained by monitoring the baths to which the system is coupled, provided that coupling is describable by a Markovian master equation. Using a family of master equations for the (atom) laser, we solve for the physically realizable (PR) ensembles. We find that for any finite self-energy chi of the bosons in the laser mode, the coherent-state ensemble is not PR; the closest one can come to it is an ensemble of squeezed states. This is particularly relevant for atom lasers, where the self-energy arising from elastic collisions is expected to be large. By contrast, the number-state ensemble is always PR. As the self-energy chi increases, the states in the PR ensemble closest to the coherent-state ensemble become increasingly squeezed. Nevertheless, there are values of chi for which states with well-defined coherent amplitudes are PR, even though the atom laser is not coherent (in the sense of having a Bose-degenerate output). We discuss the physical significance of this anomaly in terms of conditional coherence (and hence conditional Bose degeneracy).

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.

Limits to squeezing in the degenerate optical parametric oscillator

We develop a systematic theory of quantum fluctuations in the driven optical parametric oscillator, including the region near threshold. This allows us to treat the limits imposed by nonlinearities to quantum squeezing and noise reduction in this nonequilibrium quantum phase transition. In particular, we compute the squeezing spectrum near threshold and calculate the optimum value. We find that the optimal noise reduction occurs at different driving fields, depending on the ratio of damping rates. The largest spectral noise reductions are predicted to occur with a very high-Q second-harmonic cavity. 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.