Stabilizing Canonical- Ensemble Calculations in the Auxiliary-Field Monte Carlo Method
|
|
- Lorena Gaines
- 6 years ago
- Views:
Transcription
1 Stabilizing Canonical- Ensemble Calculations in the Auxiliary-Field Monte Carlo Method C. N. Gilbreth In collaboration with Y. Alhassid Yale University Advances in quantum Monte Carlo techniques for non-relativistic many-body systems Institute for Nuclear Theory July 5, 2013
2 Outline Auxiliary-field quantum Monte Carlo (AFMC) in the configuration-interaction framework Exact particle-number projection Numerical stabilization More efficient method for particle-number projection Application: Thermodynamic properties of a trapped finite cold atomic Fermi gas Heat capacity, condensate fraction, pairing gap Spatial density Density due to an extra particle 2
3 Auxiliary-field Quantum Monte Carlo (AFMC) A method for studying highly-correlated systems which is free of systematic errors Benefits: Permits finite-temperature calculations No fixed-node approximations (for good sign interactions) Useful in different contexts (electronic structure, nuclear physics, atomic physics, chemistry) Allows calculation of any one- or two-body observable Challenges: Sign problem -- for repulsive interactions and certain projections Scaling is Ns 3 N t or Ns 4 N t, depending on the application Numerical stability at low temperatures / large model spaces 3
4 AFMC in Configuration (Fock) Space Goal: compute the thermal expectation hôi = Tr(Ôe Ĥ ) Tr(e Ĥ ) for a collection of fermions. Steps to formulate AFMC: (a) Formulation of the Configuration-Interaction (CI) Hamiltonian as a sum of quadratic one-body operators (b) Trotter decomposition (imaginary time discretization) (c) Path-integral representation (d) Calculation of observables for a given set of fields (e) Monte Carlo evaluation 4
5 (a) Hamiltonian and (b)trotter decomposition Single-particle basis: ii,i=1,...,n s Model space: set of all Slater determinants from the s.p. basis Hamiltonian: Ĥ = X i,j Rewrite in terms of one-body densities a i a k, a j a l and diagonalize the matrix Ĥ = Ĥ h ij a i a j X Ô2 X i,j,k,l v ijkl v ijkl a i a j a la k to obtain where and are one-body operators. Ĥ 0 Ô Trotter decomposition: and e Ĥ = e Ĥ 0 Y e Ĥ =(e Ĥ ) N t, /N t e Ô2 /2 + O(( ) 2 ) 5
6 Discretization error in observables condensate fraction T/T F = N max =11 fit ( β ( >= apple 1/32) β The error due to the Trotter decomposition scales as We perform a linear fit and extrapolate to! 0 6
7 { { { (c) Path-Integral Formulation Hubbard-Stratonovich (HS) transformation: e Ô 2 /2 = Z 1 This linearizes the two-body part of the Hamiltonian into a one-body operator, with the addition of an auxiliary field. End result of the HS transformation: 1 d e 2 /2 e s Ô, s = ( 1, < 0 i, > 0 e Ĥ = Z Time-dependent auxiliary fields (many) D[ ]G Û( ) Integration measure Non-interacting propagator Gaussian weight A path integral of a non-interacting propagator with respect to fluctuating time-dependent auxiliary fields. 7
8 (d) Observables in the Canonical Ensemble Want to compute: Z D[ ]G Tr[Û( )Ô] Û( ) is a non-interacting propagator and Ô is a one-or two- body operator. To compute traces for fixed particle number, we use exact particle-number projection: Number of single-particle states Tr N [ÔÛ( )] = 1 XN s N s m=1 e i' mn Tr GC [ÔÛ( )ei' m ˆN ]. Canonical trace N = number of particles Can be derived by writing ' m 2 m/n s Tr GC Grand-canonical trace as a sum of canonical traces [W. E. Ormand, et al., Phys. Rev. C 49, 1422 (1994)] 8
9 (d) Canonical ensemble cont. To compute grand-canonical quantities, can use matrix algebra in the single particle space (of order ~ ): Tr GC (Û( )ei' m ˆN )=det(1+ue i' m ) Tr GC [a i ajû( )ei' m ˆN ] Tr GC [Û( )ei' m ˆN ] = 1 1+U 1 e i' m ji where U is the matrix representing space. (e) Monte Carlo evaluation Û in the single-particle Z Discretize the integrals over the auxiliary fields: D[ ]G Û( )! X w( )Û( ) Apply the Metropolis algorithm. [Koonin, et al., Phys. Rep. 278, 1 (1997)] 9
10 Accuracy of AFMC: Three cold atoms E E diagonalization β β Ground-state energy AFMC AFMC vs. diagonalization for a fixed interaction with three cold atoms (""#) in a harmonic trap 10
11 Numerical Stabilization
12 Numerical Stabilization Need to compute, for each m, Tr GC [Û( )ei' m ˆN ]=det(1+ue i' m ) where U = matrix representing Û in the single-particle space. At long imaginary times, the propagator U becomes unstable. (1) Compute U: (a) Unstabilized method (b) Stabilized method (2) Compute determinant: (i) Unstabilized method (ii) Standard stabilized method (iii) New stabilized method (faster) 12
13 Unstabilized calculation of U As N t grows large, the matrix product U = U( Nt ) U( 2 )U( 1 ) becomes ill-conditioned: ratio of largest to smallest eigenvalues becomes very large. The matrix elements of U, however, all become large, and it is impossible to extract smaller scales (those near the Fermi surface) accurately. E Energy vs. inverse temperature for 20 cold atoms in the unitary limit. At low temperatures, calculation becomes unreliable β Unstabilized calculation is only good for high temperatures (small ) 13
14 Stabilized Calculation of U Stabilized calculation of U : Use method of [E. Y. Loh Jr. and J. E. Gubernatis, in Electronic Phase Transitions, 1992] Instead of computing U, compute a decomposition U = ADB of U which explicitly displays all numerical scales. A and B are well-conditioned, and D is diagonal with positive entries: U = ADB = x x x x x x x x x E.g., singular value decomposition or QR decomposition X X X x x x x x x x x x 14
15 To compute Stabilized Matrix Multiplication decomposition of each new factor:, begin with an initial and carefully update when multiplying in U = U Nt U 2 U 1 U 1 U 1 = A 1 D 1 B 1. U n+1 (A n D n B n )=((U n+1 A n )D n )B n =(A n+1 D n+1 B 0 )B n = A n+1 D n+1 (B 0 B n ) (decompose) (group terms) (decompose) (group terms) = A n+1 D n+1 B n+1. U = A Nt D Nt B Nt The intermediate matrix (U n+1 A n )D n therefore can be decomposed stably is column-stratified and 15
16 Stabilized Matrix Multiplication cont. A column-stratified matrix displays its scales in the columns: U n+1 A n D n = x x x x x x x x x x x x x x x x x x X X X = X X X X X X X X X Similarly, a row-stratified matrix displays its scales in the rows: DBA = X X X 1 0 x x x x x x x x x 1 0 x x x x x x x x x Stratified matrices are much more stable for decomposition and diagonalization than general matrices with widely varying scales. 1 A = X X X X X X X X X 1 A 16
17 (2) Calculation of the Determinant Wish to compute det(1 + Ue i' m ) (i) Unstabilized calculation: for each m = 1, 2,..., N s Compute U = U Nt U 1 Diagonalize U and compute Only for high temperatures. Scales as using standard matrix multiplication (diagonalization) (ii) Standard stabilized calculation: [Alhassid et al., 101, (2008)] Cannot multiply out ADB to diagonalize, as this would destroy information. Instead: 1+ADBe i' m = A(A 1 B 1 + De i' m )B Requires decomposition for each as O(Ns 4 ) det(1 + Ue i' m )= O(N 3 s ) YN s k=1 (1 + k e i' m ) = A(A 0 D 0 B 0 )B (decompose) m =1,...,N s, so scales 17
18 Stabilized Matrix Diagonalization [Gilbreth & Alhassid, arxiv: ] (iii) New method: Transform ADBx = x Let x = Ay, DBAy = y The matrix DBA is row-stratified and can be multiplied out & diagonalized stably Eigenvalues are identical to those of ADB. Reduces time back to O(Ns 3 ) Test 1: Compute eigenvectors and eigenvalues of is a random complex matrix, and n =1,2,3,... (e C ) n, where C Test 2: Compare to standard stabilization method 18
19 Accuracy of Diagonalizing an ill-conditioned Matrix max relative error unstabilized eigenvalues unstabilized eigenvectors stabilized eigenvalues stabilized eigenvectors condition number (e C ) n Error in stabilized and unstabilized eigenvalues and eigenvectors of the matrix (e C ) n of 19
20 Numerical Stabilization -- Accuracy and Timing E unstabilized stabilized time per sample (s) 1.5x x x x10 0 New method Old Method β Number of single-particle states Accuracy Timing 20
21 Summary AFMC method in the configuration-interaction framework Strongly-interacting systems with arbitrary goodsign interactions can be studied Calculations in the canonical ensemble A new stabilization method in the canonical ensemble Much faster, scales as O(N 3 s ) O(N 4 s ) instead of 21
22 Thermodynamics of a finite trapped cold atomic Fermi gas
23 Advances with Atomic Fermi Gases 1999: First realization of an ideal degenerate Fermi gas ( 40 K). Evaporative cooling to T = 0.5 TF (JILA, CO) 2002: Feshbach resonance allows tuning to stronglyinteracting regime ( 6 Li, Duke Univ.) 2005: Observation of vortex lattice after stirring confirms superfluidity (MIT) 2012: Measurement of lambda peak in heat capacity (MIT) [Zwierlein, M. W. et al., Nature 435, 1047 (2005)] [Ku, et al., Science 335, 563 (2012)] B C V /Nk B C T/TF
24 Clean, strongly interacting, tunable systems Two species ( " and # ) of Fermions interact at very short range in a harmonic trap. Zero-range (s-wave) interactions s-wave scattering length a controllable via external magnetic field Unitary limit saturates (a!1) s-wave scattering cross section. Strongly interacting, nonperturbative system when a is large [Zwierlein, M. W. et al., Nature 435, 1047 (2005)] 24
25 Unitary, Finite, Trapped Fermi Gas Goal: Describe thermodynamics of trapped finite-size (~20 particle) Fermi gas in the unitary limit. Heat capacity Pairing gap Condensate fraction. Particle density Questions: Is the superfluid phase transition visible in a system of this size? If so, is there a pseudogap effect? [Q. Chen et al., Physics Reports 412 (2005) 1-88] 25
26 AFMC for Trapped Cold Atoms Hamiltonian: H = X N 1 +N 2 i=1 Single-particle basis of harmonic-oscillator states: Interaction: 1 2 r2 i m!2 r 2 i V (r i r j )=V 0 (r i r j ) Shell-model decomposition: Our calculations are done in the canonical ensemble. + X i<j V (r i r j ) ii = nlmi = R nl (r)y m l (, ), 2n + l apple N max H = XN s i,j=1 h i,j a i a j X K,, Cutoff parameter (V0 renormalized for each N max to reproduce two-particle ground-state energy) K, ( ) KX M=0 ( ˆQ K M ( ) 2 + ˆP K M ( ) 2 ) ˆQ, ˆP one-body operators 26
27 Monte Carlo Sign of the Contact Interaction It is well-known that the contact interaction has good Monte Carlo sign in coordinate space. V (r) =V 0 (r) It is not as clear in the Configuration-Interaction formalism. The proof proceeds as follows: The Hamiltonian takes the form H = nlmσ ε nl a nlmσ a nlmσ (ab ˆV cd)a aσ a bσ a dσ a cσ, and in a time-reversed density decomposition : Ĥ = Ĥ (a b ˆV c d)( ) m b ( ) m d a aσa cσ ā bσ ā dσ, One can show that the interaction matrix ā a, = time reverse V ac,bd =(a b ˆV c d)( ) m a ( ) m b has only nonpositive eigenvalues when V 0 < 0. 27
28 Sign of the Contact Interaction (cont.) We show this as follows: ac,bd v ac(a b ˆV c d)( ) m b+m d v bd This implies e = V 0 = V 0 = V 0 Ô 2 /2 = s =1 ac,bd Z 1 d 3 r d 3 r v acϕ a(r)ϕ c (r)ϕ b(r)( ) m b ϕ d(r)( ) m d v bd ( ac ac v acϕ a(r)ϕ c (r) )( d 3 r vacϕ a(r)ϕ c (r) 2 0. in the Hubbard-Stratonovich transformation: 1 d e 2 /2 e Ô, bd ϕ b (r)ϕ d(r)v bd ) So the auxiliary-field Hamiltonian is complex-conjugation invariant. 28
29 Sign of the Contact Interaction (cont.) The propagator factorizes: Û = Û"Û# Û ", Û# are invariant under complex conjugation: K 0 Û ",# K 0 = Û",# Û ",# vi = vi K 0 Û ",# vi = Û",#K 0 vi = K 0 vi ) Therefore, complex eigenvalues come in complex-conjugate pairs So Tr N",# (Û",#) = X i (N ",# ) i = real When N " = N #, Tr N",N # (Û) =Tr N " (Û")Tr N# (Û#) =[Tr N" (Û")] 2 0 so the calculation has good sign. 29
30 Accuracy of Renormalized Contact Interaction Exact 3-particle energy at unitarity 12 E 10 8 Fermi gas Contact interaction, Nmax = T/T F Tuning V0 to reproduce two-particle ground-state energy provides an accurate interaction for the three-particle system. [Exact results based on S. Tan and Werner & Castin] 30
31 C = de(t ) dt Heat capacity When computing the heat capacity, the statistical error in the energy can be greatly magnified. We use a method of correlated errors introduced in [S. Liu and Y. Alhassid, PRL 87, (2001)] to avoid this. Sample fields at a single temperature Compute energies at T + T and T T for each sample: E(T ± T )= Statistical error: Var(C) = [E(T + T ) E(T T )]/2 T R D[ ]G (T ± T )Tr[ Û(,T ± T )] R D[ ]G (T ± T )Tr Û(,T ± T ) T 1 h (2 T ) 2 Var(hĤ(T + T )i) + Var(hĤ(T T )i) i 2 Cov(hĤ(T + T )i, hĥ(t T )i) Covariance greatly reduces statistical error in heat capacity 31
32 15 10 C 5 54 Fe 54 Fe T (MeV) T (MeV) [S. Liu and Y. Alhassid, PRL 87, (2001)] 32
33 Condensate Fraction Condensate fraction does not have a standard meaning in finite-size systems We define a condensate fraction from the two-body density matrix (TBDM). C. N. Yang (1962): Off-diagonal long-range order (ODLRO) is equivalent to the existence of a large eigenvalue in the TBDM Calculation: 1. Compute C(ij, kl) ha i," a j,# a l,"a k,# i 2. Diagonalize to obtain a pair wavefunction corresponding to the largest eigenvalue of C. B X ij ' ij a i" a j# 3. The largest eigenvalue satisfies (for N " = N # = N/2) and = hb Bi 0 apple apple N/2 so n /(N/2) defines a condensate fraction. 33
34 Energy-Staggering Pairing gap Energy-staggering pairing gap where gap [E(N ",N # + 1) + E(N " +1,N # ) E(N " +1,N # + 1) E(N ",N # )]/2 E(N ",N # )= energy of a system with spin-up particles and spin-down particles. N # gap systems and a systems with an unpaired particle. measures the difference in energy between fully paired Calculation: (for N " = N # = N/2) 1. Sample auxiliary fields for the system. 2. Compute E(N ",N # ), E(N " +1,N # )=E(N ",N # + 1), and E(N " +1,N # + 1) for each sample 3. Compute variance of using correlated errors. gap N " (N ",N # + 1) 34
35 Signatures of the Phase Transition: 20 Atoms [Gilbreth and Alhassid, arxiv: ] gap /ε F condensate fraction C Fermi gas Fermi gas T/T F gap =[E(N ",N # + 1) + E(N " +1,N # ) E(N ",N # ) E(N " +1,N # + 1)]/2 C L (ab, cd) =ha LM (ab)a LM(cd)i n = max /(N/2) 35
36 Particle Density & Odd Particle Effect ρ(r) (units of a -3 osc ) T/T F = T/T F = T/T F = T/T F = ρ (units of a -3 osc ) T/T F = r/a osc r/a osc (r) = P (r) (r) (r) =[ (r) (9,10) + (r) (10,9) ]/2 [ (r) (10,10) + (r) (9,9) ]/2 Cloud spreads out at higher temperatures Unpaired particle prefers the edge of the trap 36
37 Summary: AFMC in the CI Framework for Cold Atoms We devised a more efficient numerical stabilization method in the canonical ensemble which allows calculations in much larger model spaces. First ab initio calculations of heat capacity and energy-staggering pairing gap across the superfluid phase transition in any system of cold atoms Condensate fraction and particle density for a finite trapped system of cold atoms Clear signatures of a superfluid phase transition The addition of an extra particle to the spin-balanced system produces extra density at the edge of the trap. 37
38 Questions Questions: Does gap show a pseudogap effect in the trapped system for larger numbers of particles at unitarity? Does gap (as opposed to the spectral function) show a pseudogap effect in the uniform system? How can we extend AFMC in the CI framework to larger model spaces? - Represent the interaction in coordinate space? - Optimize matrix exponential, diagonalization, etc. methods? - Optimize single-particle basis functions? - Alternatives to Metropolis? 38
39 Acknowledgements -- Thank You! Yoram Alhassid Abhishek Mukherjee Cem Özen Yale HPC Department of Physics Konstantin Nesterov Marco Bonet-Matiz 39
The 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 informationAuxiliary-field Monte Carlo methods in Fock space: sign problems and methods to circumvent them
Auxiliary-field Monte Carlo methods in Fock space: sign problems and methods to circumvent them Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo methods in Fock
More informationAuxiliary-field quantum Monte Carlo methods for nuclei and cold atoms
Introduction Auxiliary-field quantum Monte Carlo methods for nuclei and cold atoms Yoram Alhassid (Yale University) Auxiliary-field Monte Carlo (AFMC) methods at finite temperature Sign problem and good-sign
More informationStatistical properties of nuclei by the shell model Monte Carlo method
Statistical properties of nuclei by the shell model Monte Carlo method Introduction Yoram Alhassid (Yale University) Shell model Monte Carlo (SMMC) method Circumventing the odd particle-number sign problem
More informationAuxiliary-field quantum Monte Carlo methods in heavy nuclei
Mika Mustonen and Yoram Alhassid (Yale University) Introduction Auxiliary-field quantum Monte Carlo methods in heavy nuclei Auxiliary-field Monte Carlo (AFMC) methods at finite temperature Sign problem
More informationThe shell model Monte Carlo approach to level densities: recent developments and perspectives
The shell model Monte Carlo approach to level densities: recent developments and perspectives Yoram Alhassid (Yale University) Introduction: the shell model Monte Carlo (SMMC) approach Level density in
More informationShell model Monte Carlo level density calculations in the rare-earth region
Shell model Monte Carlo level density calculations in the rare-earth region Kadir Has University Istanbul, Turkey Workshop on Gamma Strength and Level Density in Nuclear Physics and Nuclear Technology
More informationarxiv: v1 [nucl-th] 18 Jan 2018
Nuclear deformation in the configuration-interaction shell model arxiv:181.6175v1 [nucl-th] 18 Jan 218 Y. Alhassid, 1 G.F. Bertsch 2,3 C.N. Gilbreth, 2 and M.T. Mustonen 1 1 Center for Theoretical Physics,
More informationGround-state properties, excitations, and response of the 2D Fermi gas
Ground-state properties, excitations, and response of the 2D Fermi gas Introduction: 2D FG and a condensed matter perspective Auxiliary-field quantum Monte Carlo calculations - exact* here Results on spin-balanced
More informationIs a system of fermions in the crossover BCS-BEC. BEC regime a new type of superfluid?
Is a system of fermions in the crossover BCS-BEC BEC regime a new type of superfluid? Finite temperature properties of a Fermi gas in the unitary regime Aurel Bulgac,, Joaquin E. Drut, Piotr Magierski
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 informationFermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations
Fermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations Aurel Bulgac,, Joaquin E. Drut and Piotr Magierski University of Washington, Seattle, WA
More informationF. Chevy Seattle May 2011
THERMODYNAMICS OF ULTRACOLD GASES F. Chevy Seattle May 2011 ENS FERMION GROUPS Li S. Nascimbène Li/K N. Navon L. Tarruell K. Magalhaes FC C. Salomon S. Chaudhuri A. Ridinger T. Salez D. Wilkowski U. Eismann
More informationThermodynamics of the polarized unitary Fermi gas from complex Langevin. Joaquín E. Drut University of North Carolina at Chapel Hill
Thermodynamics of the polarized unitary Fermi gas from complex Langevin Joaquín E. Drut University of North Carolina at Chapel Hill INT, July 2018 Acknowledgements Organizers Group at UNC-CH (esp. Andrew
More informationShell Model Monte Carlo Methods and their Applications to Nuclear Level Densities
Shell Model Monte Carlo Methods and their Applications to Nuclear Level Densities H. Nakada (Chiba U., Japan) @ RIKEN (Feb. 14, 2013) Contents : I. Introduction II. Shell model Monte Carlo methods III.
More informationThe shape distribution of nuclear level densities in the shell model Monte Carlo method
The shape distribution of nuclear level densities in the shell model Monte Carlo method Introduction Yoram Alhassid (Yale University) Shell model Monte Carlo (SMMC) method and level densities Nuclear deformation
More informationPath Integral (Auxiliary Field) Monte Carlo approach to ultracold atomic gases. Piotr Magierski Warsaw University of Technology
Path Integral (Auxiliary Field) Monte Carlo approach to ultracold atomic gases Piotr Magierski Warsaw University of Technology Collaborators: A. Bulgac - University of Washington J.E. Drut - University
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 informationTackling the Sign Problem of Ultracold Fermi Gases with Mass-Imbalance
Tackling the Sign Problem of Ultracold Fermi Gases with Mass-Imbalance Dietrich Roscher [D. Roscher, J. Braun, J.-W. Chen, J.E. Drut arxiv:1306.0798] Advances in quantum Monte Carlo techniques for non-relativistic
More informationParticle-number projection in finite-temperature mean-field approximations to level densities
Particle-number projection in finite-temperature mean-field approximations to level densities Paul Fanto (Yale University) Motivation Finite-temperature mean-field theory for level densities Particle-number
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 informationIntroduction to Cold Atoms and Bose-Einstein Condensation. Randy Hulet
Introduction to Cold Atoms and Bose-Einstein Condensation Randy Hulet Outline Introduction to methods and concepts of cold atom physics Interactions Feshbach resonances Quantum Gases Quantum regime nλ
More informationNuclear structure III: Nuclear and neutron matter. National Nuclear Physics Summer School Massachusetts Institute of Technology (MIT) July 18-29, 2016
Nuclear structure III: Nuclear and neutron matter Stefano Gandolfi Los Alamos National Laboratory (LANL) National Nuclear Physics Summer School Massachusetts Institute of Technology (MIT) July 18-29, 2016
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 informationLow- and High-Energy Excitations in the Unitary Fermi Gas
Low- and High-Energy Excitations in the Unitary Fermi Gas Introduction / Motivation Homogeneous Gas Momentum Distribution Quasi-Particle Spectrum Low Energy Excitations and Static Structure Function Inhomogeneous
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 informationThermodynamics, pairing properties of a unitary Fermi gas
Thermodynamics, pairing properties of a unitary Fermi gas Piotr Magierski (Warsaw University of Technology/ University of Washington, Seattle) Collaborators: Aurel Bulgac (Seattle) Joaquin E. Drut (LANL)
More informationSpecific heat of a fermionic atomic cloud in the bulk and in traps
Specific heat of a fermionic atomic cloud in the bulk and in traps Aurel Bulgac,, Joaquin E. Drut, Piotr Magierski University of Washington, Seattle, WA Also in Warsaw Outline Some general remarks Path
More informationEquilibrium and nonequilibrium properties of unitary Fermi gas from Quantum Monte Carlo
Equilibrium and nonequilibrium properties of unitary ermi gas from Quantum Monte Carlo Piotr Magierski Warsaw University of Technology Collaborators: A. Bulgac - University of Washington J.E. Drut - University
More informationQuantum Quantum Optics Optics VII, VII, Zakopane Zakopane, 11 June 09, 11
Quantum Optics VII, Zakopane, 11 June 09 Strongly interacting Fermi gases Rudolf Grimm Center for Quantum Optics in Innsbruck University of Innsbruck Austrian Academy of Sciences ultracold fermions: species
More informationBCS-BEC BEC Crossover at Finite Temperature in Cold Gases and Condensed Matter KITP
BCS-BEC BEC Crossover at Finite Temperature in Cold Gases and Condensed Matter KITP May 2007 Cold Atom Collaborators: Qijin Chen J. Stajic (U Chicago; LANL) Yan He (U. Chicago) ChihChun Chien (U. Chicago)
More informationIntroduction to Bose-Einstein condensation 4. STRONGLY INTERACTING ATOMIC FERMI GASES
1 INTERNATIONAL SCHOOL OF PHYSICS "ENRICO FERMI" Varenna, July 1st - July 11 th 2008 " QUANTUM COHERENCE IN SOLID STATE SYSTEMS " Introduction to Bose-Einstein condensation 4. STRONGLY INTERACTING ATOMIC
More informationCentral density. Consider nuclear charge density. Frois & Papanicolas, Ann. Rev. Nucl. Part. Sci. 37, 133 (1987) QMPT 540
Central density Consider nuclear charge density Frois & Papanicolas, Ann. Rev. Nucl. Part. Sci. 37, 133 (1987) Central density (A/Z* charge density) about the same for nuclei heavier than 16 O, corresponding
More informationNuclear Structure and Reactions using Lattice Effective Field Theory
Nuclear Structure and Reactions using Lattice Effective Field Theory Dean Lee North Carolina State University Nuclear Lattice EFT Collaboration Frontiers of Nuclear Physics Kavli Institute for Theoretical
More informationThermodynamic Measurements in a Strongly Interacting Fermi Gas
J Low Temp Phys (2009) 154: 1 29 DOI 10.1007/s10909-008-9850-2 Thermodynamic Measurements in a Strongly Interacting Fermi Gas Le Luo J.E. Thomas Received: 25 July 2008 / Accepted: 12 October 2008 / Published
More informationRecent results in lattice EFT for nuclei
Recent results in lattice EFT for nuclei Dean Lee (NC State) Nuclear Lattice EFT Collaboration Centro de Ciencias de Benasque Pedro Pascua Bound states and resonances in EFT and Lattice QCD calculations
More informationThe phases of matter familiar for us from everyday life are: solid, liquid, gas and plasma (e.f. flames of fire). There are, however, many other
1 The phases of matter familiar for us from everyday life are: solid, liquid, gas and plasma (e.f. flames of fire). There are, however, many other phases of matter that have been experimentally observed,
More 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 informationnuclear level densities from exactly solvable pairing models
nuclear level densities from exactly solvable pairing models Stefan Rombouts In collaboration with: Lode Pollet, Kris Van Houcke, Dimitri Van Neck, Kris Heyde, Jorge Dukelsky, Gerardo Ortiz Ghent University
More informationStrongly paired fermions
Strongly paired fermions Alexandros Gezerlis TALENT/INT Course on Nuclear forces and their impact on structure, reactions and astrophysics July 4, 2013 Strongly paired fermions Neutron matter & cold atoms
More information2D Bose and Non-Fermi Liquid Metals
2D Bose and Non-Fermi Liquid Metals MPA Fisher, with O. Motrunich, D. Sheng, E. Gull, S. Trebst, A. Feiguin KITP Cold Atoms Workshop 10/5/2010 Interest: A class of exotic gapless 2D Many-Body States a)
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
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 informationHigh-order Chin actions in path integral Monte Carlo
High-order Chin actions in path integral Monte Carlo High-order actions and their applications, Barcelona 2009 Jordi Boronat Departament de Física i Enginyeria Nuclear Universitat Politècnica de Catalunya
More informationPROGRESS IN UNDERSTANDING THE PROPERTIED OF MANY-BODY SYSTEMS BY QUANTUM MONTE CARLO SIMULATIONS
PROGRESS IN UNDERSTANDING THE PROPERTIED OF MANY-BODY SYSTEMS BY QUANTUM MONTE CARLO SIMULATIONS Francesco Pederiva! Physics Department - University of Trento INFN - TIFPA, Trento Institute for Fundamental
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 informationThe Cooper Problem. Problem : A pair of electrons with an attractive interaction on top of an inert Fermi sea c c FS,
Jorge Duelsy Brief History Cooper pair and BCS Theory (1956-57) Richardson exact solution (1963). Gaudin magnet (1976). Proof of Integrability. CRS (1997). Recovery of the exact solution in applications
More informationThermodynamics of the unitary Fermi gas
Ludwig-Maximilians-University, Faculty of Physics, Chair of Theoretical Solid State Physics, Theresienstr. 37, 80333 Munich, Germany E-mail: O.Goulko@physik.uni-muenchen.de M. Wingate DAMTP, University
More informationStrongly Interacting Fermi Gases: Universal Thermodynamics, Spin Transport, and Dimensional Crossover
NewSpin, College Station, 1/16/011 Strongly Interacting ermi Gases: Universal hermodynamics, Spin ransport, and Dimensional Crossover Martin Zwierlein Massachusetts Institute of echnology Center for Ultracold
More informationStrongly correlated Cooper pair insulators and superfluids
Strongly correlated Cooper pair insulators and superfluids Predrag Nikolić George Mason University Acknowledgments Collaborators Subir Sachdev Eun-Gook Moon Anton Burkov Arun Paramekanti Affiliations and
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 informationEffective Field Theory and. the Nuclear Many-Body Problem
Effective Field Theory and the Nuclear Many-Body Problem Thomas Schaefer North Carolina State University 1 Nuclear Effective Field Theory Low Energy Nucleons: Nucleons are point particles Interactions
More information(Biased) Theory Overview: Few-Body Physics
Image: Peter Engels group at WSU (Biased) Theory Overview: Few-Body Physics Doerte Blume Dept. of Physics and Astronomy, Washington State University, Pullman. With graduate students Kevin M. Daily and
More informationMomentum-space and Hybrid Real- Momentum Space DMRG applied to the Hubbard Model
Momentum-space and Hybrid Real- Momentum Space DMRG applied to the Hubbard Model Örs Legeza Reinhard M. Noack Collaborators Georg Ehlers Jeno Sólyom Gergely Barcza Steven R. White Collaborators Georg Ehlers
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 information4 Matrix product states
Physics 3b Lecture 5 Caltech, 05//7 4 Matrix product states Matrix product state (MPS) is a highly useful tool in the study of interacting quantum systems in one dimension, both analytically and numerically.
More informationCompensating strong coupling with large charge
Compensating strong coupling with large charge Orestis Loukas Institute for Theoretical Physics and Albert Einstein Center for Fundamental Physics, Bern based on arxiv: 1610.04495, 1612.08985 together
More informationSpin liquid phases in strongly correlated lattice models
Spin liquid phases in strongly correlated lattice models Sandro Sorella Wenjun Hu, F. Becca SISSA, IOM DEMOCRITOS, Trieste Seiji Yunoki, Y. Otsuka Riken, Kobe, Japan (K-computer) Williamsburg, 14 June
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 informationThe Fermion Bag Approach
The Fermion Bag Approach Anyi Li Duke University In collaboration with Shailesh Chandrasekharan 1 Motivation Monte Carlo simulation Sign problem Fermion sign problem Solutions to the sign problem Fermion
More informationA Mixture of Bose and Fermi Superfluids. C. Salomon
A Mixture of Bose and Fermi Superfluids C. Salomon Enrico Fermi School Quantum Matter at Ultralow Temperatures Varenna, July 8, 2014 The ENS Fermi Gas Team F. Chevy, Y. Castin, F. Werner, C.S. Lithium
More 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 informationarxiv: 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 informationQuantum gases in the unitary limit and...
Quantum gases in the unitary limit and... Andre LeClair Cornell university Benasque July 2 2010 Outline The unitary limit of quantum gases S-matrix based approach to thermodynamics Application to the unitary
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationStrongly Correlated Physics With Ultra-Cold Atoms
Strongly Correlated Physics With Ultra-Cold Atoms Predrag Nikolić Rice University Acknowledgments Collaborators Subir Sachdev Eun-Gook Moon Anton Burkov Arun Paramekanti Sponsors W.M.Keck Program in Quantum
More informationEffective Field Theory for light nuclear systems
Effective Field Theory for light nuclear systems Jimmy Rotureau Chalmers University of Technology, Göteborg, Sweden B. Barrett, University of Arizona, Tucson I. Stetcu, University of Washington, Seattle
More informationR. Citro. In collaboration with: A. Minguzzi (LPMMC, Grenoble, France) E. Orignac (ENS, Lyon, France), X. Deng & L. Santos (MP, Hannover, Germany)
Phase Diagram of interacting Bose gases in one-dimensional disordered optical lattices R. Citro In collaboration with: A. Minguzzi (LPMMC, Grenoble, France) E. Orignac (ENS, Lyon, France), X. Deng & L.
More informationThe 2010 US National Nuclear Physics Summer School and the TRIUMF Summer Institute, NNPSS-TSI June 21 July 02, 2010, Vancouver, BC, Canada
TU DARMSTADT The 2010 US National Nuclear Physics Summer School and the TRIUMF Summer Institute, NNPSS-TSI June 21 July 02, 2010, Vancouver, BC, Canada Achim Richter ECT* Trento/Italy and TU Darmstadt/Germany
More information2 The Density Operator
In this chapter we introduce the density operator, which provides an alternative way to describe the state of a quantum mechanical system. So far we have only dealt with situations where the state of a
More information3 Symmetry Protected Topological Phase
Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and
More informationPairing properties, pseudogap phase and dynamics of vortices in a unitary Fermi gas
Pairing properties, pseudogap phase and dynamics of vortices in a unitary Fermi gas Piotr Magierski (Warsaw University of Technology/ University of Washington, Seattle) Collaborators: Aurel Bulgac (Seattle)
More informationReview of lattice EFT methods and connections to lattice QCD
Review of lattice EFT methods and connections to lattice QCD Dean Lee Michigan State University uclear Lattice EFT Collaboration Multi-Hadron Systems from Lattice QCD Institute for uclear Theory Feburary
More informationQuantum quenches in the thermodynamic limit
Quantum quenches in the thermodynamic limit Marcos Rigol Department of Physics The Pennsylvania State University Correlations, criticality, and coherence in quantum systems Evora, Portugal October 9, 204
More informationMatrix Representation
Matrix Representation Matrix Rep. Same basics as introduced already. Convenient method of working with vectors. Superposition Complete set of vectors can be used to express any other vector. Complete set
More informationMany-body physics 2: Homework 8
Last update: 215.1.31 Many-body physics 2: Homework 8 1. (1 pts) Ideal quantum gases (a)foranidealquantumgas,showthatthegrandpartitionfunctionz G = Tre β(ĥ µ ˆN) is given by { [ ] 1 Z G = i=1 for bosons,
More informationNuclear physics: a laboratory for many-particle quantum mechanics or From model to theory in nuclear structure physics
Nuclear physics: a laboratory for many-particle quantum mechanics or From model to theory in nuclear structure physics G.F. Bertsch University of Washington Stockholm University and the Royal Institute
More informationCold atoms and AdS/CFT
Cold atoms and AdS/CFT D. T. Son Institute for Nuclear Theory, University of Washington Cold atoms and AdS/CFT p.1/27 History/motivation BCS/BEC crossover Unitarity regime Schrödinger symmetry Plan Geometric
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 informationProbing the Optical Conductivity of Harmonically-confined Quantum Gases!
Probing the Optical Conductivity of Harmonically-confined Quantum Gases! Eugene Zaremba Queen s University, Kingston, Ontario, Canada Financial support from NSERC Work done in collaboration with Ed Taylor
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 informationThe Nuclear Many-Body Problem. Lecture 2
The Nuclear Many-Body Problem Lecture 2 How do we describe nuclei? Shell structure in nuclei and the phenomenological shell model approach to nuclear structure. Ab-initio approach to nuclear structure.
More informationApplications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices
Applications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices Vahid Dehdari and Clayton V. Deutsch Geostatistical modeling involves many variables and many locations.
More informationSecond quantization. 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 informationMultipartite entanglement in fermionic systems via a geometric
Multipartite entanglement in fermionic systems via a geometric measure Department of Physics University of Pune Pune - 411007 International Workshop on Quantum Information HRI Allahabad February 2012 In
More informationDimensional BCS-BEC crossover
Dimensional BCS-BEC crossover Igor Boettcher Institute for Theoretical Physics, Heidelberg University ERG 2014 Outline Ultracold atoms BCS-BEC Crossover 2D Experiments with... Theory Experiment C. Wetterich
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 informationChapter 2 Quantum chemistry using auxiliary field Monte Carlo
Chapter 2 Quantum chemistry using auxiliary field Monte Carlo 1. The Hubbard-Stratonovich Transformation 2. Neuhauser s shifted contour 3. Calculation of forces and PESs 4. Multireference AFMC 5. Examples
More informationDiagrammatic Green s Functions Approach to the Bose-Hubbard Model
Diagrammatic Green s Functions Approach to the Bose-Hubbard Model Matthias Ohliger Institut für Theoretische Physik Freie Universität Berlin 22nd of January 2008 Content OVERVIEW CONSIDERED SYSTEM BASIC
More informationApproximation Methods in QM
Chapter 3 Approximation Methods in QM Contents 3.1 Time independent PT (nondegenerate)............... 5 3. Degenerate perturbation theory (PT)................. 59 3.3 Time dependent PT and Fermi s golden
More informationThe interacting boson model
The interacting boson model P. Van Isacker, GANIL, France Dynamical symmetries of the IBM Neutrons, protons and F-spin (IBM-2) T=0 and T=1 bosons: IBM-3 and IBM-4 The interacting boson model Nuclear collective
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 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 informationA Mixture of Bose and Fermi Superfluids. C. Salomon
A Mixture of Bose and Fermi Superfluids C. Salomon INT workshop Frontiers in quantum simulation with cold atoms University of Washington, April 2, 2015 The ENS Fermi Gas Team F. Chevy, Y. Castin, F. Werner,
More informationComputational Challenges in Nuclear and Many-Body Physics. Summary
Computational Challenges in Nuclear and Many-Body Physics Summary Yoram Alhassid (Yale University) Physical systems of interest Challenges in many-body physics Computational methods and their formulation
More informationDense Matter and Neutrinos. J. Carlson - LANL
Dense Matter and Neutrinos J. Carlson - LANL Neutron Stars and QCD phase diagram Nuclear Interactions Quantum Monte Carlo Low-Density Equation of State High-Density Equation of State Neutron Star Matter
More informationSuper Efimov effect. Sergej Moroz University of Washington. together with Yusuke Nishida and Dam Thanh Son. Tuesday, April 1, 14
Super Efimov effect together with Yusuke Nishida and Dam Thanh Son Sergej Moroz University of Washington Few-body problems They are challenging but useful: Newton gravity Quantum atoms Quantum molecules
More informationUniversal Behavior of Small Two-Component Fermi Gases with Equal Masses
Image: Peter Engels group at WSU Universal Behavior of Small Two-Component Fermi Gases with Equal Masses Doerte Blume Dept. of Physics and Astronomy, Washington State University, Pullman Collaborators:
More informationFinding eigenvalues for matrices acting on subspaces
Finding eigenvalues for matrices acting on subspaces Jakeniah Christiansen Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty advisor: Prof Todd Kapitula Department
More informationEquation of state of the unitary Fermi gas
Equation of state of the unitary Fermi gas Igor Boettcher Institute for Theoretical Physics, University of Heidelberg with S. Diehl, J. M. Pawlowski, and C. Wetterich C o ld atom s Δ13, 11. 1. 2013 tio
More information