Search

We consider the topological properties of the so-called central spin model that consists of $N$ identical spins coupled to a single distinguishable central spin which arises in physical systems such as circuit-QED and bosonic Josephson junctions coupled to an impurity atom. The model closely corresponds to the Su-Schrieffer-Heeger (SSH) model except that the chain of sites in the SSH model is replaced by a chain of states in Fock space specifying the magnetization. We find that the model accommodates topologically protected eigenstates that are `Bell-cat' states consisting of a Schr\"{o}dinger cat state of the $N$ spins that is maximally entangled with the central spin, and show how this state can be adiabatically created and moved along the chain by driving the central spin. We visualize the Bell-cat states by plotting their Wigner function and also show their robustness against random noise by solving the master equation for the density matrix., Comment: 10 pages, 6 figures

The question of the correct formulation for the momentum of light in a dielectric medium is typically referred to as the ``Abraham-Minkowski controversy". Experiments conducted to elucidate the issue have primarily focused on measuring forces and momentum transfers. In this work, we propose an interferometric approach using matter waves to measure the light-induced version of the He-McKellar-Wilkens (optical HMW) phase for a neutral atomic dipole in dynamical electromagnetic fields. Beginning from the action principle, we show that this geometric phase is directly related via the Euler-Lagrange equations of motion to the Abraham force lying at the heart of the controversy., Comment: 6 pages, 2 figures

Recent theoretical studies have predicted the existence of caustics in many-body quantum dynamics, where they manifest as extended regions of enhanced probability density that obey temporal and spatial scaling relations. Focusing on the transverse-field Ising model, we investigate the dynamics initiated by a local quench in a spin chain, resulting in outward-propagating excitations that create a distinct caustic pattern. We calculate the scaling of the first two maxima of the interference fringes dressing the caustic, finding a universal exponent of 2/3, associated with an Airy function catastrophe. We demonstrate that this property is universal in the entire paramagnetic phase of the model, and starts varying at the quantum phase transition (QPT). This robust scaling persists even under perturbations that break the integrability of the model. We additionally explore the effect of boundary conditions and find that open boundaries introduce significant edge effects, leading to complex interference patterns. Despite these edge-induced dynamics, the overall power-law scaling exponent remains robust. These findings highlight the potential of quantum caustics as a powerful diagnostic tool for QPTs, demonstrating resilience against integrability-breaking perturbations and boundary condition variations., Comment: 9 pages, 5 figures

We provide a new perspective on quantum dynamical phase transitions (DPTs) by explaining their origin in terms of caustics that form in the Fock space representation of the many-body state over time, using the fully connected transverse field Ising model as an example. In this way we establish a connection between DPTs in a quantum spin system and an everyday natural phenomenon: The dark band between the primary and seconday bows (caustics) in rainbows known as Alexander's dark band. The DPT occurs when the Loschmidt echo crosses the switching line between the evanescent tails of two back-to-back Airy functions that dress neighbouring fold caustics in Fock space and is the time-dependent analogue of what is seen as a function of angle in the sky. The structural stability and universal properties of caustics, as described mathematically by catastrophe theory, explains the generic occurrence of DPTs in the model and suggests that our analysis has wide applicability. Based on our thorough analytical understanding we propose a protocol which can be used to verify the existence of a DPT in a finite system experiment.

We perform a theoretical study of the coupled dynamics of two species of Bose-Einstein condensates (BECs) in a double well potential where both the tunneling and the interatomic interactions are driven periodically in time. The population difference between the wells of each species gives rise to a two dimensional lattice in Fock space with dimensions given by the number of atoms in each BEC. We use a Floquet analysis to derive an effective Hamiltonian that acts in this Fock space and find that it contains an artificial gauge field. This system simulates noninteracting particles in a tight binding lattice subject to an additional harmonic potential and vector potential. When the intra-species interactions are attractive there is a critical value at which the ground state of the Floquet operator undergoes a transition from a Gaussian state to a quantized vortex state in Fock space. The transition can be quantified in terms of the angular momentum as well as the entanglement entropy of the ground state with both showing sudden jumps as the intra-species interactions become stronger. The stability of the vortex state vanishes in the thermodynamic limit., Comment: 12 pages, 5 figures

In this paper we show how the quantum mechanics of the inverted harmonic oscillator can be mapped to the quantum mechanics of a particle in a super-critical inverse square potential. We demonstrate this by relating both of these systems to the Berry-Keating system with hamiltonian $H=(xp+px)/2$. It has long been appreciated that the quantum mechanics of the inverse square potential has an ambiguity in choosing a boundary condition near the origin and we show how this ambiguity is mapped to the inverted harmonic oscillator system. Imposing a boundary condition requires specifying a distance scale where it is applied and changes to this scale come with a renormalization group (RG) evolution of the boundary condition that ensures observables do not directly depend on the scale (which is arbitrary). Physical scales instead emerge as RG invariants of this evolution. The RG flow for the inverse square potential is known to follow limit cycles describing the discrete breaking of classical scale invariance in a simple example of a quantum anomaly, and we find that limit cycles also occur for the inverted harmonic oscillator. However, unlike the inverse square potential where the continuous scaling symmetry is explicit, in the case of the inverted harmonic oscillator it is hidden and occurs because the hamiltonian is part of a larger su(1,1) spectrum generating algebra. Our map does not require the boundary condition to be self-adjoint, as can be appropriate for systems that involve the absorption or emission of particles., Comment: 29 pages, 7 figures

Caustics are singularities that occur naturally in optical, hydrodynamic and quantum waves, giving rise to high amplitude patterns that can be described using catastrophe theory. In this paper we study caustics in a statistical field theory setting in the form of the sine-Gordon model that describes a variety of physical systems including coupled 1D superfluids. Specifically, we use classical field simulations to study the dynamics of two ultracold 1D Bose gases (quasi-condensates) that are suddenly coupled to each other and find that the resulting non-equilibrium dynamics are dominated by caustics. Thermal noise is included by sampling the initial states from a Boltzmann distribution for phononic excitations. We find that caustics pile up over time in both the number and phase difference observables leading to a characteristic non-thermal `circus tent' shaped probability distribution at long times., Comment: 30 pages, 15 figures. In this second version we have added a figure (Fig. 15) giving the probability distribution at long times after a quench for the phase difference variable

