Rapporto di ricerca Research report
|
|
- Norma Small
- 5 years ago
- Views:
Transcription
1 Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report September /2009 Spectral methods for dissipative nonlinear Schrödinger equations Laura M. Morato Marco Caliari Simone Zuccher Questo rapporto è disponibile su Web all indirizzo: This report is available on the web at the address:
2 Abstract In this technical report we describe an application of spectral methods to numerically solve some nonlinear Schrödinger equations with dissipative terms and carefully study the problem of vortex formation. Keywords: Bose Einstein condensation, dissipation, vortex formation,spectral methods
3 1 Introduction A Bose Einstein condensate (BEC) is a state of matter of a dilute gas of interacting bosons confined in an external potential and cooled to temperatures very near to absolute zero. Under such conditions, a large fraction of the bosons occupies the same quantum state and quantum effects become apparent on a macroscopic scale. A spectacular quantum phenomenon is the formation of quantized vortices when the condensate is stirred into rotation. The mechanism is not yet completely understood and represents an active research field both in physics and in mathematics. A BEC trapped in an external potential is usually described by a macroscopic wave function Ψ obeying the Gross Pitaevskii (GP) equation. In the 3D frame rotating with the frequency Ω around the Z axis the GP equation reads i Ψ ) ( t = 2 2m 2 + V trap + V rot + g Ψ 2 ΩL Z Ψ. (1) Here g = 4π 2 a/m (coupling constant) represents the strength of interactions characterized by the s-wave scattering length a > 0 and L Z = i (X Y Y X ) the angular momentum. The wave function is normalized by the total particle number. We assume that the external harmonic trapping potential has the form V trap (X,Y,Z) = 1 2 m [ ω 2 (X 2 + Y 2 ) + ω 2 ZZ 2], and that the potential that drives the rotation is V rot (X,Y,Z) = 1 2 mω2 (ε X X 2 + ε Y Y 2 ). with the anisotropy parameters ε X ε Y. Under a proper reduction to 2D and normalization, equation (1) can be rewritten (in the laboratory frame) i ψ ( t = 12 ) 2 + V + θ ψ 2 ψ. (2) Here V is the trapping and rotating potential defined by V (x,y,t) = (1 + ε x) 2 (x cos Ωt + y sin Ωt) 2 + (1 + ε y) ( x sin Ωt + y cos Ωt)
4 2 GPE with extrinsic dissipation In order to study the vortex formation by numerical experiments, in [6] the following dissipative generalization of (1) was proposed: (i γ) ψ [ t = 1 ] V + θ ψ 2 µ + γ(x y y x ) ψ. (3) This equation somehow models friction between the condensate and the normal cloud: quoting [6, 11], the term with γ introduces the dissipation. Although the detailed mechanism of the dissipation is yet to be understood the observation of the vortex lattices implies the presence of dissipation. When the condensate is grown from a rotating vapor cloud the growth process itself is dissipative, and in the experiments [...] dissipation can arise by transfer of atoms between the thermal cloud and the condensate. [...] A dissipative term arises from collisions between atoms in the thermal cloud trapped by the same potential as the condensate, in which one of the colliding atoms enters the condensate after the collision. One can see that a pure advection term appears. A difficulty with this equation is that, in order to preserve the total number of particles, it is necessary to introduce the ψ-dependent chemical potential µ. A further problem is to understand whether this picture is valid at zero temperature or not. Equation (3), rewritten in the rotating frame, has been discretized in space by a spectral method, namely Galerkin s method applied to the decomposition into Hermite functions (cf. [1, 3, 5]). Hermite functions naturally vanish at infinity and, thus, no artificial boundary condition is imposed as usually done in the literature. For the time discretization, an explicit Runge Kutta method of order four has been used [5]. Since the chemical potential µ depends on ψ, given the solution ψ at a certain time, µ was computed and kept constant for the next time step. The ground state solution (cf. [4]) of the pure Gross Pitaevskii Equation (i.e., γ = µ = 0) was chosen as initial condition. We note that, if the initial condition satisfies the symmetry condition ψ 0 (x,y) = ψ 0 ( x, y), then the same is true for the solution ψ(x,y,t), as it can be shown by substituting x with x and y with y in equation (3). According to [6], the values for the parameters were ε x = 0.03, ε y = 0.09 and Ω = 0.7. The dissipation constant γ was ranging between 0 and 0.03 and the coupling constant θ between 0 and 500/ 2. From our simulations, we can draw the following conclusions: If γ = 0 and the space discretization is fine enough, i.e., a sufficiently large number N (N 32) of Hermite functions is used, no vortex appears during the simulation. On the contrary, if a poor space resolution 2
5 Figure 1: ψ 2 at time t = 244 with N = 64 Hermite functions (left) and with N = 16 (right). 10 time t = time t = Figure 2: ψ 2 at time t = 360 with the whole set of Hermite basis functions (left) and with only even Hermite basis functions (right). is used (e.g., N = 16), then some artificial (and unstable) vortices appear during the simulation (see Figure 1). This confirms experiments shown in [12]. If θ = 0 (no cubic term in the equation) no vortex appears during the simulation. For γ = 0.03 and θ = 500/ 2 the simulations show results in agreement with those reported in [6], with a final stable lattice of vortices, in which the symmetry ψ(x,y,t) = ψ( x, y,t) is broken, due to numerical effects, see Figure 2 (left). If only the even Hermite functions are used, than the result is a stable and symmetric lattice of vortices, see Figure 2 (right), with an energy slightly higher than the energy corresponding to the unsymmetric lattice. 3
6 3 Schrödinger equation with vorticity In this section we consider a new dissipative generalization of the GP equation, recently proposed in [2]: i ψ [ ] 1 t = 2 (i + A)2 + V + θ ψ 2 ψ A (b t = 12 ) ( A) where A is a real vector field, b is defined as b = S A 1 2 ln ρ, (4) being S and ρ the phase and the absolute value of ψ, i.e., ψ = ρ 1 2e is. If θ = 0, this is the generalization of the Schrödinger equation which arises in Stochastic Quantization with the Stochastic Lagrangian Variational principle as introduced in [7, 9]. Indeed, for θ = 0 and A = 0, one gets the usual Schrödinger equation. Moreover, if the Hamiltonian operator H = V is bounded from below in L 2 (dx), then for generic initial data the solution asymptotically relaxes towards a solution where A is equal to zero almost surely with respect to the Lebesgue measure. Thus the set of the solutions to the canonical Schrödinger equation acts as an attracting set. Indeed, denoting by E(t) the energy functional E(t) = R 3 [ ] 1 2 i ψ + Aψ 2 + V ψ 2 dx the Energy theorem proved in [7] claims that, under mild regularity assumptions, the following equality holds d dt E(t) = 1 A 2 ψ 2 dx 0. (5) 2 R 3 If θ 0, equation (4) arises naturally in the stochastic quantization of a system of trapped pair interacting bosons [8]. Notably, the energy theorem still holds in this case, with E(t) replaced by E θ (t) = R 3 [ 1 2 i ψ + Aψ 2 + θ ] 2 ψ 2 + V ψ 2 dx while the right side of (5) remains unchanged (see [2]). 4
7 Hereafter we assume θ = 0. Since the spectral method based on Hermite functions appeared to be a good candidate for the simulation of vortex dynamics without causing spurious effects at the boundary, we chose to apply it to equation (4). The energy theorem proved in [7, 9] assures that E(t) = R 3 [ ] 1 2 i ψ + Aψ 2 + V ψ 2 dx < +, therefore Aψ = Γ L 2 (R 3 C) and Γ can be naturally expanded into Hermite functions. Equation (4) in the unknowns ψ and Γ reads ψ ψ ψ t = i 2 ψ( ψ 2 ψ + Γ 2 + i ψ Γ + i Γ ψ) iψ ψv ψ Γ ψ ψ t = i 2 ψ ( ψ Γ) [ ψ( ψ Γ ψ Γ)]+ [( ) ] 1 ψ 2 ψ Γ ( Γ) + i 2 (ψ ψ ψ ψ) ( Γ)+ ] [ 1 (i 1) ψ + Γ 2 ( ψ Γ) ψ + Γ ψ t, (6) where the first equation was multiplied by iψ with the aim of having the same mass matrix for both the equations when Galerkin s method is applied. As in the previous case, we considered the 2D version of the model. Integration by parts was applied to each second order term in the equations. Once the problem was discretized in space, the explicit Runge Kutta Fehlberg 4(5) (RKF45) method with variable step size was employed for time integration. It should be noticed that, although the method is explicit, a linear system has to be solved at each stage because of the mass matrix. If a linear-gaussian initial condition with radial symmetry, namely ψ 0 (r) = σ0 π e σ 0r2 /2 A 0 (r) = a 0 rˆr + α 0 rˆθ (7) is taken for equation (4), then the global solution is still linear-gaussian with radial symmetry (cf. [10]) σ(t) /2 ψ(r,t) = π e σ(t)r2 A(r,t) = a(t)rˆr + α(t)rˆθ 5
8 where σ(t), a(t) and α(t) satisfy σ(t) = 2a(t)σ(t) σ(0) = σ 0 ȧ(t) = σ 2 (t) + α 2 (t) a 2 (t) 1 a(0) = a 0 (8) α(t) = 2(a(t) + σ(t))α(t) α(0) = α 0 and the corresponding energy is E(t) = 1 2σ(t) (σ2 (t) + a 2 (t) + α 2 (t) + 1). In order to verify the correct implementation, we chose a linear-gaussian initial solution (7), with Γ 0 (r) = A 0 (r)ψ 0 (r), for the PDEs system (6) and compared the solution with that of the much simpler ODEs system (8). In our first experiments, we successfully reproduced the breathing phenomenon [5], i.e., we chose σ 0 = 2, a 0 = 0 and α 0 = 0, so that the first equation in system (6) is the linear Schrödinger equation multiplied by iψ ψ and the second one has the trivial solution Γ(t) Γ 0 = A 0 ψ = 0 for t 0. We observed the conservation of the energy and a perfect match between the PDEs and ODEs results. Then we tested a more general linear-gaussian solution by taking the initial solution (7) with σ 0 = 2, a 0 = 1 and α 0 = 2. The 2.5 PDEs ODEs energy time Figure 3: PDEs and ODEs computed energies for σ 0 = 2, a 0 = 0 and α 0 = 0 6
9 energy comparison is reported in Figure 3. We observe a very good agreement of the computed energies up to about t = 0.5. Then, the solution of the PDEs system departures from the correct solution, since the corresponding energy increases instead of decreasing (see (5)). From the numerical point of view, we can see that the RKF45 scheme, in the effort of limiting the error, keeps reducing the time step size, up to useless values. In order to exclude stiffness issues, we also implemented the backward (implicit) Euler method, based on a Newton-like method with a modified Jacobian for the resulting non-linear system, obtaining the same behaviour as for the explicit method. A further investigation of the problem led to the conclusion that the mass matrix is very ill-conditioned (condition numbers of order ). However, when a 0 = α 0 = 0 and thus Γ(t) 0, despite of condition numbers of order 10 18, the correct solution is very well captured and maintained in time. This suggests that the errors due to the ill-conditioned mass matrix have a fatal effect only when the second equation plays a non-trivial role in system (6). Future work may be twofold: either ad hoc numerical methods for improving the numerical solution of ill-conditioned systems should be devised, or completely new approaches, such as Monte Carlo methods, should be employed. References [1] W. Bao and J. Shen, A fourth-order time-splitting Laguerre Hermite pseudospectral method for Bose Einstein condensates, SIAM J. Sci. Comput. 26(6), (2005). [2] M. Caliari, M. I. Loffredo, L. M. Morato, and S. Zuccher, Cubic Nonlinear Schrödinger Equation with vorticity, New J. Phys. 10, (2008). [3] M. Caliari, Ch. Neuhauser, and M. Thalhammer, High-order timesplitting Hermite and Fourier spectral methods for the Gross Pitaevskii equation, J. Comput. Phys. 228(3), (2009). [4] M. Caliari, A. Ostermann, S. Rainer, and M. Thalhammer, A minimisation approach for computing the ground state of Gross Pitaevskii systems, J. Comput. Phys. 228(2), (2009). [5] C. M. Dion and E. Cancès, Spectral method for the time-dependent Gross Pitaevskii equation with a harmonic trap, Phys. Rev. E 67, (2003). 7
10 [6] K. Kasamatsu, M. Tsubota, and M. Ueda, Nonlinear dynamics of vortex lattice formation in a rotating Bose Einstein condensate, Phys. Rev. A 67, (2003). [7] M. Loffredo and L. M. Morato, Lagrangian variational principle in stochastic mechanics: Gauge structure and stability, J. Math. Phys. 30, (1989). [8] M. Loffredo and L. M. Morato, Stochastic quantization for a system of N identical interacting Bose particles, J. Phys. A: Math. Theor., 40(30), (2007). [9] L. M. Morato, Path-wise stochastic calculus of variations and quantum systems, Phys. Rev. D 31(8), (1985). [10] L. M. Morato and S. Ugolini, Gaussian solutions to the Lagrangian variational problem in Stochastic Mechanics, Ann. Inst. H. Poincaré Phys. Théor., 60(3), (1994). [11] A. A. Penkwitt and R. J. Ballagh, Nucleation, growth, and stabilization of Bose Einstein condensate vortex lattices, Phys. Rev. Lett. 89, (2002) [12] M. Tsubota, Quantum turbulence: From superfluid helium to atomic Bose Einstein condensates, invited lecture at Mathematical Models of Quantum Fluids, September 14 17, Verona (2009). 8
11
12 University of Verona Department of Computer Science Strada Le Grazie, 15 I Verona Italy
High-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More informationLecture 2: Weak Interactions and BEC
Lecture 2: Weak Interactions and BEC Previous lecture: Ideal gas model gives a fair intuition for occurrence of BEC but is unphysical (infinite compressibility, shape of condensate...) Order parameter
More informationQuantized Vortex Stability and Dynamics in Superfluidity and Superconductivity
Quantized Vortex Stability and Dynamics in Superfluidity and Superconductivity Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
More informationModeling, Analysis and Simulation for Degenerate Dipolar Quantum Gas
Modeling, Analysis and Simulation for Degenerate Dipolar Quantum Gas Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
More informationThird Critical Speed for Rotating Bose-Einstein Condensates. Mathematical Foundations of Physics Daniele Dimonte, SISSA
Third Critical Speed for Rotating Bose-Einstein Condensates Mathematical Foundations of Physics Daniele Dimonte, SISSA 1 st November 2016 this presentation available on daniele.dimonte.it based on a joint
More informationNumerical methods for computing vortex states in rotating Bose-Einstein condensates. Ionut Danaila
Numerical methods for computing vortex states in rotating Bose-Einstein condensates Ionut Danaila Laboratoire de mathématiques Raphaël Salem Université de Rouen www.univ-rouen.fr/lmrs/persopage/danaila
More informationRapidly Rotating Bose-Einstein Condensates in Strongly Anharmonic Traps. Michele Correggi. T. Rindler-Daller, J. Yngvason math-ph/
Rapidly Rotating Bose-Einstein Condensates in Strongly Anharmonic Traps Michele Correggi Erwin Schrödinger Institute, Vienna T. Rindler-Daller, J. Yngvason math-ph/0606058 in collaboration with preprint
More informationSolitons and vortices in Bose-Einstein condensates with finite-range interaction
Solitons and vortices in Bose-Einstein condensates with finite-range interaction Luca Salasnich Dipartimento di Fisica e Astronomia Galileo Galilei and CNISM, Università di Padova INO-CNR, Research Unit
More informationPAPER 84 QUANTUM FLUIDS
MATHEMATICAL TRIPOS Part III Wednesday 6 June 2007 9.00 to 11.00 PAPER 84 QUANTUM FLUIDS Attempt TWO questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY REQUIREMENTS
More informationVortices in Bose-Einstein condensates. Ionut Danaila
Vortices in Bose-Einstein condensates 3D numerical simulations Ionut Danaila Laboratoire Jacques Louis Lions Université Pierre et Marie Curie (Paris 6) http://www.ann.jussieu.fr/ danaila October 16, 2008
More informationNumerical simulation of Bose-Einstein Condensates based on Gross-Pitaevskii Equations
1/38 Numerical simulation of Bose-Einstein Condensates based on Gross-Pitaevskii Equations Xavier ANTOINE Institut Elie Cartan de Lorraine (IECL) & Inria-SPHINX Team Université de Lorraine, France Funded
More informationBose-Einstein condensates under rotation: The structures within
Bose-Einstein condensates under rotation: The structures within Peter Mason Centre de Mathématiques Appliquées, Ecole Polytechnique Paris, France In collaboration with Amandine Aftalion and Thierry Jolicoeur:
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 informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationThe Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs
The Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs RHI seminar Pascal Büscher i ( t Φ (r, t) = 2 2 ) 2m + V ext(r) + g Φ (r, t) 2 Φ (r, t) 27 Nov 2008 RHI seminar Pascal Büscher 1 (Stamper-Kurn
More informationThermalisation and vortex formation in a mechanically perturbed condensate. Vortex Lattice Formation. This Talk. R.J. Ballagh and Tod Wright
Thermalisation and vortex formation in a mechanically perturbed condensate Jack Dodd Centre ACQAO R.J. Ballagh and Tod Wright Jack Dodd centre for Photonics and Ultra-cold Atoms Department of Physics University
More informationDynamics and Statistics of Quantum Turbulence in Quantum Fluid
Dynamics and Statistics of Quantum Turbulence in Quantum Fluid Faculty of Science, Osaka City University Michikazu Kobayashi May 25, 2006, Kansai Seminar House Contents 1. 2. 3. 4. 5. Introduction - history
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 informationQuantum vortex reconnections
Quantum vortex reconnections A.W. Baggaley 1,2, S. Zuccher 4, Carlo F Barenghi 2, 3, A.J. Youd 2 1 University of Glasgow 2 Joint Quantum Centre Durham-Newcastle 3 Newcastle University 4 University of Verona
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationNumerical simulation of the dynamics of the rotating dipolar Bose-Einstein condensates
Numerical simulation of the dynamics of the rotating dipolar Bose-Einstein condensates Qinglin TANG Inria, University of Lorraine Joint work with: Prof. Weizhu BAO, Daniel MARAHRENS, Yanzhi ZHANG and Yong
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 informationVortices and other topological defects in ultracold atomic gases
Vortices and other topological defects in ultracold atomic gases Michikazu Kobayashi (Kyoto Univ.) 1. Introduction of topological defects in ultracold atoms 2. Kosterlitz-Thouless transition in spinor
More informationSuperfluidity. v s. E. V. Thuneberg Department of Physical Sciences, P.O.Box 3000, FIN University of Oulu, Finland (Dated: June 8, 2012)
Superfluidity E. V. Thuneberg Department of Physical Sciences, P.O.Box 3000, FIN-90014 University of Oulu, Finland (Dated: June 8, 01) PACS numbers: 67.40.-w, 67.57.-z, 74., 03.75.-b I. INTRODUCTION Fluids
More information1 Superfluidity and Bose Einstein Condensate
Physics 223b Lecture 4 Caltech, 04/11/18 1 Superfluidity and Bose Einstein Condensate 1.6 Superfluid phase: topological defect Besides such smooth gapless excitations, superfluid can also support a very
More informationPeriodic oscillations in the Gross-Pitaevskii equation with a parabolic potential
Periodic oscillations in the Gross-Pitaevskii equation with a parabolic potential Dmitry Pelinovsky 1 and Panos Kevrekidis 2 1 Department of Mathematics, McMaster University, Hamilton, Ontario, Canada
More informationPhysics 598 ESM Term Paper Giant vortices in rapidly rotating Bose-Einstein condensates
Physics 598 ESM Term Paper Giant vortices in rapidly rotating Bose-Einstein condensates Kuei Sun May 4, 2006 kueisun2@uiuc.edu Department of Physics, University of Illinois at Urbana- Champaign, 1110 W.
More informationVortices in Rotating Bose-Einstein Condensates A Review of (Recent) Mathematical Results. Michele Correggi
Vortices in Rotating Bose-Einstein Condensates A Review of (Recent) Mathematical Results Michele Correggi CIRM FBK, Trento 15/09/2009 Mathematical Models of Quantum Fluids Università di Verona M. Correggi
More informationLOCALIZATION OF RELATIVE ENTROPY IN BOSE-EINSTEIN CONDENSATION OF TRAPPED INTERACTING BOSONS
LOCALIZATION OF RELATIVE ENTROPY IN BOSE-EINSTEIN CONDENSATION OF TRAPPED INTERACTING BOSONS LAURA M. MORATO 1 AND STEFANIA UGOLINI 2 1 Facoltà di Scienze, Università di Verona, Strada le Grazie, 37134
More informationbeen succeeded in 1997 Rb, 23 Na, 7 Li, 1 H, 85 Rb, 41 K, 4 He, 133 Cs, 174 Yb, 52 Cr, 40 Ca, 84 Sr, 164 Dy Laser cooling Trap of atoms 87
Non-Abelian Vortices and Their Non-equilibrium Michikazu Kobayashi a University of Tokyo November 18th, 2011 at Keio University 2 nd Workshop on Quarks and Hadrons under Extreme Conditions - Lattice QCD,
More informationExcitations and dynamics of a two-component Bose-Einstein condensate in 1D
Author: Navarro Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain. Advisor: Bruno Juliá Díaz Abstract: We study different solutions and their stability for a two component
More informationQinglin TANG. INRIA & IECL, University of Lorraine, France
Numerical methods on simulating dynamics of the nonlinear Schrödinger equation with rotation and/or nonlocal interactions Qinglin TANG INRIA & IECL, University of Lorraine, France Collaborators: Xavier
More informationSuperfluidity in bosonic systems
Superfluidity in bosonic systems Rico Pires PI Uni Heidelberg Outline Strongly coupled quantum fluids 2.1 Dilute Bose gases 2.2 Liquid Helium Wieman/Cornell A. Leitner, from wikimedia When are quantum
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 informationBose-Einstein condensates in optical lattices: mathematical analysis and analytical approximate formulas
0.5 setgray0 0.5 setgray1 Bose-Einstein condensates in optical lattices: mathematical analysis and analytical approximate formulas IV EBED João Pessoa - 2011 Rolci Cipolatti Instituto de Matemática - UFRJ
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationThe XY model, the Bose Einstein Condensation and Superfluidity in 2d (I)
The XY model, the Bose Einstein Condensation and Superfluidity in 2d (I) B.V. COSTA UFMG BRAZIL LABORATORY FOR SIMULATION IN PHYSICS A Guide to Monte Carlo Simulations in Statistical Physics by Landau
More informationWe can then linearize the Heisenberg equation for in the small quantity obtaining a set of linear coupled equations for and :
Wednesday, April 23, 2014 9:37 PM Excitations in a Bose condensate So far: basic understanding of the ground state wavefunction for a Bose-Einstein condensate; We need to know: elementary excitations in
More informationEffective Dynamics of Solitons I M Sigal
Effective Dynamics of Solitons I M Sigal Porquerolles October, 008 Joint work with Jürg Fröhlich, Stephen Gustafson, Lars Jonsson, Gang Zhou and Walid Abou Salem Bose-Einstein Condensation Consider a system
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 informationA simple Dufort-Frankel type scheme for the Gross-Pitaevskii equation of Bose-Einstein condensates on different geometries
A simple Dufort-Frankel type scheme for the Gross-Pitaevskii equation of Bose-Einstein condensates on different geometries Ming-Chih Lai Chung-Yin Huang Te-Sheng Lin Department of Applied Mathematics,
More informationSpontaneous Symmetry Breaking in Bose-Einstein Condensates
The 10th US-Japan Joint Seminar Spontaneous Symmetry Breaking in Bose-Einstein Condensates Masahito UEDA Tokyo Institute of Technology, ERATO, JST collaborators Yuki Kawaguchi (Tokyo Institute of Technology)
More informationLecture I: numerical computation of the dynamics of Schrödinger equations. Extension to the GPE for BECs
/67 Lecture I: numerical computation of the dynamics of Schrödinger equations. Extension to the GPE for BECs Xavier Antoine Institut Elie Cartan de Lorraine (IECL) Université de Lorraine INRIA-CORIDA Team
More informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
More information(1) BEC (2) BEC (1) (BEC) BEC BEC (5) BEC (LT) (QFS) BEC (3) BEC. 3 He. 4 He. 4 He 3 He
22 (BEC) 4 He BEC We study theoretically and experimentally quantum hydrodynamics in quantum condensed phases at low temperature, namely superfluid helium and atomic Bose-Einstein condensates (BECs). Quantum
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationQuantized vortex stability and interaction in the nonlinear wave equation
Physica D 237 (2008) 2391 2410 www.elsevier.com/locate/physd Quantized vortex stability and interaction in the nonlinear wave equation Weizhu Bao a,b,, Rong Zeng c, Yanzhi Zhang a,1 a Department of Mathematics,
More informationComputer Physics Communications. Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations
Computer Physics Communications (13) 1 3 Computer Physics Communications Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations Xavier Antoine a,b, Weizhu Bao 1c,d,
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 informationPattern Formation in the Fractional Quantum Hall Effect
Journal of the Physical Society of Japan 72, Supplement C (2003) 18-23 Pattern Formation in the Fractional Quantum Hall Effect Pierre Gaspard Center for Nonlinear Phenomena and Complex Systems, Université
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 informationCold atoms. 1: Bose-Einstein Condensation. Emil Lundh. April 13, Department of Physics Umeå University
1: Bose-Einstein Condensation Department of Physics Umeå University lundh@tp.umu.se April 13, 2011 Umeå 114 000 inhabitants Average age 37.9 years Cultural capital of Europe 2014 400 km ski tracks 180
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 informationA uniformly convergent numerical method for singularly perturbed nonlinear eigenvalue problems
A uniformly convergent numerical method for singularly perturbed nonlinear eigenvalue problems Weizhu Bao and Ming-Huang Chai Abstract In this paper we propose a uniformly convergent numerical method for
More informationCitation PHYSICAL REVIEW A (2004), 70(4) RightCopyright 2004 American Physical So
Title Splitting instability of a multiply Einstein condensate Author(s) Kawaguchi, Y; Ohmi, T Citation PHYSICAL REVIEW A (2004), 70(4) Issue Date 2004-10 URL http://hdl.handle.net/2433/49950 RightCopyright
More informationInteraction between atoms
Interaction between atoms MICHA SCHILLING HAUPTSEMINAR: PHYSIK DER KALTEN GASE INSTITUT FÜR THEORETISCHE PHYSIK III UNIVERSITÄT STUTTGART 23.04.2013 Outline 2 Scattering theory slow particles / s-wave
More information16. GAUGE THEORY AND THE CREATION OF PHOTONS
6. GAUGE THEORY AD THE CREATIO OF PHOTOS In the previous chapter the existence of a gauge theory allowed the electromagnetic field to be described in an invariant manner. Although the existence of this
More informationVortex dynamics in finite temperature two-dimensional superfluid turbulence. Andrew Lucas
Vortex dynamics in finite temperature two-dimensional superfluid turbulence Andrew Lucas Harvard Physics King s College London, Condensed Matter Theory Special Seminar August 15, 2014 Collaborators 2 Paul
More informationQuantum Lattice Gas Algorithm for Quantum Turbulence and Vortex Reconnection in the Gross-Pitaevskii Equation
Quantum Lattice Gas Algorithm for Quantum Turbulence and Vortex Reconnection in the Gross-Pitaevskii Equation George Vahala 1, Jeffrey Yepez 1,2 and Linda Vahala 3 1 Dept. of Physics, William & Mary, Williamsburg,
More informationSuperfluidity and Supersolidity in 4 He
Superfluidity and Supersolidity in 4 He Author: Lars Bonnes Supervisor: Lode Pollet Proseminar Theoretische Physik: Phase Transitions SS 07 18.06.2007 Superfluid Liquid Helium Motivation Two-fluid Model
More informationStationary States of Bose Einstein Condensates in Single- and Multi-Well Trapping Potentials
Laser Physics, Vol., No.,, pp. 37 4. Original Tet Copyright by Astro, Ltd. Copyright by MAIK Nauka /Interperiodica (Russia). ORIGINAL PAPERS Stationary States of Bose Einstein Condensates in Single- and
More informationTurbulence in quantum hydrodynamics
Turbulence in quantum hydrodynamics Michikazu Kobayashi Department Physics, Kyoto University By numerically solving the Gross-Piteavskii equation with the energy injection and dissipation, we study quantum
More informationAdiabatic trap deformation for preparing Quantum Hall states
Marco Roncaglia, Matteo Rizzi, and Jean Dalibard Adiabatic trap deformation for preparing Quantum Hall states Max-Planck Institut für Quantenoptik, München, Germany Dipartimento di Fisica del Politecnico,
More informationSYNTHETIC GAUGE FIELDS IN ULTRACOLD ATOMIC GASES
Congresso Nazionale della Società Italiana di Fisica Università della Calabria 17/21 Settembre 2018 SYNTHETIC GAUGE FIELDS IN ULTRACOLD ATOMIC GASES Sandro Stringari Università di Trento CNR-INO - Bose-Einstein
More informationSuperfluidity & Bogoliubov Theory: Rigorous Results
Superfluidity & Bogoliubov Theory: Rigorous Results Mathieu LEWIN mathieu.lewin@math.cnrs.fr (CNRS & Université Paris-Dauphine) collaborations with P.T. Nam (Vienna), N. Rougerie (Grenoble), B. Schlein
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 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 informationArtificial Gauge Fields for Neutral Atoms
Artificial Gauge Fields for Neutral Atoms Simon Ristok University of Stuttgart 07/16/2013, Hauptseminar Physik der kalten Gase 1 / 29 Outline 1 2 3 4 5 2 / 29 Outline 1 2 3 4 5 3 / 29 What are artificial
More informationarxiv: v1 [cond-mat.quant-gas] 8 Sep 2010
arxiv:1009.1574v1 [cond-mat.quant-gas] 8 Sep 2010 Theoretical confirmation of Feynman s hypothesis on the creation of circular vortices in Bose Einstein condensates: III E Infeld, A Senatorski and A A
More informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
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 informationLecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations
Lecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations Usama Al Khawaja, Department of Physics, UAE University, 24 Jan. 2012 First International Winter School on Quantum Gases
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationHylomorphic solitons and their dynamics
Hylomorphic solitons and their dynamics Vieri Benci Dipartimento di Matematica Applicata U. Dini Università di Pisa 18th May 2009 Vieri Benci (DMA-Pisa) Hylomorphic solitons 18th May 2009 1 / 50 Types
More informationDesign and realization of exotic quantum phases in atomic gases
Design and realization of exotic quantum phases in atomic gases H.P. Büchler and P. Zoller Theoretische Physik, Universität Innsbruck, Austria Institut für Quantenoptik und Quanteninformation der Österreichischen
More informationQuantised Vortices in an Exciton- Polariton Condensate
4 th International Conference on Spontaneous Coherence in Excitonic Systems Quantised Vortices in an Exciton- Polariton Condensate Konstantinos G. Lagoudakis 1, Michiel Wouters 2, Maxime Richard 1, Augustin
More informationMax Lewandowski. Axel Pelster. March 12, 2012
Primordial Models for Dissipative Bose-Einstein Condensates Max Lewandowski Institut für Physik und Astronomie, Universität Potsdam Axel Pelster Hanse-Wissenschaftskolleg, Delmenhorst March 12, 2012 Experiment
More informationSUPERFLUIDTY IN ULTRACOLD ATOMIC GASES
College de France, May 14, 2013 SUPERFLUIDTY IN ULTRACOLD ATOMIC GASES Sandro Stringari Università di Trento CNR-INFM PLAN OF THE LECTURES Lecture 1. Superfluidity in ultra cold atomic gases: examples
More information2 Canonical quantization
Phys540.nb 7 Canonical quantization.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system?.1.1.lagrangian Lagrangian mechanics is a reformulation of classical mechanics.
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 informationMichikazu Kobayashi. Kyoto University
Topological Excitations and Dynamical Behavior in Bose-Einstein Condensates and Other Systems Michikazu Kobayashi Kyoto University Oct. 24th, 2013 in Okinawa International Workshop for Young Researchers
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 informationSuperfluidity and superconductivity. IHP, Paris, May 7 and 9, 2007
Superfluidity and superconductivity. IHP, Paris, May 7 and 9, 2007 L.P. Pitaevskii Dipartimento di Fisica, Universita di Trento, INFM BEC CNR,Trento, Italy; Kapitza Institute for Physical Problems, ul.
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 informationSuperfluidity. Krzysztof Myśliwy. October 30, Theoretical Physics Proseminar
Superfluidity Krzysztof Myśliwy Theoretical Physics Proseminar October 30, 2017 Outline The λ transition Phenomenology of He-II Landau theory- a semi-phenomenological approach Feynman s explanation- from
More informationOn the Torus Quantization of Two Anyons with Coulomb Interaction in a Magnetic Field
Preprint DFPD/97/TH/15 On the Torus Quantization of Two Anyons with Coulomb Interaction in a Magnetic Field Luca Salasnich 1 Dipartimento di Matematica Pura ed Applicata Università di Padova, Via Belzoni
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 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 informationDilute Bose-Einstein condensates and their role in the study of quantum fluids: rotations, vortices, and coherence.
Dilute Bose-Einstein condensates and their role in the study of quantum fluids: rotations, vortices, and coherence. Kirk Madison Vincent Bretin Frédéric Chevy Peter Rosenbuch Wendel Wohlleben Jean Dalibard
More informationSuperfluidity of a 2D Bose gas (arxiv: v1)
Superfluidity of a 2D Bose gas (arxiv:1205.4536v1) Christof Weitenberg, Rémi Desbuquois, Lauriane Chomaz, Tarik Yefsah, Julian Leonard, Jérôme Beugnon, Jean Dalibard Trieste 18.07.2012 Phase transitions
More informationarxiv: v1 [cond-mat.quant-gas] 15 Jun 2010
The approach to vortex reconnection Richard Tebbs 1, Anthony J. Youd 1, and Carlo F. Barenghi 1 1 School of Mathematics and Statistics, Newcastle University, arxiv:1006.3004v1 [cond-mat.quant-gas] 15 Jun
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 informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationEffects of spin-orbit coupling on the BKT transition and the vortexantivortex structure in 2D Fermi Gases
Effects of spin-orbit coupling on the BKT transition and the vortexantivortex structure in D Fermi Gases Carlos A. R. Sa de Melo Georgia Institute of Technology QMath13 Mathematical Results in Quantum
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationINTERACTING BOSE GAS AND QUANTUM DEPLETION
922 INTERACTING BOSE GAS AND QUANTUM DEPLETION Chelagat, I., *Tanui, P.K., Khanna, K.M.,Tonui, J.K., Murunga G.S.W., Chelimo L.S.,Sirma K. K., Cheruiyot W.K. &Masinde F. W. Department of Physics, University
More informationQuantum Theory of Matter
Quantum Theory of Matter Revision Lecture Derek Lee Imperial College London May 2006 Outline 1 Exam and Revision 2 Quantum Theory of Matter Microscopic theory 3 Summary Outline 1 Exam and Revision 2 Quantum
More informationQuantum mechanics. Effects on the equations of motion of the fractal structures of the geodesics of a nondifferentiable space
http ://luth.obspm.fr/~luthier/nottale/ Quantum mechanics Effects on the equations of motion of the fractal structures of the geodesics of a nondifferentiable space Laurent Nottale CNRS LUTH, Paris-Meudon
More information1 Quantum Theory of Matter
Quantum Theory of Matter: Superfluids & Superconductors Lecturer: Derek Lee Condensed Matter Theory Blackett 809 Tel: 020 7594 7602 dkk.lee@imperial.ac.uk Level 4 course: PT4.5 (Theory Option) http://www.cmth.ph.ic.ac.uk/people/dkk.lee/teach/qtm
More informationarxiv:cond-mat/ v1 13 Mar 1998
Stability of Solution of the Nonlinear Schrödinger Equation for the Bose-Einstein Condensation arxiv:cond-mat/9803174v1 13 Mar 1998 Yeong E. Kim and Alexander L. Zubarev Department of Physics, Purdue University
More information