Non-diffusive phase spreading of a Bose-Einstein condensate at finite temperature

We show that the phase of a condensate in a finite temperature gas spreads linearly in time at long times rather than in a diffusive way. This result is supported by classical field simulations, and analytical calculations which are generalized to the quantum case under the assumption of quantum ergodicity in the system. This super-diffusive behavior is intimately related to conservation of energy during the free evolution of the system and to fluctuations of energy in the prepared initial state.Comment: 16 pages, 7 figure

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

Flow graphs: interweaving dynamics and structure

The behavior of complex systems is determined not only by the topological organization of their interconnections but also by the dynamical processes taking place among their constituents. A faithful modeling of the dynamics is essential because different dynamical processes may be affected very differently by network topology. A full characterization of such systems thus requires a formalization that encompasses both aspects simultaneously, rather than relying only on the topological adjacency matrix. To achieve this, we introduce the concept of flow graphs, namely weighted networks where dynamical flows are embedded into the link weights. Flow graphs provide an integrated representation of the structure and dynamics of the system, which can then be analyzed with standard tools from network theory. Conversely, a structural network feature of our choice can also be used as the basis for the construction of a flow graph that will then encompass a dynamics biased by such a feature. We illustrate the ideas by focusing on the mathematical properties of generic linear processes on complex networks that can be represented as biased random walks and also explore their dual consensus dynamics.Comment: 4 pages, 1 figur

Simulations of thermal Bose fields in the classical limit

We demonstrate that the time-dependent projected Gross-Pitaevskii equation derived earlier [Davis, et al., J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. We find that this equation will evolve randomised initial wave functions to equilibrium, and compare our numerical data to the predictions of a gapless, second-order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)]. We find that we can determine the temperature of the equilibrium state when this theory is valid. Outside the range of perturbation theory we describe how to measure the temperature of our simulations. We also determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. As the Gross-Pitaevskii equation is non-perturbative, we expect that it can describe the correct thermal behaviour of a Bose gas as long as all relevant modes are highly occupied.Comment: 15 pages, 12 figures, revtex4, follow up to Phys. Rev. Lett. 87 160402 (2001). v2: Modified after referee comments. Extra data added to two figures, section on temperature determination expande

The stochastic Gross-Pitaevskii equation II