Waves propagating near an event horizon display interesting features including logarithmic phase singularities and caustics. We consider an acoustic horizon in a flowing Bose-Einstein condensate where the elementary excitations obey the Bogoliubov dispersion relation. In the hamiltonian ray theory the solutions undergo a broken pitchfork bifurcation near the horizon and one might therefore expect the associated wave structure to be given by a Pearcey function, this being the universal wave function that dresses catastrophes with two control parameters. However, the wave function is in fact an Airy-type function supplemented by a logarithmic phase term, a novel type of wave catastrophe. Similar wave functions arise in aeroacoustic flows from jet engines and also gravitational horizons if dispersion which violates Lorentz symmetry in the UV is included. The approach we take differs from previous authors in that we analyze the behaviour of the integral representation of the wave function using exponential coordinates. This allows for a different treatment of the branches that gives rise to an analysis based purely on saddlepoint expansions, which resolve the multiple real and complex waves that interact at the horizon and its companion caustic. We find that the horizon is a physical manifestation of a Stokes surface, marking the place where a wave is born, and that the horizon and the caustic do not in general coincide: the finite spatial region between them delineates a broadened horizon., Comment: 34 pages, 12 figures

Out of time ordered correlators (OTOCs) are useful tools for investigating foundational questions such as thermalization in closed quantum systems because they can potentially distinguish between integrable and nonintegrable dynamics. Here we discuss the properties of wavefronts of OTOCs by focusing on the region around the main wavefront at $x=v_{B}t$, where $v_{B}$ is the butterfly velocity. Using a Heisenberg spin model as an example, we find that the leading edge of a propagating Gaussian with the argument $-m(x)\left( x-v_B t \right)^2 +b(x)t$ gives an excellent fit to the region around $x=v_{B}t$ for both the free and chaotic cases. However, the scaling in these two regimes is very different: in the free case the coefficients $m(x)$ and $b(x)$ have an inverse power law dependence on $x$ whereas in the chaotic case they decay exponentially. We conjecture that this result is universal by using catastrophe theory to show that, on the one hand, the wavefront in the free case has to take the form of an Airy function and its local expansion shows that the power law scaling seen in the numerics holds rigorously, and on the other hand an exponential scaling of the OTOC wavefront must be a signature of nonintegrable dynamics. We find that the crossover between the two regimes is smooth and characterized by an S-shaped curve giving the lifting of Airy nodes as a function of a chaos parameter. This shows that the Airy form is qualitatively stable against weak chaos and consistent with the concept of a quantum Kolmogorov-Arnold-Moser theory., Comment: 9 pages (including references), 3 figures. This version includes new results on the stability of the Airy function wavefront as chaos is switched on gradually from zero to fully chaotic

We demonstrate that two-time correlation functions, which are generalizations of out-of-time-ordered correlators (OTOCs), can show 'false-flags' of chaos by exhibiting behaviour predicted by random matrix theory even in a system with classically regular dynamics. In particular, we analyze a system of bosons trapped in a double-well potential and probed by a quantum dot which is coupled to the bosons dispersively. This is an integrable system (considered both as separate parts and in total). Despite the continuous time evolution generated by the actual Hamiltonian, we find that the n-fold two-time correlation function for the probe describes an effective stroboscopic or Floquet dynamics whereby the bosons appear to be alternately driven by two different non-commuting Hamiltonians in a manner reminiscent of the Trotterized time evolution that occurs in digital quantum simulation. The classical limit of this effective dynamics can have a nonzero Lyapunov exponent, while the effective level statistics and return probability show traditional signatures of chaotic behaviour. In line with several other recent studies, this work highlights the fact that the behavior of OTOCs and their generalizations must be interpreted with some care., Comment: 14 pages, 7 figures

