Coherence time of a Bose-Einstein condensate

Temporal coherence is a fundamental property of macroscopic quantum systems, such as lasers in optics and Bose-Einstein condensates in atomic gases and it is a crucial issue for interferometry applications with light or matter waves. Whereas the laser is an "open" quantum system, ultracold atomic gases are weakly coupled to the environment and may be considered as isolated. The coherence time of a condensate is then intrinsic to the system and its derivation is out of the frame of laser theory. Using quantum kinetic theory, we predict that the interaction with non-condensed modes gradually smears out the condensate phase, with a variance growing as A t^2+B t+C at long times t, and we give a quantitative prediction for A, B and C. Whereas the coefficient A vanishes for vanishing energy fluctuations in the initial state, the coefficients B and C are remarkably insensitive to these fluctuations. The coefficient B describes a diffusive motion of the condensate phase that sets the ultimate limit to the condensate coherence time. We briefly discuss the possibility to observe the predicted phase spreading, also including the effect of particle losses.Comment: 17 pages, 8 figures; typos correcte

On the Inelastic Collapse of a Ball Bouncing on a Randomly Vibrating Platform

We study analytically the dynamics of a ball bouncing inelastically on a randomly vibrating platform, as a simple toy model of inelastic collapse. Of principal interest are the distributions of the number of flights n_f till the collapse and the total time \tau_c elapsed before the collapse. In the strictly elastic case, both distributions have power law tails characterised by exponents which are universal, i.e., independent of the details of the platform noise distribution. In the inelastic case, both distributions have exponential tails: P(n_f) ~ exp[-\theta_1 n_f] and P(\tau_c) ~ exp[-\theta_2 \tau_c]. The decay exponents \theta_1 and \theta_2 depend continuously on the coefficient of restitution and are nonuniversal; however as one approches the elastic limit, they vanish in a universal manner that we compute exactly. An explicit expression for \theta_1 is provided for a particular case of the platform noise distribution.Comment: 32 page

Exact and approximate many-body dynamics with stochastic one-body density matrix evolution

We show that the dynamics of interacting fermions can be exactly replaced by a quantum jump theory in the many-body density matrix space. In this theory, jumps occur between densities formed of pairs of Slater determinants, $D_{ab}=| \Phi_a > < \Phi_b |$, where each state evolves according to the Stochastic Schr\"odinger Equation (SSE) given in ref. \cite{Jul02}. A stochastic Liouville-von Neumann equation is derived as well as the associated Bogolyubov-Born-Green-Kirwood-Yvon (BBGKY) hierarchy. Due to the specific form of the many-body density along the path, the presented theory is equivalent to a stochastic theory in one-body density matrix space, in which each density matrix evolves according to its own mean field augmented by a one-body noise. Guided by the exact reformulation, a stochastic mean field dynamics valid in the weak coupling approximation is proposed. This theory leads to an approximate treatment of two-body effects similar to the extended Time-Dependent Hartree-Fock (Extended TDHF) scheme. In this stochastic mean field dynamics, statistical mixing can be directly considered and jumps occur on a coarse-grained time scale. Accordingly, numerical effort is expected to be significantly reduced for applications.Comment: 12 pages, 1 figur

Hybrid phase-space simulation method for interacting Bose fields

We introduce an approximate phase-space technique to simulate the quantum dynamics of interacting bosons. With the future goal of treating Bose-Einstein condensate systems, the method is designed for systems with a natural separation into highly occupied (condensed) modes and lightly occupied modes. The method self-consistently uses the Wigner representation to treat highly occupied modes and the positive-P representation for lightly occupied modes. In this method, truncation of higher-derivative terms from the Fokker-Planck equation is usually necessary. However, at least in the cases investigated here, the resulting systematic error, over a finite time, vanishes in the limit of large Wigner occupation numbers. We tested the method on a system of two interacting anharmonic oscillators, with high and low occupations, respectively. The Hybrid method successfully predicted atomic quadratures to a useful simulation time 60 times longer than that of the positive-P method. The truncated Wigner method also performed well in this test. For the prediction of the correlation in a quantum nondemolition measurement scheme, for this same system, the Hybrid method gave excellent agreement with the exact result, while the truncated Wigner method showed a large systematic error.Comment: 13 pages; 6 figures; references added; figures correcte

Analogue model of a FRW universe in Bose-Einstein condensates: Application of the classical field method

Analogue models of gravity have been motivated by the possibility of investigating phenomena not readily accessible in their cosmological counterparts. In this paper, we investigate the analogue of cosmological particle creation in a Friedmann-Robertson-Walker universe by numerically simulating a Bose-Einstein condensate with a time-dependent scattering length. In particular, we focus on a two-dimensional homogeneous condensate using the classical field method via the truncated Wigner approximation. We show that for various forms of the scaling function the particle production is consistent with the underlying theory in the long wavelength limit. In this context, we further discuss the implications of modified dispersion relations that arise from the microscopic theory of a weakly interacting Bose gas.Comment: 26 pages, 8 figure

Robust unravelings for resonance fluorescence

