Dynamics of Trapped Fermion Clouds
|
|
- Anna Melton
- 6 years ago
- Views:
Transcription
1 Laser Physics, Vol. 14, No. 2, 2004, pp Original Text Copyright 2004 by Astro, Ltd. Copyright 2004 by MAIK Nauka /Interperiodica (Russia). PHYSICS OF COLD TRAPPED ATOMS Dynamics of Trapped Fermion Clouds F. Toschi 1, 2, P. Vignolo 3, P. Capuzzi 3, S. Succi 1, 3, and M. P. Tosi 3 1 Istituto per le Applicazioni del Calcolo, CNR, Viale del Policlinico 137, Roma, Italy 2 INFM, Sezione di Roma Tor Vergata, Via della Ricerca Scientifica 1, Roma, Italy 3 NEST-INFM and Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, Italy vignolo@sns.it Received July 31, 2003 Abstract Experiments at JILA on the collisions between two oscillating spin-polarized components of a Fermi gas of 40 K atoms have shown that this setup is an important tool for studying the dynamics of quantum gases. We have developed advanced numerical methods to solve the Vlasov Landau equations for a binary mixture of atomic Fermi gases inside a harmonic trap at very low temperatures. In particular, a locally adaptive importance-sampling technique allows us to handle collisional interactions several orders of magnitude more efficiently than in standard Monte Carlo calculations. These methods have been applied to study the dynamical and transport properties of a fermion mixture in the presence of mean-field interactions and of collisions between its components for system parameters ranging from the collisionless to the hydrodynamic regime. 1. INTRODUCTION Since the achievement of Bose Einstein condensation [1], a main goal of experimentalists has been to cool a gas of fermionic atoms down to an expected superfluid regime. The cooling process is an evaporative process where at each step the most energetic atoms are expelled from the trap and the others rethermalize through collisions. However, at very low temperatures, the main collisional events should occur in the s-wave channel but are inhibited by the Pauli principle in a spin-polarized Fermi gas. A way to circumvent this problem is to add a second fermionic component, which is polarized in another isospin state inside the magnetic trap. Experiments at JILA [2] have studied the collisions between two oscillating spin-polarized components of a Fermi gas of 40 K atoms and have given evidence for a transition from collisionless (zerosound) to hydrodynamic (first-sound) behavior through measuring the damping time τ. This becomes very long at both small and large values of the estimated collision rate. Following earlier numerical studies of the dynamics of various types of atomic clouds [3 5], Geist and Kennedy [6] described the evaporative cooling process in a two-component Fermi gas from semiclassical Boltzmann equations as functions only of the energy. A more comprehensive numerical study of the dynamics of such a system in coordinate and momentum space has been realized by Toschi et al. [7] by solving the Vlasov Landau equations (VLE) for the fermionic Wigner distributions. We have examined the collisional properties of the two-component Fermi gas as functions of both temperature T and rate Γ q of quantum collisions, which are counted step by step in the simulation. In this paper, we present in detail the strategy used in the numerical solution of the VLE and report some illustrative results. We treat the quantum fluid using a particle-dynamics approach [8] accounting for both mean-field interactions and collisions between the two components. Upon inclusion of collisional events at very low temperature, the calculation would grind to a halt, basically through a saturation of phase space resulting in a vanishing efficiency of the Monte Carlo sampling. We have developed a locally adaptive importance-sampling technique, which allows us to handle collisional interactions several orders of magnitude faster than in standard Monte Carlo techniques. Moreover, since we focus on the JILA setup [2], where axial symmetry is maintained during the experiment, we can use a code in which the angular degree of freedom is taken into account via an effective weight. We find that, even though most collisions become forbidden both classically and by the Pauli principle as temperature is lowered, the few collisions that can occur involve particles increasingly clustered around the Fermi level and suffice to drive the gas from the collisionless to the hydrodynamic regime. 2. THE MODEL AND ITS NUMERICAL SOLUTION The purpose of our numerical code is to solve the VLE for two interacting Fermi gases in the presence of collisions at finite temperature T down to T ~ 0.2T F, where T F is the Fermi degeneracy temperature VLE The two fermionic components in external potentials (r) are described by distribution functions V ext 302
2 DYNAMICS OF TRAPPED FERMION CLOUDS 303 f (j) (r, p, t) with j = 1 or 2. These obey the kinetic equations t f + p m --- r f r U p f = C 12 [ f ], (1) where the Hartree Fock (HF) effective potential is U (j) (r, t) V ext (r) + gn (r, t), with j denoting the species different from j. Here, we have set = 1, g = 2πa/m r, with a being the s-wave scattering length between two atoms and m r being their reduced mass, and n (j) (r, t) is the spatial density given by integration of f (j) (r, p, t) over momentum degrees of freedom. Collisions between atoms of the same spin can be neglected at low temperature, so that in Eq. (1) the term C 12 involves only collisions between particles with different polarizations: C 12 [ f ] = 22π ( ) 4 g 2 /V 3 Δ p Δ ε f ( j f 2 f ) 3 f 4 f ( j f ) [ 2 f 3 f 4 ], p 2, p 3, p 4 f with f (j) f (j) (r, p, t), 1 f (j), f (j) (r, p i, t), and 1. V is the volume occupied by the gas, and the factors Δ p and Δ ε are the usual delta functions accounting for conservation of momentum and energy, with the energies given by p 2 + U (j) i /2m j (r, t) Numerical Solution The numerical procedure by which the VLE are solved in time consists of three basic steps: (i) initialization of the fermionic distributions, (ii) propagation in phase space driven by external and mean-field forces, and (iii) propagation in phase space through collisional events Initializing the VLE. The initial phase-space distributions at equilibrium inside harmonic traps are generated from the HF expression f eq ( rp, ) = β p2 + U ( r) μ 1 exp + 1 2m j (2) where β = 1/k B T and μ (j) is the chemical potential of species j [9]. The particle densities entering U (j) (r) are determined self-consistently together with the chemical potentials. This is done by means of a standard rootfinding subroutine. We exploit the axial symmetry of the system by introducing the angularly integrated particle densities (j) (r, z) = 2πrn (j) (r, z) defined on a {r, z} grid in cylindrical coordinates. To evaluate them numerically, we, move to a particle-in-cell description by locating a number (j) (r, z)δrδz of fermions inside each cell of volume ΔrΔz [5]. Low statistical noise is achieved by representing each fermion by means of a number N q of computational particles ( quarks ). To generate the initial momentum distribution, the ith quark is located at point {p ir, p iθ, p iz } in momentum space by using a Monte Carlo sampling and making sure that each cell of volume h 3 is occupied by no more than N q quarks. This control in 3D phase space is transferred to 2D by imposing a maximum number wn q of quarks in the 2D cell {Δr, Δz, Δp r, Δp z } of volume h 2. Here, the weight w = int(2p F r/ ) takes into account that the number of available cells in 2D depends on the radial position and on the number of particles through the Fermi momentum p F = 2mk B T F. The use of N q quarks to represent a single atom is the analog of the splitting techniques widely used in kinetic Monte Carlo studies [10]. In our work on quantum gases, it has a further beneficial effect on the acceptance rate of the Monte Carlo procedure used to implement the Pauli exclusion principle Propagating the VLE. Taking our first experiment as a specific example, in the propagation step, the two clouds are rigidly displaced from the center of their traps along the z direction and start evolving in time in the {r, z} plane by performing oscillations at their respective frequencies. The VLE are advanced in time by a standard particle-in-cell (PIC) method [11] in a modified Verlet timemarching scheme on a grid with mesh spacing dx > vdt, where v is a typical particle velocity and dt is the time step. This inequality is dictated by the accuracy and stability of the propagation step [5]. We obtain the following set of discrete algebraic equations: r i ( t + dt) = r i () t + v ri ()dt t + z i ( t + dt) = z i () t + v zi ()dt t + v ri ( t + dt) = v ri () t + a ri ()dt t v zi ( t + dt) = v zi () t + a zi ()dt t a ri ()dt t 2 /2 a zi ()dt t 2 /2 (3) for i = 1 to N. Here, a r, a z and v r, v z are the two components of the acceleration and of the velocity along the {r, z} coordinates. The algorithm is standard except for (i) the decoupling of the variable p θ from the equations of motion, exploiting the fact that, in terms of symmetry, the average value of p θ at each value of r is zero and does not change in time, and (ii) the specification of the self-consistent coupling arising from the forces due to the density gradients. The grid-forces which cause the accelerations are evaluated at the discrete particle locations by means of a bilinear cloud-in-cell interpolator [5]. With the force/acceleration fields transferred to the par-
3 304 TOSCHI et al. ticle locations, only the collisions need setting in order to march the VLE in time The Collision step. Collisions are tracked on a grid of mesh spacing λ B (the de Broglie wavelength) both along r and z. λ B is smaller than the particle mean free path l and larger than the spacing dl of the propagation mesh (l > λ B > dl). The first inequality enhances the statistical accuracy of the collision step, whereas coarse graining with respect to the propagation step avoids the need for Pauli-principle constraints in the Lagrangian evolution, at least down to 0.2T F. In this step, again we exploit axial symmetry: we make the approximation that the angular momentum of each quark is left unchanged by collisional events and take into account the third dimension in the exclusion principle by suppressing collisions whose final states are occupied by wn q quarks. The number of classically probable collisions between all possible pairs of quarks belonging to the two species in each cell of volume dv = 2πrdrdz is evaluated at each computational step as dn coll dt = v dv ij σ ij. (4) Here, v ij is the magnitude of the relative speed of quark i of species 1 and quark j of species 2 and σ ij is the corresponding cross section σ ij = σ/n q, with σ being the atom atom cross section σ = πa 2. If dn coll < 1, the collision probability is accumulated over the subsequent time steps until an integer number di coll int(dn coll ) of collisions occurs. The collision probability becomes smaller than the classical one after multiplication by the quantum suppression factor 1 N(r, z, p r, p z )/(wn q ) due to the occupancy of the final state, and the efficiency of the Monte Carlo sampling drops with decreasing temperature. To enhance the acceptance rate by two to three orders of magnitude at the lowest temperatures, we proceed as follows. First, within each spatial cell, the pairs of particles are ordered according to the value of their classical transition probability and pairs with classical probability below a given threshold are filtered out. The threshold is dynamically adjusted cell-by-cell in order to guarantee a correct supply of di coll collisions at each time step. Each particle is then allowed to collide only with the partner which maximizes the product v ij σ ij, so that the original pool of N 1 N 2 collisions is cut down to N 1 collisions. This strategy has allowed us to circumvent the numerical attrition problems due to Pauli blocking down to 0.2T F. 3. RESULTS AND CONCLUDING REMARKS As a first application of the numerical method, we have considered a system of 200 magnetically trapped ij 40 K atoms represented by a total number of quarks [7]. As in the JILA experiments [2], the atoms are equally shared among two hyperfine states (m f = 9/2 and m f = 7/2) inside harmonic traps with somewhat different longitudinal frequencies (ω 9/2 = 2π 19.8 s 1 and ω 7/2 = 2π s 1 ). The two clouds after initialization are rigidly displaced from the trap centers and in the absence onteractions would keep oscillating at their respective trap frequencies without damping. Mean-field interactions would only modify these bare frequencies. This picture is profoundly changed when collisions are allowed to occur as the two clouds meet in their trajectories. We vary the collision rate Γ q at each given temperature by changing the scattering length a, thus mimicking the exploitation of a Feshbach resonance [12]. As in the JILA experiments [2], with increasing Γ q there is an increasing transfer of relative momentum from one cloud to the other, up to the point where the relative motion of the two clouds is suppressed and they start oscillating together at the same frequency. This transition from collisionless-ballistic motions (for Knudsen number Kn ) to a collisional-hydrodynamic regime (for Kn 0), as driven by varying the scattering length at each temperature, is displayed both in the plot of the dipole-mode frequencies of the two clouds in Fig. 1a and in the plot of the damping rate γ = 1/τ of the relative motion of the centers of mass z cm (t) in Fig. 1b. It should be noticed that, since the number of atoms used in the simulation (N = 200) is much smaller than in the experiments (N exp 10 6 ), we have upscaled the scattering cross section σ by a factor of N exp /N ~ 10 4 in order to keep the same Knudsen numbers. In fact, the effects of statistical errors due to the limited number of particles in the simulation are best controlled by letting the simulated and physical system share the same Knudsen numbers, namely, the ratio of the scattering length to the macroscopic size of the experimental sample. We can argue that physical results are not affected by such rescaling as long as the diluteness condition, nσ 3/2 1, remains fulfilled. In Fig. 1a, the oscillation frequencies ω j of the two clouds at two values of T/T F, with k B T F = ω 9/2 (6N 9/2 ) 1/3, have been obtained by fitting z cm (t) with the functions cos(ω j t)e γt. At very low Γ q, the dipole-mode frequencies are given by the corresponding trap frequencies with a small shift due to mean-field interactions. At intermediate values of Γ q, the data points show large fluctuations due to the fact that in this region just a few collisions can drastically alter the motion of the clouds. The trend towards a locking of the two dipole modes at large Γ q is very clear, and the location of the locking is identified with reasonable accuracy.
4 DYNAMICS OF TRAPPED FERMION CLOUDS 305 ω, s (a) γ, s 1 10 (b) 122 T/T F = 0.3 T/T F = 1.0 T/T F = 1.0 T/T F = Fig. 1. (a) The oscillation frequencies ω (in units of s 1 ) and (b) the damping coefficient γ (in units of s 1 ) as functions of the quantum collision rate Γ q (in units of s 1 ) for the two components of the gas at various temperatures. The horizontal dashed lines in (a) show the bare trap frequencies. The damping rate γ is shown in Fig. 1b at two values of T/T F and has been obtained from the correlation functions φ (j) (t) = z cm ( t' ) z cm ( t' + t) between the magnitude of the turning points z cm, which decay exponentially as exp( γt) at long times. The rate is essentially the same for the two components and increases linearly with Γ q in the collisionless regime, 1 while in the collisional one it scales as Γ q. It is also seen from Fig. 1 that the locking transition shifts to lower Γ q as temperature is lowered (see Fig. 2). This effect is a consequence of Fermi statistics: at lower temperature, the collisions involve particles in a narrower region around the Fermi level and a smaller number of collisions is needed to produce locking of the two species. The transition to hydrodynamic behavior can thus be driven by purely thermal means, taking advantage of the increasing importance of a decreasing set of T/T F Fig. 2. Log log plot of the collision rate Γ q at mode locking (in units of s 1 ) as a function of temperature T (in units of T F ). The fitting curve (dashed line) is (T/T F ) 0.40 ± strategic collisions between particles in states clustered around the Fermi level. To quantify this effect, we have performed a fit to our data and obtained a power law behavior of the form Γ q (T/T F ) 0.40 ± This is shown in log log scale in Fig. 2. The role of trap anisotropy in the transition from the collisionless to the hydrodynamic regime continues to be explored both experimentally and numerically. Moreover, the effects of collisions during expansion provide decisive information about the nature of the trapped state. In view of these facts, we are at present developing a fully 3D concurrent code to solve the Vlasov Landau equations. This code is allowing us to explore regimes with complete anisotropy of the trapping potential and thus to remove all symmetry constraints during, e.g., the collisional step. In addition, preliminary tests suggest that larger numbers of particles (N = 10 5 ) can be handled [13]. These further studies should allow us to investigate strategies by which one may circumvent Pauli blocking in reaching down to lower temperatures under general trapping situations. ACKNOWLEDGMENTS This work was supported by INFM through the PRA-Photonmatter Initiative. M.P.T. acknowledges the hospitality of the Abdus Salam International Center for Theoretical Physics during the preparation of the manuscript. REFERENCES 1. M. H. Anderson, J. R. Ensher, M. R. Matthews, et al., Science 269, 198 (1995); K. B. Davis, M.-O. Mewes, M. R. Andrews, et al., Phys. Rev. Lett. 75, 3969 (1995); C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995).
5 306 TOSCHI et al. 2. S. D. Gensemer and D. S. Jin, Phys. Rev. Lett. 87, (2001); B. DeMarco and D. S. Jin, Phys. Rev. Lett. 88, (2002). 3. H. Wu and C. J. Foot, J. Phys. B 29, L321 (1996). 4. T. Nikuni, E. Zaremba, and A. Griffin, Phys. Rev. Lett. 83, 10 (1999). 5. P. Vignolo, M. L. Chiofalo, S. Succi, and M. P. Tosi, J. Comput. Phys. 182, 368 (2002). 6. W. Geist and T. A. B. Kennedy, Phys. Rev. A 65, (2002). 7. F. Toschi, P. Vignolo, S. Succi, and M. P. Tosi, Phys. Rev. A 67, (R) (2003). 8. G. A. Bird, Molecular Gas Dynamics and the Direct Simulation of Gas Flows (Oxford Univ. Press, Oxford, 1994). 9. M. Amoruso, I. Meccoli, A. Minguzzi, and M. P. Tosi, Eur. Phys. J. D 7, 441 (1999). 10. P. K. MacKeown, Stochastic Simulation in Physics (Springer, Singapore, 1997). 11. R. Hockney and J. Eastwood, Computer Simulation Using Particles (McGraw-Hill, New York, 1981). 12. T. Loftus, C. A. Regal, C. Ticknor, et al., Phys. Rev. Lett. 88, (2002). 13. F. Toschi, P. Capuzzi, P. Vignolo, et al. (in press).
Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)
Dynamics of trapped two-component Fermi gas : temperature dependence of the transition from the collisionless to the collisional regime Toschi, F.; Vignolo, P.; Succi, S.; Tosi, M.P. Published in: Physical
More informationarxiv: v1 [cond-mat.stat-mech] 30 Aug 2007
Journal of Low Temperature Physics manuscript No. (will be inserted by the editor) arxiv:0708.414v1 [cond-mat.stat-mech] 30 Aug 007 Luca Salasnich 1 and Flavio Toigo Shell Effects in the First Sound Velocity
More informationIntroduction to Bose-Einstein condensation 4. STRONGLY INTERACTING ATOMIC FERMI GASES
1 INTERNATIONAL SCHOOL OF PHYSICS "ENRICO FERMI" Varenna, July 1st - July 11 th 2008 " QUANTUM COHERENCE IN SOLID STATE SYSTEMS " Introduction to Bose-Einstein condensation 4. STRONGLY INTERACTING ATOMIC
More informationDynamics of a Trapped Bose-Condensed Gas at Finite Temperatures
Dynamics of a Trapped Bose-Condensed Gas at Finite Temperatures Tetsuro Nikuni, Tokyo University of Science Collaborators: Allan Griffin, University of Toronto (Canada) Eugene Zaremba, Queen s University
More informationUltracold Fermi and Bose Gases and Spinless Bose Charged Sound Particles
October, 011 PROGRESS IN PHYSICS olume 4 Ultracold Fermi Bose Gases Spinless Bose Charged Sound Particles ahan N. Minasyan alentin N. Samoylov Scientific Center of Applied Research, JINR, Dubna, 141980,
More informationarxiv:cond-mat/ v1 28 Jan 2003
Three-Fluid Description of the Sympathetic Cooling of a Boson-Fermion Mixture M. Wouters, J. Tempere, J. T. Devreese Departement Natuurkunde, Universiteit Antwerpen, Universiteitsplein, B260 arxiv:cond-mat/030544v
More informationBose-Einstein condensation of lithium molecules and studies of a strongly interacting Fermi gas
Bose-Einstein condensation of lithium molecules and studies of a strongly interacting Fermi gas Wolfgang Ketterle Massachusetts Institute of Technology MIT-Harvard Center for Ultracold Atoms 3/4/04 Workshop
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationOptimization of evaporative cooling
PHYSICAL REVIEW A VOLUME 55, NUMBER 5 MAY 1997 Optimization of evaporative cooling C. A. Sackett, C. C. Bradley, and R. G. Hulet Physics Department and Rice Quantum Institute, Rice University, Houston,
More informationMonte Carlo Simulation of Bose Einstein Condensation in Traps
Monte Carlo Simulation of Bose Einstein Condensation in Traps J. L. DuBois, H. R. Glyde Department of Physics and Astronomy, University of Delaware Newark, Delaware 19716, USA 1. INTRODUCTION In this paper
More informationThe phases of matter familiar for us from everyday life are: solid, liquid, gas and plasma (e.f. flames of fire). There are, however, many other
1 The phases of matter familiar for us from everyday life are: solid, liquid, gas and plasma (e.f. flames of fire). There are, however, many other phases of matter that have been experimentally observed,
More informationIntroduction to Cold Atoms and Bose-Einstein Condensation. Randy Hulet
Introduction to Cold Atoms and Bose-Einstein Condensation Randy Hulet Outline Introduction to methods and concepts of cold atom physics Interactions Feshbach resonances Quantum Gases Quantum regime nλ
More informationBEC of 6 Li 2 molecules: Exploring the BEC-BCS crossover
Institut für Experimentalphysik Universität Innsbruck Dresden, 12.10. 2004 BEC of 6 Li 2 molecules: Exploring the BEC-BCS crossover Johannes Hecker Denschlag The lithium team Selim Jochim Markus Bartenstein
More informationQuantum Quantum Optics Optics VII, VII, Zakopane Zakopane, 11 June 09, 11
Quantum Optics VII, Zakopane, 11 June 09 Strongly interacting Fermi gases Rudolf Grimm Center for Quantum Optics in Innsbruck University of Innsbruck Austrian Academy of Sciences ultracold fermions: species
More informationFrom BEC to BCS. Molecular BECs and Fermionic Condensates of Cooper Pairs. Preseminar Extreme Matter Institute EMMI. and
From BEC to BCS Molecular BECs and Fermionic Condensates of Cooper Pairs Preseminar Extreme Matter Institute EMMI Andre Wenz Max-Planck-Institute for Nuclear Physics and Matthias Kronenwett Institute for
More informationFermi Condensates ULTRACOLD QUANTUM GASES
Fermi Condensates Markus Greiner, Cindy A. Regal, and Deborah S. Jin JILA, National Institute of Standards and Technology and University of Colorado, and Department of Physics, University of Colorado,
More informationFrom laser cooling to BEC First experiments of superfluid hydrodynamics
From laser cooling to BEC First experiments of superfluid hydrodynamics Alice Sinatra Quantum Fluids course - Complement 1 2013-2014 Plan 1 COOLING AND TRAPPING 2 CONDENSATION 3 NON-LINEAR PHYSICS AND
More informationReference for most of this talk:
Cold fermions Reference for most of this talk: W. Ketterle and M. W. Zwierlein: Making, probing and understanding ultracold Fermi gases. in Ultracold Fermi Gases, Proceedings of the International School
More informationarxiv:quant-ph/ v2 5 Feb 2001
Understanding the production of dual BEC with sympathetic cooling G. Delannoy, S. G. Murdoch, V. Boyer, V. Josse, P. Bouyer and A. Aspect Groupe d Optique Atomique Laboratoire Charles Fabry de l Institut
More informationIntroduction to cold atoms and Bose-Einstein condensation (II)
Introduction to cold atoms and Bose-Einstein condensation (II) Wolfgang Ketterle Massachusetts Institute of Technology MIT-Harvard Center for Ultracold Atoms 7/7/04 Boulder Summer School * 1925 History
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 19 May 2004
arxiv:cond-mat/0405440v1 [cond-mat.soft] 19 May 2004 Condensate localization in a quasi-periodic structure Y. Eksioglu, P. Vignolo and M.P. Tosi NEST-INFM and Scuola Normale Superiore, Piazza dei Cavalieri
More information84 Quantum Theory of Many-Particle Systems ics [Landau and Lifshitz (198)] then yields the thermodynamic potential so that one can rewrite the statist
Chapter 6 Bosons and fermions at finite temperature After reviewing some statistical mechanics in Sec. 6.1, the occupation number representation is employed to derive some standard results for noninteracting
More informationBose-Einstein Condensate: A New state of matter
Bose-Einstein Condensate: A New state of matter KISHORE T. KAPALE June 24, 2003 BOSE-EINSTEIN CONDENSATE: A NEW STATE OF MATTER 1 Outline Introductory Concepts Bosons and Fermions Classical and Quantum
More informationA Mixture of Bose and Fermi Superfluids. C. Salomon
A Mixture of Bose and Fermi Superfluids C. Salomon Enrico Fermi School Quantum Matter at Ultralow Temperatures Varenna, July 8, 2014 The ENS Fermi Gas Team F. Chevy, Y. Castin, F. Werner, C.S. Lithium
More informationLecture 3. Bose-Einstein condensation Ultracold molecules
Lecture 3 Bose-Einstein condensation Ultracold molecules 66 Bose-Einstein condensation Bose 1924, Einstein 1925: macroscopic occupation of the lowest energy level db h 2 mk De Broglie wavelength d 1/3
More informationLandau damping of transverse quadrupole oscillations of an elongated Bose-Einstein condensate
PHYSICAL REVIEW A 67, 053607 2003 Landau damping of transverse quadrupole oscillations of an elongated Bose-Einstein condensate M. Guilleumas 1 and L. P. Pitaevskii 2,3 1 Departament d Estructura i Constituents
More informationThe Gross-Pitaevskii Equation A Non-Linear Schrödinger Equation
The Gross-Pitaevskii Equation A Non-Linear Schrödinger Equation Alan Aversa December 29, 2011 Abstract Summary: The Gross-Pitaevskii equation, also called the non-linear Schrödinger equation, describes
More informationThermodynamic Measurements in a Strongly Interacting Fermi Gas
J Low Temp Phys (2009) 154: 1 29 DOI 10.1007/s10909-008-9850-2 Thermodynamic Measurements in a Strongly Interacting Fermi Gas Le Luo J.E. Thomas Received: 25 July 2008 / Accepted: 12 October 2008 / Published
More informationNanoKelvin Quantum Engineering
NanoKelvin Quantum Engineering Few x 10 5 Yb atoms 250mm 400 nk 250 nk < 200 nk Control of atomic c.m. position and momentum. Today: Bose-Fermi double superfluid Precision BEC interferometry Ultracold
More informationAnisotropic Spin Diffusion in Trapped Boltzmann Gases
Anisotropic Spin Diffusion in Trapped Boltzmann Gases W. J. Mullin 1 and R. J. Ragan 2 1 Physics Department, Hasbrouck Laboratory University of Massachusetts, Amherst, MA 01003 USA 2 Physics Department,
More informationBEC Vortex Matter. Aaron Sup October 6, Advisor: Dr. Charles Hanna, Department of Physics, Boise State University
BEC Vortex Matter Aaron Sup October 6, 006 Advisor: Dr. Charles Hanna, Department of Physics, Boise State University 1 Outline 1. Bosons: what are they?. Bose-Einstein Condensation (BEC) 3. Vortex Formation:
More informationLow-dimensional Bose gases Part 1: BEC and interactions
Low-dimensional Bose gases Part 1: BEC and interactions Hélène Perrin Laboratoire de physique des lasers, CNRS-Université Paris Nord Photonic, Atomic and Solid State Quantum Systems Vienna, 2009 Introduction
More informationAchieving steady-state Bose-Einstein condensation
PHYSICAL REVIEW A VOLUME 57, NUMBER 3 MARCH 1998 Achieving steady-state Bose-Einstein condensation J. Williams, R. Walser, C. Wieman, J. Cooper, and M. Holland JILA and Department of Physics, University
More informationUltra-cold gases. Alessio Recati. CNR INFM BEC Center/ Dip. Fisica, Univ. di Trento (I) & Dep. Physik, TUM (D) TRENTO
Ultra-cold gases Alessio Recati CNR INFM BEC Center/ Dip. Fisica, Univ. di Trento (I) & Dep. Physik, TUM (D) TRENTO Lectures L. 1) Introduction to ultracold gases Bosonic atoms: - From weak to strong interacting
More informationQuantum Properties of Two-dimensional Helium Systems
Quantum Properties of Two-dimensional Helium Systems Hiroshi Fukuyama Department of Physics, Univ. of Tokyo 1. Quantum Gases and Liquids 2. Bose-Einstein Condensation 3. Superfluidity of Liquid 4 He 4.
More informationUltracold Fermi Gases with unbalanced spin populations
7 Li Bose-Einstein Condensate 6 Li Fermi sea Ultracold Fermi Gases with unbalanced spin populations Nir Navon Fermix 2009 Meeting Trento, Italy 3 June 2009 Outline Introduction Concepts in imbalanced Fermi
More informationNon equilibrium Ferromagnetism and Stoner transition in an ultracold Fermi gas
Non equilibrium Ferromagnetism and Stoner transition in an ultracold Fermi gas Gareth Conduit, Ehud Altman Weizmann Institute of Science See: Phys. Rev. A 82, 043603 (2010) and arxiv: 0911.2839 Disentangling
More informationBEC AND MATTER WAVES an overview Allan Griffin, University of Toronto
BEC AND MATTER WAVES an overview Allan Griffin, University of Toronto The discovery of Bose-Einstein condensation ( BEC ) in 1995 in dilute, ultracold trapped atomic gases is one of the most exciting developments
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationStrongly correlated systems in atomic and condensed matter physics. Lecture notes for Physics 284 by Eugene Demler Harvard University
Strongly correlated systems in atomic and condensed matter physics Lecture notes for Physics 284 by Eugene Demler Harvard University January 25, 2011 2 Chapter 12 Collective modes in interacting Fermi
More informationA study of the BEC-BCS crossover region with Lithium 6
A study of the BEC-BCS crossover region with Lithium 6 T.Bourdel, L. Khaykovich, J. Cubizolles, J. Zhang, F. Chevy, M. Teichmann, L. Tarruell, S. Kokkelmans, Christophe Salomon Theory: D. Petrov, G. Shlyapnikov,
More informationShock waves in the unitary Fermi gas
Shock waves in the unitary Fermi gas Luca Salasnich Dipartimento di Fisica e Astronomia Galileo Galilei, Università di Padova Banff, May 205 Collaboration with: Francesco Ancilotto and Flavio Toigo Summary.
More informationBose-Bose mixtures in confined dimensions
Bose-Bose mixtures in confined dimensions Francesco Minardi Istituto Nazionale di Ottica-CNR European Laboratory for Nonlinear Spectroscopy 22nd International Conference on Atomic Physics Cairns, July
More informationBCS Pairing Dynamics. ShengQuan Zhou. Dec.10, 2006, Physics Department, University of Illinois
BCS Pairing Dynamics 1 ShengQuan Zhou Dec.10, 2006, Physics Department, University of Illinois Abstract. Experimental control over inter-atomic interactions by adjusting external parameters is discussed.
More informationCold fermions, Feshbach resonance, and molecular condensates (II)
Cold fermions, Feshbach resonance, and molecular condensates (II) D. Jin JILA, NIST and the University of Colorado I. Cold fermions II. III. Feshbach resonance BCS-BEC crossover (Experiments at JILA) $$
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationUltracold molecules - a new frontier for quantum & chemical physics
Ultracold molecules - a new frontier for quantum & chemical physics Debbie Jin Jun Ye JILA, NIST & CU, Boulder University of Virginia April 24, 2015 NIST, NSF, AFOSR, ARO Ultracold atomic matter Precise
More informationMixtures of ultracold gases: Fermi sea and Bose-Einstein condensate of Lithium isotopes
Mixtures of ultracold gases: Fermi sea and Bose-Einstein condensate of Lithium isotopes Florian Schreck To cite this version: Florian Schreck. Mixtures of ultracold gases: Fermi sea and Bose-Einstein condensate
More informationStanding Sound Waves in a Bose-Einstein Condensate
Standing Sound Waves in a Bose-Einstein Condensate G.A.R. van Dalum Debye Institute - Universiteit Utrecht Supervisors: P. van der Straten, A. Groot and P.C. Bons July 17, 2012 Abstract We study two different
More informationarxiv:cond-mat/ v1 [cond-mat.other] 19 Dec 2005
Released momentum distribution of a Fermi gas in the BCS-BEC crossover arxiv:cond-mat/5246v [cond-mat.other] 9 Dec 25 M.L. Chiofalo, S. Giorgini 2,3 and M. Holland 2 INFM and Classe di Scienze, Scuola
More informationSuperfluidity and Superconductivity
Superfluidity and Superconductivity These are related phenomena of flow without resistance, but in very different systems Superfluidity: flow of helium IV atoms in a liquid Superconductivity: flow of electron
More informationEffective Field Theory and Ultracold Atoms
Effective Field Theory and Ultracold Atoms Eric Braaten Ohio State University support Department of Energy Air Force Office of Scientific Research Army Research Office 1 Effective Field Theory and Ultracold
More informationContents Ultracold Fermi Gases: Properties and Techniques Index
V Contents 1 Ultracold Fermi Gases: Properties and Techniques 1 Selim Jochim 1.1 Introduction 1 1.2 Ultracold Fermions in a Trap 2 1.2.1 Ideal Fermi Gas 3 1.3 Preparing an Ultracold Fermi Gas 6 1.4 Interactions
More informationBEC in one dimension
BEC in one dimension Tilmann John 11. Juni 2013 Outline 1 one-dimensional BEC 2 theoretical description Tonks-Girardeau gas Interaction exact solution (Lieb and Liniger) 3 experimental realization 4 conclusion
More informationThe non-interacting Bose gas
Chapter The non-interacting Bose gas Learning goals What is a Bose-Einstein condensate and why does it form? What determines the critical temperature and the condensate fraction? What changes for trapped
More informationarxiv: v1 [cond-mat.quant-gas] 18 Sep 2009
Sound propagation in a Bose-Einstein condensate at finite temperatures arxiv:0909.3455v1 [cond-mat.quant-gas] 18 Sep 2009 R. Meppelink, S. B. Koller, and P. van der Straten 1 1 Atom Optics and Ultrafast
More informationBCS-BEC Crossover. Hauptseminar: Physik der kalten Gase Robin Wanke
BCS-BEC Crossover Hauptseminar: Physik der kalten Gase Robin Wanke Outline Motivation Cold fermions BCS-Theory Gap equation Feshbach resonance Pairing BEC of molecules BCS-BEC-crossover Conclusion 2 Motivation
More informationStrongly paired fermions
Strongly paired fermions Alexandros Gezerlis TALENT/INT Course on Nuclear forces and their impact on structure, reactions and astrophysics July 4, 2013 Strongly paired fermions Neutron matter & cold atoms
More informationICAP Summer School, Paris, Three lectures on quantum gases. Wolfgang Ketterle, MIT
ICAP Summer School, Paris, 2012 Three lectures on quantum gases Wolfgang Ketterle, MIT Cold fermions Reference for most of this talk: W. Ketterle and M. W. Zwierlein: Making, probing and understanding
More informationQuantum Mechanics: Fundamentals
Kurt Gottfried Tung-Mow Yan Quantum Mechanics: Fundamentals Second Edition With 75 Figures Springer Preface vii Fundamental Concepts 1 1.1 Complementarity and Uncertainty 1 (a) Complementarity 2 (b) The
More informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationThe Phase of a Bose-Einstein Condensate by the Interference of Matter Waves. W. H. Kuan and T. F. Jiang
CHINESE JOURNAL OF PHYSICS VOL. 43, NO. 5 OCTOBER 2005 The Phase of a Bose-Einstein Condensate by the Interference of Matter Waves W. H. Kuan and T. F. Jiang Institute of Physics, National Chiao Tung University,
More informationBroad and Narrow Fano-Feshbach Resonances: Condensate Fraction in the BCS-BEC Crossover
Broad and Narrow Fano-Feshbach Resonances: Condensate Fraction in the BCS-BEC Crossover Luca Salasnich Dipartimento di Fisica e Astronomia Galileo Galilei and CNISM, Università di Padova INO-CNR, Research
More informationSwinburne Research Bank
Swinburne Research Bank http://researchbank.swinburne.edu.au Deuar, P., & Drummond, P. D. (001). Stochastic gauges in quantum dynamics for many-body simulations. Originally published in Computer Physics
More informationEntropy and irreversibility in gas dynamics. Joint work with T. Bodineau, I. Gallagher and S. Simonella
Entropy and irreversibility in gas dynamics Joint work with T. Bodineau, I. Gallagher and S. Simonella Kinetic description for a gas of hard spheres Hard sphere dynamics The system evolves under the combined
More information5. Gross-Pitaevskii theory
5. Gross-Pitaevskii theory Outline N noninteracting bosons N interacting bosons, many-body Hamiltonien Mean-field approximation, order parameter Gross-Pitaevskii equation Collapse for attractive interaction
More informationA Mixture of Bose and Fermi Superfluids. C. Salomon
A Mixture of Bose and Fermi Superfluids C. Salomon INT workshop Frontiers in quantum simulation with cold atoms University of Washington, April 2, 2015 The ENS Fermi Gas Team F. Chevy, Y. Castin, F. Werner,
More informationCHAPTER 9 Statistical Physics
CHAPTER 9 Statistical Physics 9.1 9.2 9.3 9.4 9.5 9.6 9.7 Historical Overview Maxwell Velocity Distribution Equipartition Theorem Maxwell Speed Distribution Classical and Quantum Statistics Fermi-Dirac
More informationDynamics of interacting vortices on trapped Bose-Einstein condensates. Pedro J. Torres University of Granada
Dynamics of interacting vortices on trapped Bose-Einstein condensates Pedro J. Torres University of Granada Joint work with: P.G. Kevrekidis (University of Massachusetts, USA) Ricardo Carretero-González
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationCondensate fraction for a polarized three-dimensional Fermi gas
Condensate fraction for a polarized three-dimensional Fermi gas Luca Salasnich Dipartimento di Fisica e Astronomia Galileo Galilei, Università di Padova, Italy Camerino, June 26, 2014 Collaboration with:
More informationNon-Linear Plasma Wave Decay to Longer Wavelength
Non-Linear Plasma Wave Decay to Longer Wavelength F. Anderegg 1, a), M. Affolter 1, A. Ashourvan 1, D.H.E. Dubin 1, F. Valentini 2 and C.F. Driscoll 1 1 University of California San Diego Physics Department
More informationCold atoms in the presence of disorder and interactions
Cold atoms in the presence of disorder and interactions Collaboration: A. Minguzzi, S. Skipetrov, B. van-tiggelen (Grenoble), P. Henseler (Bonn), J. Chalker (Oxford), L. Beilin, E. Gurevich (Technion).
More informationCOPYRIGHTED MATERIAL. Index
347 Index a AC fields 81 119 electric 81, 109 116 laser 81, 136 magnetic 112 microwave 107 109 AC field traps see Traps AC Stark effect 82, 84, 90, 96, 97 101, 104 109 Adiabatic approximation 3, 10, 32
More informationFluids with dipolar coupling
Fluids with dipolar coupling Rosensweig instability M. D. Cowley and R. E. Rosensweig, J. Fluid Mech. 30, 671 (1967) CO.CO.MAT SFB/TRR21 STUTTGART, ULM, TÜBINGEN FerMix 2009 Meeting, Trento A Quantum Ferrofluid
More informationSymmetry of the Dielectric Tensor
Symmetry of the Dielectric Tensor Curtis R. Menyuk June 11, 2010 In this note, I derive the symmetry of the dielectric tensor in two ways. The derivations are taken from Landau and Lifshitz s Statistical
More informationPlan of the lectures
Plan of the lectures 1. Introductory remarks on metallic nanostructures Relevant quantities and typical physical parameters Applications. Linear electron response: Mie theory and generalizations 3. Nonlinear
More informationPHYSICS DEPARTMENT OF THE ÉCOLE NORMALE SUPÉRIEURE LABORATOIRE KASTLER BROSSEL
PHYSICS DEPARTMENT OF THE ÉCOLE NORMALE SUPÉRIEURE LABORATOIRE KASTLER BROSSEL DOCTORAL THESIS OF THE UNIVERSITÉ PARIS VI speciality : quantum physics presented by Florian SCHRECK to obtain the title Docteur
More informationControlled collisions of a single atom and an ion guided by movable trapping potentials
Controlled collisions of a single atom and an ion guided by movable trapping potentials Zbigniew Idziaszek CNR-INFM BEC Center, I-38050 Povo (TN), Italy and Center for Theoretical Physics, Polish Academy
More informationF. Chevy Seattle May 2011
THERMODYNAMICS OF ULTRACOLD GASES F. Chevy Seattle May 2011 ENS FERMION GROUPS Li S. Nascimbène Li/K N. Navon L. Tarruell K. Magalhaes FC C. Salomon S. Chaudhuri A. Ridinger T. Salez D. Wilkowski U. Eismann
More informationThe Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs
The Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs i ( ) t Φ (r, t) = 2 2 2m + V ext(r) + g Φ (r, t) 2 Φ (r, t) (Mewes et al., 1996) 26/11/2009 Stefano Carignano 1 Contents 1 Introduction
More informationRoton Mode in Dipolar Bose-Einstein Condensates
Roton Mode in Dipolar Bose-Einstein Condensates Sandeep Indian Institute of Science Department of Physics, Bangalore March 14, 2013 BECs vs Dipolar Bose-Einstein Condensates Although quantum gases are
More informationRevolution in Physics. What is the second quantum revolution? Think different from Particle-Wave Duality
PHYS 34 Modern Physics Ultracold Atoms and Trappe Ions Today and Mar.3 Contents: a) Revolution in physics nd Quantum revolution b) Quantum simulation, measurement, and information c) Atomic ensemble and
More informationLow dimensional quantum gases, rotation and vortices
Goal of these lectures Low dimensional quantum gases, rotation and vortices Discuss some aspect of the physics of quantum low dimensional systems Planar fluids Quantum wells and MOS structures High T c
More informationSTIMULATED RAMAN ATOM-MOLECULE CONVERSION IN A BOSE-EINSTEIN CONDENSATE. Chisinau, Republic of Moldova. (Received 15 December 2006) 1.
STIMULATED RAMAN ATOM-MOLECULE CONVERSION IN A BOSE-EINSTEIN CONDENSATE P.I. Khadzhi D.V. Tkachenko Institute of Applied Physics Academy of Sciences of Moldova 5 Academiei str. MD-8 Chisinau Republic of
More informationSignatures of Superfluidity in Dilute Fermi Gases near a Feshbach Resonance
Signatures of Superfluidity in Dilute ermi Gases near a eshbach Resonance A. Bulgac (Seattle), Y. Yu (Seattle Lund) P.. Bedaque (Berkeley), G.. Bertsch (Seattle), R.A. Broglia (Milan), A.C. onseca (Lisbon)
More informationStructure and stability of trapped atomic boson-fermion mixtures
Structure and stability of trapped atomic boson-fermion mixtures Robert Roth Clarendon Laboratory, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom (Dated: May 8, 22) The structure of trapped
More informationExperiments with an Ultracold Three-Component Fermi Gas
Experiments with an Ultracold Three-Component Fermi Gas The Pennsylvania State University Ken O Hara Jason Williams Eric Hazlett Ronald Stites John Huckans Overview New Physics with Three Component Fermi
More informationHigh-Temperature Superfluidity
High-Temperature Superfluidity Tomoki Ozawa December 10, 2007 Abstract With the recent advancement of the technique of cooling atomic gases, it is now possible to make fermionic atom gases into superfluid
More informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
More informationSupported by NIST, the Packard Foundation, the NSF, ARO. Penn State
Measuring the electron edm using Cs and Rb atoms in optical lattices (and other experiments) Fang Fang Osama Kassis Xiao Li Dr. Karl Nelson Trevor Wenger Josh Albert Dr. Toshiya Kinoshita DSW Penn State
More informationDissipative nuclear dynamics
Dissipative nuclear dynamics Curso de Reacciones Nucleares Programa Inter universitario de Fisica Nuclear Universidad de Santiago de Compostela March 2009 Karl Heinz Schmidt Collective dynamical properties
More information1 Fluctuations of the number of particles in a Bose-Einstein condensate
Exam of Quantum Fluids M1 ICFP 217-218 Alice Sinatra and Alexander Evrard The exam consists of two independant exercises. The duration is 3 hours. 1 Fluctuations of the number of particles in a Bose-Einstein
More information( ) /, so that we can ignore all
Physics 531: Atomic Physics Problem Set #5 Due Wednesday, November 2, 2011 Problem 1: The ac-stark effect Suppose an atom is perturbed by a monochromatic electric field oscillating at frequency ω L E(t)
More informationCold Metastable Neon Atoms Towards Degenerated Ne*- Ensembles
Cold Metastable Neon Atoms Towards Degenerated Ne*- Ensembles Supported by the DFG Schwerpunktprogramm SPP 1116 and the European Research Training Network Cold Quantum Gases Peter Spoden, Martin Zinner,
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More informationWe already came across a form of indistinguishably in the canonical partition function: V N Q =
Bosons en fermions Indistinguishability We already came across a form of indistinguishably in the canonical partition function: for distinguishable particles Q = Λ 3N βe p r, r 2,..., r N ))dτ dτ 2...
More informationNonlinear BEC Dynamics by Harmonic Modulation of s-wave Scattering Length
Nonlinear BEC Dynamics by Harmonic Modulation of s-wave Scattering Length I. Vidanović, A. Balaž, H. Al-Jibbouri 2, A. Pelster 3 Scientific Computing Laboratory, Institute of Physics Belgrade, Serbia 2
More informationSuperfluidity and Condensation
Christian Veit 4th of June, 2013 2 / 29 The discovery of superfluidity Early 1930 s: Peculiar things happen in 4 He below the λ-temperature T λ = 2.17 K 1938: Kapitza, Allen & Misener measure resistance
More information