The attractive inverse square potential arises in a number of physical problems such as a dipole interacting with a charged wire, the Efimov effect, the Calgero-Sutherland model, near-horizon black hole physics and the optics of Maxwell fisheye lenses. Proper formulation of the inverse-square problem requires specification of a boundary condition (regulator) at the origin representing short-range physics not included in the inverse square potential and this generically breaks the Hamiltonian's continuous scale invariance in an elementary example of a quantum anomaly. The system's spectrum qualitatively changes at a critical value of the inverse-square coupling, and we here point out that the transition at this critical potential strength can be regarded as an example of a $\mathcal{PT}$ symmetry breaking transition. In particular, we use point particle effective field theory (PPEFT), as developed by Burgess et al [J. High Energy Phys., 2017(4):106, 2017], to characterize the renormalization group (RG) evolution of the boundary coupling under rescalings. While many studies choose boundary conditions to ensure the system is unitary, these RG methods allow us to systematically handle the richer case of nonunitary physics describing a source or sink at the origin (such as is appropriate for the charged wire or black hole applications). From this point of view the RG flow changes character at the critical inverse-square coupling, transitioning from a sub-critical regime with evolution between two real, unitary fixed points ($\mathcal{PT}$ symmetric phase) to a super-critical regime with imaginary, dissipative fixed points ($\mathcal{PT}$ symmetry broken phase) that represent perfect-sink and perfect-source boundary conditions, around which the flow executes limit-cycle evolution., Comment: 15 pages, 4 figures. This article is a contribution to the proceedings for the online seminar series on Pseudo-Hermitian Hamiltonians in Quantum Physics

By studying the effects of quadratic anisotropy and quartic perturbations on the two-dimensional harmonic oscillator, one arrives at a simple model termed here the Ince oscillator, whose analytic solutions are given in terms of Ince polynomials. This one model unifies diverse physical systems, including aberrated optical cavities that are shown to support Ince-Gauss beams as their modes, and the two-mode Bose-Hubbard dimer describing two coupled superfluids. The Ince oscillator model describes a topological transition which can have very different origins: in the optical case, which is fundamentally linear, it is driven by the ratio of astigmatic to spherical mirror aberrations, whereas in the superfluid case it is driven by the ratio of particle tunneling to interparticle interactions and corresponds to macroscopic quantum self trapping.

We describe a new class of nonequilibrium quantum many-body phenomena in the form of networks of caustics that dominate the many-body wavefunction in the semiclassical regime following a sudden quench. It includes the light cone-like propagation of correlations as a particular case. Caustics are singularities formed by the birth and death of waves and form a hierarchy of universal patterns whose natural mathematical description is via catastrophe theory. Examples in classical waves range from rainbows and gravitational lensing in optics to tidal bores and rogue waves in hydrodynamics. Quantum many-body caustics are discretized by second-quantization (``quantum catastrophes'') and live in Fock space which can potentially have many dimensions. We illustrate these ideas using the Bose Hubbard dimer and trimer models which are simple enough that the caustic structure can be elucidated from first principles and yet run the full range from integrable to nonintegrable dynamics. The dimer gives rise to discretized versions of fold and cusp catastrophes whereas the trimer allows for higher catastrophes including the codimension-3 hyperbolic and elliptic umbilics which are organized by, and projections of, an 8-dimensional corank-2 catastrophe known as $X_9$. These results describe a hitherto unrecognized form of universality in quantum dynamics organized by singularities that manifest as strong fluctuations in mode population probabilities., Comment: 30 pages, 20 figures

We calculate the out-of-time-ordered correlation function (OTOC) of a single impurity qubit coupled to fully a connected many-particle system such as a bosonic Josephson junction or spins with long-range interactions. In these systems the qubit OTOC can be used to detect both ground state and excited state quantum phase transitions (QPTs), making it a robust order parameter that is considerably more sensitive than the standard one-body correlation function. Finite size scaling exponents for an $N$ body system can also be accurately extracted from the long-time OTOC dynamics, however, for short times there is a discrepancy due to the fact that the qubit has not had enough time to couple to the larger system. Our results show that the OTOC of even the smallest probe is enough to diagnose a QPT in fully connected models but, like a continuous measurement, can still cause a backaction effect which leads to weakly chaotic dynamics and gradual information scrambling., Comment: 11 pages, 6 figures

We carry out a realistic, yet simple, calculation of the Casimir-Polder interaction in the presence of a metallic shield in order to aid the design of experiments to test non-Newtonian gravity. In particular, we consider a rubidium atom near a movable silicon slab with a gold film in between. We show that by moving the slab to various distances and making precise measurements of the force exerted on the atom, one could in principle discern the existence of short-range modifications to Newtonian gravity. This avoids the need for a patterned surface where calculations are much harder and for which the probe must be moved laterally at a fixed distance. We also briefly discuss the case where an atomic cloud undergoes Bloch oscillations within an optical lattice created by reflecting a laser off the shield. We find that our scheme has the potential to improve current constraints if relatively modest improvements in atom localisation in optical lattices are made.

The Hamiltonian Mean Field (HMF) model describes particles on a ring interacting via a cosine interaction, or equivalently, rotors coupled by infinite-range XY interactions. Conceived as a generic statistical mechanical model for long-range interactions such as gravity (of which the cosine is the first Fourier component), it has recently been used to account for self-organization in experiments on cold atoms with long-range optically mediated interactions. The significance of the HMF model lies in its ability to capture the universal effects of long-range interactions and yet be exactly solvable in the canonical ensemble. In this work we consider the quantum version of the HMF model in 1D and provide a classification of all possible stationary solutions of its generalized Gross-Pitaevskii equation (GGPE) which is both nonlinear and nonlocal. The exact solutions are Mathieu functions that obey a nonlinear relation between the wavefunction and the depth of the meanfield potential, and we identify them as bright solitons. Using a Galilean transformation these solutions can be boosted to finite velocity and are increasingly localized as the meanfield potential becomes deeper. In contrast to the usual local GPE, the HMF case features a tower of solitons, each with a different number of nodes. Our results suggest that long-range interactions support solitary waves in a novel manner relative to the short-range case., Comment: 10 pages, 4 figures, 7 pages of appendices

