Variational and Diffusion Monte Carlo in the Continuum
|
|
- Colleen Higgins
- 6 years ago
- Views:
Transcription
1 Variational and Diffusion Monte Carlo in the Continuum Describe VMC and DMC Trial functions Application to the homogeneous electron gas Lab will use QMCPACK to illustrate a simple calculation David Ceperley lecturer Jeremy McMinis Lab instructor University of Illinois Urbana-Champaign
2 Variational Monte Carlo Historically first quantum simulation method Slater-Jastrow trial function Calculations of properties: n(k). Examples: liquid helium and electron gas. Wavefunctions for Quantum solids Ewald Sums for Charged systems WaveFunctions beyond Slater-Jastrow: back flow and 3-body Twist Averaged Boundary Conditions
3 First Major QMC Calculation PhD thesis of W. McMillan (1964) University of Illinois. VMC calculation of ground state of liquid helium 4. Applied MC techniques from classical liquid theory. Ceperley, Chester and Kalos (1976) generalized to fermions. Zero temperature (single state) method VMC can be generalized to finite temperature by using trial density matrix instead of trial wavefunction.
4 Notation Individual coordinate of a particle r i All 3N coordinates R= (r 1,r,. r N ) Total potential energy = V(R) Kinetic energy : Hamiltonian : N λ i where λ i=1 Hˆ = Tˆ+ Vˆ m
5 Liquid helium the prototypic quantum fluid Interatomic potential is known more accurately than any other atom because electronic excitations are so high. A helium atom is an elementary particle. A weakly interacting hard sphere. Two isotopes: 3 He (fermion: antisymmetric trial function, spin 1/) 4 He(boson: symmetric trial function, spin zero) vlj 1 6 σ σ () r 4ε = r r
6 Helium phase diagram Because interaction is so weak helium does not crystallize at low temperatures. Quantum exchange effects are important Both isotopes are quantum fluids and become superfluids below a critical temperature. One of the goals of computer simulation is to understand these states, and see how they differ from classical liquids starting from nonrelativistic Hamiltonian: H = i m i λ m i i +V (R)
7 Variational Monte Carlo (VMC) Variational Principle. Given an appropriate trial function: Continuous Proper symmetry Normalizable Finite variance Quantum chemistry uses a product of single particle functions With MC we can use any computable function. Sample R from Ψ using MCMC. Take average of local energy: Optimize Ψ to get the best upper bound Error in energy is nd order Better wavefunction, lower variance! (non-classical) zero variance principle. E σ V dr ψ H ψ = dr ψψ V E dr ψ H ψ = E dr ψψ 1 EL( R) =R ψ ( R) Hψ( R) E = E ( R) E V L ψ 0 0
8 How do we treat spin in QMC? Spin & real vs. complex For extended systems we use the S z representation. We have a fixed number of up and down electrons and we antisymmetrize among electrons with the same spin. This leads to Slater determinants. For a given trial function, its real part is also a trial function (but it may have different symmetries) ( ikr ikr e, e ) or ( cos( kr),sin( kr) ) For the ground state, without magnetic fields or spin-orbit interaction, we can always work with real functions. However, it may be better to work with complex functions.
9 Trial function for 4 He: Jastrow or pair product We want finite variance of the local energy. Whenever atoms get close together wavefunction should vanish. The pseudopotential u(r) is similar to classical potential Local energy has the form: G is the pseudoforce: If v(r) diverges as r -n how should u(r) diverge? Assume: u(r)=r -m Keep N-1 atoms fixed and let 1 atom approach another and analyze the singular parts of the local energy. Gives a condition on u at small r. For Lennard-Jones 6-1 potential, Jastrow, u ~r -5 ψ ( R) = i< j e u( r ij ) E ( R) = v( r ) λ u( r ) λ G G ψ i i ij j ij ij i i< j i = u( r ) m 1 ( ) n εr = λ αmr for n > n m = 1 1 α = m ε λ
10 Fermions: antisymmetric trial function At mean field level the wavefunction is a Slater determinant. Orbitals for homogenous systems are a filled set of plane waves. We can compute this energy analytically (HF). To include correlation we multiply by a pseudopotential. We need MC to evaluate properties. New feature: how to compute the derivatives of a deteminant and sample the determinant. Use tricks from linear algebra. Reduces complexity to O(N ). det Ψ ( R) = Det e { ik ( )} i r jη σ s i j PBC: k L= π n+ { θ} ikirj i j Ψ SJ ( R) = Det{ e } e < Slater-Jastrow trial function. T ( ( )) = det ( ) ( M ) T ( ) ( ) u( r ) φ r φ r φ r M k j k j k j 1 k, i k 1 det 1 M = Tr M det( M) a a ij
11
12 Optimization of trial function Try to optimize u(r) using reweighting (correlated sampling) Sample R using P(R)=Ψ (R,a 0 ) Now find minima of the analytic function E v (a) Or minimize the variance (more stable but wavefunctions less accurate). Statistical accuracy declines away from a 0. New methods allow many more parameters E V ( a) = wr (, a) i = ψ( ah ) ψ( a) ψ ( a) 1 (, ) = ψ (, ) ψ(, ) eff i wr (, aer ) (, a) ψ ( Ra, ) = PR ( ) = k i i k i wi w wr (, a) ERa RaH Ra N i i
13 modern optimization With more computer time, we do a MC rw in both R and a (parameter space). Do usual VMC for a block and collect statistics on E, de/da, d E/(da i da j ). Special estimators for these quantities. Then make a change in a: a new =a old + c de/da+ Iterate until convergence. Lots more tricks to make it stable. Can do hundreds of parameters.
14 Scalar Properties, Static Correlations and Order Parameters What do we get out of a simulation? Energy by itself doesn t tell you very much. Other properties do NOT have an upper bound property Only first order in accuracy EXAMPLES Static properties: pressure, specific heat etc. Density Pair correlation in real space and fourier space. Order parameters and broken symmetry: How to tell a liquid from a solid Specifically quantum: the momentum distribution
15 Kinetic energy Potential energy Pair correlation function Structure function Pressure (virial relation) Other quantum properties Like properties from classical simulations No upper bound property Only first order in accuracy Momentum distribution Non-classical showing effects of bose or fermi statistics Fourier transform is the single particle off-diagonal density matrix Compute with McMillan Method. Condensate fraction ~10% nrr (, ') = dr... drψ ( rr,...) ψ( r', r...) = * N * ψ ( r', r,...) ψ ( rr,,...)
16 Momentum distribution Classically momentum distribution is always a Gaussian Non-classical showing effects of bose or fermi statistics Fourier transform is the single particle off-diagonal density matrix Compute with McMillan Method. For fermions we need to use the determinant update formulas to find the effect of the movement of 1 electron. Momentum Distribution 1 * n( r, r ') = dr... drnψ ( r, r...) ψ( r ', r...) Z = * ψ ( rr,...) ψ ( r', r...)
17 Derivation of momentum formula Want probability that a given atom has momentum hk. Find wavefunction in momentum space by FT wrt all the coordinates and integrating out all but one electron Expanding out the square and performing the integrals we get. Where: i( k11 r knrn) Pr( k,.. k ) = dr e Ψ( R) 1 n = dk... dk Pr( k, k,... k ) k N N For a homogeneous system, n(r,s)=n( r-s ) N drds dr ikr exp( ( )) (, ) () 3 3 e n r nk = ik r s n r s = ( π ) V (π ) V nrs (,) = dr... drψ rr,... r sr,... r Q ( ) ψ ( ) * N N N
18 Standard model for electrons in metals Basis of DFT. Characterized by dimensionless parameters: Density Temperature r = a/ a Γ= e / Ta What is energy? When does it freeze? What is spin polarization? What are properties? Γ< r s s classical OCP Γ= 175 classical melting 0 The electron gas log( r s ) H h = i + m log( Γ) 1 r i i< j ij
19 Charged systems How can we handle charged systems? Just treat like short-ranged potential: cutoff potential at r>l/. Problems: Effect of discontinuity never disappears: (1/r) (r ) gets bigger. Will violate Stillinger-Lovett conditions because Poisson equation is not satisfied Image potential solves this: V I = Σ v(r i -r j +nl) But summation diverges. We need to resum. This gives the ewald image potential. For one component system we have to add a background to make it neutral. Even the trial function is long ranged and needs to be resummed.
20 Ewald summation method Key idea is to split potential into k-space part and realspace part. We can do since FT is linear. V = V = i<j,l k φ(r -r +nl) k i j ( ) k N φ ρ where ρ = 1 ik gr and ϕk = dre φ( r) Ω 4π e For φ(r)=e /r ϕk = k Bare potential converges slowly at large r (in r-space) and at large k (in k-space) k i e ik gr i
21 Classic Ewald Split up using Gaussian charge distribution erfc( κ r) ( r) = decays fast at large r r φ φ k = 4π e k ( k / κ ) decays fast at large k κ = convergence parameter If we make it large enough we can use the minimum image potential in r-space. V Extra term for insulators: dipole = π µ i (ε + 1) Ω i
22 Jastrow factor for the e-gas Look at local energy either in r space or k-space: r-space: as electrons get close gives cusp condition: du/dr 0 =-1 k-space, charge-sloshing or plasmon modes. ρu k = Can combine exact properties in the Gaskell form. Write E V in terms structure factor making random phase approximation. (RPA). 1 1 V ρ uk = S + + S ideal structure factor k S k k = k k λ V k λk 1 k Optimization can hardly improve this form for the e-gas in either or 3 dimensions. RPA works better for trial function than for the energy. NEED EWALD SUMS because potential trial function is long range, it also decays as 1/r, but it is not a simple power. 1 r 3D = log( r) 1D 1/ lim r ur ( ) r D Long range properties Give rise to dielectric properties Energy is insensitive to u k at small k Those modes converge t~1/k
23 Spline Jastrow factor For the HEG, the most general Jastrow factor has the form: u(r ) = u sr (r) + u k e i k r k u sr (r) = 0 for r > L / u (r) must be continuous, with a continuous derivative. We can impose the cusp condition at r=0, and BC at r=l/. It is a smooth function: represent it as piecewise cubic polynomial in the region 0< r < L/. M Knots at b n. Total number of unknowns is M-1 Also u k k-space Jastrows. Do we use RPA values?
24 Generalized Feynman-Kacs formula gives relation between trial function and exact wavefunction! average population starting from a single point R 0 after a t imaginary time t : ( ; ) ( R ) ( R ) 0 0 ( ) 0 ( R) ( R ) ψ ( E ) P R t = dr R e R = e t H T 0 0 ψ 0 expand the density matrix in terms of exact eigenstates lim P R ; t t 0 ~ ( R) ( ) R0 α φ0( R0) ψ ( R ) ψ * P( R0; t) = dr φα ( R) φα ( R0) e ψ φ ψ e t = dt E () t L ψ ψφ 0 ( E ) t E α T dte L () t ψ
25 Wavefunctions beyond Jastrow Use method of residuals construct a sequence of increasingly better trial wave functions. Justify from the Importance sampled DMC. Zeroth order is Hartree-Fock wavefunction First order is Slater-Jastrow pair wavefunction (RPA for electrons gives an analytic formula) Second order is 3-body backflow wavefunction Three-body form is like a squared force. It is a bosonic term that does not change the nodes. exp{ ξ [ ij( rij)( ri rj)] } i j φ n+ 1 0 E ( R) φ ( R) e i j φ = e = V( R) φ = φ e n k j rj ( ) U R 1 n τ< φ Hφ > n smoothing [ ] k r ( ) ( ) E = U( R) W( R) + i Y R j j j j
26 Backflow wave function Backflow means change the coordinates to quasi- coordinates. Leads to a much improved energy and to improvement in nodal surfaces. Couples nodal surfaces together. ikr ikx i j i j Det{ e } Det{ e } xi = ri + ηij( rij)( ri rj) j 3DEG Kwon PRB 58, 6800 (1998).
27 Dependence of energy on wavefunction 3d Electron fluid at a density r s =10 Kwon, Ceperley, Martin, Phys. Rev. B58,6800, 1998 Wavefunctions Slater-Jastrow (SJ) three-body (3) backflow (BF) fixed-node (FN) Energy <f H f> converges to ground state Variance <f [H-E] f> to zero. Using 3B-BF gains a factor of 4. Using DMC gains a factor of 4. Energy FN -SJ FN-BF Variance
28 Twist averaged boundary conditions In periodic boundary conditions (Γ point), the wavefunction is periodicè Large finite size effects for metals because of shell effects. Fermi liquid theory can be used to correct the properties. In twist averaged BC we use an arbitrary phase as r è r+l If one integrates over all phases the momentum distribution changes from a lattice of k-vectors to a fermi sea. Smaller finite size effects ϕ = kl ( π ) e ikr = πn + θ k x iθ Ψ ( x + L) = e Ψ( x) π 1 3 A d A 3 θ θ π = Ψ Ψ θ PBC TABC
29 Twist averaged MC Make twist vector dynamical by changing during the random walk. π < θi π i= ( 1,,3) Within GCE, change the number of electrons Within TA-VMC Initialize twist vector. Run usual VMC (with warmup) Resample twist angle within cube (iterate) Or do in parallel.
30 Grand Canonical Ensemble QMC GCE at T=0K: choose N such that E(N)-µN is minimized. According to Fermi liquid theory, interacting states are related to noninteracting states and described by k. Instead of N, we input the fermi wavevector(s) k F. Choose all states with k < k F (assuming spherical symmetry) N will depend on the twist angle = number of points inside a randomly placed sphere. k n = π L k n k F θ n + L After we average over twist we get a sphere of filled states. Is there a problem with Ewald sums as the number of electrons varies? No! average density is exactly that of the background. We only work with averaged quantities.
31 Single particle size effects Exact single particle properties with TA within HF Implies momentum distribution is a continuous curve with a sharp feature at k F. With PBC only 5 k points for k<k F Holzmann et al. PRL 107,11040 (011) No size effect within single particle theory! Kinetic energy will have much smaller size effects. T = d 3 k k m n( k)
32 Potential energy Write potential as integral over structure function: 4π k 3 ir ( 1 r) k V = d k S( k) S( k) = ρ 1 ( 1) kρk = + N e Error comes from effects. Approximating integral by sum Finite size effects in S(k) at a given k. Within HF we get exact S(k) with TABC. Discretization errors come only from non-analytic points of S(k). ( k) = 1 δ + HF q q ' k N qq, ' the absence of the k=0 term in sum. We can put it in by hand since we know the limit S(k) at small k (plasmon regime) S Remaining size effects are smaller, coming from the non-analytic behavior of S(k) at k F.
33 Brief History of Ferromagnetism in electron gas What is polarization state of fermi liquid at low density? N N ζ = N + N Bloch 199 got polarization from exchange interaction: r s > 5.4 3D r s >.0 D Stoner 1939: include electron screening: contact interaction Herring 1960 Ceperley-Alder 1980 r s >0 is partially polarized Young-Fisk experiment on doped CaB r s ~5. Ortiz-Balone 1999 : ferromagnetism of e gas at r s >0. Zong et al Redo QMC with backflow nodes and TABC.
34 Polarization of 3DEG We see second order partially polarized transition at r s =5 Is the Stoner model (replace interaction with a contact potential) appropriate? Screening kills long range interaction. Wigner Crystal at r s =105 Twist averaging makes calculation possible--much smaller size effects. Jastrow wavefunctions favor the ferromagnetic phase. Backflow 3-body wavefunctions more paramagnetic Polarization transition
35 Partially polarized phase at low density. But at lower energy and density than before. As accuracy gets higher, polarized phase shrinks Phase Diagram Real systems have different units.
36 Summary of Variational (VMC) error (au) 1.E+01 1.E+00 1.E-01 1.E-0 1.E-03 1.E-04 Better trial function Simple trial function 1.E-05 1.E+00 1.E+01 1.E+0 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 computer time (sec)
37 Summary and problems with variational methods Powerful method since you can use any trial function Scaling (computational effort vs. size) is almost classical Learn directly about what works in wavefunctions No sign problem Optimization is time consuming Energy is insensitive to order parameter Non-energetic properties are less accurate. O(1) vs. O() for energy. Difficult to find out how accurate results are. Favors simple states over more complicated states, e.g. Solid over liquid Polarized over unpolarized What goes into the trial wave function comes out! GIGO We need a more automatic method! Projector Monte Carlo
38 Projector Monte Carlo Originally suggested by Fermi and implemented in 1950 by Donsker and Kac for H atom. Practical methods and application developed by Kalos:
39 Diffusion Monte Carlo How do we analyze this operator? Expand into exact eigenstates of H. Then the evolution is simple in this basis. Long time limit is lowest energy state that overlaps with the initial state, usually the ground state. How to carry out on the computer? ψ (H E T )t ( Rt, ) = e ψ( R,0) Hφ ψ( R,0) = φ ( R) φ ψ(0) t( Eα ET ) ( Rt, ) = α( Re ) α (0) α t( E0 ET ) lim ψ( Rt, ) = φ ( Re ) φ ψ(0) 0 = t T E φ α α α α α ψ φ φ ψ 0 0 E E normalization fixed α
40 Monte Carlo process Now consider the variable t as a continuous time (it is really imaginary time). Take derivative with respect to time to get evolution. This is a diffusion + branching process. Justify in terms of Trotter s formula. Requires interpretation of the wavefunction as a probability density. But is it? Only in the boson ground state. Otherwise there are nodes. Come back to later. ψ ( Rt, ) = ( H ET ) ψ ( R, t) t h H = i + V( R) m i i ψ ( Rt, ) h = iψ ( Rt, ) t m i ψ ( Rt, ) = ( V( R) ET ) ψ ( R, t) t i
41 How do we find the solution of: The operator solution is: Trotter s formula (1959): Trotter s formula d ˆ ρ = ( A+ B) ˆ ρ dt ( A+ B) t ˆ ρ = e ˆ ρ = lim n Assumes that A,B and A+B are reasonable operators. t ˆ t ˆ t ˆ t ˆ t ˆ t ˆ A B n n n A B A B n n n n 0 n = 0 ' 1 ' n 1 ' n ' n n R e e R R e R R e R R e R R e R This means we just have to figure out what each operator does independently and then alternate their effect. This is rigorous in the limit as nè. In the DMC case A is diffusion operator, B is a branching operator. Just like molecular dynamics : at small time we evaluate each operator separately. e t n Aˆ e t n Bˆ n
42 Evaluation of kinetic density matrix r e τ T r ' = α φ α * (r)φ α (r ')e τt α In PBC eigenfunctions of ˆ T = 1 Ω e i k r and eigenvalues are λk r e τ T r ' = k 1 Ω e i convert to an integral r e τ T r ' = 1 dke i ( π ) 3 k r e i k r ' e τλk k ( r ' r ) τλk = ( 4πλτ ) 3/ e ( r r ' ) /4λτ Danger: makes assumption about boundaries and statistics. This is a diffusion process.
43 Putting this together n is number of time slices. τ is the time-step ˆρ = e β (T +V ) ˆρ = lim n τ = β / n e τ T τ V e n V is diagonal r e τ T r ' = ( 4πλτ ) 3/ e ( r r ' ) /4λτ r e τ V τv (r ) r ' = δ (r r ')e nh τ ˆ ˆ ( ˆ 1) ( ) 0 ~ τt τ V R V Rn τt n n 1 n Re R R e R e R e R e τ Error at finite n comes from commutator is roughly: Diffusion preserves normalization but potential does not! e τ ˆ, ˆ TV
44 Basic DMC algorithm Construct an ensemble (population P(0)) sampled from the trial wavefunction. {R 1,R,,R P } Go through ensemble and diffuse each one (timestep τ) R' = R + λτζ ( t) k k ndrn uprn τ ( V( R) ET ) number of copies= e + u floor function Trial energy E T adjusted to keep population fixed. E drhφ( R, t) = lim V( R) drφ ( R, t) 0 t φ ( ) Problems: 1. Branching is uncontrolled. Population unstable 3. What do we do about fermi statistics?
45 Population Bias Having the right trial energy guarantees that population will on the average be stable, but fluctuations will always cause the population to either grow too large or too small. Various ways to control the population Suppose P 0 is the desired population and P(t) is the current population. How much do we have to adjust E T to make P(t+T)=P 0? T( δ E ) Pt ( + T) = e T Pt ( ) = P Feedback procedure: δ E T ln( Pt ( ) / P0 ) = T ( ) E = E +κ ln P/ P T T0 0 There will still be a (much smaller) bias in the energy caused by a limited population. 0
46 Importance Sampling Why should we sample the wavefunction? The physically correct pdf is ϕ. Importance sample (multiply) by trial wave function. f( R, t) ψ ( R) φ( R, t) lim f( R, t) ψ ( R) φ ( R) T t T f( R,) t = ψt ( RH ) f( Rt, )/ t f( R, t) = t [ ψ ( R) ] ( ) ( 1 ln + ) T T T λ f λ f ψ R ψ Hψ ( ) f( R, t) Evolution = diffusion + drift + branching We have three terms in the evolution equation. Trotter s formula still applies. T 0 Commute Ψ through H
47 To the pure diffusion algorithm we have added a drift step that pushes the random walk in directions of increasing trial function: R' = R+ λτ ln ψ T ( R) Branching is now controlled by the local energy E L (R) E T =ψ 1 (R)H ψ (R) E T Because of zero variance principle, fluctuations are controlled. Cusp condition can limit infinities coming from singular potentials. We still determine E T by keeping asymptotic population stable. drφ( R, t) HψT ( R) E0 lim t Eψ ( R) drf ( R, t) = f ( ) Must have accurate time evolution. Adding accept/reject step is a major improvement.
48 Importanced sampled Green s function: ( R ') ( R) ψ τ H GR ( R') = Re R' ψ Exact property of DMC Green s function ( R) G( R R' ) ( R' ) G( R' R) Ψ = Ψ We enforce detailed balance to decrease time step errors. VMC satisfies detailed balance. Typically we choose time step to have 99% acceptance ratio. Method gives exact result if either time step is zero or trial function is exact. Gs ( ' s) ψ ( s') As ( s' ) = min 1, Gs ( s') ψ ( s)
49 Ensemble evolves according to Schematic of DMC Diffusion Drift branching ensemble
50
51 Mixed estimators Problem is that PMC samples the wrong distribution. OK for the energy Linear extrapolation helps correct this systematic error Other solutions: Maximum overlap Forward walking Reptation/path integrals A A A M o V 0 0 * drψ R Aφ R * drψ R φ R * drφ R Aφ R * drφ R φ R * drψ R Aψ R * drψ R ψ R ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ( φ ψ) ) ( φ ψ) ( ( φ ψ) ) A ; A A + O A M ; A A M V V M V + O A = A dr for the density minimized wrt A
52 Forward Walking Let s calculate the average population resulting from DMC starting from a single point R 0 after a time `t. ( ; ) We can estimate the correction to the mixed estimator by weighting with the number of descendants of a given configuration. 1 A = lim ( ; ) ( ) 0 t P Ri t A Ri M Problem: the fluctuations in the weights eventually diverge. Don t make t too large. ( ) ( R) ( R ) ψ ( E ) P R t dr R e R t H T 0 = 0 ψ 0 expand the density matrix in terms of exact eigenstates lim P R ; t t ( R) ( ) R0 α φ0( R0) ψ ( R ) ψ P( R0; t) = dr φα ( R) φα ( R0) e ψ = ψφ i ( E ) t H T
53 Fusion sticking coefficient Phys. Rev. A 31, 1999 (1985). Consider the 3 body system (µ d t) For the sticking coefficient, we need the exact wavefunction at the point where nuclei are at the same position. (this is a singular point) ( r r ) ψ = r, 1 3
54 Other projector functions can be used 0 ( E E ) τ T e Diffusion MC 1 GE ( ) = 1 + τ ( E ET ) Green's Function MC G 1 τ ( E E ) Power MC T GE ( ) = 1 ground state remains after many iterations T dg τ = = de time step for all 3 cases: lim GE ( ) n n = e nτ ( E E ) T E Common effect on long-time (iteration) limit. 3 rd choice generates a Krylov sequence. Only works for bounded spectra such as a lattice model.
55 Green s Function Monte Carlo Kalos, Levesque, Verlet Phys. Rev. A9, 178 (1974). It is possible to make a zero time-step-error method Works with the integral formulation of DMC τ ( 1 T ) 0 dβ GRR (, ') R 1 H E R' e β τ = + = τ 1 + H E T Sample time-step from Poisson distribution Express operator in a series expansion and sample the terms stochastically. GRR (, ') = HRR (, ') + drgrr " (, ") KR ( ", R') Recent Revival: Continuous time Monte Carlo for lattice models.
56 Fermions? How can we do fermion simulations? The initial condition can be made real but not positive (for more than 1 electron in the same spin state) In transient estimate or released-node methods one carries along the sign as a weight and samples the modulus. (Hˆ E )t T φ() t = e sign( φ( R,0)) φ( R,0) Do not forbid crossing of the nodes, but carry along sign when walks cross. What s wrong with node release: Because walks don t die at the nodes, the computational effort increases (bosonic noise) The signal is in the cancellation which dominates Monte Carlo can add but not subtract
57 ( β ) β H Ψ = e Ψ Transient Estimate Approach βh τh τh ( ) ( ) 0 p ( 0) 0 1 p 1 p ( p) Z( β) = Ψ β Ψ β = Ψe Ψ = dr... dr Ψ R R e R... R e R Ψ R ( ) HΨ( ) ( β) Ψ( β) Ψ β β β E( β) = = EL( R0 ) τ = β Ψ p ψ(β) converges to the exact ground state E is an upper bound converging to the exact answer monotonically fermi bose ( R ) σ ( R ) τh τh ( ) ( R ) σ ( R ) σ ( R ) Z( β) = dr... dr Ψ R R e R... R e R Ψ Z Z = σ 0 p p 1 p p 0 P 0 P
58 Model fermion problem: Particle in a box Symmetric potential: V(r) =V(-r) Antisymmetric state: f(r)=-f(-r) Initial (trial) state Positive walkers Final (exact) state Node Negative walkers Sign of walkers fixed by initial position. They are allowed to diffuse freely. f(r)= number of positive-negative walkers. Node is dynamically established by diffusion process. (cancellation of positive and negative walkers.) Et () σ(0) σ( tet ) ( ) = σ(0) σ( t)
59 Scaling in Released-Node Initial distribution Later distribution At any point, positive and negative walkers will tend to cancel so the signal is drown out by the fluctuations. e Signal/noise ratio is : t=projection time E F and E B are Fermion, Bose energy (proportional to N) Converges but at a slower rate. Higher accuracy, larger t. For general excited states: EF ef Exponential complexity! t[ E F EB ] CPUtime ε (1+ ) N E E Not a fermion problem but an excited state problem. Cancellation is difficult in high dimensions. g ε g
60 Exact fermion calculations Possible for the electron gas for up to 60 electrons. DEG at rs=1 N=6 Transient estimate calculation with SJ and BF-3B trial functions. Ψ T th e Ψ T
61 General statement of the fermion problem Given a system with N fermions and a known Hamiltonian and a property O. (usually the energy). How much time T will it take to estimate O to an accuracy ε? How does T scale with N and ε? If you can map the quantum system onto an equivalent problem in classical statistical mechanics then: T N ε α With 0 <a < 4 This would be a solved quantum problem! All approximations must be controlled! Algebraic scaling in N! e.g. properties of Boltzmann or Bose systems in equilibrium.
62 Fixed-node method Initial distribution is a pdf. It comes from a VMC simulation. Drift term pushes walks away from the nodes. Impose the condition: This is the fixed-node BC Will give an upper bound to the exact energy, the best upper bound consistent with the FNBC. f( R,0) = ψ ( R) φ( R) = 0 when ψ ( R) = 0. E FN 0 T E = E if φ ( R) ψ( R) 0 all R FN E 0 0 f(r,t) has a discontinuous gradient at the nodal location. Accurate method because Bose correlations are done exactly. Scales well, like the VMC method, as N 3. Classical complexity. Can be generalized from the continuum to lattice finite temperature, magnetic fields, One needs trial functions with accurate nodes. T
63 Proof of fixed-node theorem Suppose we solve S.E. in a subvolume V determined by the nodes of an antisymetric trial function. Ĥφ FN = E FN φ FN inside V Extend the solution to all space with the permutation operator. ( ) ˆφ FN (R) 1 1 P φ N! FN PR P ( ) Inside a given sub-volume only permutations of a given sign (±) contribute. Hence the extended solution is non-zero. Evaluate the variational energy the extended trial function. ( 1) P+P' * drφ FN ( PR) Ĥφ FN ( P' R ) PP' E 0 = E ( 1) P+P' * FN E VMC drφ FN ( PR)φ FN ( P' R) PP' Edges of volumes do not contribute to the integral since the solution vanishes
64 Nodal Properties If we know the sign of the exact wavefunction (the nodes), we can solve the fermion problem with the fixed-node method. If f(r) is real, nodes are f(r)=0 where R is the 3N dimensional vector. Nodes are a 3N-1 dimensional surface. (Do not confuse with single particle orbital nodes!) Coincidence points r i = r j are 3N-3 dimensional hyper-planes In 1 spatial dimension these points exhaust the nodes: fermion problem is easy to solve in 1D with the no crossing rule. Coincidence points (and other symmetries) only constrain nodes in higher dimensions, they do not determine them. The nodal surfaces define nodal volumes. How many nodal volumes are there? Conjecture: there are typically only different volumes (+ and -) except in 1D. (but only demonstrated for free particles.)
65 Free electron Other electrons Nodal Picture: d slice thru 3d space Nodes pass thru their positions Divides space into regions Wavelength given by interparticle spacing
66 Fixed-Phase method Generalize the FN method to complex trial functions: Since the Hamiltonian is Hermitian, the variational energy is real: E V [ ] + λ[ I ] RU( R) dr e V R + λ U R λ R U R UR = ( ) ( ) ( ) ( ) dr e RU( R) ( ) We see only one place where the energy depends on the phase of the wavefunction. If we fix the phase, then we add this term to the potential energy. In a magnetic field we get also the vector potential. ( ) effective potential= V( R) + λi Ari +I iur ( ) We can now do VMC or DMC and get upper bounds as before. The imaginary part of the local energy will not be zero unless the right phase is used. Used for twisted boundary conditions, magnetic fields, vortices, phonons, spin states, i ( ) Ψ R = e U R
67 Problem with core electrons Bad scaling in both VMC and DMC In VMC, energy fluctuations from core dominate the calculation In DMC, time step will be controlled by core dynamics Solution is to eliminate core states by a pseudopotential Conventional solution: semi-local form rvˆ e core r' = vlocal ( r) δ ( r r') + vl ( r) Pl (cos( r r')) Ensures that valence electrons go into well defined valence states with the wavefunction and energy for each angular momentum state prescribed. PP is non-local: OK for VMC. Leads to an extra MC integral. But DMC uses a locality approximation and good trial functions. Extra approximation. l
68 Summary of T=0 methods: Variational(VMC), Fixed-node(FN), Released-node(RN) error (au) 1.E+01 1.E+00 1.E-01 1.E-0 1.E-03 Better trial function Simple trial function VMC FN 1.E-04 RN 1.E-05 1.E+00 1.E+01 1.E+0 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 computer time (sec)
69 Problems with projector methods Fixed-node is a super-variational method DMC dynamics is determined by Hamiltonian Zero-variance principle allows very accurate calculation of ground state energy if trial function is good. Projector methods need a trial wavefunction for accuracy. They are essentially methods that perturb from the trial function to the exact function. (Note: if you don t use a trial function, you are perturbing from the ideal gas) Difficulty calculating properties other than energy. We must use extrapolated estimators or forward walking. f( R, ) =φ ( R) ψ ( R) not φ ( R) 0 T 0 Bad for phase transitions and finite temperature, complex systems. Path Integral MC (reptation MC) solves some of these problems.
Quantum Monte Carlo. QMC methods in the continuum
Quantum Monte Carlo Premise: need to use simulation techniques to solve manybody quantum problems just as you need them classically. Both the wavefunction and expectation values are determined by the simulations.
More informationthe abdus salam international centre for theoretical physics
united nations educational, scientific and cultural organization the abdus salam international centre for theoretical physics international atomic energy agency SMR 1595-7 Joint DEMOCRIOS - ICP School
More informationIntroduction to Quantum Monte Carlo Methods Applied to the Electron Gas
Introduction to Quantum Monte Carlo Methods Applied to the Electron Gas D. M. Ceperley Physics Department and NCSA, University of Illinois Urbana-Champaign, Urbana, IL, USA 1. Introduction In these lectures,
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 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 informationI. QUANTUM MONTE CARLO METHODS: INTRODUCTION AND BASICS
I. QUANTUM MONTE CARLO METHODS: INTRODUCTION AND BASICS Markus Holzmann LPMMC, UJF, Grenoble, and LPTMC, UPMC, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus (Dated: January 24, 2012)
More informationMany-body wavefunctions for normal liquid 3 He
Many-body wavefunctions for normal liquid 3 He Markus Holzmann, 1 Bernard Bernu, 1 and D. M. Ceperley 2 1 LPTMC, UMR 7600 of CNRS, Université Pierre et Marie Curie, Paris, France 2 Department of Physics
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 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 informationMethods of Fermion Monte Carlo
Methods of Fermion Monte Carlo Malvin H. Kalos kalos@llnl.gov Institute for Nuclear Theory University of Washington Seattle, USA July, 2013 LLNL-PRES-XXXXXX This work was performed under the auspices of
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationDiffusion Monte Carlo
Diffusion Monte Carlo Notes for Boulder Summer School 2010 Bryan Clark July 22, 2010 Diffusion Monte Carlo The big idea: VMC is a useful technique, but often we want to sample observables of the true ground
More informationTwist-averaged boundary conditions in continuum quantum Monte Carlo algorithms
PHYSICAL REVIEW E, VOLUME 64, 016702 Twist-averaged boundary conditions in continuum quantum Monte Carlo algorithms C. Lin, F. H. Zong, and D. M. Ceperley Department of Physics and NCSA, University of
More informationMomentum Distribution and Effective Mass of Jellium and Simple Metals. Markus Holzmann. LPTMC, CNRS, Universite Pierre et Marie Curie, Paris
Momentum Distribution and Effective Mass of Jellium and Simple Metals Markus Holzmann LPTMC, CNRS, Universite Pierre et Marie Curie, Paris LPMMC, CNRS, Université Joseph Fourier, Grenoble, France QMC:
More informationFermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations
Fermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations Aurel Bulgac,, Joaquin E. Drut and Piotr Magierski University of Washington, Seattle, WA
More 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 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 informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationQuantum Monte Carlo methods
Quantum Monte Carlo methods Lubos Mitas North Carolina State University Urbana, August 2006 Lubos_Mitas@ncsu.edu H= 1 2 i i 2 i, I Z I r ii i j 1 r ij E ion ion H r 1, r 2,... =E r 1, r 2,... - ground
More informationQuantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid
Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge
More informationWave Function Optimization and VMC
Wave Function Optimization and VMC Jeremy McMinis The University of Illinois July 26, 2012 Outline Motivation History of Wave Function Optimization Optimization in QMCPACK Multideterminant Example Motivation:
More informationAn Overview of Quantum Monte Carlo Methods. David M. Ceperley
An Overview of Quantum Monte Carlo Methods David M. Ceperley Department of Physics and National Center for Supercomputing Applications University of Illinois Urbana-Champaign Urbana, Illinois 61801 In
More informationNeighbor Tables Long-Range Potentials
Neighbor Tables Long-Range Potentials Today we learn how we can handle long range potentials. Neighbor tables Long-range potential Ewald sums MSE485/PHY466/CSE485 1 Periodic distances Minimum Image Convention:
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationIs a system of fermions in the crossover BCS-BEC. BEC regime a new type of superfluid?
Is a system of fermions in the crossover BCS-BEC BEC regime a new type of superfluid? Finite temperature properties of a Fermi gas in the unitary regime Aurel Bulgac,, Joaquin E. Drut, Piotr Magierski
More informationObservations on variational and projector Monte Carlo Methods
Observations on variational and projector Monte Carlo Methods C. J. Umrigar Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853, USA (Dated: September 29, 2015) Variational
More informationDepartment of Physics and NCSA University of Illinois, Urbana-Champaign, IL, 61801, USA
UNDERSTANDING ELECTRONIC WAVE FUNCTIONS D. M. Ceperley Department of Physics and NCSA University of Illinois, Urbana-Champaign, IL, 61801, USA INTRODUCTION In this article I discuss some aspects of what
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 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 information7.4. Why we have two different types of materials: conductors and insulators?
Phys463.nb 55 7.3.5. Folding, Reduced Brillouin zone and extended Brillouin zone for free particles without lattices In the presence of a lattice, we can also unfold the extended Brillouin zone to get
More informationIs the homogeneous electron gas homogeneous?
Is the homogeneous electron gas homogeneous? Electron gas (jellium): simplest way to view a metal homogeneous and normal Hartree-Fock: simplest method for many-electron systems a single Slater determinant
More informationQMC dissociation energy of the water dimer: Time step errors and backflow calculations
QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.
More informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationElectronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node
Electronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node approximations Cody Melton, M. Chandler Bennett, L. Mitas, with A. Ambrosetti, F. Pederiva North Carolina State
More informationCondensed matter theory Lecture notes and problem sets 2012/2013
Condensed matter theory Lecture notes and problem sets 2012/2013 Dmitri Ivanov Recommended books and lecture notes: [AM] N. W. Ashcroft and N. D. Mermin, Solid State Physics. [Mar] M. P. Marder, Condensed
More informationAb-initio molecular dynamics for High pressure Hydrogen
Ab-initio molecular dynamics for High pressure Hydrogen Claudio Attaccalite Institut d'electronique, Microélectronique et Nanotechnologie (IEMN), Lille Outline A brief introduction to Quantum Monte Carlo
More informationPhysics 127c: Statistical Mechanics. Weakly Interacting Fermi Gas. The Electron Gas
Physics 7c: Statistical Mechanics Wealy Interacting Fermi Gas Unlie the Boson case, there is usually no ualitative change in behavior going from the noninteracting to the wealy interacting Fermi gas for
More informationPhase Diagram of the Two-Dimensional Homogeneous Electron Gas
Phase Diagram of the Two-Dimensional Homogeneous Electron Gas Neil D. Drummond and Richard J. Needs TCM Group, Cavendish Laboratory, University of Cambridge ESDG Meeting Wednesday 13th February, 2008 Two-Dimensional
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 informationVan der Waals Interactions Between Thin Metallic Wires and Layers
Van der Waals Interactions Between Thin Metallic Wires and Layers N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, United Kingdom Quantum Monte
More informationLectures on Quantum Monte Carlo. 1 Random Walks 3. 2 Variational Monte Carlo Projector Monte Carlo Path Integral Monte Carlo 27
Lectures on Quantum Monte Carlo David M. Ceperley May 1, 1996 Contents 1 Random Walks 3 1.1 Metropolis : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 4 1.2 Dynamical
More informationContinuum variational and diffusion quantum Monte Carlo calculations
Continuum variational and diffusion quantum Monte Carlo calculations R J Needs, M D Towler, N D Drummond and P López Ríos Theory of Condensed Matter Group, Cavendish Laboratory, Cambridge CB3 0HE, UK Abstract.
More informationBEC-BCS Crossover in Cold Atoms
BEC-BCS Crossover in Cold Atoms (2 years later...) Andrew Morris Pablo López Ríos Richard Needs Theory of Condensed Matter Cavendish Laboratory University of Cambridge TTI 31 st July 2009 Outline Theory
More informationLattice Vibrations. Chris J. Pickard. ω (cm -1 ) 200 W L Γ X W K K W
Lattice Vibrations Chris J. Pickard 500 400 300 ω (cm -1 ) 200 100 L K W X 0 W L Γ X W K The Breakdown of the Static Lattice Model The free electron model was refined by introducing a crystalline external
More informationFermi liquid theory Can we apply the free fermion approximation to a real metal? Phys540.nb Strong interaction vs.
Phys540.nb 7 Fermi liquid theory.. Can we apply the free fermion approximation to a real metal? Can we ignore interactions in a real metal? Experiments says yes (free fermion models work very well), but
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 informationSupplementary Material: Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas
Supplementary Material: Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas (Dated: February 26, 2013) In the main text, we describe our calculations of the 3D finite-temperature
More informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
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 informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationIdentical Particles. Bosons and Fermions
Identical Particles Bosons and Fermions In Quantum Mechanics there is no difference between particles and fields. The objects which we refer to as fields in classical physics (electromagnetic field, field
More informationSpecific heat of a fermionic atomic cloud in the bulk and in traps
Specific heat of a fermionic atomic cloud in the bulk and in traps Aurel Bulgac,, Joaquin E. Drut, Piotr Magierski University of Washington, Seattle, WA Also in Warsaw Outline Some general remarks Path
More informationProject: Vibration of Diatomic Molecules
Project: Vibration of Diatomic Molecules Objective: Find the vibration energies and the corresponding vibrational wavefunctions of diatomic molecules H 2 and I 2 using the Morse potential. equired Numerical
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 informationarxiv: v1 [cond-mat.mes-hall] 10 Feb 2010
Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationEXCHANGE IN QUANTUM CRYSTALS: MAGNETISM AND MELTING OF THE LOW DENSITY 2D WIGNER CRYSTAL. David M. CEPERLEY
united nations educational, scientific and cultural organization the abdus salam international centre for theoretical physics international atomic energy agency SMR 1595-28 Joint DEMOCRITOS - ICTP School
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 informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationOptimization of quantum Monte Carlo wave functions by energy minimization
Optimization of quantum Monte Carlo wave functions by energy minimization Julien Toulouse, Roland Assaraf, Cyrus J. Umrigar Laboratoire de Chimie Théorique, Université Pierre et Marie Curie and CNRS, Paris,
More informationElectrons in Crystals. Chris J. Pickard
Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity
More information2D Bose and Non-Fermi Liquid Metals
2D Bose and Non-Fermi Liquid Metals MPA Fisher, with O. Motrunich, D. Sheng, E. Gull, S. Trebst, A. Feiguin KITP Cold Atoms Workshop 10/5/2010 Interest: A class of exotic gapless 2D Many-Body States a)
More informationInteracting Fermi Gases
Interacting Fermi Gases Mike Hermele (Dated: February 11, 010) Notes on Interacting Fermi Gas for Physics 7450, Spring 010 I. FERMI GAS WITH DELTA-FUNCTION INTERACTION Since it is easier to illustrate
More information(1) Consider a lattice of noninteracting spin dimers, where the dimer Hamiltonian is
PHYSICS 21 : SISICL PHYSICS FINL EXMINION 1 Consider a lattice of noninteracting spin dimers where the dimer Hamiltonian is Ĥ = H H τ τ Kτ where H and H τ are magnetic fields acting on the and τ spins
More informationLuigi Paolasini
Luigi Paolasini paolasini@esrf.fr LECTURE 7: Magnetic excitations - Phase transitions and the Landau mean-field theory. - Heisenberg and Ising models. - Magnetic excitations. External parameter, as for
More informationQuantum Monte Carlo simulation of spin-polarized tritium
Higher-order actions and their applications in many-body, few-body, classical problems Quantum Monte Carlo simulation of spin-polarized tritium I. Bešlić, L. Vranješ Markić, University of Split, Croatia
More informationQuantum Monte Carlo wave functions and their optimization for quantum chemistry
Quantum Monte Carlo wave functions and their optimization for quantum chemistry Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France CEA Saclay, SPhN Orme des Merisiers April 2015 Outline
More informationChapter 29. Quantum Chaos
Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical
More informationV. ELECTRON GAS, RANDOM PHASE APPROXIMATION
V. ELECTRON GAS, RANDOM PHASE APPROXIMATION Marus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris marus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/marus/cours.html (Dated:
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 informationVariational Monte Carlo Notes for Boulder Summer School Bryan Clark. July 21, 2010
Variational Monte Carlo Notes for Boulder Summer School 010 Bryan Clark July 1, 010 VMC The Variational Monte Carlo method helps to accomplish two things: Given a wave function, T,extractinformationoutofit(i.e.
More informationDynamics of Quantum Many-Body Systems from Monte Carlo simulations
s of I.N.F.N. and Physics Department, University of Trento - ITALY INT - Seattle - October 24, 2012 Collaborators s of Alessandro Roggero (Trento) Giuseppina Orlandini (Trento) Outline s of : transform
More information8.3. Feynman diagram. Phys620.nb multi-particle Green s functions the equations of motion of the Green s functions
Phys60.nb 0 8.3. Feynman diagram 8.3.. multi-particle Green s functions Single particle Green s function: G, ' ' (8.7) Two-particle Green s function: G, ; ', ' ' ' (8.73) three-particle Green s function:
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More information1 Quantum field theory and Green s function
1 Quantum field theory and Green s function Condensed matter physics studies systems with large numbers of identical particles (e.g. electrons, phonons, photons) at finite temperature. Quantum field theory
More 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 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 informationOrbital-dependent backflow transformations in quantum Monte Carlo
transformations in quantum Monte Carlo P. Seth, P. López Ríos, and R. J. Needs TCM group, Cavendish Laboratory, University of Cambridge 5 December 2012 VMC and DMC Optimization Wave functions Variational
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationTime-dependent linear-response variational Monte Carlo.
Time-dependent linear-response variational Monte Carlo. Bastien Mussard bastien.mussard@colorado.edu https://mussard.github.io/ Julien Toulouse julien.toulouse@upmc.fr Sorbonne University, Paris (web)
More informationQuantum Monte Carlo Methods for Strongly Correlated Electron Systems
This is page 1 Printer: Opaque this Quantum Monte Carlo Methods for Strongly Correlated Electron Systems Shiwei Zhang ABSTRACT We review some of the recent development in quantum Monte Carlo (QMC) methods
More informationAll-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe
All-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe A. Ma, N. D. Drummond, M. D. Towler, and R. J. Needs Theory of Condensed Matter Group, Cavendish Laboratory, University of
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 informationRandom Walks A&T and F&S 3.1.2
Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will
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 informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationLow- and High-Energy Excitations in the Unitary Fermi Gas
Low- and High-Energy Excitations in the Unitary Fermi Gas Introduction / Motivation Homogeneous Gas Momentum Distribution Quasi-Particle Spectrum Low Energy Excitations and Static Structure Function Inhomogeneous
More 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 informationQuantum Theory of Many-Particle Systems, Phys. 540
Quantum Theory of Many-Particle Systems, Phys. 540 Questions about organization Second quantization Questions about last class? Comments? Similar strategy N-particles Consider Two-body operators in Fock
More informationQUANTUM MONTE CARLO METHOD: PERFORMANCE ANALYSIS ON MODEL SYSTEMS
Vol. 79 (1991) ACTA PHYSICA POLONICA A No 6 ΟΝ THE FIXED-NODE IMPORTANCE-SAMPLING QUANTUM MONTE CARLO METHOD: PERFORMANCE ANALYSIS ON MODEL SYSTEMS A. CHŁOBOWSKI Κ. Gumiński Department of Theoretical Chemistry,
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More 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 informationIntroduction to Theory of Mesoscopic Systems
Introduction to Theory of Mesoscopic Systems Boris Altshuler Princeton University, Columbia University & NEC Laboratories America Lecture 3 Beforehand Weak Localization and Mesoscopic Fluctuations Today
More informationReleased-phase quantum Monte Carlo method
PHYSICAL REVIEW E VOLUME 55, NUMBER 5 MAY 1997 Released-phase quantum Monte Carlo method Matthew D. Jones, 1,2, * Gerardo Ortiz, 2, * and David M. Ceperley 1,2 1 National Center for Supercomputing Applications,
More informationComputational Physics Lectures: Variational Monte Carlo methods
Computational Physics Lectures: Variational Monte Carlo methods Morten Hjorth-Jensen 1,2 Department of Physics, University of Oslo 1 Department of Physics and Astronomy and National Superconducting Cyclotron
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 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 information