We provide a derivation of a more accurate version of the stochastic Gross-Pitaevskii equation, as introduced by Gardiner et al. (J. Phys. B 35,1555,(2002). The derivation does not rely on the concept of local energy and momentum conservation, and is based on a quasi-classical Wigner function representation of a "high temperature" master equation for a Bose gas, which includes only modes below an energy cutoff E_R that are sufficiently highly occupied (the condensate band). The modes above this cutoff (the non-condensate band) are treated as being essentially thermalized. The interaction between these two bands, known as growth and scattering processes, provide noise and damping terms in the equation of motion for the condensate band, which we call the stochastic Gross-Pitaevskii equation. This approach is distinguished by the control of the approximations made in its derivation, and by the feasibility of its numerical implementation.Comment: 24 pages of LaTeX, one figur

Hot Streaks in Artistic, Cultural, and Scientific Careers

The hot streak, loosely defined as winning begets more winnings, highlights a specific period during which an individual's performance is substantially higher than her typical performance. While widely debated in sports, gambling, and financial markets over the past several decades, little is known if hot streaks apply to individual careers. Here, building on rich literature on lifecycle of creativity, we collected large-scale career histories of individual artists, movie directors and scientists, tracing the artworks, movies, and scientific publications they produced. We find that, across all three domains, hit works within a career show a high degree of temporal regularity, each career being characterized by bursts of high-impact works occurring in sequence. We demonstrate that these observations can be explained by a simple hot-streak model we developed, allowing us to probe quantitatively the hot streak phenomenon governing individual careers, which we find to be remarkably universal across diverse domains we analyzed: The hot streaks are ubiquitous yet unique across different careers. While the vast majority of individuals have at least one hot streak, hot streaks are most likely to occur only once. The hot streak emerges randomly within an individual's sequence of works, is temporally localized, and is unassociated with any detectable change in productivity. We show that, since works produced during hot streaks garner significantly more impact, the uncovered hot streaks fundamentally drives the collective impact of an individual, ignoring which leads us to systematically over- or under-estimate the future impact of a career. These results not only deepen our quantitative understanding of patterns governing individual ingenuity and success, they may also have implications for decisions and policies involving predicting and nurturing individuals with lasting impact

Microcanonical temperature for a classical field: application to Bose-Einstein condensation

We show that the projected Gross-Pitaevskii equation (PGPE) can be mapped exactly onto Hamilton's equations of motion for classical position and momentum variables. Making use of this mapping, we adapt techniques developed in statistical mechanics to calculate the temperature and chemical potential of a classical Bose field in the microcanonical ensemble. We apply the method to simulations of the PGPE, which can be used to represent the highly occupied modes of Bose condensed gases at finite temperature. The method is rigorous, valid beyond the realms of perturbation theory, and agrees with an earlier method of temperature measurement for the same system. Using this method we show that the critical temperature for condensation in a homogeneous Bose gas on a lattice with a UV cutoff increases with the interaction strength. We discuss how to determine the temperature shift for the Bose gas in the continuum limit using this type of calculation, and obtain a result in agreement with more sophisticated Monte Carlo simulations. We also consider the behaviour of the specific heat.Comment: v1: 9 pages, 5 figures, revtex 4. v2: additional text in response to referee's comments, now 11 pages, to appear in Phys. Rev.

Numerical study of the coupled time-dependent Gross-Pitaevskii equation: Application to Bose-Einstein condensation

We present a numerical study of the coupled time-dependent Gross-Pitaevskii equation, which describes the Bose-Einstein condensate of several types of trapped bosons at ultralow temperature with both attractive and repulsive interatomic interactions. The same approach is used to study both stationary and time-evolution problems. We consider up to four types of atoms in the study of stationary problems. We consider the time-evolution problems where the frequencies of the traps or the atomic scattering lengths are suddenly changed in a stable preformed condensate. We also study the effect of periodically varying these frequencies or scattering lengths on a preformed condensate. These changes introduce oscillations in the condensate which are studied in detail. Good convergence is obtained in all cases studied.Comment: 9 pages, 10 figures, accepted in Physical Review

Dynamical thermalization and vortex formation in stirred 2D Bose-Einstein condensates

We present a quantum mechanical treatment of the mechanical stirring of Bose-Einstein condensates using classical field techniques. In our approach the condensate and excited modes are described using a Hamiltonian classical field method in which the atom number and (rotating frame) energy are strictly conserved. We simulate a T = 0 quasi-2D condensate perturbed by a rotating anisotropic trapping potential. Vacuum fluctuations in the initial state provide an irreducible mechanism for breaking the initial symmetries of the condensate and seeding the subsequent dynamical instability. Highly turbulent motion develops and we quantify the emergence of a rotating thermal component that provides the dissipation necessary for the nucleation and motional-damping of vortices in the condensate. Vortex lattice formation is not observed, rather the vortices assemble into a spatially disordered vortex liquid state. We discuss methods we have developed to identify the condensate in the presence of an irregular distribution of vortices, determine the thermodynamic parameters of the thermal component, and extract damping rates from the classical field trajectories.Comment: 22 pages, 15 figures. v2: Minor refinements made at suggestion of referee. Discussion of other treatments revised. To appear in Phys. Rev.

Causality and defect formation in the dynamics of an engineered quantum phase transition in a coupled binary Bose-Einstein condensate

Continuous phase transitions occur in a wide range of physical systems, and provide a context for the study of non-equilibrium dynamics and the formation of topological defects. The Kibble-Zurek (KZ) mechanism predicts the scaling of the resulting density of defects as a function of the quench rate through a critical point, and this can provide an estimate of the critical exponents of a phase transition. In this work we extend our previous study of the miscible-immiscible phase transition of a binary Bose-Einstein condensate (BEC) composed of two hyperfine states in which the spin dynamics are confined to one dimension [J. Sabbatini et al., Phys. Rev. Lett. 107, 230402 (2011)]. The transition is engineered by controlling a Hamiltonian quench of the coupling amplitude of the two hyperfine states, and results in the formation of a random pattern of spatial domains. Using the numerical truncated Wigner phase space method, we show that in a ring BEC the number of domains formed in the phase transitions scales as predicted by the KZ theory. We also consider the same experiment performed with a harmonically trapped BEC, and investigate how the density inhomogeneity modifies the dynamics of the phase transition and the KZ scaling law for the number of domains. We then make use of the symmetry between inhomogeneous phase transitions in anisotropic systems, and an inhomogeneous quench in a homogeneous system, to engineer coupling quenches that allow us to quantify several aspects of inhomogeneous phase transitions. In particular, we quantify the effect of causality in the propagation of the phase transition front on the resulting formation of domain walls, and find indications that the density of defects is determined during the impulse to adiabatic transition after the crossing of the critical point.Comment: 23 pages, 10 figures. Minor corrections, typos, additional referenc