Polarizable atoms interacting with a charged wire do so through an inverse-square potential, $V = - g/r^2$. This system is known to realize scale invariance in a nontrivial way and to be subject to ambiguities associated with the choice of boundary condition at the origin, often termed the problem of `fall to the center'. Point-particle effective field theory (PPEFT) provides a systematic framework for determining the boundary condition in terms of the properties of the source residing at the origin. We apply this formalism to the charged-wire/polarizable-atom problem, finding a result that is not a self-adjoint extension because of absorption of atoms by the wire. We explore the RG flow of the complex coupling constant for the dominant low-energy effective interactions, finding flows whose character is qualitatively different when $g$ is above or below a critical value, $g_c$. Unlike the self-adjoint case, (complex) fixed points exist when $g> g_c$, which we show correspond to perfect absorber (or perfect emitter) boundary conditions. We describe experimental consequences for wire-atom interactions and the possibility of observing the anomalous breaking of scale invariance., Comment: 31 pages, 4 figures, minor changes in ordering of sections, new references added

We show that the light cone-like structures that form in spin chains after a quench are quantum caustics. Their natural description is in terms of catastrophe theory and this implies: 1) a hierarchy of light cone structures corresponding to the different catastrophes; 2) dressing by characteristic wave functions that obey scaling laws determined by the Arnol'd and Berry indices; 3) a network of vortex-antivortex pairs in space-time inside the cone. We illustrate the theory by giving explicit calculations for the transverse field Ising model and the XY model, finding fold catastrophes dressed by Airy functions and cusp catastrophes dressed by Pearcey functions; multisite correlation functions are described by higher catastrophes such as the hyperbolic umbilic. Furthermore, we find that the vortex pairs created inside the cone are sensitive to phase transitions in these spin models with their rate of production being determined by the dynamical critical exponent. More broadly, this work illustrates how catastrophe theory can be applied to singularities in quantum fields., Comment: Changed title (old version was: Light-cones and quantum caustics in quenched spin chains), added new figures and extended the paper to include XY model. 20 pages, 10 figures

Caustics occur widely in dynamics and take on shapes classified by catastrophe theory. At finite wavelengths they produce interference patterns containing networks of vortices (phase singularities). Here we investigate caustics in quantized fields, focusing on the collective dynamics of quantum spins. We show that, following a quench, caustics are generated in the Fock space amplitudes specifying the many-body configuration and which are accessible in experiments with cold atoms, ions or photons. The granularity of quantum fields removes all singularities, including phase singularities, converting point vortices into nonlocal vortices that annihilate in pairs as the quantization scale is increased. Furthermore, the continuous scaling laws of wave catastrophes are replaced by discrete versions. Such `quantum catastrophes' are expected to be universal dynamical features of quantized fields., Comment: The original title of this paper was "Morphology of a quantum catastrophe". The title was changed to reflect its relevance to collective spin dynamics. There are 4 figures

We use a master equation to study the dynamics of two coupled macroscopic quantum systems (e.g.\ a Josephson junction made of two Bose-Einstein condensates or two spin states of an ensemble of trapped ions) subject to a weak continuous measurement. If the coupling between the two systems is suddenly switched on the resulting dynamics leads to caustics (fold and cusp catastrophes) in the number-difference probability distribution, and at the same time the measurement gradually induces a quantum-to-classical transition. Decoherence is often invoked to help resolve paradoxes associated with macroscopic quantum mechanics, but here, on the contrary, caustics are well-behaved in the quantum (many-particle) theory and divergent in the classical (mean-field) theory. Caustics thus represent a breakdown of the classical theory towards which decoherence seems to inevitably lead. We find that measurement backaction plays a crucial role in softening the resulting singularities and calculate the modification to the Arnol'd index which governs the scaling of the caustic's amplitude with the number of atoms. The Arnol'd index acts as a critical exponent for the formation of singularities during quantum dynamics and its modification by the open nature of the system is analogous to the modification of the critical exponents of phase transitions occurring in open systems., Comment: 14 pages, 8 figures. New version has a fourth author and introduces the idea of the modification of the Arnold scaling exponent by decoherence

As part of the quest to uncover universal features of quantum dynamics, we study catastrophes that form in simple many-body wave functions after a quench, focusing on two-mode systems that include the two-site Bose Hubbard model, and under some circumstances optomechanical systems and the Dicke model. When the wave function is plotted in Fock space plus time certain characteristic structures generically appear that we identify as cusp caustics. In the vicinity of such a catastrophe the wave function takes on a universal form described by the Pearcey function and obeys scaling relations which depend on the total number of particles $N$. In the thermodynamic limit ($N \rightarrow \infty$) the cusp becomes singular, but at finite $N$ it is decorated by an interference pattern. This pattern contains an intricate network of vortex-antivortex pairs, initiating a theory of topological structures in Fock space. In the case where the quench takes the form of a $\delta$-kick we show how to analytically map the wave function onto the Pearcey function and hence obtain the scaling exponents for the size and position of the cusp, as well as those for the amplitude and characteristic length scales of its interference pattern., Comment: 16 pages, 8 figures