Monitoring the fluorescent radiation of an atom unravels the master equation evolution by collapsing the atomic state into a pure state which evolves stochastically. A robust unraveling is one that gives pure states that, on average, are relatively unaffected by the master equation evolution (which applies once the monitoring ceases). The ensemble of pure states arising from the maximally robust unraveling has been suggested to be the most natural way of representing the system [H.M. Wiseman and J.A. Vaccaro, Phys. Lett. A {\bf 250}, 241 (1998)]. We find that the maximally robust unraveling of a resonantly driven atom requires an adaptive interferometric measurement proposed by Wiseman and Toombes [Phys. Rev. A {\bf 60}, 2474 (1999)]. The resultant ensemble consists of just two pure states which, in the high driving limit, are close to the eigenstates of the driving Hamiltonian $\Omega\sigma_{x}/2$. This ensemble is the closest thing to a classical limit for a strongly driven atom. We also find that it is possible to reasonably approximate this ensemble using just homodyne detection, an example of a continuous Markovian unraveling. This has implications for other systems, for which it may be necessary in practice to consider only continuous Markovian unravelings.Comment: 12 pages including 5 .eps figures, plus one .jpg figur

Exact stochastic simulation of dissipation and non-Markovian effects in open quantum systems

The exact dynamics of a system coupled to an environment can be described by an integro-differential stochastic equation of its reduced density. The influence of the environment is incorporated through a mean-field which is both stochastic and non-local in time and where the standard two-times correlation functions of the environment appear naturally. Since no approximation is made, the presented theory incorporates exactly dissipative and non-Markovian effects. Applications to the spin-boson model coupled to a heat-bath with various coupling regimes and temperature show that the presented stochastic theory can be a valuable tool to simulate exactly the dynamics of open quantum systems. Links with stochastic Schroedinger equation method and possible extensions to "imaginary time" propagation are discussed.Comment: accepted for publication in Physical Review

Hyperfine interaction induced decoherence of electron spins in quantum dots

We investigate in detail, using both analytical and numerical tools, the decoherence of electron spins in quantum dots (QDs) coupled to a bath of nuclear spins in magnetic fields or with various initial bath polarizations, focusing on the longitudinal relaxation in low and moderate field/polarization regimes. An increase of the initial polarization of nuclear spin bath has the same effect on the decoherence process as an increase of the external magnetic field, namely, the decoherence dynamics changes from smooth decay to damped oscillations. This change can be observed experimentally for a single QD and for a double-QD setup. Our results indicate that substantial increase of the decoherence time requires very large bath polarizations, and the use of other methods (dynamical decoupling or control of the nuclear spins distribution) may be more practical for suppressing decoherence of QD-based qubits.Comment: Rev. Tex, 5 pages, 3 eps color figures, submitted to Phys. Rev.

Time evolution towards q-Gaussian stationary states through unified Ito-Stratonovich stochastic equation

We consider a class of single-particle one-dimensional stochastic equations which include external field, additive and multiplicative noises. We use a parameter $\theta \in [0,1]$ which enables the unification of the traditional It\^o and Stratonovich approaches, now recovered respectively as the $\theta=0$ and $\theta=1/2$ particular cases to derive the associated Fokker-Planck equation (FPE). These FPE is a {\it linear} one, and its stationary state is given by a $q$-Gaussian distribution with $q = \frac{\tau + 2M (2 - \theta)}{\tau + 2M (1 - \theta)}<3$, where $\tau \ge 0$ characterizes the strength of the confining external field, and $M \ge 0$ is the (normalized) amplitude of the multiplicative noise. We also calculate the standard kurtosis $\kappa_1$ and the $q$-generalized kurtosis $\kappa_q$ (i.e., the standard kurtosis but using the escort distribution instead of the direct one). Through these two quantities we numerically follow the time evolution of the distributions. Finally, we exhibit how these quantities can be used as convenient calibrations for determining the index $q$ from numerical data obtained through experiments, observations or numerical computations.Comment: 9 pages, 2 figure

Metallicity and Physical Conditions in the Magellanic Bridge

We present a new analysis of the diffuse gas in the Magellanic Bridge (RA>3h) based on HST/STIS E140M and FUSE spectra of 2 early-type stars lying within the Bridge and a QSO behind it. We derive the column densities of HI (from Ly\alpha), NI, OI, ArI, SiII, SII, and FeII of the gas in the Bridge. Using the atomic species, we determine the first gas-phase metallicity of the Magellanic Bridge, [Z/H]=-1.02+/-0.07 toward one sightline, and -1.7<[Z/H]<-0.9 toward the other one, a factor 2 or more smaller than the present-day SMC metallicity. Using the metallicity and N(HI), we show that the Bridge gas along our three lines of sight is ~70-90% ionized, despite high HI columns, logN(HI)=19.6-20.1. Possible sources for the ongoing ionization are certainly the hot stars within the Bridge, hot gas (revealed by OVI absorption), and leaking photons from the SMC and LMC. From the analysis of CII*, we deduce that the overall density of the Bridge must be low (<0.03-0.1 cm^-3). We argue that our findings combined with other recent observational results should motivate new models of the evolution of the SMC-LMC-Galaxy system.Comment: Accepted for publication in the Ap