Stabilizing Canonical- Ensemble Calculations in the Auxiliary-Field Monte Carlo Method

Size: px
Start display at page:

Download "Stabilizing Canonical- Ensemble Calculations in the Auxiliary-Field Monte Carlo Method"

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 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 information

Auxiliary-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 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 information

Auxiliary-field quantum Monte Carlo methods for nuclei and cold atoms

Auxiliary-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 information

Statistical properties of nuclei by the shell model Monte Carlo method

Statistical 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 information

Auxiliary-field quantum Monte Carlo methods in heavy nuclei

Auxiliary-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 information

The 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 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 information

Shell model Monte Carlo level density calculations in the rare-earth region

Shell 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 information

arxiv: v1 [nucl-th] 18 Jan 2018

arxiv: 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 information

Ground-state properties, excitations, and response of the 2D Fermi gas

Ground-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 information

Is 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? 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 information

Small Trapped s-wave Interacting Fermi Gases: How to Quantify Correlations?

Small 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 information

Fermions 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 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 information

F. Chevy Seattle May 2011

F. 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 information

Thermodynamics 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 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 information

Shell Model Monte Carlo Methods and their Applications to Nuclear Level Densities

Shell 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 information

The 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 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 information

Path 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 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 information

Reference for most of this talk:

Reference 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 information

Tackling the Sign Problem of Ultracold Fermi Gases with Mass-Imbalance

Tackling 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 information

Particle-number projection in finite-temperature mean-field approximations to level densities

Particle-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 information

v(r i r j ) = h(r i )+ 1 N

v(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 information

Introduction to Cold Atoms and Bose-Einstein Condensation. Randy Hulet

Introduction 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 information

Nuclear 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. 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 information

From 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. 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 information

Low- and High-Energy Excitations in the Unitary Fermi Gas

Low- 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 information

Ultracold Fermi and Bose Gases and Spinless Bose Charged Sound Particles

Ultracold 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 information

Thermodynamics, pairing properties of a unitary Fermi gas

Thermodynamics, 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 information

Specific 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 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 information

Equilibrium and nonequilibrium properties of unitary Fermi gas from Quantum Monte Carlo

Equilibrium 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 information

Quantum Quantum Optics Optics VII, VII, Zakopane Zakopane, 11 June 09, 11

Quantum 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 information

BCS-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 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 information

Introduction to Bose-Einstein condensation 4. STRONGLY INTERACTING ATOMIC FERMI GASES

Introduction 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 information

Central 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) 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 information

Nuclear Structure and Reactions using Lattice Effective Field Theory

Nuclear 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 information

Thermodynamic Measurements in a Strongly Interacting Fermi Gas

Thermodynamic 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 information

Recent results in lattice EFT for nuclei

Recent 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 information

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

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 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 information

Condensate fraction for a polarized three-dimensional Fermi gas

Condensate 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 information

nuclear level densities from exactly solvable pairing models

nuclear 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 information

Strongly paired fermions

Strongly 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 information

2D Bose and Non-Fermi Liquid Metals

2D 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 information

Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory

Electronic 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 information

BCS-BEC Crossover. Hauptseminar: Physik der kalten Gase Robin Wanke

BCS-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 information

High-order Chin actions in path integral Monte Carlo

High-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 information

PROGRESS 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 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 information

ICAP Summer School, Paris, Three lectures on quantum gases. Wolfgang Ketterle, MIT

ICAP 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 information

The Cooper Problem. Problem : A pair of electrons with an attractive interaction on top of an inert Fermi sea c c FS,

The 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 information

Thermodynamics of the unitary Fermi gas

Thermodynamics 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 information

Strongly Interacting Fermi Gases: Universal Thermodynamics, Spin Transport, and Dimensional Crossover

Strongly 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 information

Strongly correlated Cooper pair insulators and superfluids

Strongly 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 information

From laser cooling to BEC First experiments of superfluid hydrodynamics