We perform a theoretical investigation into the classical and quantum dynamics of an optical field in a cavity containing a moving membrane ("membrane-in-the-middle" set-up). Our approach is based on the Maxwell wave equation, and complements previous studies based on an effective Hamiltonian. The analysis shows that for slowly moving and weakly reflective membranes the dynamics can be approximated by unitary, first-order-in-time evolution given by an effective Schr\"{o}dinger-like equation with a Hamiltonian that does not depend on the membrane speed. This approximate theory is the one typically adopted in cavity optomechanics and we develop a criterion for its validity. However, in more general situations the full second-order wave equation predicts light dynamics which do not conserve energy, giving rise to parametric amplification (or reduction) that is forbidden under first order dynamics and can be considered to be the classical counterpart of the dynamical Casimir effect. The case of a membrane moving at constant velocity can be mapped onto the Landau-Zener problem but with additional terms responsible for field amplification. Furthermore, the nature of the adiabatic regime is rather different from the ordinary Schr\"odinger case, since mode amplitudes need not be constant even when there are no transitions between them. The Landau-Zener problem for a field is therefore richer than in the standard single-particle case. We use the work-energy theorem applied to the radiation pressure on the membrane as a self-consistency check for our solutions of the wave equation and as a tool to gain an intuitive understanding of energy pumped into/out of the light field by the motion of the membrane., Comment: 27 pages, 15 figures

The experimental realization of quantum-degenerate Bose gases made of atoms with sizeable magnetic dipole moments has created a new type of fluid, known as a quantum ferrofluid, which combines the extraordinary properties of superfluidity and ferrofluidity. A hallmark of superfluids is that they are constrained to rotate through vortices with quantized circulation. In quantum ferrofluids the long-range dipolar interactions add new ingredients by inducing magnetostriction and instabilities, and also affect the structural properties of vortices and vortex lattices. Here we give a review of the theory of vortices in dipolar Bose-Einstein condensates, exploring the interplay of magnetism with vorticity and contrasting this with the established behaviour in non-dipolar condensates. We cover single vortex solutions, including structure, energy and stability, vortex pairs, including interactions and dynamics, and also vortex lattices. Our discussion is founded on the mean-field theory provided by the dipolar Gross-Pitaevskii equation, ranging from analytic treatments based on the Thomas-Fermi (hydrodynamic) and variational approaches to full numerical simulations. Routes for generating vortices in dipolar condensates are discussed, with particular attention paid to rotating condensates, where surface instabilities drive the nucleation of vortices, and lead to the emergence of rich and varied vortex lattice structures. We also present an outlook, including potential extensions to degenerate Fermi gases, quantum Hall physics, toroidal systems and the Berezinskii-Kosterlitz-Thouless transition., Comment: 42 pages with 23 Figures

We study the stability, form and interaction of single and multiple dark solitons in quasi-one-dimensional dipolar Bose-Einstein condensates. The solitons are found numerically as stationary solutions in the moving frame of a non-local Gross Pitaevskii equation, and characterized as a function of the key experimental parameters, namely the ratio of the dipolar atomic interactions to the van der Waals interactions, the polarization angle and the condensate width. The solutions and their integrals of motion are strongly affected by the phonon and roton instabilities of the system. Dipolar matter-wave dark solitons propagate without dispersion, and collide elastically away from these instabilities, with the dipolar interactions contributing an additional repulsion or attraction to the soliton-soliton interaction. However, close to the instabilities, the collisions are weakly dissipative., Comment: 13 pages, 9 figures

We analyze the optomechanics of an atomic Bose-Einstein condensate interacting with the optical lattice inside a laser-pumped optical cavity and subject to a uniform bias force such as gravity. An atomic wave packet in a tilted lattice undergoes Bloch oscillations; in a cavity the backaction of the atoms on the light leads to a time-dependent modulation of the intracavity lattice at the Bloch frequency. When the Bloch frequency is on the order of the cavity damping rate we find transport of the atoms either up or down the lattice. The transport dynamics can be interpreted as a manifestation of dynamical backaction-induced sideband damping/amplification of the optomechanical Bloch oscillator. Depending on the sign of the pump-cavity detuning, atoms are transported either with or against the bias force accompanied by an up- or down-conversion of the frequency of the pump laser light. We also evaluate the prospects for using the optomechanical Bloch oscillator to make continuous measurements of forces by reading out the Bloch frequency. In this context we establish the significant result that the optical spring effect is absent and the Bloch frequency is not modified by the backaction., Comment: 24 pages, 5 figures, published version

We study the family of static and moving dark solitons in quasi-one-dimensional dipolar Bose-Einstein condensates, exploring their modified form and interactions. The density dip of the soliton acts as a giant anti-dipole which adds a non-local contribution to the conventional local soliton-soliton interaction. We map out the stability diagram as a function of the strength and polarization direction of the atomic dipoles, identifying both roton and phonon instabilities. Away from these instabilities, the solitons collide elastically. Varying the polarization direction relative to the condensate axis enables tuning of this non-local interaction between repulsive and attractive; the latter case supports unusual dark soliton bound states. Remarkably, these bound states are themselves shown to behave like solitons, emerging unscathed from collisions with each other., Comment: 6 pages, 5 figures

