arxiv: v1 [cond-mat.quant-gas] 2 Jun 2014
|
|
- Byron Eaton
- 6 years ago
- Views:
Transcription
1 Emergence of pairing in a system of a few strongly attractive fermions in a one-dimensional harmonic trap Tomasz Sowiński,2, Mariusz Gajda,2, and Kazimierz Rz ażewski 2 Institute of Physics of the Polish Academy of Sciences, Al. Lotników 32/46, Warszawa, Poland 2 Center for Theoretical Physics of the Polish Academy of Sciences, Al. Lotników 32/46, Warszawa, Poland (Dated: June 5, 24) We study a strongly attractive system of a few spin-/2 fermions confined in a one-dimensional harmonic trap, interacting via two-body contact potential. Performing exact diagonalization of the Hamiltonian we find that for strong attractions the correlated pairs occur in the system. We find that the fraction of correlated pairs depends on temperature and we show that it has universal properties analogous to the gap function known from the theory of superconductivity. Our predictions can be verified in current experiments on ultra-cold atoms. PACS numbers: Lm arxiv:46.4v [cond-mat.quant-gas] 2 Jun 24 Introduction. Quantum engineering, i.e. coherent controlling of atomic systems, is a rapidly developing field of modern physics. Its development is important because it gives us the possibility of a deep study of quantum mechanics. Present experimental methods in quantum engineering allow one to prepare and control the system in a state with a rigorously chosen number of interacting atoms. Spectacular experiments on ultra-cold fewbody fermionic systems confined in a one-dimensional harmonic trap [ 4] have opened alternative way to study strongly correlated quantum systems. It is believed that these kind of studies can shed a new light on our understanding of the crossover from the few-body world to the situation with a macroscopic number of quantum particles [4]. They can help to answer the fundamental question on emerging of collective behavior from the two-body interactions. Recently, harmonically trapped effectively one dimensional few-body systems were intensively studied on theoretical [5 9] and experimental [ 4] level. For example, it is known that in a strictly one-dimensional geometry, in the limit of infinite repulsion (known as Tonks-Girardeau limit) the bosonic system can be mapped on the system of noninteracting fermions [ 3]. Moreover, in this limit, the ground-state of the one-dimensional system of a few spin-/2 fermions is highly degenerated [6, 9] and a combination which is fully antisymmetric with respect to the position of fermions is identical to the ground state of the system of spinless fermions [5, 4]. This degeneracy of the ground-state and its consequences for the experimental results were also studied [9]. The mapping mentioned above, in the limit of infinite repulsion was probed experimentally for two-fermion system [2], for which the exact solution is known [5 8]. In the recent paper [3] the attractive system was probed experimentally and formation of fermionic pairs was observed. This experiment opened an alternative path to the experimental exploration of the few-body systems. Therefore the quantitative predictions based on the first-principle theory of a few-body system are necessary. We should mentioned here that pairing in a system of fermionic atoms in the mean-filed Hertree-Fock approach was studied recently in [9, 2]. The Model. In this Letter we study a few fermions of mass m, confined in a one-dimensional harmonic trap of frequency Ω, mutually interacting via short-range deltalike attractive potential. In the second-quantization formalism the Hamiltonian of the model has a form: Ĥ = [ dx ˆΨ σ(x) d 2 2 dx 2 + ] 2 x2 ˆΨ σ (x) σ {, } + g dx dx ˆΨ (x) ˆΨ (x) ˆΨ (x) ˆΨ (x), () where Ψ σ (x) annihilates a fermion with spin σ at a point x, and g measures a strength of the contact interactions. For convenience we expressed all quantities in harmonic oscillator units, i.e. energies in hω, lengths in h/mω, temperatures in hω/k B etc. The dimensionless interaction constant g = (m/ h 3 Ω) /2 g D is proportional to the effective one-dimensional interaction g D between fermions of opposite spins. The Hamiltonian () commutes separately with the number of fermions in each spin component ˆN σ = dx Ψ σ(x)ψ σ (x). Here, we consider a situation where both species are distinguishable. Therefore, [ the field operators of different spins do commute, ˆΨ (x), ˆΨ ] [ (x ) = ˆΨ (x), ˆΨ (x )] =. There are no any symmetry conditions related to the exchange of particles with different spins. However, particles of the same { spin satisfy the standard anticommutation relations ˆΨσ (x), ˆΨ } σ(x ) = δ(x x ) and { ˆΨσ (x), ˆΨ σ (x )} =. We focus our attention on correlations between fermions with opposite spins in a balanced system N = N = N/2. The Method. We fix the number of fermions in each component, then we decompose the field operator Ψ σ (x) in the single-particle basis φ n (x) of harmonic oscillator eigenstates with sufficiently large cutoff, n n max, and we diagonalize the full many-body Hamiltonian () in this basis via Implicitly Restarted Arnoldi Method [2]. In this way we find N eigenstates G i and the corresponding lowest eigenenergies E i. In our case N = 6.
2 2 Properties of the ground state. To get some intuition, first we discuss properties of the system at zero temperature, i.e. we analyze the ground-state G. In the limit of vanishing interaction g, the ground state of the system is simply formed by filling up the lowest single-particle states of the harmonic potentials by pairs of opposite spin fermions. This fact is directly manifested in the spectral decomposition of the reduced two-fermion density matrix ρ (2) (x, x 2 ; x, x 2) = G ˆΨ (x ) ˆΨ (x 2) ˆΨ (x 2) ˆΨ (x ) G, (2) (normalized to Tr ρ (2) = ). Let us notice that coordinates x and x (x 2 and x 2 ) correspond to the spin-down (spin-up) fermions. From now on we will use this convention throughout our Letter. Distinction between the two species is possible only because in our Hamiltonian there are no terms responsible for spin flipping mixing different species. All nonzero eigenvalues of ρ (2) are equal to λ i = /N N. They determine the relative populations of the corresponding eigenvectors, i.e. the two-particle orbitals ψ i (x, x 2 ) describing possible pairings of fermions with opposite spins. Situation changes dramatically when attractive interaction are switched on. Then, we observe that one eigenvalue starts to dominate in the system (Fig. ). According to the Penrose-Onsager criterion, an occurrence of the dominant orbital in the spectral decomposition of the reduced density matrix is an indication of the condensation of fermionic pairs in the selected orbital. For convenience we denote the dominant orbital and its eigenvalue by ψ and λ respectively. To get better understanding of the structure of the ground state of the system it is convenient to relate the orbitals ψ i (x, x 2 ) to the orbitals ϕ i (x) the eigenvectors of the reduced single-fermion density matrices (normalized to one) ρ () (x; x ) = G ˆΨ σ(x) ˆΨ σ (x ) G, (3) where σ {, }. Since we assumed exact symmetry between up and down fermions, ρ () is exactly the same for both spin components. By projecting two-fermion orbitals of ρ (2) on the product of the single-fermion orbitals of ρ () we find that the set of all orbitals ψ i (x, x 2 ) can be divided into two classes. For a given two-fermion orbital ψ i belonging to the first class one can find exactly two single-fermion orbitals ϕ j and ϕ k (j k) such that ψ i (x, x 2 ) = [ϕ j (x )ϕ k (x 2 ) ± ϕ k (x )ϕ j (x 2 )] / 2. (4a) Note that the two-fermion orbital has no any symmetry with respect of interchanging spin up and spin down fermions, i.e. x with x 2. In the second class each twofermion orbital ψ i (x, x 2 ) has a structure of the correlated pair the superposition of states with two particles of opposite spin occupying the same singe-particle state: ψ i (x, x 2 ) = j α i (j) ϕ j (x )ϕ j (x 2 ), (4b) where j α i(j) 2 =. From our numerical analysis it follows that the dominant orbital, ψ (x, x 2 ), always belongs to the second class. Moreover it is symmetric with respect to exchange of spin-up and spin-down fermions. Moreover, the structure of the dominant orbital is unique in the sense that only in this orbital we find non vanishing expectation value of the operator  = i,j,k ˆb iˆb jˆb kˆb k, (5) where the operator ˆb σi annihilates fermion in the orbital ϕ i of the reduced single-fermion density matrix ρ (). Expectation value S = G  G is a direct coun- of the superconducting order parameter S = terpart dr ˆΨ (r)ψ (r) known from the BCS theory [22]. Operator  can be understood as describing the process of destruction of a correlated fermionic pair occupying the same one-particle state and putting the two fermions anywhere. This operator can be viewed as the number conserving destruction operator of the correlated pair. The strong correlation in the dominant orbital can be observed in the probability density of finding fermions with opposite spins (Fig. 2). In the position representation we find strong correlation between x and x 2. In contrast, in the momentum space we find strong anti-correlation between p and p 2. Moreover, for strong enough attraction we find that probability of finding a pair with opposite momenta does not depend on the value of the momentum (Fig. 3). The dominant orbital supports such correlated pairs that both particles are very likely to be found in the same position with opposite momenta. This observation is one of the fundamental features of correlated pairs known from the theory of superconductivity. Therefore, the dominant orbital ψ = j α (j) ˆb jˆb j Ω is an analog of a Cooper pair [23]. Non zero temperatures. Now, we discuss a situation when the system studied is in contact with a thermostat of given temperature T. Therefore, instead of the manybody ground state G, we assume that the state of the system is well described by the density matrix operator ˆρ T = Z exp( Ĥ/T ), where Z is a canonical partition function assuring a proper normalization Tr ˆρ T =. Numerically, the density matrix operator ˆρ T is constructed from a finite number of N = 6 eigenstates G i of the full many-body Hamiltonian as follows ˆρ T = i p i G i G i, where p i = exp( E i /T )/ i exp( E i/t ). In analogy to the zero temperature case, we define a reduced twofermion density matrix [ ρ (2) T (x, x ; x 2, x 2) = Tr ˆρ T ˆΨ (x ) ˆΨ (x 2) ˆΨ (x 2) ˆΨ ] (x ) (6) and we diagonalize it. We perform the whole spectral analysis of ρ (2) T in analogy to the zero temperature limit. First, we find temperature dependent single-particle orbitals annihilated by temperature dependent operators
3 3 Eigenvalue N = Eigenvalue Interaction g N = 8 FIG. 2: The density distribution of finding fermions with opposite spins in the dominate orbital in the position domain ψ (x, x 2) 2 (left panel) and in the momentum domain ψ (p, p 2) 2 (right panel). We observe that for stronger attractions the positions of paired fermions are strongly correlated. In contrast strong anti-correlation of momenta is present. FIG. : First eight eigenvalues λ i of the reduced two-fermion density matrix ρ (2) as a function of interaction g, calculated in the ground-state of the system obtained numerically for N = N = 3 and 4. For attractive interactions one of the eigenvalues starts to dominate in the system indicating occurrence of the fermionic pair condensation. ˆbσi. These states serve as a basis for two-fermion orbitals of ρ (2) T. We find that structural properties of the orbitals are not affected by finite temperatures, i.e. the division of all orbitals into two classes with respect to their representation by finite-temperature single-fermion orbitals holds. In particular coefficients α (j) practically do not depend on the temperature. We also checked that the dominant orbital of ρ (2) T remains the only orbital with non-vanishing expectation value of the particle conserving destruction operator  (defined with temperature dependent ˆb σi ). It means [ that] the superconducting order parameter S(T ) = Tr ˆρ (2) T  is well defined and can be directly related to the superconducting fraction represented by the dominant eigenvalue at given temperature λ (T ): S = λ (T ) λ () S. (7) With increasing temperature we observe a rapid decay of the occupation of the dominant orbital, and in consequence, the superconducting fraction tends to zero. It is worth to notice that the dominant eigenvalue λ (T ) of the reduced two-fermion density matrix for the system studied is closely related to the gap function known from the theory of superconductivity for isotropic system of interacting fermions in the thermodynamic limit. probability density position x momentum p FIG. 3: Section of the density distributions plotted in Fig.2 (red dashed line for g = and blue solid line for g = ) along the lines x = x 2 and p = p 2 for position and momentum domain respectively. We observe that for strong attraction the density distribution in momentum space becomes flat. This fact shows that probability of finding correlated fermions with opposite momentum does not depend on the value of the momentum. This is one of the features of superconducting fraction well known in the many-body theory of superconductivity. Therefore it is quite obvious that the behavior of λ (T ) is in one-to-one correspondence with the behavior of the gap function known from the BCS theory. In Fig. 4a we plot a temperature dependence of the superconducting fraction λ (T ) for N = N = 4 and different interactions g. As it is seen, this quantity crucially depends on the interaction constant g. Nevertheless, it has some fundamental universality which indicates correspondence to the BCS theory. To show this universality we define two natural quantities which depend on g, i.e. superconducting fraction at zero temperature, λ (), and the characteristic temperature T at which, by definition, the superconducting fraction drops by some chosen fac-
4 4 (a) λ (T).7.6 (b) λ (T)/λ () T (c) (d) T/T CDF Probability cutoff level n cutoff level n 2 3 T g ln[λ ()] / g FIG. 4: Universal properties of the superconducting fraction represented by dominant eigenvalue λ (T ) for N = N = 4. (a) Dominant eigenvalue λ (T ) of the two-fermion density matrix ρ (2) T as a function of temperature for different twobody interaction constant g. (b) Rescaled λ (T )/λ () as a function of rescaled temperature T/T (see the main text for definition) for different two-body interaction constant g. All numerical points collapse to the one universal curve, which confirms universal behavior of λ (T ). (c) Characteristic temperature T as a function of interaction. Numerical data points fit almost perfectly to the linear regression. (d) Zero temperature superconducting fraction λ () as a function of interaction. In chosen scaling numerical data points fit to the linear regression, which confirms universal behavior (8) predicted by the theory of superconductivity for uniform and infinite systems. tor η = λ (T )/λ (). We checked that our predictions are insensitive to any reasonable choice of the value of η and in the following we set η = 85%. Then, by plotting the normalized data, i.e. λ (T )/λ () vs. T/T, we find that all numerical points collapse to the universal curve which does not depend on model details (see Fig.4b). We checked that this universal behavior of λ (T ) is independent of the definition of the characteristic temperature T. Moreover, we find that the characteristic temperature T is a linear function of interaction g (Fig. 4c). On the other hand, the dependence of the zero temperature superconducting fraction λ () on the interaction constant g is exactly the same as the dependence of the gap-function in the BCS theory (Fig. 4d) λ () = Λe /M g. (8) Amazingly, the value of the parameter M obtained from the numerical fit to the data is equal to 2. ±.3, which is in perfect agreement with the intuitive picture of the theory of superconductivity. In the case of homogenous and infinite system, the theory predicts that FIG. 5: Two-particle cumulative distribution function (CDF) for few first orbitals ψ i at T = and g = for N = 8 (the gray scale corresponds to the orbital index i as indicated in the legend). The hight of the peak is equal to the probability of finding simultaneously two particles (described by ψ i) above chosen harmonic oscillator level n. In the inset we show the probability that the pair detected above a level n comes from the dominant strongly correlated orbital ψ. Note, that for large n the pair, if detected, is the Cooper-like pair with probability approaching one. M should be equal to the density of states at the surface of the Fermi sea. In the one-dimensional case studied here, the size of the Fermi surface does not depend on the number of particles and is simply equal to 2 obviously, the Fermi level can be occupied only by two fermions with opposite spins. Finally, let us notice that dependence of the characteristic temperature T on g is linear and different than in the BCS theory for homogeneous and 3D system. This comes from the fact that the number of particles here is small and the system is one-dimmensional. Final remarks. We believe that our predictions can be verified in current experiments with ultra-cold fermions in optical traps. By tilting the trap potential (like in experiments of Jochim group []) the particles of energies above selected harmonic oscillator level are allowed to escape and can be detected. The probability of simultaneous detection of two particles above selected energy n is related to the two-particle cumulative distribution function C(n) (CDF). In Fig. 5 the CDF for a few first orbitals ψ i of the two-fermion density matrix ρ (2) (i =,..., 3) at T = and g = for N = 8 is shown. The gray scale corresponds to the orbital index i as indicated in the legend. The hight of the peak is equal to the probability of finding simultaneously two particles (described by ψ i ) above chosen harmonic oscillator level n. In the inset we show the probability that the pair detected above a level n comes from the dominant strongly correlated orbital ψ. Note, that for large n the pair, if detected, is the Cooper-like pair with probability approaching one. The figure clearly shows that if the pair is detected above the harmonic oscillator level n = 9 (both particles must have energies larger than n) then with probability about 85% it has to be the correlated Cooper-like pair occupying initially the dominant orbital. This way the correlated
5 5 pair can be distilled from the trap. Measuring the two particle correlation function of the pair in the momentum domain (Fig. 2) is the ultimate proof of detection of the Cooper pair. We checked that this is true also at finite temperatures. Acnowledgements. This research was supported by the (Polish) National Science Center Grant No. DEC- 22/4/A/ST2/9. [] F. Serwane, G. Zürn, T. Lompe, T. B. Ottenstein, A. N. Wenz, and S. Jochim, Science 332, 336 (2). [2] G. Zürn et al., Phys. Rev. Lett. 8, 7533 (22). [3] G. Zürn et al., Phys. Rev. Lett., 7532 (23). [4] A. N. Wenz et al., Science 342, 457 (23). [5] L. Guan, S. Chen, Phys. Rev. Lett. 5, 753 (2). [6] L. Guan, S. Chen, Y. Wang, Z.-Q. Ma, Phys. Rev. Lett. 2, 642 (29). [7] X. Cui, T.-L. Ho, Phys. Rev. A 89, 236 (24). [8] S. E. Gharashi, D. Blume, Phys. Rev. Lett., 4532 (23). [9] T. Sowiński, T. Grass, O. Dutta, M. Lewenstein, Phys. Rev. A 88, 3367 (23). [] L. Tonks, Phys. Rev. 5, 955 (936). [] M. Girardeau, J. Math. Phys., 56 (96). [2] E. Lieb, W. Liniger, Phys. Rev. 3, 65 (963). [3] B. Paredes et al., Nature (London) 429, 277 (24). [4] M. D. Girardeau, Phys. Rev. A 82, 67(R) (2). [5] T. Sowiński, M.Brewczyk, M.Gajda, K. Rz ażewski, Phys. Rev. A 82, 5363 (2). [6] T. Busch, B.-G. Englert, K. Rz ażewski, M. Wilkens, Found. Phys. 28, 549 (998). [7] M. Rontani, Phys. Rev. Lett. 8, 532 (22). [8] Z. Idziaszek, T. Calarco, Phys. Rev. A 74, 2272 (26). [9] T. Karpiuk, M. Brewczyk, K. Rz ażewski, Phys. Rev. A 69, 4363 (24). [2] T. Świs locki, T. Karpiuk, M. Brewczyk, Phys. Rev. A 77, 3363 (28). [2] R. B. Lehoucq, D. C. Sorensen, C. Yang, Arpack User s Guide: Solution of Large-Scale Eigenvalue Problems With Implicityly Restorted Arnoldi Methods (Society for Industrial & Applied, 998) [22] J. Bardeen, L. N. Cooper, J. R. Schrieffer, Phys. Rev. 8, 75 (957). [23] L. N. Cooper, Phys. Rev. 4, 89 (956).
arxiv: v1 [cond-mat.quant-gas] 18 Sep 2015
Slightly imbalanced system of a few attractive fermions in a one-dimensional harmonic trap Tomasz Sowiński arxiv:1509.05515v1 [cond-mat.quant-gas] 18 Sep 2015 Institute of Physics of the Polish Academy
More informationarxiv: v1 [quant-ph] 8 Sep 2010
Few-Body Systems, (8) Few- Body Systems c by Springer-Verlag 8 Printed in Austria arxiv:9.48v [quant-ph] 8 Sep Two-boson Correlations in Various One-dimensional Traps A. Okopińska, P. Kościk Institute
More informationSecond Quantization: Quantum Fields
Second Quantization: Quantum Fields Bosons and Fermions Let X j stand for the coordinate and spin subscript (if any) of the j-th particle, so that the vector of state Ψ of N particles has the form Ψ Ψ(X
More informationarxiv: v2 [cond-mat.quant-gas] 3 Dec 2017
Coherent Oscillations in Small Fermi Polaron Systems arxiv:1706.09466v2 [cond-mat.quant-gas] 3 Dec 2017 Marek Tylutki, 1, 2, G. E. Astrakharchik, 3 and Alessio Recati 4, 5 1 INO-CNR BEC Center and Dipartimento
More informationIsotropic harmonic oscillator
Isotropic harmonic oscillator 1 Isotropic harmonic oscillator The hamiltonian of the isotropic harmonic oscillator is H = h m + 1 mω r (1) = [ h d m dρ + 1 ] m ω ρ, () ρ=x,y,z a sum of three one-dimensional
More informationLecture 6. Fermion Pairing. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 6 Fermion Pairing WS2010/11: Introduction to Nuclear and Particle Physics Experimental indications for Cooper-Pairing Solid state physics: Pairing of electrons near the Fermi surface with antiparallel
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 informationThe Dirac Equation. Topic 3 Spinors, Fermion Fields, Dirac Fields Lecture 13
The Dirac Equation Dirac s discovery of a relativistic wave equation for the electron was published in 1928 soon after the concept of intrisic spin angular momentum was proposed by Goudsmit and Uhlenbeck
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 informationarxiv:quant-ph/ v5 10 Feb 2003
Quantum entanglement of identical particles Yu Shi Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom and Theory of
More informationarxiv: v2 [cond-mat.str-el] 9 Jul 2012
Exact diagonalization of the one dimensional Bose-Hubbard model with local -body interactions Tomasz Sowiński Institute of Physics of the Polish Academy of Sciences, Aleja Lotników /46, -668 Warsaw, Poland
More informationMany Body Quantum Mechanics
Many Body Quantum Mechanics In this section, we set up the many body formalism for quantum systems. This is useful in any problem involving identical particles. For example, it automatically takes care
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 informationFor example, in one dimension if we had two particles in a one-dimensional infinite potential well described by the following two wave functions.
Identical particles In classical physics one can label particles in such a way as to leave the dynamics unaltered or follow the trajectory of the particles say by making a movie with a fast camera. Thus
More informationPHY 396 K. Problem set #5. Due October 9, 2008.
PHY 396 K. Problem set #5. Due October 9, 2008.. First, an exercise in bosonic commutation relations [â α, â β = 0, [â α, â β = 0, [â α, â β = δ αβ. ( (a Calculate the commutators [â αâ β, â γ, [â αâ β,
More information221B Lecture Notes Quantum Field Theory II (Fermi Systems)
1B Lecture Notes Quantum Field Theory II (Fermi Systems) 1 Statistical Mechanics of Fermions 1.1 Partition Function In the case of fermions, we had learnt that the field operator satisfies the anticommutation
More informationMany-Body Fermion Density Matrix: Operator-Based Truncation Scheme
Many-Body Fermion Density Matrix: Operator-Based Truncation Scheme SIEW-ANN CHEONG and C. L. HENLEY, LASSP, Cornell U March 25, 2004 Support: NSF grants DMR-9981744, DMR-0079992 The Big Picture GOAL Ground
More information6 Multi-particle Systems (2)
6 Multi-particle Systems (2) Exercise 6.1: Four particles in one dimension Four one-dimensional particles of mass m are confined to a length L with periodic boundary conditions. The Hamiltonian Ĥ for one
More informationIdentical Particles. Bosons and Fermions
Identical Particles Bosons and Fermions In Quantum Mechanics there is no difference between particles and fields. The objects which we refer to as fields in classical physics (electromagnetic field, field
More information221B Lecture Notes Quantum Field Theory II (Fermi Systems)
1B Lecture Notes Quantum Field Theory II (Fermi Systems) 1 Statistical Mechanics of Fermions 1.1 Partition Function In the case of fermions, we had learnt that the field operator satisfies the anticommutation
More informationSECOND QUANTIZATION. Lecture notes with course Quantum Theory
SECOND QUANTIZATION Lecture notes with course Quantum Theory Dr. P.J.H. Denteneer Fall 2008 2 SECOND QUANTIZATION 1. Introduction and history 3 2. The N-boson system 4 3. The many-boson system 5 4. Identical
More informationarxiv: v1 [cond-mat.other] 23 Feb 2008
Reduced density matrices and coherence of trapped interacting bosons Kaspar Sakmann, Alexej I. Streltsov, Ofir E. Alon, and Lorenz S. Cederbaum Theoretische Chemie, Universität Heidelberg, D-69120 Heidelberg,
More informationBrief review of Quantum Mechanics (QM)
Brief review of Quantum Mechanics (QM) Note: This is a collection of several formulae and facts that we will use throughout the course. It is by no means a complete discussion of QM, nor will I attempt
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 informationLecture notes for QFT I (662)
Preprint typeset in JHEP style - PAPER VERSION Lecture notes for QFT I (66) Martin Kruczenski Department of Physics, Purdue University, 55 Northwestern Avenue, W. Lafayette, IN 47907-036. E-mail: markru@purdue.edu
More informationMany-Body Problems and Quantum Field Theory
Philippe A. Martin Francois Rothen Many-Body Problems and Quantum Field Theory An Introduction Translated by Steven Goldfarb, Andrew Jordan and Samuel Leach Second Edition With 102 Figures, 7 Tables and
More informationSection 10: Many Particle Quantum Mechanics Solutions
Physics 143a: Quantum Mechanics I Section 10: Many Particle Quantum Mechanics Solutions Spring 015, Harvard Here is a summary of the most important points from this week (with a few of my own tidbits),
More informationThe nature of superfluidity in the cold atomic unitary Fermi gas
The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas
More informationStatic and Dynamic Properties of One-Dimensional Few-Atom Systems
Image: Peter Engels group at WSU Static and Dynamic Properties of One-Dimensional Few-Atom Systems Doerte Blume Ebrahim Gharashi, Qingze Guan, Xiangyu Yin, Yangqian Yan Department of Physics and Astronomy,
More informationFrom unitary dynamics to statistical mechanics in isolated quantum systems
From unitary dynamics to statistical mechanics in isolated quantum systems Marcos Rigol Department of Physics The Pennsylvania State University The Tony and Pat Houghton Conference on Non-Equilibrium Statistical
More informationGeneral Exam Part II, Fall 1998 Quantum Mechanics Solutions
General Exam Part II, Fall 1998 Quantum Mechanics Solutions Leo C. Stein Problem 1 Consider a particle of charge q and mass m confined to the x-y plane and subject to a harmonic oscillator potential V
More informationFirst Problem Set for Physics 847 (Statistical Physics II)
First Problem Set for Physics 847 (Statistical Physics II) Important dates: Feb 0 0:30am-:8pm midterm exam, Mar 6 9:30am-:8am final exam Due date: Tuesday, Jan 3. Review 0 points Let us start by reviewing
More informationFermionic Algebra and Fock Space
Fermionic Algebra and Fock Space Earlier in class we saw how harmonic-oscillator-like bosonic commutation relations [â α,â β ] = 0, [ ] â α,â β = 0, [ ] â α,â β = δ α,β (1) give rise to the bosonic Fock
More informationLecture 12. The harmonic oscillator
Lecture 12 The harmonic oscillator 107 108 LECTURE 12. THE HARMONIC OSCILLATOR 12.1 Introduction In this chapter, we are going to find explicitly the eigenfunctions and eigenvalues for the time-independent
More informationSuperfluid 3 He. Miguel A. Morales
Superfluid 3 He Miguel A. Morales Abstract In this report I will discuss the main properties of the superfluid phases of Helium 3. First, a brief description of the experimental observations and the phase
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 informationQuantum Mechanics II
Quantum Mechanics II Prof. Boris Altshuler March 8, 011 1 Lecture 19 1.1 Second Quantization Recall our results from simple harmonic oscillator. We know the Hamiltonian very well so no need to repeat here.
More informationarxiv: v2 [cond-mat.quant-gas] 19 Mar 2018
Noname manuscript No. (will be inserted by the editor) Effects of interaction imbalance in a strongly repulsive one-dimensional Bose gas arxiv:1802010v2 [cond-mat.quant-gas] 19 Mar 2018 R. E. Barfknecht
More informationSecond quantization (the occupation-number representation)
Second quantization (the occupation-number representation) February 14, 2013 1 Systems of identical particles 1.1 Particle statistics In physics we are often interested in systems consisting of many identical
More information5. Superconductivity. R(T) = 0 for T < T c, R(T) = R 0 +at 2 +bt 5, B = H+4πM = 0,
5. Superconductivity In this chapter we shall introduce the fundamental experimental facts about superconductors and present a summary of the derivation of the BSC theory (Bardeen Cooper and Schrieffer).
More informationSymmetry properties, density profiles and momentum distribution of multicomponent mixtures of strongly interacting 1D Fermi gases
Symmetry properties, density profiles and momentum distribution of multicomponent mixtures of strongly interacting 1D Fermi gases Anna Minguzzi LPMMC Université Grenoble Alpes and CNRS Multicomponent 1D
More informationQuantum Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
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 informationMagnets, 1D quantum system, and quantum Phase transitions
134 Phys620.nb 10 Magnets, 1D quantum system, and quantum Phase transitions In 1D, fermions can be mapped into bosons, and vice versa. 10.1. magnetization and frustrated magnets (in any dimensions) Consider
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationQuantum Mechanics I Physics 5701
Quantum Mechanics I Physics 5701 Z. E. Meziani 02/24//2017 Physics 5701 Lecture Commutation of Observables and First Consequences of the Postulates Outline 1 Commutation Relations 2 Uncertainty Relations
More information1 Bose condensation and Phase Rigidity
1 Bose condensation and Phase Rigidity 1.1 Overview and reference literature In this part of the course will explore various ways of describing the phenomenon of Bose-Einstein Condensation (BEC), and also
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 informationNANOSCALE SCIENCE & TECHNOLOGY
. NANOSCALE SCIENCE & TECHNOLOGY V Two-Level Quantum Systems (Qubits) Lecture notes 5 5. Qubit description Quantum bit (qubit) is an elementary unit of a quantum computer. Similar to classical computers,
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationarxiv: v2 [cond-mat.quant-gas] 12 Jan 2018
Effective three-body interactions for bosons in a double-well confinement Jacek Dobrzyniecki,, Xikun Li, Anne E. B. Nielsen,, and Tomasz Sowiński Institute of Physics, Polish Academy of Sciences, Aleja
More informationLecture notes on topological insulators
Lecture notes on topological insulators Ming-Che Chang Department of Physics, National Taiwan Normal University, Taipei, Taiwan Dated: May 8, 07 I. D p-wave SUPERCONDUCTOR Here we study p-wave SC in D
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 September 18, 2014 2 Chapter 5 Atoms in optical lattices Optical lattices
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 informationSmall Trapped s-wave Interacting Fermi Gases: How to Quantify Correlations?
Image: Peter Engels group at WSU Small Trapped s-wave Interacting Fermi Gases: How to Quantify Correlations? Doerte Blume and Kevin M. Daily Dept. of Physics and Astronomy, Washington State University,
More informationBreakdown and restoration of integrability in the Lieb-Liniger model
Breakdown and restoration of integrability in the Lieb-Liniger model Giuseppe Menegoz March 16, 2012 Giuseppe Menegoz () Breakdown and restoration of integrability in the Lieb-Liniger model 1 / 16 Outline
More informationSECOND QUANTIZATION. notes by Luca G. Molinari. (oct revised oct 2016)
SECOND QUANTIZATION notes by Luca G. Molinari (oct 2001- revised oct 2016) The appropriate formalism for the quantum description of identical particles is second quantisation. There are various equivalent
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 informationFundamentals and New Frontiers of Bose Einstein Condensation
Contents Preface v 1. Fundamentals of Bose Einstein Condensation 1 1.1 Indistinguishability of Identical Particles.......... 1 1.2 Ideal Bose Gas in a Uniform System............ 3 1.3 Off-Diagonal Long-Range
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 informationQUALIFYING EXAMINATION, Part 2. Solutions. Problem 1: Quantum Mechanics I
QUALIFYING EXAMINATION, Part Solutions Problem 1: Quantum Mechanics I (a) We may decompose the Hamiltonian into two parts: H = H 1 + H, ( ) where H j = 1 m p j + 1 mω x j = ω a j a j + 1/ with eigenenergies
More informationChapter Microscopic Theory
TT1-Chap4-1 Chapter 4 4. Microscopic Theory 4.1 Attractive Electron-Electron Interaction 4.1.1 Phonon Mediated Interaction 4.1.2 Cooper Pairs 4.1.3 Symmetry of Pair Wavefunction 4.2 BCS Groundstate 4.2.1
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 informationarxiv: v1 [cond-mat.quant-gas] 29 Mar 2016
arxiv:1603.08794v1 [cond-mat.quant-gas] 29 Mar 2016 One-dimensional multicomponent Fermi gas in a trap: quantum Monte Carlo study N. Matveeva Universite Grenoble Alpes, CNRS, LPMMC, UMR 5493, F-38042 Grenoble,
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 informationψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions.
1. Quantum Mechanics (Fall 2004) Two spin-half particles are in a state with total spin zero. Let ˆn a and ˆn b be unit vectors in two arbitrary directions. Calculate the expectation value of the product
More informationarxiv: v4 [cond-mat.quant-gas] 31 Aug 2016
epl draft arxiv:1511.136v4 [cond-mat.quant-gas] 31 Aug 16 Mariusz Gajda, Jan Mostowski, Tomasz Sowiński, and Magdalena Za luska-kotur Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 3/46,
More informationProblem Set No. 1: Quantization of Non-Relativistic Fermi Systems Due Date: September 14, Second Quantization of an Elastic Solid
Physics 56, Fall Semester 5 Professor Eduardo Fradkin Problem Set No. : Quantization of Non-Relativistic Fermi Systems Due Date: September 4, 5 Second Quantization of an Elastic Solid Consider a three-dimensional
More informationThe general solution of Schrödinger equation in three dimensions (if V does not depend on time) are solutions of time-independent Schrödinger equation
Lecture 17 Page 1 Lecture 17 L17.P1 Review Schrödinger equation The general solution of Schrödinger equation in three dimensions (if V does not depend on time) is where functions are solutions of time-independent
More information1.1 Creation and Annihilation Operators in Quantum Mechanics
Chapter Second Quantization. Creation and Annihilation Operators in Quantum Mechanics We will begin with a quick review of creation and annihilation operators in the non-relativistic linear harmonic oscillator.
More informationLS coupling. 2 2 n + H s o + H h f + H B. (1) 2m
LS coupling 1 The big picture We start from the Hamiltonian of an atomic system: H = [ ] 2 2 n Ze2 1 + 1 e 2 1 + H s o + H h f + H B. (1) 2m n e 4πɛ 0 r n 2 4πɛ 0 r nm n,m Here n runs pver the electrons,
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationIntermission: Let s review the essentials of the Helium Atom
PHYS3022 Applied Quantum Mechanics Problem Set 4 Due Date: 6 March 2018 (Tuesday) T+2 = 8 March 2018 All problem sets should be handed in not later than 5pm on the due date. Drop your assignments in the
More information129 Lecture Notes More on Dirac Equation
19 Lecture Notes More on Dirac Equation 1 Ultra-relativistic Limit We have solved the Diraction in the Lecture Notes on Relativistic Quantum Mechanics, and saw that the upper lower two components are large
More informationTwo and Three-Dimensional Systems
0 Two and Three-Dimensional Systems Separation of variables; degeneracy theorem; group of invariance of the two-dimensional isotropic oscillator. 0. Consider the Hamiltonian of a two-dimensional anisotropic
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 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 informationFermionic Algebra and Fock Space
Fermionic Algebra and Fock Space Earlier in class we saw how the harmonic-oscillator-like bosonic commutation relations [â α,â β ] = 0, [ ] â α,â β = 0, [ ] â α,â β = δ α,β (1) give rise to the bosonic
More informationarxiv: v1 [cond-mat.mes-hall] 2 Oct 2012
Quantum rings for beginners II: Bosons versus fermions M. Manninen a, S. Viefers b,, and S.M. Reimann c a Nanoscience Center, Department of Physics, FIN-40014 University of Jyväskylä, Finland b Department
More informationLECTURES ON STATISTICAL MECHANICS E. MANOUSAKIS
LECTURES ON STATISTICAL MECHANICS E. MANOUSAKIS February 18, 2011 2 Contents 1 Need for Statistical Mechanical Description 9 2 Classical Statistical Mechanics 13 2.1 Phase Space..............................
More informationThe Klein-Gordon equation
Lecture 8 The Klein-Gordon equation WS2010/11: Introduction to Nuclear and Particle Physics The bosons in field theory Bosons with spin 0 scalar (or pseudo-scalar) meson fields canonical field quantization
More informationIntroduction to Quantum Mechanics PVK - Solutions. Nicolas Lanzetti
Introduction to Quantum Mechanics PVK - Solutions Nicolas Lanzetti lnicolas@student.ethz.ch 1 Contents 1 The Wave Function and the Schrödinger Equation 3 1.1 Quick Checks......................................
More informationLecture 11: Long-wavelength expansion in the Neel state Energetic terms
Lecture 11: Long-wavelength expansion in the Neel state Energetic terms In the last class we derived the low energy effective Hamiltonian for a Mott insulator. This derivation is an example of the kind
More informationarxiv:quant-ph/ v1 15 Dec 2004
Entanglement in the XX Spin Chain with Energy Current V. Eisler, and Z. Zimborás 2, Institute for Theoretical Physics, Eötvös University, 7 Budapest, Pázmány sétány /a, Hungary 2 Research Institute for
More information4.3 Lecture 18: Quantum Mechanics
CHAPTER 4. QUANTUM SYSTEMS 73 4.3 Lecture 18: Quantum Mechanics 4.3.1 Basics Now that we have mathematical tools of linear algebra we are ready to develop a framework of quantum mechanics. The framework
More informationarxiv: v2 [cond-mat.quant-gas] 17 Jun 2014
How to experimentally detect a GGE? - Universal Spectroscopic Signatures of the GGE in the Tons gas Garry Goldstein and Natan Andrei Department of Physics, Rutgers University and Piscataway, New Jersey
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 informationQuantum Physics II (8.05) Fall 2002 Assignment 12 and Study Aid
Quantum Physics II (8.05) Fall 2002 Assignment 12 and Study Aid Announcement This handout includes 9 problems. The first 5 are the problem set due. The last 4 cover material from the final few lectures
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 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 informationP3317 HW from Lecture and Recitation 7
P3317 HW from Lecture 1+13 and Recitation 7 Due Oct 16, 018 Problem 1. Separation of variables Suppose we have two masses that can move in 1D. They are attached by a spring, yielding a Hamiltonian where
More informationIdentical Particles in Quantum Mechanics
Identical Particles in Quantum Mechanics Chapter 20 P. J. Grandinetti Chem. 4300 Nov 17, 2017 P. J. Grandinetti (Chem. 4300) Identical Particles in Quantum Mechanics Nov 17, 2017 1 / 20 Wolfgang Pauli
More informationParis Samuel Miles-Brenden. June 23, 2017
Symmetry of Covariance & Exchange: Particle, Field Theory & The Two Body Equation Paris Samuel Miles-Brenden June 23, 2017 Introduction The interior and effective representational algebra of the raising
More informationQuantum ideal gases: bosons
Quantum ideal gases: bosons Any particle with integer spin is a boson. In this notes, we will discuss the main features of the statistics of N non-interacting bosons of spin S (S =,,...). We will only
More informationarxiv: v2 [quant-ph] 12 Aug 2008
Complexity of thermal states in quantum spin chains arxiv:85.449v [quant-ph] Aug 8 Marko Žnidarič, Tomaž Prosen and Iztok Pižorn Department of physics, FMF, University of Ljubljana, Jadranska 9, SI- Ljubljana,
More informationtunneling theory of few interacting atoms in a trap
tunneling theory of few interacting atoms in a trap Massimo Rontani CNR-NANO Research Center S3, Modena, Italy www.nano.cnr.it Pino D Amico, Andrea Secchi, Elisa Molinari G. Maruccio, M. Janson, C. Meyer,
More informationEntanglement in Many-Body Fermion Systems
Entanglement in Many-Body Fermion Systems Michelle Storms 1, 2 1 Department of Physics, University of California Davis, CA 95616, USA 2 Department of Physics and Astronomy, Ohio Wesleyan University, Delaware,
More informationQuantum Physics 2006/07
Quantum Physics 6/7 Lecture 7: More on the Dirac Equation In the last lecture we showed that the Dirac equation for a free particle i h t ψr, t = i hc α + β mc ψr, t has plane wave solutions ψr, t = exp
More informationMajorana Fermions in Superconducting Chains
16 th December 2015 Majorana Fermions in Superconducting Chains Matilda Peruzzo Fermions (I) Quantum many-body theory: Fermions Bosons Fermions (II) Properties Pauli exclusion principle Fermions (II)
More information