Introduction to Path Integral Monte Carlo. Part III.
|
|
- Doreen Sutton
- 5 years ago
- Views:
Transcription
1 Introduction to Path Integral Monte Carlo. Part III. Alexey Filinov, Michael Bonitz Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D Kiel, Germany November 19, 2008
2 Outline 1 Calculation of superfluid fraction: finite and macroscopic systems 2 Winding number: ergodicity. Worm algorithm. Grand Canonical ensemble 3 Fermion sign problem 4 Comparison of Fermi/Bose statistics: superfluidity 5 Numerical issues of PIMC Improved high-temperature action Effective pair potential Thermodynamic averages 6 Summary
3 Outline 1 Calculation of superfluid fraction: finite and macroscopic systems 2 Winding number: ergodicity. Worm algorithm. Grand Canonical ensemble 3 Fermion sign problem 4 Comparison of Fermi/Bose statistics: superfluidity 5 Numerical issues of PIMC Improved high-temperature action Effective pair potential Thermodynamic averages 6 Summary
4 Examples: Superfluidity Superfluidity: loss of viscosity of interacting bosons below critical temperature. Discovered in liquid 4 He (P.L.Kapitza, 1938). Rotating bucket experiment (Andronikashvili): T > T BKT : Spontaneous creation of vortices by thermal excitation. Vanishing of superfluid density. Uniform 2D system: Superfluid - normal fluid phase transition at critical temperature T BKT (Berezinskii, Kosterlitz, Thouless): k B T BKT = ρ s π 2 2m 2
5 Computation of superfluid fraction γ sf Two-fluid model (Landau): Only normal fluid component of a liquid responds to slow rotation of the container walls.
6 Computation of superfluid fraction γ sf Quantum mechanical moment of inertia I qm deviates from classical expectation value I class non-classical rotational inertia (NCRI) γ sf = 1 Iqm I class, I qm = d ˆL z dω, I class = Hamiltonian in the rotating frame: NX m i ri 2 i=1 Ĥ ω = Ĥ0 ωˆl z, ˆL z = Tr[ˆL z e βĥω ], ˆLz = NX (x i p iy y i p ix ) i=1
7 Computation of superfluid fraction γ sf Quantum mechanical moment of inertia I qm deviates from classical expectation value I class non-classical rotational inertia (NCRI) γ sf = 1 Iqm I class, I qm = d ˆL z dω, I class = Hamiltonian in the rotating frame: NX m i ri 2 i=1 Ĥ ω = Ĥ0 ωˆl z, ˆL z = Tr[ˆL z e βĥω ], ˆLz = NX (x i p iy y i p ix ) i=1 Derivative of the exponential operator de βĥ ω dω = MX k=1 e (k 1)τĤω de τĥ ω dω e (M k)τĥω = Z β 0 dte tĥω dĥ ω dω dte (β t)ĥω
8 Computation of superfluid fraction γ sf Quantum mechanical moment of inertia I qm deviates from classical expectation value I class non-classical rotational inertia (NCRI) γ sf = 1 Iqm I class, I qm = d ˆL z dω, I class = Hamiltonian in the rotating frame: Ĥ ω = Ĥ 0 ωˆl z, ˆL z = Tr[ˆL z e βĥ ω ], ˆL z = Derivative of the exponential operator NX m i ri 2 i=1 NX (x i p iy y i p ix ) i=1 de βĥ ω dω = MX k=1 e (k 1)τĤω de τĥ ω d ˆL z dω ω 0 = dω * Zβ 0 e (M k)τĥω = Z β dtˆl z e tĥ 0 ˆL z e (β t)ĥ 0 + = I qm 0 tĥω dĥω dte dω dte (β t)ĥω
9 Computation of superfluid fraction γ sf Quantum mechanical moment of inertia I qm deviates from classical expectation value I class non-classical rotational inertia (NCRI) γ sf = 1 Iqm, I qm = d ˆL z NX I class dω, I class = m i ri 2 Hamiltonian in the rotating frame: Ĥ ω = Ĥ0 ωˆl z, ˆL z = Tr[ˆL z e βĥω ], ˆLz = Derivative of the exponential operator de βĥ ω dω = MX k=1 e (k 1)τĤω de τĥ ω dω e (M k)τĥω = Z β 0 i=1 NX (x i p iy y i p ix ) Superfluid fraction (in linear response ω 0) * Zβ + γ sf = ρs ρ = 1 1 dtˆl z e tĥ 0 ˆL z e (β t)ĥ 0 I class 0 i=1 tĥω dĥω dte dω dte (β t)ĥω
10 Computation of superfluid fraction γ sf d L z dω = τ Tr " ˆL 2 ze MτĤ 0 + # MX ˆL ze (k 1)τĤ 0 ˆLze (M (k 1))τĤ 0 k=2 Consider any term in the second sum. Angular momentum operates only on the kinetic energy part of the action and commutes with the internal potential energy. In the coordinate representation we obtain Z NX «drdr 1... dr M 1 ( i ) x i y i R e τĥ 0 R 1... = y i x i i=1 Z NX «2π drdr 1... dr M 1 ( i ) [r 1i r i ] z R e τĥ 0 R 1... i=1 Area of the path segment NX A 1z = [r 1i r i ] z, A kz = i=1 λ 2 τ NX [r (k+1)i r ki ] z Hence, the second sum can be written as.. = Tr[..ˆρ] = R dr... R e τĥ 0 R 1... «2 ( i ) 2 2π MX «2 2π D E! A λ 2 1z Akz A 2 τ λ 2 1z τ i=1
11 Computation of superfluid fraction γ sf Similar for the first term (operator ˆL z operates twice on one link) we obtain «2 «ˆL 2 z R e τĥ0 R 1 = ( i ) 2 2π A 2 1z R e τĥ0 R 1 +( i ) 2 2π X N (x λ 2 τ λ 2 i x 1i +y 1i y i ) τ i=1 Now we Combine the first and the second term I qm = I c I q where I c is the classical part of the responce and the quantum part NX I c = m (x i x 1i + y 1i y i ) i=1 I q = m2 2 τ A1zAz, Az = M X k=1 A kz In the end do symmetrization: A 1zA z = 1 M MP A kz A z = 1 P M A M kz A z = A 2 z. k=1 k=1
12 Superfluid fraction: finite systems Final result. Superfluid fraction in a finite system: particles are placed in a external field (e.g. a rotating cylinder around ω) γ sf = 4m2 A 2 z 2 βi class, A z = A ω ω [P.Sindzingre, M.Klein, D.Ceperley, Phys.Rev. Lett. 63, 1601 (1981)] Area enclosed by the paths A = 1 2 NX M 1 X r ki r (k+1)i i=1 k=
13 Superfluid fraction: finite systems Final result. Superfluid fraction in a finite system: particles are placed in a external field (e.g. a rotating cylinder around ω) γ sf = 4m2 A 2 z 2 βi class, A z = A ω ω [P.Sindzingre, M.Klein, D.Ceperley, Phys.Rev. Lett. 63, 1601 (1981)] Area enclosed by the paths A = NX M 1 X r ki r (k+1)i i=1 k=0 + Superfluid fraction for N = 5 charged bosons in a two-dimensional trap. λ = (e 2 /ɛl 0)/( ω) is the coupling constant.
14 Superfluid fraction: macroscopic systems Instead of a filled cylinder (with N particles) we consider two cylinders with the radius R and spacing d, with d R. Such a torus is topologically equivalent to the usual periodic boundary conditions: I class = mnr 2, A z = N round πr 2 = W z R/2 R d W z = 2πR N round Winding number W : total length of the paths along the torus W = Superfluid fraction: NX [r i (β) r i (0)] i=1 z γ sf = 4m2 A 2 z m W 2 2 βi class 2 β N Simulation box with periodic boundary conditions
15 Outline 1 Calculation of superfluid fraction: finite and macroscopic systems 2 Winding number: ergodicity. Worm algorithm. Grand Canonical ensemble 3 Fermion sign problem 4 Comparison of Fermi/Bose statistics: superfluidity 5 Numerical issues of PIMC Improved high-temperature action Effective pair potential Thermodynamic averages 6 Summary
16 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. We start from identity permutations
17 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Two particle exchanges are very probable: r λ D (T )
18 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Three particle exchanges can happen but very infrequent: 2 r > λ D (T ). Probability of longer permutations is exponentially suppressed! We are trapped in 2-3-particle permutation sector Zero winding numbers and zero superfluidity!
19 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Idea: expand the configuration space to the offdiagonal sector G and sample offdiagonal density matrix ρ(r, ˆPR ; β) Worm algorithm [Prokof ev, Svistunov and Tupitsyn (1997); N.Prokof ev, B.Svistunov, Boninsegni, Phys.Rev.Lett 96, (2006)]
20 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Continue to work in the offdiagonal sector G to sample two-particle exchanges
21 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Enjoy stable high acceptance rate for any permutation length!
22 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Enjoy stable high acceptance rate for any permutation length!
23 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Occasionally close the trajectory to return back in the diagonal sector Z. Measure thermodynamic observables related to the diagonal density matrix ρ(r, ˆPR; β) or partition function Z.
24 Winding number: ergodicity Changing a winding number requires a global update in the permutation space. Conclusion: Accurate statistics on W requires permutation lengths N 1/d. With typical few-particle update sampling (diagonal sector) we stay in one permutation sector non-ergodic sampling of W. More advance algorithms are required: Ising model cluster algorithm; PIMC worm algorithm (expanded conf. space Z and G).
25 Worm algorithm Key features: All updates in open-path-configurations are performed exclusively through the end-points of the disconnected paths Local sampling with high acceptance rates! No global updates with exponentially low acceptance! Topological classes are sampled efficiently. No critical slowing down in most cases. Open paths are related to important physics, i.e correlation functions, and are not merely an algorithm trick. Usual PIMC in Canonical ensemble can now be easily generalized to Grand Canonical ensemble New tool for corresponding experimental systems.
26 Why Grand Canonical ensemble is better Z GCE = X e N=0 βr 0 µ N(τ) dτ ZCE (N, V, β) Advantages of the simulations in grand canonical ensemble (GCE): Off-diagonal one-particle d. matrix: D E n(r, r, t t) = Ψ(r, t ) Ψ (r, t) Condensate fraction: n(r, r ) =N 0φ 0 (r)φ 0(r )+ + X i 0 n i φ i (r)φ i (r ) Canonical (CE) Grand Canonical (GCE) Ψ(r, t ) Ψ (r, t) n(r, r ) r r N 0
27 Why Grand Canonical ensemble is better Z GCE = X e N=0 βr 0 µ N(τ) dτ ZCE (N, V, β) Advantages of the simulations in grand canonical ensemble (GCE): Off-diagonal one-particle d. matrix: D E n(r, r, t t) = Ψ(r, t ) Ψ (r, t) Condensate fraction: n(r, r ) =N 0φ 0 (r)φ 0(r )+ + X i 0 n i φ i (r)φ i (r ) Prokof ev, Svistunov, Boninsegni, PRE 74, (2007) n(r, r ) r r N 0
28 Why Grand Canonical ensemble is better Z GCE = X e N=0 βr 0 µ N(τ) dτ ZCE (N, V, β) Advantages of the simulations in grand canonical ensemble (GCE): µ is an input parameter and N µ is a simple diagonal property. Compressibility, kvt = (N N ) 2 µ. N µ Average particle number N 10 Average particle number for chemical potential (T=0.5 K) λ Chemical potential µ in K µ-dependence of the partice number N in 2D He 4 clusters in a parabolic trap ω 1/λ 2. µ
29 Why Grand Canonical ensemble is better Z GCE = X e N=0 βr 0 µ N(τ) dτ ZCE (N, V, β) Advantages of the simulations in grand canonical ensemble (GCE): e 2 R p(x)dx R 1 R 2 Solve ergodicity issue for disorder problems Allows for efficient sampling of exponentially rare event
30 Outline 1 Calculation of superfluid fraction: finite and macroscopic systems 2 Winding number: ergodicity. Worm algorithm. Grand Canonical ensemble 3 Fermion sign problem 4 Comparison of Fermi/Bose statistics: superfluidity 5 Numerical issues of PIMC Improved high-temperature action Effective pair potential Thermodynamic averages 6 Summary
31 Sign problem Let us consider a standard Monte Carlo problem: Integration over space of states ν. Each configuration has a weight factor: Eν /T W ν > 0, e Expectation values: P A = ν Aν Wν P Wν
32 Sign problem Let us consider a standard Monte Carlo problem: Integration over space of states ν. Each configuration has a weight factor: Eν /T W ν > 0, e Expectation values: P A = ν Aν Wν P Wν Frequently for quantum mechanical systems W ν is not-positive definite function. Then P ν Aν sign(wν) Wν A = P sign(wν) Wν ν Now we can proceed with the standart MC using W ν as a sampling probability P ν A A = ν sign(w ν ) A sign P = ν sign(w ν ) sign
33 Sign problem The trouble comes in cases: sign 0. Both A sign and sign have finite errorbars which fluctuate 1 + or A + δa = A sign + δ AS sign + δ S A sign sign δa A δ AS A sign 0 + δ S sign 0 δ AS A sign + δ S sign «
34 Sign problem The trouble comes in cases: sign 0. Both A sign and sign have finite errorbars which fluctuate 1 + or A + δa = A sign + δ AS sign + δ S A sign sign δa A δ AS A sign 0 + δ S sign 0 δ AS A sign + δ S sign There is no generic solution of the sign-problem. This prevents MC methods from studies: interacting fermions magnetic systems real time dynamics, etc. But the sign-problem can be reduced or eliminated by a proper choice of the basis set. Example: ĤΨ ν = E νψ ν, Z = X ν e βeν Ψ-eigenfunctions - all terms are positive «
35 Fermion sign problem in PIMC Metropolis algorithm gives the same distribution of permutations for both Fermi and Bose systems. The reason is that for sampling permutations we use the modulus of the off-diagonal density matrix ρ S/A (R, R; β) = 1 X (±1) P ρ(r, N! ˆPR; β) = P = 1 X Z (±1) P dr 1... dr M 1 ρ(r, R 1; β)... ρ(r M 1, N! ˆPR; β) P Bosons: all permutations contribute with the same (positive) sign Fermions: positive and negative terms (corresponding to even and odd permutations) are close in their absolute value and cancel each other. Accurate calculation of this small difference is drastically hampered with the increase of quantum degeneracy (low T, high density).
36 Fermion problem: partial solutions for Quantum MC 1 Fixed-node (fixed-phase) approximation Use restricted (reduced) area of PIMC integration which contains only even permutations. Most of the area with the cancellation of even and odd permutations are excluded using an approximate trial ansatz for the N-particle fermion density matrix. Requires knowledge of nodes of DM. References: D.M.Ceperley, Fermion Nodes, J. Stat. Phys. 63, 1237 (1991); D.M.Ceperley, Path Integral Calculations of Normal Liquid 3He, Phys. Rev. Lett. 69, 331 (1992).
37 Fermion problem: partial solutions for Quantum MC 2 Direct PIMC Do not sample individual permutations in the sum. Instead use the full expression presented in a form of an determinant. In this case the absolute value of the determinant is used in the sampling probabilities. Its value becomes close to zero in the regions of equal contributions of even and odd permutations and Monte Carlo sampling successfully avoids such regions. References: V.S.Filinov, M.Bonitz, W.Ebeling, and V.E.Fortov, Thermodynamics of hot dense H-plasmas: Path integral Monte Carlo simulations and analytical approximations, Plasma Physics and Controlled Fusion 43, 743 (2001).
38 Fermion problem: partial solutions for Quantum MC 3 Multilevel-blocking PIMC Trace the cancellations of permutations by grouping the path coordinates into blocks (levels). Use numerical integration to get good estimation of the fermion density matrix at high temperature. Further use it in the sampling probabilities of path coordinates on the next level (corresponding to the density matrix at lower temperature). Most of the sign fluctuations are already excluded at higher levels and sampling at low levels (lower temperatures) becomes more efficient. References: R.Egger, W.Hausler, C.H.Mak, and H.Grabert, Crossover from Fermi Liquid to Wigner Molecule Behavior in Quantum Dots, Phys. Rev. Lett. 82, 3320 (1999).
39 Comparison of Fermi/Bose statistics for excitons Exact treatment of excitons as composite particles consisting of two fermions: N e!n h! permutations, total density matrix is antisymmetric.
40 Comparison of Fermi/Bose statistics for excitons Exact treatment of excitons as composite particles consisting of two fermions: N e!n h! permutations, total density matrix is antisymmetric. Exact fermionic calculations clearly show the effect of Fermi repulsion on the radial distributions. This favors localization of individual particles and reduces the density in between the central particle and the shell. Although the effect is small it has a large effect on the superfluid fraction.
41 Outline 1 Calculation of superfluid fraction: finite and macroscopic systems 2 Winding number: ergodicity. Worm algorithm. Grand Canonical ensemble 3 Fermion sign problem 4 Comparison of Fermi/Bose statistics: superfluidity 5 Numerical issues of PIMC Improved high-temperature action Effective pair potential Thermodynamic averages 6 Summary
42 Comparison of Fermi/Bose statistics: superfluidity Simulation: Spin polarized electron-hole bilayer (two coupled ZnSe quantum wells). External confinement: parabolic trap. Parameters: T = 312mK, d = 20.1nm. A.Filinov, M.Bonitz, P.Ludwig, and Yu.E.Lozovik, phys. stat. sol. (c) 3, No. 7, 2457 (2006) Applicability of bosonic model: Critical density: ρ cm 2 ρ ρ : both models agree. ρ > ρ : fermionic calculations indicate drop in the superfluid fraction.
43 Outline 1 Calculation of superfluid fraction: finite and macroscopic systems 2 Winding number: ergodicity. Worm algorithm. Grand Canonical ensemble 3 Fermion sign problem 4 Comparison of Fermi/Bose statistics: superfluidity 5 Numerical issues of PIMC Improved high-temperature action Effective pair potential Thermodynamic averages 6 Summary
44 Improved high-temperature action We need to improve the simple factorization formula for unbound potentials, e.g. V (r) = 1/r 2 3 NY ρ(r, R, τ) ρ F (r i, r i, τ) exp 4 τ Y V (r jk ) 5 j<k i=1 not be normalized (due to singularity). Alternative: take into account two-body correlations exactly ρ(r, R ; τ) NY i=1 ρ F (r i, r i ; τ) Y j<k ρ [2] (r j, r k, r j, r k; τ) ρ F (r i, r i ; τ)ρ F (r k, r k ; τ) + O(ρ[3] ), Three-body and higher order terms become negligible by decreasing τ = β/m.
45 Two-body density matrix Effective pair potential U pair ρ [2] (r 1, r 2, r 1, r 2; τ) ρ F (r 1, r 1 ; τ)ρ F (r 2, r 2 ; τ) = e τu U pair (r 1 r 2, r 1 r 2; τ) is temperature-dependent and finite at r 12 = 0. Two-body density matrix ρ [2] can be obtained by solving two-particle problem. Replace singular potentials with bounded effective potentials defined as h i ρ(r 12, r 12; τ) = ρ F (r 12, r 12; β) exp τu pair (r 12, r 12; τ) pair 12
46 Advantages of effective potentials Exact treatment of pair correlations allows: drastically reduce number of factorization factors reduce dimension of integrals simplification of path integral sampling. Approaches to computation of the pair density matrix: 1 Cumulant approximation (Feynman-Kacs) [R.Feynman and A.R.Hibbs Quantum Mechanics and Path Integral] 2 Solution of two-particle Bloch equation (matrix squaring technique) [Klemm and Storer (1974), D.Ceperley] 3 Perturbation or semi-classical approximation [Kelbg, Ebeling, Deutsch, Feynman, Kleinert ( )]
47 Two-body density matrix: eigenstates and Feynman-Kacs Sum over eigenstates of the Hamiltonian ρ(r, r, τ) = X i e τe i Ψ i (r) Ψ i (r ) Can be used if all eigenstates are known analytically, e.g. for Coulomb potential (Pollock Comm. Phys. 52, 49 (1988)). Feynman-Kac formula: ρ(r, r ; τ) = 2 Z r(τ)=r Z τ dt 3 * = Dr(t) exp 4 mṙ 2 (t)/2 + V 12(r(t)) 5 = r(0)=r 0 τr e 0 + V 12 (r(t))dt The average can be calculated by Monte Carlo sampling of all Gaussian random walks from r to r. ρ F
48 Two-body density matrix: matrix squaring Matrix squaring technique: Factorization into a center-of-mass, relative coordinates, ρ(r i, r j, r i, r j; τ) = ρ cm(r, R ; τ)ρ(r, r ; τ), and expansion in partial waves: ρ 2D (r, r ; τ) = ρ 3D (r, r ; τ) = 1 2π r r + X l= ρ l (r, r ; τ)e i lθ, 1 + X (2l + 1) ρ 4πr r l (r, r ; τ) P l (cos Θ) l=0 Convolution equation: k-iterations raise the temperature by 2 k : τ/2 k τ Z 0 dr ρ l (r, r ; τ 2 m+1 ) ρ l(r, r ; τ 2 m+1 ) = ρ l(r, r ; τ ), m = k 1,..., 0 2m Semi-classical approximation: the start for the matrix-squaring iterations Z! r ρ l (r, r ; τ/2 k ) = exp τ/2k V (x)dx r r r
49 Two-body density matrix: perturbative solution First order perturbation solution for the two-particle Bloch equation: ρ ij = (m im j ) 3/2 (2π τ) 3 Solution for Coulomb interaction: τ ρ(r i, r j, r i, r j; τ) = Ĥ ρ(r i, r j, r i, r j; τ) h exp m i 2 2 τ (r i r j) 2i h exp m j 2 2 τ (r i r j) 2i exp[ τφ ij ], Φ ij (r ij, r ij; τ) e i e j Z 1 0 dα d ij (α) erf d ij (α)/λ ij 2 p α(1 α)!, where d ij (α) = αr ij + (1 α)r ij, erf(x) is the error function, R erf(x) = 2 x π 0 dte t2, and λ 2 ij = 2 τ 2µ ij with µ 1 ij = m 1 i + m 1 j. The diagonal element (r ij = r ij ) is called the Kelbg potential (DKP) 2 Φ (r ij, τ) = q iq j 41 e r 2 ij λ 2 ij + π r» «3 ij r ij 1 erf γ ij 5. r ij λ ij γ ij λ ij
50 Effective pair potential Solid lines: electron electron (e e) and electron proton (e p) potential U pair at T = 10 6 K. Dotted lines: Coulomb interaction; dashed lines: Deutsch potential; open circles: U F variational perturbative potential (Feynman and Kleinert). Distances are in a B = 2 /m ee 2 and potentials are in Ha = e 2 /a B units.
51 Effective pair potential: temperature dependence 0.0 a) U(β) 0.0 b) U(β)+β U(β)/ β K K K K K K K U pair Φ Kelbg W ω, x m Φ 0 Kelbg r/a B r/a B [Filinov,Golubnychiy,Bonitz,Ebeling,Dufty, PRE 70, (2004)] K U pair Φ Kelbg W ω, x m Φ 0 Kelbg (a): Effective electron-proton potential (in units of Ha): the DKP Φ 0 (r; β), the improved DKP Φ(r; β), variational potential W Ω,xm 1, pair potential U p corresponding to the exact density matrix. Temperatures 5 000, , and K. (b): Contribution to the potential energy: f (β) = U(β) + β U(β), Ep = Tr[f (β)ˆρ]/z. β
52 Estimators for thermodynamic averages Quantities of interest: energy, pressure (equation of state), specific heat, structure factor, pair distribution functions, condensate or superfluid fraction, etc. All quantities can be obtained by averaging with the thermal density matrix or as derivatives of the partition function Z.
53 Estimators for thermodynamic averages Quantities of interest: energy, pressure (equation of state), specific heat, structure factor, pair distribution functions, condensate or superfluid fraction, etc. All quantities can be obtained by averaging with the thermal density matrix or as derivatives of the partition function Z. 1 We usually calculate only ratios of integrals. Free energy and entropy require special techniques. Examples: Total energy E = 1 Z Z β = 1 Z dr Z Z ρ(r, R; β) = E = dmn 2β ρ(r, R; β), β dr 1... dr M 1 (e S kin S V )/λ dmn τ * * M 1 1 X β i=0 + π (R i R i+1 ) 2 + λ 2 τ Kinetic energy and pressure E kin = m βz Z m, P = F V = 1 3V * 1 M M 1 X i=0 2E kin X i<j + d dτ (τu(r i, τ)). + d(τu(r i, τ)) r ij dr ij
54 Estimators for thermodynamic averages Quantities of interest: energy, pressure (equation of state), specific heat, structure factor, pair distribution functions, condensate or superfluid fraction, etc. All quantities can be obtained by averaging with the thermal density matrix or as derivatives of the partition function Z. 1 We usually calculate only ratios of integrals. Free energy and entropy require special techniques. MP 2 The variance of some estimators can be large: A = 1 A(R M i ) ± σ A / M. Use virial estimators. Introduce temperature-dependent measure i=1
55 Estimators for thermodynamic averages Quantities of interest: energy, pressure (equation of state), specific heat, structure factor, pair distribution functions, condensate or superfluid fraction, etc. All quantities can be obtained by averaging with the thermal density matrix or as derivatives of the partition function Z. 1 We usually calculate only ratios of integrals. Free energy and entropy require special techniques. MP 2 The variance of some estimators can be large: A = 1 A(R M i ) ± σ A / M. Use virial estimators. Introduce temperature-dependent measure ix R i (τ) = R 0 + λ τ ξ m, i = 1,..., M 1. S kin = M 1 X i=0 m=1 M 1 π X (R i R i+1 ) 2 = πξi+1 2 f (β), λ 2 τ Z ρ(r, R; β) = E = dn 2β + i=0 Z dr1..dr M 1 (..) = λ dmn τ dξ1..dξ M 1 λ dmn τ * X U(Ri (τ)) + β U(R i(τ) R i (τ) i=1 λ d(m 1)N τ (..) 1 R i (τ) β +. λ dn τ
56 Estimators for thermodynamic averages Quantities of interest: energy, pressure (equation of state), specific heat, structure factor, pair distribution functions, condensate or superfluid fraction, etc. All quantities can be obtained by averaging with the thermal density matrix or as derivatives of the partition function Z. 1 We usually calculate only ratios of integrals. Free energy and entropy require special techniques. MP 2 The variance of some estimators can be large: A = 1 A(R M i ) ± σ A / M. Use virial estimators. Introduce temperature-dependent measure 3 Take care for other sources of errors: systematic (Trotter formula), statistical (autocorrelation times) and finite-size (scaling with system size). i=1
57 Outline 1 Calculation of superfluid fraction: finite and macroscopic systems 2 Winding number: ergodicity. Worm algorithm. Grand Canonical ensemble 3 Fermion sign problem 4 Comparison of Fermi/Bose statistics: superfluidity 5 Numerical issues of PIMC Improved high-temperature action Effective pair potential Thermodynamic averages 6 Summary
58 Summary PIMC allows for first-principle simulations of interacting quantum particles. Many-body correlations can be treated at any accuracy. Problems remain to be solved: Fermion sign problem Limitations to several hundred particles due to large computational demands Efficient simulations of spin and magnetic effects
59 References R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals, McGraw Hill, New York H. Kleinert, Path Integrals in Quantum Mechanics, Statistics and Polymer Physics, World Scientific, Second edition, D.M. Ceperley, Path Integrals in the Theory of Condensed Helium, Rev. Mod. Phys. 67, 279 (1995). D.M. Ceperley, Microscopic simulations in physics, Rev. Mod. Phys. 71, 438 (1999) David P. Landau and K. Binder, A Guide to Monte Carlo Simulations in Statistical Physics. Cambridge, Cambridge University Press, A. Filinov and M. Bonitz, First-principle Approaches to Correlated Quantum Systems, p in Introduction to Computational Methods in Many-Body Physics, M. Bonitz and D. Semkat (Eds.), Rinton Press, Princeton, 2006
Clas s ical and path integral Monte Carlo s imulation of charged particles in traps
Clas s ical and path integral Monte Carlo s imulation of charged particles in traps Alexei Filinov, Michael B onitz, in c o lla b o r a tio n w ith Vladimir Filinov, Patrick Ludwig, Jens B öning, Henning
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 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 informationWorm Algorithm PIMC Application to Liquid and Solid 4 He
Worm Algorithm PIMC Application to Liquid and Solid 4 He KITP - Santa Barbara, 2/8/06 Massimo Boninsegni University of Alberta Nikolay Prokof ev and Boris Svistunov University of Massachusetts Outline
More informationImproved Kelbg potential for correlated Coulomb systems
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSA: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 36 (2003) 5957 5962 PII: S0305-4470(03)55108-3 Improved Kelbg potential for correlated Coulomb systems
More informationBosonic Path Integrals
Bosonic Path Integrals 1. Overview of effect of bose statistics 2. Permutation sampling considerations 3. Calculation of superfluid density and momentum distribution. 4. Applications of PIMC to liquid
More informationI. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS
I. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS Marus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris marus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/marus (Dated:
More informationPhysics 127c: Statistical Mechanics. Application of Path Integrals to Superfluidity in He 4
Physics 17c: Statistical Mechanics Application of Path Integrals to Superfluidity in He 4 The path integral method, and its recent implementation using quantum Monte Carlo methods, provides both an intuitive
More informationLECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS
LECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS LECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS Path-integral for lattice bosons*: oriented closed loops, of course LECTURE 3 WORM ALGORITHM
More informationSuperfluidity and Condensation
Christian Veit 4th of June, 2013 2 / 29 The discovery of superfluidity Early 1930 s: Peculiar things happen in 4 He below the λ-temperature T λ = 2.17 K 1938: Kapitza, Allen & Misener measure resistance
More informationHexatic and microemulsion phases in a 2D quantum Coulomb gas
Hexatic and microemulsion phases in a 2D quantum Coulomb gas Bryan K Clark (University of Illinois at Urbana Champaign) Michele Casula (Ecole Polytechnique, Paris) David M Ceperley (University of Illinois
More informationPath Integral Methods for Continuum Quantum Systems
Path Integral Methods for Continuum Quantum Systems D. Ceperley University of Illinois Urbana-Champaign Imaginary time path integrals Path integrals for bosons Restricted path integral method for fermions
More informationWigner crystallization in mesoscopic semiconductor structures
Wigner crystallization in mesoscopic semiconductor structures Michael Bonitz Rostock University In collaboration with: V. Golubnichy, A. Filinov, Yu. Lozovik (Moscow) Talk given at ISG, Forschungszentrum
More informationPath Integral Monte Carlo and the Worm Algorithm in the Spatial Continuum
Path Integral Monte Carlo and the Worm Algorithm in the Spatial Continuum Adrian Del Maestro Department of Physics, University of Vermont, Burlington, VT 545, USA (Dated: June 2, 214) The goal of these
More informationQuantum Monte Carlo Simulations of Exciton Condensates
Quantum Monte Carlo Simulations of Exciton Condensates J. Shumway a and D. M. Ceperley b a Dept. of Physics and Astronomy, Arizona State University, Tempe, AZ 8583 b Dept. of Physics, University of Illinois,
More informationCPP. ContributionstoPlasmaPhysics REPRINT. Editors W. Ebeling G. Fußmann T. Klinger K.-H. Spatschek
CPP ContributionstoPlasmaPhysics Editors W. Ebeling G. Fußmann T. Klinger K.-H. Spatschek Coordinating Editors M. Dewitz C. Wilke REPRINT Contrib. Plasma Phys. 51, No. 8, 687 697 (2011) / DOI 10.1002/ctpp.201100012
More informationSuperfluidity in Hydrogen-Deuterium Mixed Clusters
Journal of Low Temperature Physics - QFS2009 manuscript No. (will be inserted by the editor) Superfluidity in Hydrogen-Deuterium Mixed Clusters Soomin Shim Yongkyung Kwon Received: date / Accepted: date
More informationI. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS
I. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS Markus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus/cours.html
More informationLECTURE 4 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS II
LECTURE 4 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS II LECTURE 4 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS II Path-integral for continuous systems: oriented closed loops LECTURE 4 WORM ALGORITHM
More informationBosonic Path Integrals
Bosonic Path Integrals 1. Overview of effect of bose statistics. Permutation sampling considerations 3. Calculation of superfluid density and momentum distribution. 4. Applications of PIMC to liquid helium
More information5 Topological defects and textures in ordered media
5 Topological defects and textures in ordered media In this chapter we consider how to classify topological defects and textures in ordered media. We give here only a very short account of the method following
More informationQuantum Momentum Distributions
Journal of Low Temperature Physics, Vol. 147, Nos. 5/6, June 2007 ( 2007) DOI: 10.1007/s10909-007-9344-7 Quantum Momentum Distributions Benjamin Withers and Henry R. Glyde Department of Physics and Astronomy,
More information1 Superfluidity and Bose Einstein Condensate
Physics 223b Lecture 4 Caltech, 04/11/18 1 Superfluidity and Bose Einstein Condensate 1.6 Superfluid phase: topological defect Besides such smooth gapless excitations, superfluid can also support a very
More informationChemistry 120A 2nd Midterm. 1. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (1-electron):
April 6th, 24 Chemistry 2A 2nd Midterm. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (-electron): E n = m e Z 2 e 4 /2 2 n 2 = E Z 2 /n 2, n =, 2, 3,... where Ze is
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationComputational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen
! Computational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen WHO DID THE WORK? Miguel Morales, Livermore Carlo Pierleoni: L Aquila, Italy Jeff McMahon
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 informationUnusual criticality in a generalized 2D XY model
Unusual criticality in a generalized 2D XY model Yifei Shi, Austen Lamacraft, Paul Fendley 22 October 2011 1 Vortices and half vortices 2 Generalized XY model 3 Villain Model 4 Numerical simulation Vortices
More informationChallenges in Path Integral Monte Carlo. David Ceperley Physics UIUC
Challenges in Path Integral Monte Carlo David Ceperley Physics UIUC Imaginary time path integrals Exchange in quantum crystals The fermion sign problem Restricted paths The Penalty method Quantum dynamics?
More informationα Particle Condensation in Nuclear systems
Particle Condensation in Nuclear systems A. Tohsaki, H. Horiuchi, G. Röpke, P. Sch. T. Yamada and Y. Funaki -condensation in matter 8 Be and Hoyle state in 12 C -condensate wave function Effective GPE
More informationin-medium pair wave functions the Cooper pair wave function the superconducting order parameter anomalous averages of the field operators
(by A. A. Shanenko) in-medium wave functions in-medium pair-wave functions and spatial pair particle correlations momentum condensation and ODLRO (off-diagonal long range order) U(1) symmetry breaking
More information763620SS STATISTICAL PHYSICS Solutions 8 Autumn 2015
76360SS STATISTICAL PHYSICS Solutions 8 Autumn 05. Thermodynamics of a Black Hole The energy and entropy for a non-rotating black hole of mass M can be defined as E = Mc S = k BA, 8πl P where A = 4πR is
More informationarxiv: v1 [cond-mat.stat-mech] 1 Oct 2009
Shell Structure of Confined Charges at Strong Coupling J. Wrighton, J. Dufty, M. Bonitz, 2 and H. Kählert 2 Department of Physics, University of Florida, Gainesville, FL 326 arxiv:9.76v [cond-mat.stat-mech]
More informationFunctional renormalization for ultracold quantum gases
Functional renormalization for ultracold quantum gases Stefan Floerchinger (Heidelberg) S. Floerchinger and C. Wetterich, Phys. Rev. A 77, 053603 (2008); S. Floerchinger and C. Wetterich, arxiv:0805.2571.
More information8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles
8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles 1. Surfactant condensation: N surfactant molecules are added to the surface of water over an area
More informationLow dimensional quantum gases, rotation and vortices
Goal of these lectures Low dimensional quantum gases, rotation and vortices Discuss some aspect of the physics of quantum low dimensional systems Planar fluids Quantum wells and MOS structures High T c
More informationMotivation for Path Integral MC
Motivation for Path Integral MC We need to go beyond classical simulations since the microscopic world is quantum mechanical. PIMC makes nice generalization of classical simulations related to polymers,.
More informationInternational Summer School on Numerical Methods for Strongly Correlated Systems in Condensed Matter. Sherbrooke, May 26 -June
International Summer School on Numerical Methods for Strongly Correlated Systems in Condensed Matter Sherbrooke, May 26 -June 7 2008 The Worm Algorithm M. Boninsegni University of Alberta Credit Nikolay
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 information1 Interaction of Quantum Fields with Classical Sources
1 Interaction of Quantum Fields with Classical Sources A source is a given external function on spacetime t, x that can couple to a dynamical variable like a quantum field. Sources are fundamental in the
More informationDesign and realization of exotic quantum phases in atomic gases
Design and realization of exotic quantum phases in atomic gases H.P. Büchler and P. Zoller Theoretische Physik, Universität Innsbruck, Austria Institut für Quantenoptik und Quanteninformation der Österreichischen
More informationThermodynamics of finite quantum systems
Thermodynamics of finite quantum systems Shell structure in finite quantum systems Erice summer school, July 25-30, 2010 Klavs Hansen Department of Physics, University of Gothenburg Topics Metal clusters
More informationIntroduction to Quantum Mechanics PVK - Solutions. Nicolas Lanzetti
Introduction to Quantum Mechanics PVK - Solutions Nicolas Lanzetti lnicolas@student.ethz.ch 1 Contents 1 The Wave Function and the Schrödinger Equation 3 1.1 Quick Checks......................................
More informationStatistical Mechanics
Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2
More informationThe nature of superfluidity in the cold atomic unitary Fermi gas
The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas
More informationPath Integral Monte Carlo Simulations on the Blue Waters System. Burkhard Militzer University of California, Berkeley
Path Integral Monte Carlo Simulations on the Blue Waters System Burkhard Militzer University of California, Berkeley http://militzer.berkeley.edu Outline 1. Path integral Monte Carlo simulation method
More informationComputation of the High Temperature Coulomb Density Matrix in Periodic Boundary Conditions
Computation of the High Temperature Coulomb Density Matrix in Periodic Boundary Conditions B. Militzer Department of Earth and Planetary Science and Department of Astronomy, University of California, Berkeley,
More informationSimulating Quantum High Energy Density Plasmas: Approximations According to the Time-Dependent Variational Principle
Simulating Quantum High Energy Density Plasmas: Approximations According to the Time-Dependent Variational Principle (LANL) with contributions from Andreas Markmann (Yale), Mike Surh (LLNL), Michael S.
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
More informationFermi gas model. Introduction to Nuclear Science. Simon Fraser University Spring NUCS 342 February 2, 2011
Fermi gas model Introduction to Nuclear Science Simon Fraser University Spring 2011 NUCS 342 February 2, 2011 NUCS 342 (Lecture 9) February 2, 2011 1 / 34 Outline 1 Bosons and fermions NUCS 342 (Lecture
More informationRenormalization of microscopic Hamiltonians. Renormalization Group without Field Theory
Renormalization of microscopic Hamiltonians Renormalization Group without Field Theory Alberto Parola Università dell Insubria (Como - Italy) Renormalization Group Universality Only dimensionality and
More informationTwo-Dimensional Spin-Polarized Hydrogen at Zero
J Low Temp Phys (2013) 171:685 692 DOI 10.1007/s10909-012-0756-7 Two-Dimensional Spin-Polarized Hydrogen at Zero Temperature L. Vranješ Markić J. Boronat Received: 13 July 2012 / Accepted: 16 September
More informationPath Integral Calculations of exchange in solid 4 He
Path Integral Calculations of exchange in solid 4 He B. Bernu LPTL, UMR 7600 of CNRS, Université P. et M. Curie, Paris, France D. M. Ceperley Dept. of Physics and NCSA, University of Illinois at Urbana-Champaign,
More informationPhys 622 Problems Chapter 5
1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit
More informationLandau s Fermi Liquid Theory
Thors Hans Hansson Stockholm University Outline 1 Fermi Liquids Why, What, and How? Why Fermi liquids? What is a Fermi liquids? Fermi Liquids How? 2 Landau s Phenomenological Approach The free Fermi gas
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationSuperflow in Solid Helium
Superflow in Solid Helium Abhishek Roy December 16, 2004 Abstract The small mass and significant zero point motion of solid 4 He make it very different from its conventional counterparts. In the seventies
More informationBosonic Path Integrals
Bosonic Path Integrals 1. Overview of effect of bose statistics. Permutation sampling considerations 3. Calculation of superfluid density and momentum distribution. 4. Applications of PIMC to liquid helium
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 15 Jul 2005
Correlation functions of one-dimensional Bose-Fermi mixtures Holger Frahm and Guillaume Palacios Institut für Theoretische Physik, Universität Hannover, Appelstr. 2, 30167 Hannover, Germany (Dated: July
More informationSupersolidity of excitons
Supersolidity of excitons Michał Matuszewski Institute of Physics, Polish Academy of Sciences, Warsaw Thomas R. Taylor and Alexey V. Kavokin University of Southampton, UK ISNP 2012, Phuket Outline 1. What
More informationFirst Problem Set for Physics 847 (Statistical Physics II)
First Problem Set for Physics 847 (Statistical Physics II) Important dates: Feb 0 0:30am-:8pm midterm exam, Mar 6 9:30am-:8am final exam Due date: Tuesday, Jan 3. Review 0 points Let us start by reviewing
More informationPhysics 127b: Statistical Mechanics. Lecture 2: Dense Gas and the Liquid State. Mayer Cluster Expansion
Physics 27b: Statistical Mechanics Lecture 2: Dense Gas and the Liquid State Mayer Cluster Expansion This is a method to calculate the higher order terms in the virial expansion. It introduces some general
More informationWang-Landau sampling for Quantum Monte Carlo. Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart
Wang-Landau sampling for Quantum Monte Carlo Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart Overview Classical Monte Carlo First order phase transitions Classical Wang-Landau
More informationContents. 1.1 Prerequisites and textbooks Physical phenomena and theoretical tools The path integrals... 9
Preface v Chapter 1 Introduction 1 1.1 Prerequisites and textbooks......................... 1 1.2 Physical phenomena and theoretical tools................. 5 1.3 The path integrals..............................
More informationTrapping Centers at the Superfluid-Mott-Insulator Criticality: Transition between Charge-quantized States
Trapping Centers at the Superfluid-Mott-Insulator Criticality: Transition between Charge-quantized States Boris Svistunov University of Massachusetts, Amherst DIMOCA 2017, Mainz Institute for Theoretical
More informationMonatomic ideal gas: partition functions and equation of state.
Monatomic ideal gas: partition functions and equation of state. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Statistical Thermodynamics, MC260P105, Lecture 3,
More informationElements of Statistical Mechanics
Elements of Statistical Mechanics Thermodynamics describes the properties of macroscopic bodies. Statistical mechanics allows us to obtain the laws of thermodynamics from the laws of mechanics, classical
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationSuperfluidity. Krzysztof Myśliwy. October 30, Theoretical Physics Proseminar
Superfluidity Krzysztof Myśliwy Theoretical Physics Proseminar October 30, 2017 Outline The λ transition Phenomenology of He-II Landau theory- a semi-phenomenological approach Feynman s explanation- from
More informationLecture 10. Central potential
Lecture 10 Central potential 89 90 LECTURE 10. CENTRAL POTENTIAL 10.1 Introduction We are now ready to study a generic class of three-dimensional physical systems. They are the systems that have a central
More informationCircumventing the pathological behavior of path-integral Monte Carlo for systems with Coulomb potentials
Circumventing the pathological behavior of path-integral Monte Carlo for systems with Coulomb potentials M. H. Müser and B. J. Berne Department of Chemistry, Columbia University, New York, New York 10027
More informationINTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition
INTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition Kerson Huang CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Group an Informa
More informationQuantum Mechanics (Draft 2010 Nov.)
Quantum Mechanics (Draft 00 Nov) For a -dimensional simple harmonic quantum oscillator, V (x) = mω x, it is more convenient to describe the dynamics by dimensionless position parameter ρ = x/a (a = h )
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationSuperfluidity in bosonic systems
Superfluidity in bosonic systems Rico Pires PI Uni Heidelberg Outline Strongly coupled quantum fluids 2.1 Dilute Bose gases 2.2 Liquid Helium Wieman/Cornell A. Leitner, from wikimedia When are quantum
More informationClassical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More 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 informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
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 informationEvaporation/Condensation of Ising Droplets
, Elmar Bittner and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11, D-04109 Leipzig, Germany E-mail: andreas.nussbaumer@itp.uni-leipzig.de Recently Biskup et
More informationFilippo Tramonto. Miniworkshop talk: Quantum Monte Carlo simula9ons of low temperature many- body systems
Miniworkshop talk: Quantum Monte Carlo simulations of low temperature many-body systems Physics, Astrophysics and Applied Physics Phd school Supervisor: Dott. Davide E. Galli Outline Interests in quantum
More informationBose-Einstein condensation in Quantum Glasses
Bose-Einstein condensation in Quantum Glasses Giuseppe Carleo, Marco Tarzia, and Francesco Zamponi Phys. Rev. Lett. 103, 215302 (2009) Collaborators: Florent Krzakala, Laura Foini, Alberto Rosso, Guilhem
More informationLecture 6 Photons, electrons and other quanta. EECS Winter 2006 Nanophotonics and Nano-scale Fabrication P.C.Ku
Lecture 6 Photons, electrons and other quanta EECS 598-002 Winter 2006 Nanophotonics and Nano-scale Fabrication P.C.Ku From classical to quantum theory In the beginning of the 20 th century, experiments
More informationINTRODUCTION TO QUANTUM ELECTRODYNAMICS by Lawrence R. Mead, Prof. Physics, USM
INTRODUCTION TO QUANTUM ELECTRODYNAMICS by Lawrence R. Mead, Prof. Physics, USM I. The interaction of electromagnetic fields with matter. The Lagrangian for the charge q in electromagnetic potentials V
More informationTensor network renormalization
Coogee'15 Sydney Quantum Information Theory Workshop Tensor network renormalization Guifre Vidal In collaboration with GLEN EVENBLY IQIM Caltech UC Irvine Quantum Mechanics 1920-1930 Niels Bohr Albert
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 informationSimulations of Dense Atomic Hydrogen in the Wigner Crystal Phase
Published in Journal of Physics and Chemistry of Solids, 67 (2006) 2136 Simulations of Dense Atomic Hydrogen in the Wigner Crystal Phase Burkhard Militzer Geophysical Laboratory, Carnegie Institution of
More informationφ α (R) φ α ψ 0 e nτ(e α E T ). (2)
Atomic Scale Simulations Projector Quantum Monte Carlo Projector Monte Carlo We now turn to a potentially more powerful method where a function of the Hamiltonian projects out the the ground state, hence
More informationQUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES
QUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES 1 Tapio T. Rantala, Department of Physics, Tampere University of Technology http://www.tut.fi/semiphys CONTENTS MOTIVATION PATH INTEGRAL APPROACH
More informationQuantum impurities in a bosonic bath
Ralf Bulla Institut für Theoretische Physik Universität zu Köln 27.11.2008 contents introduction quantum impurity systems numerical renormalization group bosonic NRG spin-boson model bosonic single-impurity
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationA. Field theoretical methods for the partition function, perturbative expansions
IV. QUANTUM FIELD THEORY, GREEN S FUNCTION, AND PERTURBATION THEORY Markus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus
More informationSuperfluid 3 He. Miguel A. Morales
Superfluid 3 He Miguel A. Morales Abstract In this report I will discuss the main properties of the superfluid phases of Helium 3. First, a brief description of the experimental observations and the phase
More 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 informationEffects of spin-orbit coupling on the BKT transition and the vortexantivortex structure in 2D Fermi Gases
Effects of spin-orbit coupling on the BKT transition and the vortexantivortex structure in D Fermi Gases Carlos A. R. Sa de Melo Georgia Institute of Technology QMath13 Mathematical Results in Quantum
More informationPhysics 127a: Class Notes
Physics 127a: Class Notes Lecture 15: Statistical Mechanics of Superfluidity Elementary excitations/quasiparticles In general, it is hard to list the energy eigenstates, needed to calculate the statistical
More informationThe Quantum Hall Effect
The Quantum Hall Effect David Tong (And why these three guys won last week s Nobel prize) Trinity Mathematical Society, October 2016 Electron in a Magnetic Field B mẍ = eẋ B x = v cos!t! y = v sin!t!!
More information