We use fidelity susceptibility to calculate quantum critical scaling exponents for a system consisting of $N$ identical bosons interacting with a single impurity atom in a double well potential (bosonic Josephson junction). Above a critical value of the boson-impurity interaction energy there is a spontaneous breaking of $\mathbb{Z}_2$ symmetry corresponding to a second order quantum phase transition from a balanced to an imbalanced number of particles in either the left or right hand well. We show that the exponents match those in the Lipkin-Meshkov-Glick and Dicke models suggesting that the impurity model is in the same universality class. The phase transition can be interpreted as a measurement of the position of the impurity by the bosons., Comment: 9 pages, 5 figures

We consider the "membrane in the middle" optomechanical model consisting of a laser pumped cavity which is divided in two by a flexible membrane that is partially transmissive to light and subject to radiation pressure. Steady state solutions at the mean-field level reveal that there is a critical strength of the light-membrane coupling above which there is a symmetry breaking bifurcation where the membrane spontaneously acquires a displacement either to the left or the right. This bifurcation bears many of the signatures of a second order phase transition and we compare and contrast it with that found in the Dicke model. In particular, by studying limiting cases and deriving dynamical critical exponents using the fidelity susceptibility method, we argue that the two models share very similar critical behaviour. For example, the obtained critical exponents indicate that they fall within the same universality class. Away from the critical regime we identify, however, some discrepancies between the two models. Our results are discussed in terms of experimentally relevant parameters and we evaluate the prospects for realizing Dicke-type physics in these systems., Comment: 14 pages, 6 figures

We predict that an atomic Bose-Einstein condensate strongly coupled to an intracavity optical lattice can undergo resonant tunneling and directed transport when a constant and uniform bias force is applied. The bias force induces Bloch oscillations, causing amplitude and phase modulation of the lattice which resonantly modifies the site-to-site tunneling. For the right choice of parameters a net atomic current is generated. The transport velocity can be oriented oppositely to the bias force, with its amplitude and direction controlled by the detuning between the pump laser and the cavity. The transport can also be enhanced through imbalanced pumping of the two counter-propagating running wave cavity modes. Our results add to the cold atoms quantum simulation toolbox, with implications for quantum sensing and metrology., Comment: Published version: 5 pages, 4 figures; Supplementary Material included

