Clas s ical and path integral Monte Carlo s imulation of charged particles in traps
|
|
- Henry Ford
- 5 years ago
- Views:
Transcription
1 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 B aumg artner Electron energy Christian-Albrechts-Universitaet zu Kiel, Institut für Theoretische Physik und Astrophysik, Institute of Physics, Rostock University Institute for High Energy Density RAS, Russia Plasma crystal (classical particles) CB Ee ehω Eg VB Eh h+ L z Electrons in nanostructures (quantum particles)
2 Contents Introduction: C las s ical MC for trapped particles G eneral idea of PIMC s imulations Quantum s tatis tics. PIMC for bos ons and fermions Practical is s ues. Applications Examples : Many-body correlation effects of excitonic s ys tems in s emiconductors
3 Idea of classical thermodynamic Monte Carlo Goal: - Obtain exact thermodynamic equilibrium configuration R of interacting particles at given temperature, particle number, external fields etc. R= r 1, r 2,..., r N Coordinates of all particles (microstate) - Evaluate measurable quantities, such as energy, potential energy V, pressure, pair distribution function etc. 1 V, N = dr V R e V R, Z = 1, k BT - Averaging with canonical probability distribution (Boltzmann factor) Realization: - try all possible configurations ( move particles) and choose the most probable one - Metropolis Monte Carlo procedure (Markov process) - compute averages from fluctuating microstates --> See examples
4 Observables for Quantum Systems Quantum mechanical average using N-particle wave function (pure ensemble) V = 1 dr V R R R, Z Direct generalization to coherent superposition of N-particle wave functions Finite temp. (T>0, mixed ensemble), use density matrix instead of wave function V = 1 dr V R R, R ;, Z R, R ; Diagonal N-particle density matrix. Monte Carlo calculations like in classical case - are possible since the density matrix is positive on the diagonal (true probability density) But: for interacting systems DM ρ is unknown. --> We can sample particle's coordinates R only after we find an analytical representation for the DM.
5 Side remark: representations of quantum mechanics There exist several equivalent formulation of quantum mechanics: Schrödinger picture: time-dependent wave functions Heisenberg picture: time-dependent operators (matrices) Feynman: path integral representation - describes time evolution (wave function) or thermodynamic state ( imaginary time ) via density operator - analytical solutions only for few models - numerical solutions: computationally very demanding breakthrough in connection with Monte Carlo methods since the 1970s
6 Properties of the density matrix (1) Operator notation (thermodynamic equilibrium) ] =exp [ H t =exp[ i H t /ℏ] Original idea (Feynman): solution of Schrödinger equation U via time evolution operator --> path integral representation of quantum mechanics System Hamiltonian includes N-particle kinetic energy and potential energy operators (not commuting!) = K V, H [ K, V ] 0 Coordinate representation: superposition of wavefunctions weighted with probability density (with N-particle energy eigenvalues) R R ' = R, R' ; = R R ' e E Density operator known only for few cases, Example: high temperature -> classical system with weak interaction Idea (Feynman): express Density operator of interacting low temperature system by known high-temperature limit --> gives rise to path integral
7 Properties of the density matrix (2) Convolution property of density matrix: low-temperature DM related to DM at higher temperature H R ', R, R' ; 1 2 = R e 1 dr1 R e 1 H 2 H R1 R1 e 2 1 2= 1 k B T 1 T 2 R' = dr1 R, R1 ; 1 R1, R ' ; 2 1 (complete set) Generalization: use convolution property M times R R' = R exp[ H ]... exp[ H ] R ', = / M ]... exp[ H ] R ', R R' = R exp[ H R, R' ; = dr1 dr2...dr M 1 R, R1 ; R1, R2 ;... R M 2, R M 1 ; R M 1, R ' ; M factors connected by M-1 intermediate integrations. Note: total dimension of integral: (M-1) 3 N, may be very large number success of method relies on highly efficient Monte Carlo integration
8 Properties of the density matrix (3) Visualization of diagonal elements of density matrix: R, R ; = dr1 dr 2... dr M 1 R, R1 ; R1, R 2 ;... R M 2, R M 1 ; R M 1, R ; 2 arbitrary particles (e, h) represented by set ( chain ) of points ( beads ) this chain is closed. It is called path (therefore: path integral ) r2 λe rm-1 r r2 λh r1 r1 rm-1 re rh r All intermediate points do occur, We have to integrate over all possible intermediate points (paths)
9 High temperature approximation Obtain high-temperature expression of the DM, use Trotter's theorem (1959) T V T V M =lim M [e =e, [e T V M e ] e 2 M [T, V ]/2 ] O e [e T V M e ] O 1/ M Coordinate representation of the high temperature density matrix R i, R i 1 ; = Ri e T e V Ri 1 =exp 2 T = ℏ 2m V diagonal (local), [ 2 R R i 1 i V Ri 1 ; 2 DeBroglie wave length 2 r2 e λe rm-1 2 e, h =2 ℏ / me, h r2 λ r1 = r + λ ξ 1 r Resulting Particle extension~λ(β) λh V(R) re r1 rm-1 r rh h ]
10 Analytical form of the N-particle density matrix Density matrix has analogy to a classical system of interacting polymers (use action S) [ M ] R, R ; = V dr1 dr 2... dr M 1 exp S i, S i= 2 i=1 Ri 1 Ri 2 U Ri ''Spring'' term holding polymers (paths) together Full interaction energy of polymers (Example: electron-hole system) e h e h U Ri =U ee Ri ; U hh Ri ; U eh Ri, Ri ; Classical Coulomb interaction or effective regularized potential Averaging over many configurations yields smooth probability densities of interacting quantum systems which become exact with M -->
11 Example: two electrons in 2D parabolic trap r2 r1 rm-1 r Classical particles low density PIMC simulations of Jens Böning and Alexei Filinov medium density high density --> PIMC simulations yield probability density of correlated particles
12 Basic numerical issues of PIMC r2 How to sample the paths. It is necessary to explore the whole coordinate space for each intermediate point.this is very time consuming. To speed up convergence: move several slices (points of path) at once r1 rm-1 r Choose action as accurate as possible e.g. use effective interaction potentials which take into account two, three and higher order correlation effects. More accurate actions help to reduce the number of time slices by a factor of 10 or more. How to calculate physical properties. There are different approaches for calculating expectation values of physical observables, such as the energy, momentum distribution, etc. They are called estimators. In each particular case convergence can be improved by choosing the proper estimator.
13 General idea. Metropolis method In statistical mechanics thermodynamic averages are calculated as R' R ' R ] dr dr ' R O [ O =Tr O = [Tr ] R dr R O = dr dr ' O R, R' p R, R ' sign[ R, R ' ; ] p(r,r') normalized probability density (recover same form like in classical case) For diagonal operators: O R, R ' =O R R R' probability density given by p R ; = 1 R, R ; Z With DM obtained from high-temperature decomposition: R, R' ; = dr 1 dr 2... dr M 1 R, R1 ; R1, R 2 ;... R M 2, R M 1 ; R M 1, R ' ;
14 Numerical realization of PIMC (neglecting exchange) Sampling of configurations with probability p: use Metropolis method Configuration space consists of positions and spin projections of all N particles: Old configuration: New configuration: s= R, R1, R2,..., R M 1 ; 1,..., N s ' = R ', R' 1, R' 2,..., R ' M 1 ; ' 1,..., ' N Particle spin projections Construct sequence of uncorrelated configurations (Markov chain) Generated states (configurations) should have the correct probability distribution given by: [ s1, s 2,..., s L ] p R, = 1 R, R ; Z Realization in thermodynamic equilibrium (Metropolis): use ratio of probabilities: M 1 p s= exp[ S i Ri ; ], Z i=1 M 1 p s' = exp[ S i R ' i ; ], Z i=1 M ps ' T s, s ' = =exp[ S i s' ; S i s ; ] ps i=1 T(s',s) : transition probability used to construct Markov chain. Final result for physical observables in PIMC: average over all configurations L O =lim L 1 O si 1 / T si si 1 L i=1 L 1 1/T si si 1 L i=1, L: length of Markov chain. Averages of all quantities of interest are computed in one simulation
15 Generic canonical PIMC algorithm (1) Specify total number of time slices M. This affects the error of the high-temperature representation. (Try different M, check convergence with respect to M) (2) Randomly choose positions of N particles. (3) Set maximum number L of Monte Carlo steps. Initialize variables for thermodynamic averages and MC step counter (set to zero). (4) Set a ''menu'' of different Monte Carlo moves in the Metropolis algorithm. (5) Choose a particle i at random among N particles. (6) Proceed with the sampling of the path coordinates depending on the chosen possibility from the ''menu'' of Monte Carlo moves, e.g. (a ) Whole particle displacement (like for classical particle). (b ) Multilevel (bisection) deformation of particle trajectory. [c) Choose a new permutation for given particle, realize it with the bisection algorithm.] [d) possible additional moves] (7) Accept or reject the move using the acceptance probability. Change the system state if the move was accepted, otherwise stay in the old configuration. (8) Accumulate thermodynamic averages. (9) Increase MC step counter by one. If the maximum number is not reached proceed to (5). (10) Else: finish simulations. Output of final estimators for thermodynamic averages.
16 How to sample paths. Metropolis Monte Carlo that moves a single coordinate of the path is too slow, --> move more time slices together (multilevel method). but: path deformation of N quantum particles not independent due to quantum statistics (particles are either fermions or bosons) R t /2 Path deformation should simultaneously realize the N-particle exchange properties (symmetry/antisymmetry of N-particle wave function or density operator) Therefore: consider first extension of PIMC to bosons/fermions We will return to path sampling later
17 Quantum exchange. PIMC for fermions/bosons (1) 1. Properties of the system of N electrons and holes at a finite temperature T are determined by the density operator =exp [ H / k B T ] 2. Due to the Fermi/Bose statistics the total density matrix should be (anti)symmetric under arbitrary exchange of identical particles (e.g. electrons, holes, e with same spin projection etc.). A/ S for fermions/bosons We have to replace 3. Construct particles. A/S : as superposition of all N! permutations of N identical Example: two types (e,h) of particles of with numbers Ne, Nh A/S R e, R h, Re, Rh ; = 1 P P 1 1 Re, Rh, P e Re, P h R h ; P P N e! N h! e e h h P parity of permutation (number of equivalent pair transpositions) bosons: all terms have positive sign, fermions: alternating sign
18 Quantum exchange. PIMC for fermions/bosons (2) Permutation operators --> exchange path ends Illustration: 2 electrons and 2 holes (fermions) A R e, Rh, Re, R h ; = 1 P P e Re, P h R h ; 1 1 R, R, P e h 2! 2! P P e e System of quantum particles h h Classical system of interacting polymers Exact treatment of many-body interactions PIMC: mapping here all paths are closed --> identity permutation must add 2-particle exchanges N e! N h!=4 Number of possible permutations: e 1 e 1 Ph h Pe h
19 Quantum exchange. PIMC for fermions/bosons (3) Initial configuration (a): Particle 1 and 2 are in identy permutations; 3,4,5 form a permutation of length l=3 (cycle 4-3-5) Final configuration (b): Particle 2 is in identy permutation; Particles 1 and 4 are exchanged; 1,3,4,5 form a new exchange cycle of length l=4 (cycle ).
20 How to sample paths. Bisection algorithm. Metropolis Monte Carlo that moves a single coordinate of the path is too slow, --> move more time slices together (multilevel method). R /2 Key point: sample a path using mid-points R / 2 Guiding rule to sample mid-points Rt / 2 : P R /2 = R e H /2 R /2 R /2 e H /2 R ' R e H R ' [ exp R R 0 /2 /2 2 d /2 2 /2 ] Gaussian with R R' R0 =, 2 Exact for ideal systems ℏ2 = 2m Construct entire path by recursion from this formula (bisection). Possible improved sampling: use optimized mean R 0 and variance which account for local interaction strength (nearest neighbor interaction). U R0 R R ' R0=, 2 R ℏ ℏ 2 = U R0 2m m 2
21 Construction of new paths
22 How to sample paths. Bisection Construction of new pathsalgorithm. Select time slices. 4
23 How to sample paths. Bisection Construction of new pathsalgorithm. Select time slices. Select permutation from possible ones (pairs, triplets, etc.), using ' ; 4 p R, PR 4
24 How to sample paths. Bisection Construction of new pathsalgorithm. Select time slices. Select permutation from possible ones (pairs, triplets, etc.), using ' ; 4 p R, PR Sample mid-points. 4 2
25 How to sample paths. Bisection Construction of new pathsalgorithm. Select time slices. Select permutation from possible ones (pairs, triplets, etc.), using ' ; 4 p R, PR Sample mid-points. Bisect again, up to lowest level. 4
26 How to sample paths. Bisection Construction of new pathsalgorithm. Select time slices. Select permutation from possible ones (pairs, triplets, etc.), using ' ; 4 p R, PR 4 Sample mid-points. Bisect again, up to lowest level. Accept or reject entire move.
27 How to sample paths. algorithm. Advantages of thebisection bisection method Detailed balance property is satisfied at each level. We do not waste time on moves for which paths come close and the poten tial energy strongly increases (for repulsive interaction). Such configurati ons are rejected already on early steps. Computer time is spent more efficiently because we consider mainly configurations where the acceptance rate is high. Sampling of particle permutations is easy to perform in the bisection method. Old configuration of 5 particles New configuration with 2 particles in a pair exchange
28 Indirect excitons in an electrostatic trap T w o c o u p le d e -h la y e r s Ez exciton electric field e h ZnMgSe/ZnSe-quantum well Dipole moment in L=30nm er= 8.7 wide QW mel=0.15m0, at field E=20 kv / cm m =0.86m, hh d =20.4 nm 6.65 a ZnSe B Stronger dipole interaction of excitons 0 mhh=0.37m0 1aB= 3.07nm 1Ha=53.9meV GaAs/AlGaAS-quantum well N = N e= N h We consider the system with equal number electron-hole pair (optically excited) N =N e= N h er=12.58 mel=0.067m0 mh=0.34m0 1aB = 9.98nm 1Ha=11.5meV The density is controlled by the trap frequency Single particle Hamiltonian m e 2e =m h 2h and characterized by coupling parameter = e 2 / l 0 /ℏ e At low temperature and small d, electrons and holes form bound state: exciton (boson)
29 Quantum crystal of Excitons in an electrostatic trap e 2 / ε l0 λ =, l0 = ω mrω Vary trap frequency and Coupling parameter by changing field strength Structural transitions: Classical exciton liquid Exciton crystal Exciton Bose fluid PIMC simulations of Patrick Ludwig and Alexei Filinov
30 Quantum crystal of Excitons in an electrostatic trap e 2 / ε l0 λ =, l0 = ω mrω Vary trap frequency and Coupling parameter by changing field strength Structural transitions: Classical exciton liquid Exciton crystal Exciton Bose fluid PIMC simulations of Patrick Ludwig and Alexei Filinov
31 Permutation probability. Bose condensation of excitons 1. Spin statistics: Fully spin polarized electrons and holes (to enhance quantum statistical effects in small systems) 2. Classical system: only identity permutation 3. Ideal quantum gas at T=0: equal probability of all permutations =1 5 NB Using a threshold value P cr =0.8/ N x we obtain the Bose condensate fraction according to B =N B / N x At low temperatures reaches 80%. N u m b e r o f p a r tic le s in th e p e r m u ta tio n S n a p s h o ts o f p e r m u ta tio n s in th e e le c tr o n la y e r (3 0 p a r tic le s ) in 2 D h a r m o n ic tr a p T = 208 mk T = K T = 1.25 K
32 Permutation probability. Influence of particle number and interaction. C o u lo m b c o u p lin g p a r a m e t e r = e / l / ℏ 2 0 e
33 Superfluidity Loss of viscosity below critical temperature - discovered in liquid He by P.L. Kapitza L.D. Landau: two-fluid (normal+suprafluid) model Andronikashvili experiment Superfluid density (fraction) from Classical/quantum moment of inertia I I quant ρs = 1 ρ I class First principle result for interacting bosons from path area A: Quantum Monte Carlo simulations (D. Ceperley, 1995 ) Moment of inertia computed in PIMC fromarea enclosed by pathe
34 Superfluid fraction of excitons Path integral Monte Carlo results for strongly correlated excitons ZnSe 30nm quantum well, E=20kV/cm Coupling parameter e 2 / ε l0 λ =, l0 = ω mrω A.Filinov, M. Bonitz, P. Ludwig, and Yu. Lozovik, phys. stat. sol. (c) 3, 2457 (2006)
35 General aspects of PIMC Alternative strategies Limitations of the simulations Scope of applications
36 Thermodynamic averages. Estimators. Quantities of interest: Energy, pressure (equation of state), specific heat, fluc tuations, condensate or superfluid fraction, pair distribution function... In principle: all quantities follow from density matrix or partition function Z Limitations Example: thermodynamic estimator for the internal energy 1. We usually calculate only ratios of integrals, e.g., free energy and entropy require special techniques. Direct estimator 2. The variance of some estim ators can be too high. 3. Many quantities are defined as dynamical quantities, but we are limited to only imaginary time (static quantities). One of the forms of virial estimator
37 Fermion sign problem. Concepts of fermionic PIMC (1) 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 R ; ±1 P R, P N! P Bose systems: all permutations contribute with the same (positive) sign Fermions: essential cancellation of positive and negative terms (corresponding to even and odd permutations), both are close in their absolute value. Accurate calculation of this small difference is drastically hampered with the in crease of quantum degeneracy (low T, high density).
38 Fermion sign problem. Concepts of fermionic PIMC (2) A R, R ; = 1 R ; 1 P R, P N! P Fermions: essential cancellation of positive and negative terms (corresponding to even and odd permutations), both are close in their absolute value. Accurate calculation of this small difference is drastically hampered with the in crease of quantum degeneracy (low T, high density). Used approaches to overcome Fermion problem in PIMC: (a ) Fixed-node (fixed-phase) approximation. Idea: 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).
39 Fermion sign problem. Concepts of fermionic PIMC (3) (b ) Direct PIMC. Idea: Do not sample individual permutations in the sum. Instead: use the full ex pression presented in a form of an N N determinant. In this case the absoute 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). (c ) Multilevel-blocking PIMC. Idea: 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 lower temperature. Further use it in the sampling probabilities of path coordinates on the next level (corresponding to the density matrix at even 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).
40 How to sampleof paths. Bisection algorithm. Application PIMC to quantum charged particles General fields of use low temperature systems (relevance of quantum effects) Small dimensions (system size comparable to DeBroglie wavelength λ high density: 2-particle distance comparable to / smaller than λ 1. Single particle in complicated potential e.g. disorder effects --> effective solution of Schrödinger's equation (ground state, T=0) or finite temperature extension (density matrix) 2. Two-particle interaction improvement of pair potential at small distances (include quantum effects) 3. Finite number of particles in traps atoms, ions at ultralow temperature (Bose condensates etc.) electrons, holes in quantum dots 4. Macroscopic quantum systems electrons and ions in astrophysics: planet cores, dwarf stars, highly excited solids (many electrons, holes in nanostructures)
41 How to sampleof paths. Bisection algorithm. Application PIMC to quantum charged particles (1) low temperature systems (relevance of quantum effects) Small dimensions (system size comparable to DeBroglie wavelength λ high density: 2-particle distance comparable to / smaller than λ V Ti / Al AlGaAs X X + h - GaAs - e AlGaAs 1. Single particle or bound complexes in complicated potential (e.g. disorder effects) --> effective solution of Schrödinger's equation (ground state, T=0) or finite temperature extension (density matrix) Example:1 realistic semiconductor heterostructure (quantum well) with well fluctuations n-gaas 0V Example:2 electrons, holes, excitons in strong external E-field
42 Indirect excitons in a Quantum well in electric field Ez CB Ee ehω Eg Eh VB h+ L z Electron energy Electron energy CB e- Ee Eg VB hω Eh h+ L z Shift of energy bands in electric field: E z =e E z z Change of quantum probability distribution of electrons and holes
43 Exact Quantum pair potentials 1. Exact pair potential from exact 2-particle density matrix (numer.) U ab pair (r ) = k BT ln ρ ab (r, r ) 2. Derive analytical potential (1 fit parameter) Result: Drastic improvement of previous potentials (Kelbg, Deutsch etc.) applicable to strong coupling including bound states A.Filinov, M. Bonitz, W. Ebeling, J. Phys. A 36, 5899 (2003) U, units of Ryd T=30,000K Hydrogen Yields exact thermodynamic properties (on 2-particle level)
44 Summary We have presented an introduction to first principle path integral Monte Carlo simulations (PIMC) and applications Virtually exact method for thermodynamics of quantum systems at any interaction strength. Limited only by fermion sign problem (min temperature) There are many technical and simulation details: they will be discussed in the lab course on Thursday by Alexei Filinov Proceedings: basics of PIMC simulation of trapped particles (brief) A detailed description (110 pages) is given in our book Introduction to Computational methods for Many-body physics Rinton Press, Princeton 2006 (on the table) it includes chapters on - Particle in cell simulations, by H. Ruhl - Density functional theory, by G. Bertsch - Quantum Kinetic equations, by M. Bonitz and D. Semkat - Classcial and Path integral Monte Carlo, by A. Filinov and M. Bonitz - Quantum MD, by A. Filinov, V. Filinov, Yu. Lozovik and M. Bonitz
Introduction to Path Integral Monte Carlo. Part III.
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-24098 Kiel, Germany November
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 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 informationClassical Monte-Carlo simulations
Classical Monte-Carlo simulations Graduate Summer Institute Complex Plasmas at the Stevens Insitute of Technology Henning Baumgartner, A. Filinov, H. Kählert, P. Ludwig and M. Bonitz Christian-Albrechts-University
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 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 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 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 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 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 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. Permutation sampling considerations 3. Calculation of superfluid density and momentum distribution. 4. Applications of PIMC to liquid helium
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 informationLandau Bogolubov Energy Spectrum of Superconductors
Landau Bogolubov Energy Spectrum of Superconductors L.N. Tsintsadze 1 and N.L. Tsintsadze 1,2 1. Department of Plasma Physics, E. Andronikashvili Institute of Physics, Tbilisi 0128, Georgia 2. Faculty
More informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
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 informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
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 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: 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 informationStructure Formation in Strongly Correlated Coulomb Systems
Structure Formation in Strongly Correlated Coulomb Systems M. Bonitz 1, V. Filinov, A. Filinov 1, V. Golubnychiy 1, P. Ludwig 1,3, H. Baumgartner 1, P.R. Levashov, V.E. Fortov, and H. Fehske 4 1 Institut
More informationarxiv: v1 [cond-mat.quant-gas] 10 Jan 2011
Crystallization of an exciton superfluid J. öning, 1 A. Filinov, 1 and M. onitz 1, 1 Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität, Leibnizstr. 15, D-24098 Kiel, Germany
More informationBCS-BEC Crossover. Hauptseminar: Physik der kalten Gase Robin Wanke
BCS-BEC Crossover Hauptseminar: Physik der kalten Gase Robin Wanke Outline Motivation Cold fermions BCS-Theory Gap equation Feshbach resonance Pairing BEC of molecules BCS-BEC-crossover Conclusion 2 Motivation
More 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 informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
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 informationThe Uniform Electron Gas at Warm Dense Matter Conditions
The Uniform Electron Gas at Warm Dense Matter Conditions Tobias Dornheim, Simon Groth, and Michael Bonitz, Physics Reports 744, 1-86 (2018) Institute of Theoretical Physics and Astrophysics Kiel University
More informationVariational wave function for a two-electron quantum dot
Physica B 255 (1998) 145 149 Variational wave function for a two-electron quantum dot A. Harju *, V.A. Sverdlov, B. Barbiellini, R.M. Nieminen Laboratory of Physics, Helsinki University of Technology,
More informationIn-class exercises. Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 8 Exercises due Mon March 19 Last correction at March 5, 2018, 8:48 am c 2017, James Sethna,
More informationSupersolids. Bose-Einstein Condensation in Quantum Solids Does it really exist?? W. J. Mullin
Supersolids Bose-Einstein Condensation in Quantum Solids Does it really exist?? W. J. Mullin This is a lively controversy in condensed matter physics. Experiment says yes. Theory says no, or at best maybe.
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 information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
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 information1) K. Huang, Introduction to Statistical Physics, CRC Press, 2001.
Chapter 1 Introduction 1.1 Literature 1) K. Huang, Introduction to Statistical Physics, CRC Press, 2001. 2) E. M. Lifschitz and L. P. Pitajewski, Statistical Physics, London, Landau Lifschitz Band 5. 3)
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
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 informationElectron spins in nonmagnetic semiconductors
Electron spins in nonmagnetic semiconductors Yuichiro K. Kato Institute of Engineering Innovation, The University of Tokyo Physics of non-interacting spins Optical spin injection and detection Spin manipulation
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 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 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 informationComplexity of the quantum adiabatic algorithm
Complexity of the quantum adiabatic algorithm Peter Young e-mail:peter@physics.ucsc.edu Collaborators: S. Knysh and V. N. Smelyanskiy Colloquium at Princeton, September 24, 2009 p.1 Introduction What is
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 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 informationQuantum Properties of Two-dimensional Helium Systems
Quantum Properties of Two-dimensional Helium Systems Hiroshi Fukuyama Department of Physics, Univ. of Tokyo 1. Quantum Gases and Liquids 2. Bose-Einstein Condensation 3. Superfluidity of Liquid 4 He 4.
More informationComputational Physics (6810): Session 13
Computational Physics (6810): Session 13 Dick Furnstahl Nuclear Theory Group OSU Physics Department April 14, 2017 6810 Endgame Various recaps and followups Random stuff (like RNGs :) Session 13 stuff
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 informationUltracold Fermi and Bose Gases and Spinless Bose Charged Sound Particles
October, 011 PROGRESS IN PHYSICS olume 4 Ultracold Fermi Bose Gases Spinless Bose Charged Sound Particles ahan N. Minasyan alentin N. Samoylov Scientific Center of Applied Research, JINR, Dubna, 141980,
More informationMesoscopic Nano-Electro-Mechanics of Shuttle Systems
* Mesoscopic Nano-Electro-Mechanics of Shuttle Systems Robert Shekhter University of Gothenburg, Sweden Lecture1: Mechanically assisted single-electronics Lecture2: Quantum coherent nano-electro-mechanics
More informationAyan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur
Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur Under the supervision of: Dr. Marcel Nooijen Associate Professor Department of Chemistry University of Waterloo
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
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 informationThermodynamics of nuclei in thermal contact
Thermodynamics of nuclei in thermal contact Karl-Heinz Schmidt, Beatriz Jurado CENBG, CNRS/IN2P3, Chemin du Solarium B.P. 120, 33175 Gradignan, France Abstract: The behaviour of a di-nuclear system in
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 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 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 informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More informationLab 70 in TFFM08. Curie & Ising
IFM The Department of Physics, Chemistry and Biology Lab 70 in TFFM08 Curie & Ising NAME PERS. -NUMBER DATE APPROVED Rev Aug 09 Agne 1 Introduction Magnetic materials are all around us, and understanding
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
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 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 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 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 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 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 informationThermodynamic and Transport Properties of Dense Quantum Coulomb Systems
John von Neumann Institute for Computing Thermodynamic and Transport Properties of Dense Quantum Coulomb Systems M. Bonitz, V.S. Filinov published in NIC Symposium 2004, Proceedings, Dietrich Wolf, Gernot
More informationPrinciples of Equilibrium Statistical Mechanics
Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC
More informationPhase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality. Hans-Henning Klauss. Institut für Festkörperphysik TU Dresden
Phase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality Hans-Henning Klauss Institut für Festkörperphysik TU Dresden 1 References [1] Stephen Blundell, Magnetism in Condensed
More informationJ ij S i S j B i S i (1)
LECTURE 18 The Ising Model (References: Kerson Huang, Statistical Mechanics, Wiley and Sons (1963) and Colin Thompson, Mathematical Statistical Mechanics, Princeton Univ. Press (1972)). One of the simplest
More informationQuantum Fluids and Solids. Christian Enss (4He) & Henri Godfrin (3He)
Quantum Fluids and Solids Christian Enss (4He) & Henri Godfrin (3He) What Do We Consider as Quantum Matter? φ Lennard Jones Potential σ ε r De Boer quantum parameter: De Broglie wavelength: Introduction
More informationinterband transitions in semiconductors M. Fox, Optical Properties of Solids, Oxford Master Series in Condensed Matter Physics
interband transitions in semiconductors M. Fox, Optical Properties of Solids, Oxford Master Series in Condensed Matter Physics interband transitions in quantum wells Atomic wavefunction of carriers in
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 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 informationBound states of two particles confined to parallel two-dimensional layers and interacting via dipole-dipole or dipole-charge laws
PHYSICAL REVIEW B VOLUME 55, NUMBER 8 15 FEBRUARY 1997-II Bound states of two particles confined to parallel two-dimensional layers and interacting via dipole-dipole or dipole-charge laws V. I. Yudson
More informationInternational Physics Course Entrance Examination Questions
International Physics Course Entrance Examination Questions (May 2010) Please answer the four questions from Problem 1 to Problem 4. You can use as many answer sheets you need. Your name, question numbers
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationIntroduction to Heisenberg model. Javier Junquera
Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics
More informationQuantum Theory of Matter
Quantum Theory of Matter Revision Lecture Derek Lee Imperial College London May 2006 Outline 1 Exam and Revision 2 Quantum Theory of Matter Microscopic theory 3 Summary Outline 1 Exam and Revision 2 Quantum
More 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 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 informationChapter 2 Ensemble Theory in Statistical Physics: Free Energy Potential
Chapter Ensemble Theory in Statistical Physics: Free Energy Potential Abstract In this chapter, we discuss the basic formalism of statistical physics Also, we consider in detail the concept of the free
More information9.1 System in contact with a heat reservoir
Chapter 9 Canonical ensemble 9. System in contact with a heat reservoir We consider a small system A characterized by E, V and N in thermal interaction with a heat reservoir A 2 characterized by E 2, V
More informationTheory and simulation of strong correlations in quantum Coulomb systems
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSA: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 36 (2003) 5921 5930 PII: S0305-4470(03)55274-X Theory and simulation of strong correlations in quantum
More informationH.O. [202] 3 2 (2) (2) H.O. 4.0 [200] 1 2 [202] 5 2 (2) (4) (2) 3.5 [211] 1 2 (2) (6) [211] 3 2 (2) 3.0 (2) [220] ε
E/ħω H r 0 r Y0 0 l s l l N + l + l s [0] 3 H.O. ε = 0.75 4.0 H.O. ε = 0 + l s + l [00] n z = 0 d 3/ 4 [0] 5 3.5 N = s / N n z d 5/ 6 [] n z = N lj [] 3 3.0.5 0.0 0.5 ε 0.5 0.75 [0] n z = interaction of
More informationSpontaneous Spin Polarization in Quantum Wires
Spontaneous Spin Polarization in Quantum Wires Julia S. Meyer The Ohio State University with A.D. Klironomos K.A. Matveev 1 Why ask this question at all GaAs/AlGaAs heterostucture 2D electron gas Quantum
More informationto satisfy the large number approximations, W W sys can be small.
Chapter 12. The canonical ensemble To discuss systems at constant T, we need to embed them with a diathermal wall in a heat bath. Note that only the system and bath need to be large for W tot and W bath
More informationSuperfluidity and Superconductivity
Superfluidity and Superconductivity These are related phenomena of flow without resistance, but in very different systems Superfluidity: flow of helium IV atoms in a liquid Superconductivity: flow of electron
More informationORIGINS. E.P. Wigner, Conference on Neutron Physics by Time of Flight, November 1956
ORIGINS E.P. Wigner, Conference on Neutron Physics by Time of Flight, November 1956 P.W. Anderson, Absence of Diffusion in Certain Random Lattices ; Phys.Rev., 1958, v.109, p.1492 L.D. Landau, Fermi-Liquid
More informationStatistical Physics of The Symmetric Group. Mobolaji Williams Harvard Physics Oral Qualifying Exam Dec. 12, 2016
Statistical Physics of The Symmetric Group Mobolaji Williams Harvard Physics Oral Qualifying Exam Dec. 12, 2016 1 Theoretical Physics of Living Systems Physics Particle Physics Condensed Matter Astrophysics
More informationThe Cooper Problem. Problem : A pair of electrons with an attractive interaction on top of an inert Fermi sea c c FS,
Jorge Duelsy Brief History Cooper pair and BCS Theory (1956-57) Richardson exact solution (1963). Gaudin magnet (1976). Proof of Integrability. CRS (1997). Recovery of the exact solution in applications
More 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 informationThe Nuclear Equation of State
The Nuclear Equation of State Abhishek Mukherjee University of Illinois at Urbana-Champaign Work done with : Vijay Pandharipande, Gordon Baym, Geoff Ravenhall, Jaime Morales and Bob Wiringa National Nuclear
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 informationCopyright 2001 University of Cambridge. Not to be quoted or copied without permission.
Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions.
1. Quantum Mechanics (Fall 2004) Two spin-half particles are in a state with total spin zero. Let ˆn a and ˆn b be unit vectors in two arbitrary directions. Calculate the expectation value of the product
More informationQuantum 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 informationAn Effective Potential Approach to Modelling 25nm MOSFET Devices S. Ahmed 1 C. Ringhofer 2 D. Vasileska 1
An Effective Potential Approach to Modelling 25nm MOSFET Devices S. Ahmed 1 C. Ringhofer 2 D. Vasileska 1 1 Department of Electrical Engineering, 2 Department of Mathematics, Arizona State University,
More informationSimulations with MM Force Fields. Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi
Simulations with MM Force Fields Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi Some slides taken with permission from Howard R. Mayne Department of Chemistry University of New Hampshire Walking
More informationI. Collective Behavior, From Particles to Fields
I. Collective Behavior, From Particles to Fields I.A Introduction The object of the first part of this course was to introduce the principles of statistical mechanics which provide a bridge between the
More information