arxiv:cond-mat/ v1 [cond-mat.other] 24 Mar 2005
|
|
- Mabel Owen
- 5 years ago
- Views:
Transcription
1 arxiv:cond-mat/ v1 [cond-mat.other] 4 Mar 005 Abstract Vortices as Yukawa bosons in the two-dimensional Harmonic confinement K. K. Rajagopal A system of N vortices in a two-dimensional harmonic trap is considered. The vortices are treated as Yukawa bosons. The problem is solved using the Density functional theory within Kohn-Sham scheme. Key words: Bose-Einstein condensation; Vortices and DFT. PACS: Fi, Jp, 3.80.Pj 1 Introduction Nelson et al. [1] have shown that the statistical mechanics of the flux-line lattice (FLL) of high-t c superconductors can be studied through an appropriate mapping onto the D Yukawa boson system interacting with potential V Y (r) = ǫk 0 (r/σ) where ǫ is an energy scale and σ is a length scale having the meaning of a screening length. The function K 0 (x) being the modified Bessel function and the coefficient ǫ attached to it is the coupling strength that scales the energy of interaction. The Hamiltonian of N bosons describing this system reads H = i Λ i + i<j K 0 (r ij ). (1) Here i, j are particle indices and the de Boer dimensionless parameter is defined by Λ = h /mσ ǫ. Meanwhile the reduced density of the system is denoted by ρ = N/Aσ, where A is the surface area. These two dimensionless parameter characterize the system. Corresponding author, bose@gmx.co.uk Preprint submitted to Elsevier Science 6 March 008
2 Liquid Liquid 0.07 Solid Λ ρ 1 10 Fig. 1. The phase diagram of Yukawa bosons reproduced from Magro and Ceperley []. The pluses are transition points computed with DMC. The dashed lines at high densities is the scaling law: Λ 0.04/ ρ. On the other hand, the Hamiltonian for a system of N vortices mapped as a system of N boson interacting with Yukawa potential is written as H = i (k B T) i ǫ + i<j φ 0 8π λ K 0(r ij /λ). () The system described by the Hamiltonian () is equal to one described by Hamiltonian (1) if we set σ = λ and Λ = (πk B T) )/ ǫφ 0. Following this idea, Magro and Ceperley [] have performed a Diffusion Monte Carlo (DMC) and Variational Monte Carlo (VMC) numerical calculation to calibrate the ground state vortex properties for a homogeneous D Yukawa boson system. Of particular interest to us is the phase diagram obtained for the data of (ρ, Λ) at transition points displayed in Fig (1). They observed that for a particular region Λ > 0.09 the system is dominated by kinetic energy and does not crystallize. Below this threshold however a peculiar behaviour of reentrant liquid has been observed. Which means that system is in the liquid phase at very low density as well at high density. The crystal melts on compression and expansion. Nevertheless this reentrant liquid behaviour has not been fully understood. Inspired by their work, I will investigate the ground state property of a system of N vortices mapped as N bosons interacting via D Yukawa potential V (x) = ǫ K 0 (x) confined in a D harmonic planar geometry 1 mω r. As a
3 matter of fact, the Yukawa potential has an interesting combination of short range with soft core. For small x, K 0 (x) diverges like ln(x) while for large x it decays as exp( x)/ x. The logaritmic term represents the inter-vortex repulsion and the exponential decay signifies that vortices far apart do not feel each other. 1.1 Density Functional Theory With the Yukawa potential we cannot use the similar approach used for the contact delta potential V 0 δ(r r ) in deriving a dynamical equation of motion as described in chapter 1. Here we use the Density functional theory (DFT) which is originally based on the notion that for a many-electron system there is a one-to-one mapping between the external potential and the electron density: v ext (r) ρ(r). In other words, the density is uniquely determined given a potential, and vice versa. All properties are therefore a functional of the density, because the density determines the potential, which determines the Hamiltonian, which determines the energy and the wave function. Following this train of reasoning, the inhomogeneous dilute system of N interacting bosons can be described within the second quantization language as Ĥ = Ĥ0 + drψ (r)v ext (r)ψ(r) + 1 dr dr ψ (r)ψ (r )V ( r r )ψ(r )ψ(r) = Ĥ0 + ˆV ext + ˆV int (3) where Ĥ0 = drψ (r) [ h m µ ] ψ(r) and V ( r r ) is the inter-atomic interaction potential. The annihilation and creation field operators are denoted by ψ (r) and ψ(r ) respectively and obey Bose-Einstein commutation relations: [ψ(r), ψ (r )] = δ(r r ) [ψ(r), ψ(r )] = [ψ (r), ψ (r )] = 0. (4) Let us denote the ground state of the system as g so the ground state energy is defined as E 0 = g Ĥ g and the density as n(r) = g ψ ψ g. The Hohenberg-Kohn (HK) theorem [6] guarantees that there exists a unique functional of the density, F[n(r)] = g Ĥ0 + ˆV int g, (5) irrespective of the external potential. The theorem was originally proved for fermions but its generalization also covers bosons. Following HK, we can write 3
4 the total energy functional of the system as following, E[n(r)] = F[n(r)] + drv ext (r)n(r) (6) Determination of the ground state energy follows by imposing the stationary conditions δe[n(r)] = 0 (7) δn(r) In general, the ground state cannot be determined exactly, but one has to resort to the Kohn-Sham procedure [5] to introduce an accurate approximation. The idea is to map the interacting system of interest to an auxiliary system [7]. In order to describe the auxiliary system, let us first decompose the Bose quantum field operators ψ(r) as a sum of spatially varying condensate Φ(r) and a fluctuation field operator ψ : ψ(r) = Φ(r) + ψ(r), ψ (r) = Φ (r) + ψ (r) (8) The wave function Φ(r) is also known as the order parameter and is defined as the statistical average of the particle field operator Φ(r) =< ψ >. On the other hand, the condensate density is defined in term of order parameter as n c (r) = Φ(r). In the well-known Bogoliubov approximation [8] for a dilute weakly interacting gas, in which almost all the atoms are Bose condensed, one keeps only the term up to quadratic in the non-condensate field operators ψ and ψ. When the fluctuation field is comparably small to the spatial order parameter (in an average sense, ψ Φ) one can comfortably utilise a Mean-field approach. Using the decomposition in Eq.(8) above and applying the Bogoliubov approximation, one can replace the product of field operators ψ (r)ψ (r )ψ(r )ψ(r) in the interaction energy term of Eq.(3) by ψ (r)ψ (r )ψ(r )ψ(r)= Φ(r) Φ(r ) + Φ(r) [ Φ(r ) ψ (r ) + Φ (r ) ψ(r ) ] + Φ(r ) ψ (r) ψ(r) + Φ (r )Φ(r) ψ (r) ψ (r ) + 1 Φ (r )Φ (r) ψ(r) ψ(r ) + 1 Φ(r )Φ(r) ψ (r) ψ (r ). (9) So now we are in a position to define the Hamiltonian of the auxiliary system that includes the inter-particle potential defined above as Ĥ s = dr ψ (r) ] [ h m µ ψ(r) + ˆV s + dr ψ (r)vext s (r)ψ(r). (10) where our definition of auxiliary interparticle potential reads 4
5 ˆV s = dr dr V ( r r ) Φ(r) [ Φ(r ) ψ (r ) + Φ (r ) ψ(r ) ] + 1 dr dr V ( r r )[Φ (r )Φ(r) ψ (r) ψ(r ) + Φ(r )Φ(r) ψ (r) ψ (r ) + Φ (r )Φ (r) ψ(r) ψ(r )]. (11) Thus for the ground state g s we can define the unique ground state energy functional of the auxiliary system as F s [n s (r)] = g s drψ (r) ] [ h m µ ψ(r) g s + g s ˆV s g s, (1) which does not depend on the auxiliary external potential V s ext. Hence, the total energy functional of the auxiliary system can be written as E s [n s (r)] = F s [n s (r)] + dr n(r)v s ext(r) (13) in which E s [n s (r)] can be approximated to E[n(r)] or in other word the density of the auxiliary system n s (r) is identical to the real system n(r) by choosing a proper choice of auxiliary external potential Vext. s Since the non-interacting part of Eq. (1) is equal to the non-interacting part of Eq. (5) and comparing terms in Eq.(11) with Eq. (9) we obtain the following functional relation: F[n(r)] = F s [n(r)] + dr dr V (r r )n(r)n(r ) + F xc [n(r)]. (14) The second term ˆV H is called the Hartree-energy defined as ˆV H = 1 dr dr V ( r r )n(r)n(r ) (15) in which the total density n(r) includes a non-condensate local density defined as ñ(r) < ψ (r) ψ(r) >= n(r) Φ(r). (16) The last term in Eq. (14) represents the exchange-correlation energy E xc [n(r)] and includes all the contributions to the interaction energy beyond mean field. Calculating the variational derivatives in Eq. (7) using Eq.(14), one finds δf s [n(r)] δn(r) + V H (r) + δf xc[n(r)] δn(r) + V ext (r) = 0 (17) where the Hartree field read V H (r) = dr V ( r r )n(r ). (18) 5
6 Performing a similar variational calculation on Eq. (13) we deduce that the density of the auxiliary system is identical to the actual system if Vext(r) s = V ext (r) + V H (r) + δf xc[n(r)]. (19) δn(r) 1. Interacting auxiliary system : Bogoliubov approach In general, if an interacting auxiliary system is chosen, one needs to do further analysis by using a Bogoliubov approach. We rewrite the auxiliary Hamiltonian Eq. (10) using the decomposition in Eq. (8) by Ĥ s = + drφ (r)lφ(r) + dr ψ (r)lφ(r) + drφ (r)l ψ(r) dr ψ (r)l ψ(r) + 1 dr dr V Y ( r r )[Φ (r )Φ(r) ψ (r) ψ(r ) + Φ(r )Φ(r) ψ (r) ψ (r ) + Φ (r )Φ (r) ψ(r) ψ(r )] (0) where the operator ˆL is defined by ˆL ( ) m + V ext s µ. (1) In order to diagonalize Eq. (0), we first eliminate the terms linear in ψ and ψ by requiring that Φ(r) satisfy the equation Hence our auxiliary Hamiltonian Eq. (0)reduces to ˆLΦ(r) = 0. () Ĥ s = dr ψ (r)l ψ(r) + dr dr V ( r r )Φ (r )Φ(r) ψ (r) ψ(r ) + 1 dr dr V ( r r )[Φ(r )Φ(r) ψ (r) ψ (r ) + Φ (r )Φ (r) ψ(r) ψ(r )]. (3) The quadratic expression given by this equation can be diagonalized by the usual Bogoliubov transformation [8] ψ(r) = i ψ (r) = i [u j (r)α j v j (r)α j] [u j (r)α j v j (r)α j ], (4) 6
7 where the quasiparticle operators α i and α i satisfy Bose commutation relations. The values of u i and v i are obtained by solving the following coupled Bogoliubov equations: ˆLu i + + dr V ( r r )Φ (r )Φ(r)u i (r ) dr V ( r r )Φ(r )Φ(r)v i (r ) = E i u i (r) (5) and ˆLv i + + dr V ( r r )Φ (r)φ(r )v i (r ) dr V ( r r )Φ (r)φ (r)u i (r ) = E i v i (r). (6) Hence with the solution of u i and v i from the above coupled equations, the Hamiltonian in Eq. (0) can be shown to reduce to Ĥ s new = i E i dr v j (r) + i E i α iα i, (7) which describes a non-interacting gas of quasiparticles of energy E j. The ground state expectation value of Eq. (7) is given by E s 0 = i E i dr v j (r), (8) taking into consideration that the ground state g s obeys α g s = 0. Besides, using Eq. (??), the local non-condensate density defined in Eq. (16) can be written as ñ(r ) g s ψ (r) ψ(r) g s = v i (r). (9) i Thus performing a little algebra, the new auxiliary functional which includes depletion of the condensate is defined by F s [n(r)] = ( ) drφ (r) m µ Φ(r) i dr v i (r) [E j + V s (r)]. (30) The key element in the above derivation is that one can take into account the depletion of the condensate as shown in Eq. (9). Using the the auxiliary system defined by the Hamiltonian Eq. (10) and the relation in Eq (19), both the condensate density n c Φ(r) and also the total density n(r) can be computed. The difficulties in finding the solution for Eq. (), Eqs (5) and (6) lies in the fact that the exact functional dependence of the exchange-correlation 7
8 potential v xc [n(r)] = δf xc [n(r)]/δn(r) is unknown for most of the system of interest and thus one has to resort to approximations such as the Local Density Approximation (LDA). 1.3 Application of DFT to N vortex problem We are now in a position to analyse our problem of N vortices treated as N Yukawa bosons via the one-to-one mapping of Nelson and Seung [1] with the DFT tool summarized in the previous section within our approriate scaling. By choosing an auxiliary system of N non-interacting bosons (referring to Eq. (1) and Eq. (19)), Eq. () then can be written more explicitly as satisfying the following equation: ( h m + V ext (r) + V H (r) + δf ) xc[n(r)] Φ(r) = µφ(r). (31) δn(r) This equation is a sort of generalization of Gross-Pitaevskii equation for a dilute inhomogeneous gas at T = 0. Since we have chosen a non-interacting auxillary system, we expect that at T = 0 all the atoms are Bose condensed as their total density n(r) is approximated by condensate density n c (r) = Φ(r). Setting F xc [n(r)] = 0 one recovers a closed non-linear Schrödinger equation (NLSE). The available numerical data of Magro and Ceperley [] and Strepparola et. al [3] permit one to obtain information on the homogeneous excess free energy f ex [n] rather than the homogeneous exchange correlation energy f xc [n]. Thus it is much convenient to work with excess free energy defined as [9] F ex [n(r)] = ˆV H [n(r)] + F xc [n(r)]. (3) In general, the excess correlation functional energy F ex [n(r)] in a DFT calculation is not known exactly. One can resort to approximations such as the Local density approximation (LDA) which reads, F ex [n(r)] drf hom ex [n] n n(r) (33) where Fex hom [n] = nfex hom [n]. The functional derivatives (excess-correlation potential) of the above relation can be written as [10] V ex (r, n(r)) = δf ex[n(r)] δn(r) = (nf ex[n]) n n(r). (34) n Information of the homogeneous excess correlation energy f ex [n] can be obtained by subtracting the kinetic energy from the total ground state energy of 8
9 a homogeneous system with N boson. A plot of the DMC data of Magro and Ceperley [] and the STLS data of Strepparola et. al. [3] are depicted along with the following fit f ex [ρ] = ρ ln(1/ρ), (35) in Fig (). To be consistent with the scaling of Magro and Ceperley [] I will scale all lengths by harmonic oscillator length σ = h/mω and the energy by ǫ = hω /. Eq. (31) is solved numerically by using the excess correlation potentials of Eq.(34) and Eq. (35) self-consistently with the condition that the areal integral of the condensate density Φ(r) is equal to the total number N of bosons. 9
10 fex[ρ] ρ 1 Fig.. Homogeneous excess-correlation energy f ex (ρ) versus the dimensionless homogenous density ρ gas of the VMC data of Magro and Ceperley [] (pluses) and STLS data of Strepparola et. al. [3] (crosses) compared to the fit function Eq. (35) (dashed-lines) N = 50 n(r) a ho /N N = r/a ho Fig. 3. Density profiles n(r) ( in units a ho /N, with a ho = h/mω ) versus the radial distance r (in units a ho ) obtained using Eq.(31) for N = 50 and N = 100 bosons 10
11 Acknowledgement I would like to acknowledge SNS Pisa for the hospitability and financial support during early part of this work. 11
12 References [1] D. R. Nelson and S. Seung. Phys. Rev. B, 39:9153, [] W. R. Magro and D. M. Ceperley. Phys. Rev. B, 48:411, [3] E. Strepparola and R. Nifosi and M. P. Tosi. J. Phys: Condens. Matter, 10:11645, [4] A. Görlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, T. L. Gustavson, J. R. Abo-Shaeer, A. P. Chikkatur, S. Gupta, S. Inouye, T. Rosenband, and W. Ketterle. Phys. Rev. Lett., 87:13040, 001. [5] W. Kohn and L. J. Sham. Phys. Rev., 140:1133, [6] P. Hohenberg and W. Kohn. Phys. Rev., B136:864, [7] A. Griffin cond-matt/ , [8] N. Bogoliubov J. Phys. USSR, 11:3, [9] S. Moroni and G. Senatore. Phys. Rev. B., 44:9864, [10] A. P. Albus, F. Illuminati and M. Wilkins. Phys. Rev. A., 67:063606,
arxiv:cond-mat/ v1 [cond-mat.stat-mech] 24 Jul 2001
arxiv:cond-mat/010751v1 [cond-mat.stat-mech] 4 Jul 001 Beyond the Thomas-Fermi Approximation for Nonlinear Dynamics of Trapped Bose-Condensed Gases Alexander L. Zubarev and Yeong E. Kim Department of Physics,
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 informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationarxiv:cond-mat/ v1 2 Apr 1997
Thermodynamics of a trapped Bose-condensed gas S. Giorgini 1,2, L. P. Pitaevskii 2,3,4 and S. Stringari 2 1 ECT, Villa Tambosi, Strada delle Tabarelle 286, I-38050 Villazzano, Italy 2 Dipartimento di Fisica,
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 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 informationDensity matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
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 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 information2.5 Time dependent density functional theory
.5 Time dependent density functional theory The main theorems of Kohn-Sham DFT state that: 1. Every observable is a functional of the density (Hohenger-Kohn theorem).. The density can be obtained within
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More information1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More informationVortices in Atomic Bose-Einstein Condensates in the large gas parameter region. Abstract
Vortices in Atomic Bose-Einstein Condensates in the large gas parameter region J. K. Nilsen 1), J. Mur-Petit 2), M. Guilleumas 2), M. Hjorth-Jensen 1), and A.Polls 2) 1 ) Department of Physics, University
More informationIntroduction to density-functional theory. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
More informationDensity functional theory
Density functional theory Seminar Ib - 1. letnik, II. stopnja Author: Matic Pečovnik Mentor: Doc. Dr. Tomaž Rejec Faculty of Mathematics and Physics, University in Ljubljana Abstract Many-body problems
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 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 informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 2 Aug 2004
Ground state energy of a homogeneous Bose-Einstein condensate beyond Bogoliubov Christoph Weiss and André Eckardt Institut für Physik, Carl von Ossietzky Universität, D-6 Oldenburg, Germany (Dated: November
More informationInauguration Meeting & Celebration of Lev Pitaevskii s 70 th Birthday. Bogoliubov excitations. with and without an optical lattice.
Inauguration Meeting & Celebration of Lev Pitaevskii s 7 th Birthday Bogoliubov excitations with and without an optical lattice Chiara Menotti OUTLINE OF THE TALK Bogoliubov theory: uniform system harmonic
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
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 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 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 informationarxiv:cond-mat/ v2 [cond-mat.soft] 7 Mar 2005
Center of mass rotation and vortices in an attractive Bose gas A. Collin 1, E. Lundh,1, K.-A. Suominen 3 1 Helsinki Institute of Physics, PL 64, arxiv:cond-mat/040336v [cond-mat.soft] 7 Mar 005 FIN-00014
More informationˆΨ( x) = 1 V. a k exp(i k x), (1) d 3 xd 3 x ˆΨ ( x) ˆΨ( x )V ( x x ) ˆΨ( x ) ˆΨ( x). (2) [a k, a k ] = δ k, k (3)
Dilute Boson Gas In this section we consider a gas of spinless bosons of mass m interacting through a static potential. The treatment given here is not realistic enough to quantitatively describe liquid
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 informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
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 informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 17 Feb 2003
Beyond Fermi pseudopotential: a modified GP equation arxiv:cond-mat/030331v1 [cond-mat.stat-mech 17 Feb 003 Haixiang Fu, 1 Yuzhu Wang, 1 and Bo Gao 1, 1 Shanghai Institute of Optics and Fine Mechanics,
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationQuantum Monte Carlo methods
Quantum Monte Carlo methods Lubos Mitas North Carolina State University Urbana, August 2006 Lubos_Mitas@ncsu.edu H= 1 2 i i 2 i, I Z I r ii i j 1 r ij E ion ion H r 1, r 2,... =E r 1, r 2,... - ground
More informationA Hartree - Fock study of a Bose condensed gas
Home Search Collections Journals About Contact us My IOPscience A Hartree - Fock study of a Bose condensed gas This article has been downloaded from IOPscience. Please scroll down to see the full text
More informationarxiv: v2 [cond-mat.mes-hall] 11 Sep 2007
Vortex pump for dilute Bose-Einstein condensates Mikko Möttönen, 1,2 Ville Pietilä, 1 and Sami M. M. Virtanen 1 1 Laboratory of Physics, Helsinki University of Technology P. O. Box 4100, FI-02015 TKK,
More informationarxiv: v1 [cond-mat.str-el] 18 Jul 2007
arxiv:0707.2704v1 [cond-mat.str-el] 18 Jul 2007 Determination of the Mott insulating transition by the multi-reference density functional theory 1. Introduction K. Kusakabe Graduate School of Engineering
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 informationIntroduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionals by Michel Côté Université de Montréal Département de physique Eamples Carbazole molecule Inside of diamant Réf: Jean-François Brière http://www.phys.umontreal.ca/~michel_
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
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 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 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 informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationHexatic and microemulsion phases in a 2D quantum Coulomb gas
Hexatic and microemulsion phases in a 2D quantum Coulomb gas Bryan K Clark (University of Illinois at Urbana Champaign) Michele Casula (Ecole Polytechnique, Paris) David M Ceperley (University of Illinois
More informationFirst-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design
First-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design Peter Kratzer 5/2/2018 Peter Kratzer Abeokuta School 5/2/2018 1 / 34 Outline
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 informationQuantum Monte Carlo Simulations of Exciton Condensates
Quantum Monte Carlo Simulations of Exciton Condensates J. Shumway a and D. M. Ceperley b a Dept. of Physics and Astronomy, Arizona State University, Tempe, AZ 8583 b Dept. of Physics, University of Illinois,
More informationBose-Einstein Condensation
Bose-Einstein Condensation Kim-Louis Simmoteit June 2, 28 Contents Introduction 2 Condensation of Trapped Ideal Bose Gas 2 2. Trapped Bose Gas........................ 2 2.2 Phase Transition.........................
More informationAdvanced Solid State Theory SS Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt
Advanced Solid State Theory SS 2010 Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt i 0. Literatur R. M. Martin, Electronic Structure: Basic Theory and
More information2 Quantization of the Electromagnetic Field
2 Quantization of the Electromagnetic Field 2.1 Basics Starting point of the quantization of the electromagnetic field are Maxwell s equations in the vacuum (source free): where B = µ 0 H, D = ε 0 E, µ
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
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 informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 12 Jan 2007
Stability and dynamics of vortex clusters in nonrotated Bose-Einstein condensates arxiv:cond-mat/0701285v1 [cond-mat.stat-mech] 12 Jan 2007 V. Pietilä, 1 M. Möttönen, 1,2 T. Isoshima, 3 J. A. M. Huhtamäki,
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 informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
More informationCLASSICAL REPRESENTATION OF QUANTUM SYSTEMS AT EQUILIBRIUM
CLASSICAL REPRESENTATION OF QUANTUM SYSTEMS AT EQUILIBRIUM By SANDIPAN DUTTA A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR
More informationSpacetime analogue of Bose-Einstein condensates
Spacetime analogue of Bose-Einstein condensates Bogoliubov-de Gennes formulation Hideki ISHIHARA Osaka City Univ., JAPAN Y.Kurita, M.Kobayashi, T.Morinari, M.Tsubota, and H.I., Phys. Rev. A79, 043616 (2009)
More informationDipolar Bose-Einstein condensates at finite temperature
PHYSICAL REVIEW A 76, 4367 7 Dipolar Bose-Einstein condensates at finite temperature Shai Ronen JILA and Department of Physics, University of Colorado, Boulder, Colorado 839-44, USA John L. Bohn* JILA,
More informationBogoliubov theory of the weakly interacting Bose gas
Chapter 4 Bogoliubov theory of the wealy interacting Bose gas In the presence of BEC, the ideal Bose gas has a constant pressure against variation of volume, so that the system features infinite compressibility.
More informationQuantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid
Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge
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 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 informationGreen s Function Approach to B. E. Condensation in the alkali atoms
arxiv:quant-ph/010048v1 8 Oct 00 Green s Function Approach to B. E. Condensation in the alkali atoms Mahendra Sinha Roy October 9, 018 Department of Physics, Presidency College, Calcutta 700 073,India.
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 informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
More informationSupplementary Figure 3: Interaction effects in the proposed state preparation with Bloch oscillations. The numerical results are obtained by
Supplementary Figure : Bandstructure of the spin-dependent hexagonal lattice. The lattice depth used here is V 0 = E rec, E rec the single photon recoil energy. In a and b, we choose the spin dependence
More informationarxiv: v1 [cond-mat.quant-gas] 11 Jul 2015
Testing the nonlocal kinetic energy functional of an inhomogeneous, two-dimensional degenerate Fermi gas within the average density approximation J. Towers Jack Dodd Centre for Quantum Technology, Department
More information1 Back to the ground-state: electron gas
1 Back to the ground-state: electron gas Manfred Lein and Eberhard K. U. Gross 1.1 Introduction In this chapter, we explore how concepts of time-dependent density functional theory can be useful in the
More informationWill be published: Phys. Rev. Lett. 96, (2006)
Will be published: Phys. Rev. Lett. 96, 230402 (2006) Vortex-lattice melting in a one-dimensional optical lattice Michiel Snoek and H. T. C. Stoof Institute for Theoretical Physics, Utrecht University,
More informationSqueezing and superposing many-body states of Bose gases in confining potentials
Squeezing and superposing many-body states of Bose gases in confining potentials K. B. Whaley Department of Chemistry, Kenneth S. Pitzer Center for Theoretical Chemistry, Berkeley Quantum Information and
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 informationTwo-Dimensional Spin-Polarized Hydrogen at Zero
J Low Temp Phys (2013) 171:685 692 DOI 10.1007/s10909-012-0756-7 Two-Dimensional Spin-Polarized Hydrogen at Zero Temperature L. Vranješ Markić J. Boronat Received: 13 July 2012 / Accepted: 16 September
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationIntroduction to Path Integral Monte Carlo. Part I.
Introduction to Path Integral Monte Carlo. Part I. Alexey Filinov, Jens Böning, Michael Bonitz Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Sep 2005
1 arxiv:cond-mat/0509103v1 [cond-mat.stat-mech] 5 Sep 2005 Normal and Anomalous Averages for Systems with Bose-Einstein Condensate V.I. Yukalov 1,2 and E.P. Yukalova 1,3 1 Institut für Theoretische Physik,
More informationAll electron optimized effective potential method for solids
All electron optimized effective potential method for solids Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. 22
More informationBardeen Bardeen, Cooper Cooper and Schrieffer and Schrieffer 1957
Unexpected aspects of large amplitude nuclear collective motion Aurel Bulgac University of Washington Collaborators: Sukjin YOON (UW) Kenneth J. ROCHE (ORNL) Yongle YU (now at Wuhan Institute of Physics
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 informationarxiv:cond-mat/ v1 28 Feb 1996
Bose-Einstein Condensation of Atoms in a Trap T.T. Chou, Chen Ning Yang 2 and L.H. Yu 3 arxiv:cond-mat/960253v 28 Feb 996 Department of Physics, University of Georgia, Athens, Georgia 30602 2 Institute
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 3 Aug 2004
APS/1-QED Population and phase dynamics of F = 1 spinor condensates in an external magnetic field arxiv:cond-mat/4869v1 [cond-mat.stat-mech] Aug 4 D. R. Romano and E. J. V. de Passos Instituto de Física,
More information1 Quantum field theory and Green s function
1 Quantum field theory and Green s function Condensed matter physics studies systems with large numbers of identical particles (e.g. electrons, phonons, photons) at finite temperature. Quantum field theory
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 informationMechanism of Tunneling in Interacting Open Ultracold Few-Boson Systems. Abstract
Mechanism of Tunneling in Interacting Open Ultracold Few-Boson Systems Axel U. J. Lode, Alexej I. Streltsov, Ofir E. Alon, and Lorenz S. Cederbaum arxiv:1005.0093v1 [cond-mat.quant-gas] 1 May 2010 Theoretische
More informationWe can instead solve the problem algebraically by introducing up and down ladder operators b + and b
Physics 17c: Statistical Mechanics Second Quantization Ladder Operators in the SHO It is useful to first review the use of ladder operators in the simple harmonic oscillator. Here I present the bare bones
More informationBose Gases, Bose Einstein Condensation, and the Bogoliubov Approximation
Bose Gases, Bose Einstein Condensation, and the Bogoliubov Approximation Robert Seiringer IST Austria Mathematical Horizons for Quantum Physics IMS Singapore, September 18, 2013 R. Seiringer Bose Gases,
More informationThe phonon dispersion relation of a Bose-Einstein condensate
The phonon dispersion relation of a Bose-Einstein condensate I. Shammass, 1 S. Rinott, 2 A. Berkovitz, 2 R. Schley, 2 and J. Steinhauer 2 1 Department of Condensed Matter Physics, Weizmann Institute of
More informationTime Dependent Density Functional Theory. Francesco Sottile
An Introduction Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) Belfast, 29 Jun 2007 Outline 1 Introduction: why TD-DFT? 2 (Just)
More informationPseudopotential generation and test by the ld1.x atomic code: an introduction
and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)
More informationElectronic Structure Calculations and Density Functional Theory
Electronic Structure Calculations and Density Functional Theory Rodolphe Vuilleumier Pôle de chimie théorique Département de chimie de l ENS CNRS Ecole normale supérieure UPMC Formation ModPhyChem Lyon,
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 informationRabi oscillations within TDDFT: the example of the 2 site Hubbard model
Rabi oscillations within TDDFT: the example of the 2 site Hubbard model Johanna I. Fuks, H. Appel, Mehdi,I.V. Tokatly, A. Rubio Donostia- San Sebastian University of Basque Country (UPV) Outline Rabi oscillations
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 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 informationFinite-Temperature Hartree-Fock Exchange and Exchange- Correlation Free Energy Functionals. Travis Sjostrom. IPAM 2012 Workshop IV
1 of 45 Finite-Temperature Hartree-Fock Exchange and Exchange- Correlation Free Energy Functionals Travis Sjostrom Quantum Theory Project Depts. of Physics and Chemistry IPAM 2012 Workshop IV 2012 2 of
More informationin-medium pair wave functions the Cooper pair wave function the superconducting order parameter anomalous averages of the field operators
(by A. A. Shanenko) in-medium wave functions in-medium pair-wave functions and spatial pair particle correlations momentum condensation and ODLRO (off-diagonal long range order) U(1) symmetry breaking
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 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 condensation in trapped bosons: A variational Monte Carlo analysis
PHYSICAL REVIEW A, VOLUME 63, 023602 Bose-Einstein condensation in trapped bosons: A variational Monte Carlo analysis J. L. DuBois and H. R. Glyde Department of Physics and Astronomy, University of Delaware,
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 information