From 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 information

Effective Field Theory and. the Nuclear Many-Body Problem

Effective 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

(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 information

Momentum-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 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 information

MP463 QUANTUM MECHANICS

MP463 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 information

4 Matrix product states

4 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 information

Compensating strong coupling with large charge

Compensating 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 information

Spin liquid phases in strongly correlated lattice models

Spin 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 information

Introduction to Path Integral Monte Carlo. Part I.

Introduction 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 information

The Fermion Bag Approach

The 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 information

A Mixture of Bose and Fermi Superfluids. C. Salomon

A 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 information

Many-Body Problems and Quantum Field Theory

Many-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 information

arxiv: v1 [cond-mat.quant-gas] 18 Sep 2015

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 information

Quantum gases in the unitary limit and...

Quantum 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 information

Linear Algebra: Matrix Eigenvalue Problems

Linear 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 information

Strongly Correlated Physics With Ultra-Cold Atoms

Strongly 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 information

Effective Field Theory for light nuclear systems

Effective 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 information

R. Citro. In collaboration with: A. Minguzzi (LPMMC, Grenoble, France) E. Orignac (ENS, Lyon, France), X. Deng & L. Santos (MP, Hannover, Germany)

R. 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 information

The 2010 US National Nuclear Physics Summer School and the TRIUMF Summer Institute, NNPSS-TSI June 21 July 02, 2010, Vancouver, BC, Canada

The 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 information

2 The Density Operator

2 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 information

3 Symmetry Protected Topological Phase

3 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 information

Pairing 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 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 information

Review of lattice EFT methods and connections to lattice QCD

Review 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 information

Quantum quenches in the thermodynamic limit

Quantum 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 information

Matrix Representation

Matrix 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 information

Many-body physics 2: Homework 8

Many-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 information

Nuclear 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 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 information

Cold atoms and AdS/CFT

Cold 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 information

Squeezing and superposing many-body states of Bose gases in confining potentials

Squeezing 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

Probing the Optical Conductivity of Harmonically-confined Quantum Gases!

Probing 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 information

Thermalisation 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. 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 information

The Nuclear Many-Body Problem. Lecture 2

The 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 information

Applications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices

Applications 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 information

Second quantization. Emmanuel Fromager

Second 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 information

Multipartite entanglement in fermionic systems via a geometric

Multipartite 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 information

Dimensional BCS-BEC crossover

Dimensional 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 information

SECOND QUANTIZATION. Lecture notes with course Quantum Theory

SECOND 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 information

Chapter 2 Quantum chemistry using auxiliary field Monte Carlo

Chapter 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 information

Diagrammatic Green s Functions Approach to the Bose-Hubbard Model

Diagrammatic 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 information

Approximation Methods in QM

Approximation 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 information

The interacting boson model

The 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 information

1 Fluctuations of the number of particles in a Bose-Einstein condensate

1 Fluctuations of the number of particles in a Bose-Einstein condensate Exam of Quantum Fluids M1 ICFP 217-218 Alice Sinatra and Alexander Evrard The exam consists of two independant exercises. The duration is 3 hours. 1 Fluctuations of the number of particles in a Bose-Einstein

More information

Low-dimensional Bose gases Part 1: BEC and interactions

Low-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 information

A Mixture of Bose and Fermi Superfluids. C. Salomon

A 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 information

Computational Challenges in Nuclear and Many-Body Physics. Summary

Computational 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 information

Dense Matter and Neutrinos. J. Carlson - LANL

Dense 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 information

Super Efimov effect. Sergej Moroz University of Washington. together with Yusuke Nishida and Dam Thanh Son. Tuesday, April 1, 14

Super 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 information

Universal Behavior of Small Two-Component Fermi Gases with Equal Masses

Universal 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 information

Finding eigenvalues for matrices acting on subspaces

Finding 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 information

Equation of state of the unitary Fermi gas

Equation 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