We present vortex solutions for the homogeneous two-dimensional Bose-Einstein condensate featuring dipolar atomic interactions, mapped out as a function of the dipolar interaction strength (relative to the contact interactions) and polarization direction. Stable vortex solutions arise in the regimes where the fully homogeneous system is stable to the phonon or roton instabilities. Close to these instabilities, the vortex profile differs significantly from that of a vortex in a nondipolar quantum gas, developing, for example, density ripples and an anisotropic core. Meanwhile, the vortex itself generates a mesoscopic dipolar potential which, at distance, scales as 1/r^2 and has an angular dependence which mimics the microscopic dipolar interaction., Comment: 11 pages, 5 figures. Proceedings of the 22nd International Laser Physics Workshop (Prague, July 2013

We perform a theoretical study into how dipole-dipole interactions modify the properties of superfluid vortices within the context of a two-dimensional atomic Bose gas of co-oriented dipoles. The reduced density at a vortex acts like a giant anti-dipole, changing the density profile and generating an effective dipolar potential centred at the vortex core whose most slowly decaying terms go as $1/\rho^2$ and $\ln(\rho)/\rho^3$. These effects modify the vortex-vortex interaction which, in particular, becomes anisotropic for dipoles polarized in the plane. Striking modifications to vortex-vortex dynamics are demonstrated, i.e. anisotropic co-rotation dynamics and the suppression of vortex annihilation., Comment: PRL accepted, 6 pages, 5 figures

We study a model describing $N$ identical bosonic atoms trapped in a double-well potential together with a single impurity atom, comparing and contrasting it throughout with the Dicke model. As the boson-impurity coupling strength is varied, there is a symmetry-breaking pitchfork bifurcation which is analogous to the quantum phase transition occurring in the Dicke model. Through stability analysis around the bifurcation point, we show that the critical value of the coupling strength has the same dependence on the parameters as the critical coupling value in the Dicke model. We also show that, like the Dicke model, the mean-field dynamics go from being regular to chaotic above the bifurcation and macroscopic excitations of the bosons are observed. Overall, the boson-impurity system behaves like a poor man's version of the Dicke model., Comment: 17 pages, 16 figures

We extend our theory of Bloch oscillations of cold atoms inside an optical cavity [ Venkatesh et al. Phys. Rev. A 80 063834 (2009)] to include the effects of quantum noise arising from coupling to external modes. The noise acts as a form of quantum measurement backaction by perturbing the coupled dynamics of the atoms and the light. We take it into account by solving the Heisenberg-Langevin equations for linearized fluctuations about the atomic and optical mean fields and examine how this influences the signal-to-noise ratio of a measurement of external forces using this system. In particular, we investigate the effects of changing the number of atoms, the intracavity lattice depth, and the atom-light coupling strength, and show how resonances between the Bloch oscillation dynamics and the quasiparticle spectrum have a strong influence on the signal-to-noise ratio, as well as heating effects. One of the hurdles we overcome in this paper is the proper treatment of fluctuations about time-dependent mean fields in the context of cold-atom cavity QED., Comment: 25 pages, 14 figures

We study rainbow (fold) and cusp catastrophes that form in Fock space following a quench in a Bose Josephson junction. In the Gross-Pitaevskii mean-field theory the rainbows are singular caustics, but in the second-quantized theory a Poisson resummation of the wave function shows that they are described by well behaved Airy functions. The structural stability of these Fock space caustics against variations in the initial conditions and Hamiltonian evolution is guaranteed by catastrophe theory. We also show that the long-time dynamics are ergodic. Our results are relevant to the question posed by Berry [M.V. Berry, Nonlinearity 21, T19 (2008)]: are there circumstances when it is necessary to second-quantize wave theory in order to avoid singularities?, Comment: 5 pages, 3 figures. New references added plus a section on experimental realizability

We compare and contrast the mean-field and many-body properties of a Bose-Einstein condensate trapped in a double well potential with a single impurity atom. The mean-field solutions display a rich structure of bifurcations as parameters such as the boson-impurity interaction strength and the tilt between the two wells are varied. In particular, we study a pitchfork bifurcation in the lowest mean-field stationary solution which occurs when the boson-impurity interaction exceeds a critical magnitude. This bifurcation, which is present for both repulsive and attractive boson-impurity interactions, corresponds to the spontaneous formation of an imbalance in the number of particles between the two wells. If the boson-impurity interaction is large, the bifurcation is associated with the onset of a Schroedinger cat state in the many-body ground state. We calculate the coherence and number fluctuations between the two wells, and also the entanglement entropy between the bosons and the impurity. We find that the coherence can be greatly enhanced at the bifurcation., Comment: 19 pages, 17 figures. The second version contains minor corrections and some better figures (thicker lines)

We analyze the evolution of an electromagnetic field inside a double cavity when the difference in length between the two cavities is changed, e.g. by translating the common mirror. We find that this allows photons to be moved deterministically from one cavity to the other. We are able to obtain the conditions for adiabatic transfer by first mapping the Maxwell wave equation for the electric field onto a Schroedinger-like wave equation, and then using the Landau-Zener result for the transition probability at an avoided crossing. Our analysis reveals that this mapping only rigorously holds when the two cavities are weakly coupled (i.e. in the regime of a highly reflective common mirror), and that, generally speaking, care is required when attempting a hamiltonian description of cavity electrodynamics with time-dependent boundary conditions., Comment: 24 pages, 18 figures. Version 2 includes a new section (Sec. VIII) on the regimes of validity of the Schroedinger-like equations and also of the adiabatic approximation, together with a new figure (Fig. 10). The discussion section (Sec. XI) has also been enhanced

We calculate the band structure of ultracold atoms located inside a laser-driven optical cavity. For parameters where the atom-cavity system exhibits bistability, the atomic band structure develops loop structures akin to the ones predicted for Bose-Einstein condensates in ordinary (non-cavity) optical lattices. However, in our case the nonlinearity derives from the cavity back-action rather than from direct interatomic interactions. We find both bi- and tri-stable regimes associated with the lowest band, and show that the multistability we observe can be analyzed in terms of swallowtail catastrophes. Dynamic and energetic stability of the mean-field solutions is also discussed, and we show that the bistable solutions have, as expected, one unstable and two stable branches. The presence of loops in the atomic band structure has important implications for proposals concerning Bloch oscillations of atoms inside optical cavities [Peden et al., Phys. Rev. A 80, 043803 (2009), Prasanna Venkatesh et al., Phys. Rev. A 80, 063834 (2009)]., Comment: 26 pages, 22 figures

Engineering of synthetic magnetic flux in Bose-Einstein condensates [Lin et al., Nature {\bf 462}, 628 (2009)] has prospects for attaining the high vortex densities necessary to emulate the fractional quantum Hall effect. We analytically establish the hydrodynamical behaviour of a condensate in a uniform synthetic magnetic field, including its density and velocity profile. Importantly, we find that the onset of vortex nucleation observed experimentally corresponds to a dynamical instability in the hydrodynamical solutions and reveal other routes to instability and anticipated vortex nucleation., Comment: 4 pages with 2 figures

We present a general method for obtaining the exact static solutions and collective excitation frequencies of a trapped Bose-Einstein condensate (BEC) with dipolar atomic interactions in the Thomas-Fermi regime. The method incorporates analytic expressions for the dipolar potential of an arbitrary polynomial density profile, thereby reducing the problem of handling non-local dipolar interactions to the solution of algebraic equations. We comprehensively map out the static solutions and excitation modes, including non-cylindrically symmetric traps, and also the case of negative scattering length where dipolar interactions stabilize an otherwise unstable condensate. The dynamical stability of the excitation modes gives insight into the onset of collapse of a dipolar BEC. We find that global collapse is consistently mediated by an anisotropic quadrupolar collective mode, although there are two trapping regimes in which the BEC is stable against quadrupole fluctuations even as the ratio of the dipolar to s-wave interactions becomes infinite. Motivated by the possibility of fragmented BEC in a dipolar Bose gas due to the partially attractive interactions, we pay special attention to the scissors modes, which can provide a signature of superfluidity, and identify a long-range restoring force which is peculiar to dipolar systems. As part of the supporting material for this paper we provide the computer program used to make the calculations, including a graphical user interface., Comment: 23 pages, 11 figures

We compare the classical (mean-field) dynamics with the quantum dynamics of atomic Bose-Einstein condensates in double-well potentials. The quantum dynamics are computed using a simple scheme based upon the Raman-Nath equations. Two different methods for exciting a non-equilbrium state are considered: an asymmetry between the wells which is suddenly removed, and a periodic time oscillating asymmetry. The first method generates wave packets that lead to collapses and revivals of the expectation values of the macroscopic variables, and we calculate the time scale for these revivals. The second method permits the excitation of a single energy eigenstate of the many-particle system, including Schroedinger cat states. We also discuss a band theory interpretation of the energy level structure of an asymmetric double-well, thereby identifying analogies to Bloch oscillations and Bragg resonances. Both the Bloch and Bragg dynamics are purely quantum and are not contained in the mean-field treatment., Comment: 31 pages, 14 figures

We present a theoretical analysis of dilute gas Bose-Einstein condensates with dipolar atomic interactions under rotation in elliptical traps. Working in the Thomas-Fermi limit, we employ the classical hydrodynamic equations to first derive the rotating condensate solutions and then consider their response to perturbations. We thereby map out the regimes of stability and instability for rotating dipolar Bose-Einstein condensates and in the latter case, discuss the possibility of vortex lattice formation. We employ our results to propose several novel routes to induce vortex lattice formation in a dipolar condensate., Comment: 12 pages with 6 figures

Cold atoms in an optical lattice execute Bloch-Zener oscillations when they are accelerated. We have performed a theoretical investigation into the case when the optical lattice is the intra-cavity field of a driven Fabry-Perot resonator. When the atoms oscillate inside the resonator, we find that their back-action modulates the phase and intensity of the light transmitted through the cavity. We solve the coupled atom-light equations self-consistently and show that, remarkably, the Bloch period is unaffected by this back-action. The transmitted light provides a way to observe the oscillation continuously, allowing high precision measurements to be made with a small cloud of atoms., Comment: 5 pages, 2 figures. Updated version including cavity heating effects

We investigate the collapse of a trapped dipolar Bose-Einstein condensate. This is performed by numerical simulations of the Gross-Pitaevskii equation and the novel application of the Thomas-Fermi hydrodynamic equations to collapse. We observe regimes of both global collapse, where the system evolves to a highly elongated or flattened state depending on the sign of the dipolar interaction, and local collapse, which arises due to dynamically unstable phonon modes and leads to a periodic arrangement of density shells, disks or stripes. In the adiabatic regime, where ground states are followed, collapse can occur globally or locally, while in the non-adiabatic regime, where collapse is initiated suddenly, local collapse commonly occurs. We analyse the dependence on the dipolar interactions and trap geometry, the length and time scales for collapse, and relate our findings to recent experiments., Comment: In this version (the published version) we have slightly rewritten the manuscript in places and have corrected some typos. 15 pages and 13 figures

We investigate the time taken for global collapse by a dipolar Bose-Einstein condensate. Two semi-analytical approaches and exact numerical integration of the mean-field dynamics are considered. The semi-analytical approaches are based on a Gaussian ansatz and a Thomas-Fermi solution for the shape of the condensate. The regimes of validity for these two approaches are determined, and their predictions for the collapse time revealed and compared with numerical simulations. The dipolar interactions introduce anisotropy into the collapse dynamics and predominantly lead to collapse in the plane perpendicular to the axis of polarization., Comment: 5 pages with 3 figures

We derive the criteria for the Thomas-Fermi regime of a dipolar Bose-Einstein condensate in cigar, pancake and spherical geometries. This also naturally gives the criteria for the mean-field one- and two-dimensional regimes. Our predictions, including the Thomas-Fermi density profiles, are shown to be in excellent agreement with numerical solutions. Importantly, the anisotropy of the interactions has a profound effect on the Thomas-Fermi/low-dimensional criteria., Comment: 5 pages, 2 figures

We calculate the hydrodynamic solutions for a dilute Bose-Einstein condensate with long-range dipolar interactions in a rotating, elliptical harmonic trap, and analyse their dynamical stability. The static solutions and their regimes of instability vary non-trivially on the strength of the dipolar interactions. We comprehensively map out this behaviour, and in particular examine the experimental routes towards unstable dynamics, which, in analogy to conventional condensates, may lead to vortex lattice formation. Furthermore, we analyse the centre of mass and breathing modes of a rotating dipolar condensate., Comment: 4 pages, including 2 figures

We study a gaseous Bose-Einstein condensate with laser-induced dipole-dipole interactions using the Hartree-Fock-Bogoliubov theory within the Popov approximation. The dipolar interactions introduce long-range atom-atom correlations, which manifest themselves as increased depletion at momenta similar to that of the laser wavelength, as well as a "roton" dip in the excitation spectrum. Surprisingly, the roton dip and the corresponding peak in the depletion are enhanced by raising the temperature above absolute zero., Comment: 10 pages, 6 figures

A gaseous Bose-Einstein condensate (BEC) irradiated by a far off-resonance laser has long-range interatomic correlations caused by laser-induced dipole-dipole interactions. These correlations, which are tunable via the laser intensity and frequency, can produce a `roton' minimum in the excitation spectrum--behavior reminiscent of the strongly correlated superfluid liquid helium II., Comment: 6 pages, includes 3 figures