Thermal Casimir Effect for Colloids at a Fluid Interface

Size: px
Start display at page:

Download "Thermal Casimir Effect for Colloids at a Fluid Interface"

Transcription

1 Thermal Casimir Effect for Colloids at a Fluid Interface Jef Wagner Ehsan Noruzifar Roya Zandi Department of Physics and Astronomy University of California, Riverside

2 Outline

3 : Capillary Wave Hamiltonian H int = σa + U int g

4 : Capillary Wave Hamiltonian H int = σa + U int g A = d 2 x 1+ φ 2

5 : Capillary Wave Hamiltonian H int = σa + U int g A = d 2 x 1+ φ 2

6 : Capillary Wave Hamiltonian H int = σ d 2 x 1+ φ 2 + Ug int Ug int = d 2 x ρg 2 φ2

7 : Capillary Wave Hamiltonian H int = σ d 2 x 1+ φ 2 + Ug int Ug int = d 2 x ρg 2 φ2

8 : Capillary Wave Hamiltonian [ ( ) ] H int [φ] =σ d 2 x 1+ φ λ 2 R 2 2 φ2 σ λ = ρg 3cm

9 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3

10 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 ( Hpart col 1[φ] = σ d 2 x Ω 1+ φ 2 ) d 2 x Ω eq

11 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 ( Hpart col 1[φ] = σ d 2 x Ω 1+ φ 2 ) d 2 x Ω eq

12 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 H col part 1[φ] = σ [ ( ) ( d 2 x 1+ φ d 2 x Ω Ω ) ] d 2 x Ω eq

13 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 H col part 1[φ] = σ [ ( ) ( d 2 x 1+ φ d 2 x Ω Ω ) ] d 2 x Ω eq

14 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 H col part 1[φ] = σ [ ( ) ] d 2 x 1+ φ Ω[φ] Ω

15 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 U col g = d 2 x ρg (ρ I d 3 x + ρ II V I 2 φ2 + mgh d 3 x V II ) g

16 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 U col g = d 2 x ρg 2 φ2 + mgh

17 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 H col part 3[φ] =σ I A I + σ II A II

18 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 H col part 3[φ] =(σ I σ II ) A I

19 Colloid Contribution H col = σ da Ω part 1 + U col g part 2 + σ I A I + σ II A II part 3 Young s relation σ II σ I + σ cos θ =0 H col part 3[φ] =σ cos θ A I

20 Total Hamiltonian H tot [φ] =σ ( d 2 x R 2 / Ω i ( ( ) 1+ φ λ 2 + i 2 φ2 ( Ω i + cos θ A ii + m ) ) igh σ ) Expanding for spherical colloids H tot [φ] σ 2 ( d 2 x φ ( 2 + λ 2) φ + R 2 / Ω i i R s δω i dx φ 2 )

21 Proposed Experiment Silver coated hollow glass mircrospheres (10µm R s 100µm). Air Water interface. Potential decays as r 8. Potential has depth of 1k B T at sep 1µm.

22 Scattering Method for Casimir Physics (MIT) Z = Dφe 1 2 φ, ˆDφ R 2 \Ω1 Ω 2 C

23 Scattering Method for Casimir Physics (MIT) Z = Dφe 1 2 φ, ˆDφ R 2 \Ω1 Ω 2 C Delta functional constraints Dφ = C Dφδ δω1 [φ]δ δω2 [φ]

24 Scattering Method for Casimir Physics (MIT) Z = Dφδ δω1 [φ]δ δω2 [φ]e 1 2 φ, ˆDφ R 2 \Ω1 Ω 2 Fourier representation of the delta functional δ δωi [φ] = ψ i e ı ψ i,φ δωi

25 Scattering Method for Casimir Physics (MIT) Z = Dψ 1 Dψ 2 Dφ e 1 2 φ, ˆDφ R 2 \Ω1 Ω +i 2 i φ,ψ i δωi Gaussian functional integral Df e 1 2 f,af = ( det A ) 1 2

26 Scattering Method for Casimir Physics (MIT) Z = Dψ 1 Dψ 2 Dφ e 1 2 φ, ˆDφ R 2 \Ω1 Ω +i 2 i φ,ψ i δωi Gaussian functional integral Df e 1 2 f,af = ( det A ) 1 2

27 Scattering Method for Casimir Physics (MIT) Z = Dψ 1 Dψ 2 e 1 2 i,j ψ i,g 0 ψ j δωi,δω j Complete basis set ψ i = α Ψ iα φ inc iα Dφ i = α dψ iα

28 Scattering Method for Casimir Physics (MIT) { Z = dψ 1α dψ 2α exp 1 2 α ( Ψ1 Ψ 2 ) ) T (Ψ1 )} (Ĝ11 Ĝ 12 Ĝ 21 Ĝ 22 Ψ 2 Reminder: Free Energy βf = ln Z eff

29 Scattering Method for Casimir Physics (MIT) βf = 1 2 ln det ( 1 Ĝ 11 1 ) 1 Ĝ12Ĝ21 Ĝ22 Scattering Properties Ĝ 1 ii φ inc iα + β = T i T i αβ φsct iβ =0 Ĝ ij = U ij

30 Scattering Method for Casimir Physics (MIT) βf = 1 2 ln det ( 1 T 1 U 12 T 2 U 21) Questions: Why is the inverse of the ˆD operator on the space function space L 2 (R \ Ω 1 Ω 2 ) the free Green s function for R 2? How do I treat fluctuating boundaries?

31 Z = ( ) { Dφ Df i δ δωi [φ f i ] exp β 2 φ, ˆDφ β 2 i } f i, Ĥcol f i i Follow previous steps: Expand the delta functions. Integrate out the φ degree of freedom. Write the ψ i and f i fields in a complete basis set {φ iα }. Use the scattering properties.

32 Z = { DV exp 1 } 2 VT MV Vector and Matrix forms Ψ 1 V = P 1 Ψ 2 M = P 2 (T 1 ) 1 ıi U 12 0 ıi H col 0 0 U 21 0 (T 2 ) 1 ıi 0 0 ıi H col

33 Reminder: Free Energy βf = ln Z eff ( (T) 1 ıi ıi H col ) 1 ( ) U 0 = 0 0 ) ( TU 0 X 0 T = T T ( T + H col) 1 T

34 βf = 1 2 ln det ( 1 T 1 U 12 T2 U 21) The Casimir energy is still written in terms of scattering parameters. The term in the Hamiltonian due to fluctuating boundaries changes the scattering matrix.

35 Non-homogeneous Basis Functions Z = ( ) { Dφ δ Ωi [φ] exp β } φ, ˆDφ 2 i Reminder: Question Why is the inverse operator of ˆD on the function space L 2 (R 2 \ Ω 1 Ω 2 ) the free Green s function for R 2? Define the delta functional over the interior. Proceed in a similar fashion as before.

36 Non-homogeneous Basis Functions { Z = Dψ i exp β } ψ i, G 0 ψ j 2 i i,j New orthonormal basis set Let { φ iα } be a complete basis set, and partition the set into the homogeneous solutions to the operator ˆD and the rest. such that { φ iα } = { φ inc iα, ˆφ ia } ˆD φ inc iα =0 ˆD ˆφ ia 0

37 Non-homogeneous Basis Functions { Z = Dψ i exp β } ψ i, G 0 ψ j 2 i New orthonormal basis set The auxiliary field is expanded in the complete basis set i,j ψ i = α Ψ iα φ inc iα + a Ψ ia ˆφ ia Dψ i = α,a dψ iα d Ψ ia

38 Non-homogeneous Basis Functions Z = { DV exp 1 } 2 VT MV Vector and Matrix forms Ψ 1 G αβ G αb U 12 0 V = Ψ 1 Ψ 2 M = G aβ G ab 0 0 U 21 0 G αβ G αb Ψ G aβ G ab

39 Non-homogeneous Basis Functions Useful definition g(x) = d 2 x G 0 (x, x )ψ(x ) Homogeneous and particular solutions Ω ˆDg(x) = { 0 ψ(x) g(x) =g h (x)+g p (x) g h (x) = h α φ α (x) α g p (x) = g α φ α (x)+ α a ĝ a ˆφ a (x)

40 Non-homogeneous Basis Functions Matrix coefficients ˆDg p(x) =ψ(x) = g a = M 1 aβ Ψ β + M 1 Ψ ab b β b g p(x) 0 = g α = δω= S 1 αb M 1 bγ Ψ γ bγ bc Matrix Definitions S 1 αb M 1 bc M { } { α a }b = φ { α a }, ˆD ˆφ b ( M 1 M 1 ) ( ) M = I M S aβ = ˆφ a,φ β δω Ψ c

41 Non-homogeneous Basis Functions Matrix coefficients ( ) ( Gαβ G αb (T) = 1 S 1 M 1 S 1 M 1 ) G aβ G ab Block-wise Inverse ( A B ) 1 = M 1 M 1 C D X ( (A BD 1 C) 1 X ) X

42 Conclusions The thermal Casimir force appears measurable. Incorporated fluctuating boundaries into the scattering method. Is there a better proof that the non-homogenous basis functions do not contribute?

43 Conclusions The thermal Casimir force appears measurable. Incorporated fluctuating boundaries into the scattering method. Is there a better proof that the non-homogenous basis functions do not contribute?

44 Conclusions The thermal Casimir force appears measurable. Incorporated fluctuating boundaries into the scattering method. Is there a better proof that the non-homogenous basis functions do not contribute?

Chemistry 432 Problem Set 4 Spring 2018 Solutions

Chemistry 432 Problem Set 4 Spring 2018 Solutions Chemistry 4 Problem Set 4 Spring 18 Solutions 1. V I II III a b c A one-dimensional particle of mass m is confined to move under the influence of the potential x a V V (x) = a < x b b x c elsewhere and

More information

II.D Scattering and Fluctuations

II.D Scattering and Fluctuations II.D Scattering and Fluctuations In addition to bulk thermodynamic experiments, scattering measurements can be used to probe microscopic fluctuations at length scales of the order of the probe wavelength

More information

MULTIVARIABLE INTEGRATION

MULTIVARIABLE INTEGRATION MULTIVARIABLE INTEGRATION (SPHERICAL POLAR COORDINATES) Question 1 a) Determine with the aid of a diagram an expression for the volume element in r, θ, ϕ. spherical polar coordinates, ( ) [You may not

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION In the format provided by the authors and unedited. NATURE PHYSICS www.nature.com/naturephysics 1 285 x 5 15 5 25 1 1 Ω/µ 3 x C c = 4 6.5 3 µ 2 ε r =7 1.4 9 11 2 1 54 12 1 2 4 12 2 12 µ 2 12 12 12 12 NATURE

More information

Mean field theories of quantum spin glasses

Mean field theories of quantum spin glasses Mean field theories of quantum spin glasses Antoine Georges Olivier Parcollet Nick Read Subir Sachdev Jinwu Ye Talk online: Sachdev Classical Sherrington-Kirkpatrick model H = JS S i j ij i j J ij : a

More information

2 Feynman rules, decay widths and cross sections

2 Feynman rules, decay widths and cross sections 2 Feynman rules, decay widths and cross sections 2.1 Feynman rules Normalization In non-relativistic quantum mechanics, wave functions of free particles are normalized so that there is one particle in

More information

Supporting Information

Supporting Information Supporting Information A: Calculation of radial distribution functions To get an effective propagator in one dimension, we first transform 1) into spherical coordinates: x a = ρ sin θ cos φ, y = ρ sin

More information

Landau s Fermi Liquid Theory

Landau s Fermi Liquid Theory Thors Hans Hansson Stockholm University Outline 1 Fermi Liquids Why, What, and How? Why Fermi liquids? What is a Fermi liquids? Fermi Liquids How? 2 Landau s Phenomenological Approach The free Fermi gas

More information

A Numerical Study of Boson Star Binaries

A Numerical Study of Boson Star Binaries A Numerical Study of Boson Star Binaries Bruno C. Mundim with Matthew W. Choptuik (UBC) 12th Eastern Gravity Meeting Center for Computational Relativity and Gravitation Rochester Institute of Technology

More information

Introduction to Phase Transitions in Statistical Physics and Field Theory

Introduction to Phase Transitions in Statistical Physics and Field Theory Introduction to Phase Transitions in Statistical Physics and Field Theory Motivation Basic Concepts and Facts about Phase Transitions: Phase Transitions in Fluids and Magnets Thermodynamics and Statistical

More information

Investigations on the SYK Model and its Dual Space-Time

Investigations on the SYK Model and its Dual Space-Time 2nd Mandelstam Theoretical Physics Workshop Investigations on the SYK Model and its Dual Space-Time Kenta Suzuki A. Jevicki, KS, & J. Yoon; 1603.06246 [hep-th] A. Jevicki, & KS; 1608.07567 [hep-th] S.

More information

Stochastic Calculus. Kevin Sinclair. August 2, 2016

Stochastic Calculus. Kevin Sinclair. August 2, 2016 Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed

More information

Gauge-Higgs Unification and the LHC

Gauge-Higgs Unification and the LHC Gauge-Higgs Unification and the LHC If the Higgs boson is 124 or 126 or? GeV with SM couplings, Explain SM Higgs. with non-sm couplings, is not seen at LHC, Higgs is stable. Higgs does not exist. 2 If

More information

Fluctuations of cosmic parameters in the local universe

Fluctuations of cosmic parameters in the local universe Fluctuations of cosmic parameters in the local universe Alexander Wiegand Dominik Schwarz Fakultät für Physik Universität Bielefeld 6. Kosmologietag, Bielefeld 2011 A. Wiegand (Universität Bielefeld) Fluctuations

More information

MATH325 - QUANTUM MECHANICS - SOLUTION SHEET 11

MATH325 - QUANTUM MECHANICS - SOLUTION SHEET 11 MATH35 - QUANTUM MECHANICS - SOLUTION SHEET. The Hamiltonian for a particle of mass m moving in three dimensions under the influence of a three-dimensional harmonic oscillator potential is Ĥ = h m + mω

More information

Heisenberg s inequality

Heisenberg s inequality Heisenberg s inequality Michael Walter 1 Introduction 1 Definition. Let P be a probability measure on R. Then its mean or expectation is defined as E(P) := x dp(x) Its variance is the expected square deviation

More information

Surface x(u, v) and curve α(t) on it given by u(t) & v(t). Math 4140/5530: Differential Geometry

Surface x(u, v) and curve α(t) on it given by u(t) & v(t). Math 4140/5530: Differential Geometry Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt ( ds dt )2 Surface x(u, v)

More information

Statistical Mechanics of Active Matter

Statistical Mechanics of Active Matter Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017

More information

Week 7: Integration: Special Coordinates

Week 7: Integration: Special Coordinates Week 7: Integration: Special Coordinates Introduction Many problems naturally involve symmetry. One should exploit it where possible and this often means using coordinate systems other than Cartesian coordinates.

More information

Physics 215b: Problem Set 5

Physics 215b: Problem Set 5 Physics 25b: Problem Set 5 Prof. Matthew Fisher Solutions prepared by: James Sully April 3, 203 Please let me know if you encounter any typos in the solutions. Problem 20 Let us write the wavefunction

More information

8.05 Quantum Physics II, Fall 2011 FINAL EXAM Thursday December 22, 9:00 am -12:00 You have 3 hours.

8.05 Quantum Physics II, Fall 2011 FINAL EXAM Thursday December 22, 9:00 am -12:00 You have 3 hours. 8.05 Quantum Physics II, Fall 0 FINAL EXAM Thursday December, 9:00 am -:00 You have 3 hours. Answer all problems in the white books provided. Write YOUR NAME and YOUR SECTION on your white books. There

More information

arxiv:hep-ex/ v1 17 Sep 1999

arxiv:hep-ex/ v1 17 Sep 1999 Propagation of Errors for Matrix Inversion M. Lefebvre 1, R.K. Keeler 1, R. Sobie 1,2, J. White 1,3 arxiv:hep-ex/9909031v1 17 Sep 1999 Abstract A formula is given for the propagation of errors during matrix

More information

Offshore Hydromechanics Module 1

Offshore Hydromechanics Module 1 Offshore Hydromechanics Module 1 Dr. ir. Pepijn de Jong 4. Potential Flows part 2 Introduction Topics of Module 1 Problems of interest Chapter 1 Hydrostatics Chapter 2 Floating stability Chapter 2 Constant

More information

Exact Solutions of the Einstein Equations

Exact Solutions of the Einstein Equations Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.

More information

t Hooft loop path integral in N = 2 gauge theories

t Hooft loop path integral in N = 2 gauge theories t Hooft loop path integral in N = 2 gauge theories Jaume Gomis (based on work with Takuya Okuda and Vasily Pestun) Perimeter Institute December 17, 2010 Jaume Gomis (Perimeter Institute) t Hooft loop path

More information

Random Eigenvalue Problems Revisited

Random Eigenvalue Problems Revisited Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html

More information

Alternatives to conventional Monte Carlo

Alternatives to conventional Monte Carlo Alternatives to conventional Monte Carlo Recursive umerical Integration & Julia Volmer Andreas Ammon, Alan enz, Tobias Hartung, Karl Jansen, Hernan Leövey DESY Zeuthen 16. September 216 Julia Volmer (DESY

More information

Methoden moderner Röntgenphysik II: Streuung und Abbildung

Methoden moderner Röntgenphysik II: Streuung und Abbildung Methoden moderner Röntgenphysik II: Streuung und Abbildung Lecture 4 Location Vorlesung zum Haupt- oder Masterstudiengang Physik, SoSe 2015 G. Grübel, M. Martins, E. Weckert Lecture hall AP, Physics, Jungiusstraße

More information

MAE 101A. Homework 7 - Solutions 3/12/2018

MAE 101A. Homework 7 - Solutions 3/12/2018 MAE 101A Homework 7 - Solutions 3/12/2018 Munson 6.31: The stream function for a two-dimensional, nonviscous, incompressible flow field is given by the expression ψ = 2(x y) where the stream function has

More information

Surface Plasmon Polaritons on Structured Surfaces. Alexei A. Maradudin and Tamara A. Leskova

Surface Plasmon Polaritons on Structured Surfaces. Alexei A. Maradudin and Tamara A. Leskova Surface Plasmon Polaritons on Structured Surfaces Alexei A. Maradudin and Tamara A. Leskova Department of Physics and Astronomy and Institute for Surface and Interface Science, University of California,

More information

Finite-temperature Field Theory

Finite-temperature Field Theory Finite-temperature Field Theory Aleksi Vuorinen CERN Initial Conditions in Heavy Ion Collisions Goa, India, September 2008 Outline Further tools for equilibrium thermodynamics Gauge symmetry Faddeev-Popov

More information

σ 2 + π = 0 while σ satisfies a cubic equation λf 2, σ 3 +f + β = 0 the second derivatives of the potential are = λ(σ 2 f 2 )δ ij, π i π j

σ 2 + π = 0 while σ satisfies a cubic equation λf 2, σ 3 +f + β = 0 the second derivatives of the potential are = λ(σ 2 f 2 )δ ij, π i π j PHY 396 K. Solutions for problem set #4. Problem 1a: The linear sigma model has scalar potential V σ, π = λ 8 σ + π f βσ. S.1 Any local minimum of this potential satisfies and V = λ π V σ = λ σ + π f =

More information

On the hydrodynamic diffusion of rigid particles

On the hydrodynamic diffusion of rigid particles On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:

More information

ADVANCED TOPICS IN THEORETICAL PHYSICS II Tutorial problem set 2, (20 points in total) Problems are due at Monday,

ADVANCED TOPICS IN THEORETICAL PHYSICS II Tutorial problem set 2, (20 points in total) Problems are due at Monday, ADVANCED TOPICS IN THEORETICAL PHYSICS II Tutorial problem set, 15.09.014. (0 points in total) Problems are due at Monday,.09.014. PROBLEM 4 Entropy of coupled oscillators. Consider two coupled simple

More information

Front propagation directed by a line of fast diffusion : existence of travelling waves.

Front propagation directed by a line of fast diffusion : existence of travelling waves. Front propagation directed by a line of fast diffusion : existence of travelling waves. Laurent Dietrich Ph.D supervised by H. Berestycki and J.-M. Roquejoffre Institut de mathématiques de Toulouse Workshop

More information

Introduction Statistical Thermodynamics. Monday, January 6, 14

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

Y. Hosotani, SI2010, 15 August 2010, - 2

Y. Hosotani, SI2010, 15 August 2010, - 2 = Y. osotani, SI2010, 15 August 2010, - 2 Gauge-iggs Unification a la Kaluza-Klein osotani 1983, 1989 Davies, McLachlan 1988, 1989 atanaka, Inami, Lim, 1998 Gauge theory A M in higher dimensions 4-dim.

More information

Chap. 1. Some Differential Geometric Tools

Chap. 1. Some Differential Geometric Tools Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U

More information

QMI PRELIM Problem 1. All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work.

QMI PRELIM Problem 1. All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work. QMI PRELIM 013 All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work. Problem 1 L = r p, p = i h ( ) (a) Show that L z = i h y x ; (cyclic

More information

Complex functions in the theory of 2D flow

Complex functions in the theory of 2D flow Complex functions in the theory of D flow Martin Scholtz Institute of Theoretical Physics Charles University in Prague scholtz@utf.mff.cuni.cz Faculty of Transportation Sciences Czech Technical University

More information

Finite Temperature Field Theory

Finite Temperature Field Theory Finite Temperature Field Theory Dietrich Bödeker, Universität Bielefeld 1. Thermodynamics (better: thermo-statics) (a) Imaginary time formalism (b) free energy: scalar particles, resummation i. pedestrian

More information

1 Phase Spaces and the Liouville Equation

1 Phase Spaces and the Liouville Equation Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the

More information

St Hugh s 2 nd Year: Quantum Mechanics II. Reading. Topics. The following sources are recommended for this tutorial:

St Hugh s 2 nd Year: Quantum Mechanics II. Reading. Topics. The following sources are recommended for this tutorial: St Hugh s 2 nd Year: Quantum Mechanics II Reading The following sources are recommended for this tutorial: The key text (especially here in Oxford) is Molecular Quantum Mechanics, P. W. Atkins and R. S.

More information

FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017

FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked

More information

Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics

Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html

More information

Gradings on Lie Algebras of Cartan and Melikian Type

Gradings on Lie Algebras of Cartan and Melikian Type Gradings on Lie Algebras of Cartan and Melikian Type Jason McGraw Department of Mathematics and Statistics Memorial University of Newfoundland St. John's, Canada jason.mcgraw@mun.ca July 22, 2010 Jason

More information

e ik(x y) m 2 where the last equality follows from eq. (3.10a) for the sum over polarizations. = i

e ik(x y) m 2 where the last equality follows from eq. (3.10a) for the sum over polarizations. = i PHY 396 K. Solutions for homework set #4. Problem 1a): Each quantum field is a linear combination 3.11) of creation and annihilation operators. The product of two fields µ x)âν y) therefore decomposes

More information

Critical exponents for two-dimensional statistical

Critical exponents for two-dimensional statistical California State University Department of Mathematics Outline Two different type of Statistical Mechanics models in 2D: Spin Systems Coulomb Gas Spin Systems Ising Model Configuration: Place +1 or 1 at

More information

The Fermion Bag Approach

The Fermion Bag Approach The Fermion Bag Approach Anyi Li Duke University In collaboration with Shailesh Chandrasekharan 1 Motivation Monte Carlo simulation Sign problem Fermion sign problem Solutions to the sign problem Fermion

More information

the renormalization group (RG) idea

the renormalization group (RG) idea the renormalization group (RG) idea Block Spin Partition function Z =Tr s e H. block spin transformation (majority rule) T (s, if s i ; s,...,s 9 )= s i > 0; 0, otherwise. b Block Spin (block-)transformed

More information

A simple model for block-copolymer layer morphology formation

A simple model for block-copolymer layer morphology formation A simple model for block-copolymer layer morphology formation A. J. Wagner and A. Croll, Department of Physics, North Dakota State University, Fargo, ND 58102 Prato Workshop, February 12, 2016 A big thanks

More information

Dynamics of Relativistic Particles and EM Fields

Dynamics of Relativistic Particles and EM Fields October 7, 2008 1 1 J.D.Jackson, Classical Electrodynamics, 3rd Edition, Chapter 12 Lagrangian Hamiltonian for a Relativistic Charged Particle The equations of motion [ d p dt = e E + u ] c B de dt = e

More information

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n. Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier

More information

IFE for Stokes interface problem

IFE for Stokes interface problem IFE for Stokes interface problem Nabil Chaabane Slimane Adjerid, Tao Lin Virginia Tech SIAM chapter February 4, 24 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2 Problem statement

More information

Kosterlitz-Thouless Transition

Kosterlitz-Thouless Transition Heidelberg University Seminar-Lecture 5 SS 16 Menelaos Zikidis Kosterlitz-Thouless Transition 1 Motivation Back in the 70 s, the concept of a phase transition in condensed matter physics was associated

More information

Lecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations

Lecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations Lecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations Usama Al Khawaja, Department of Physics, UAE University, 24 Jan. 2012 First International Winter School on Quantum Gases

More information

2.20 Fall 2018 Math Review

2.20 Fall 2018 Math Review 2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more

More information

Quantum Physics III (8.06) Spring 2008 Final Exam Solutions

Quantum Physics III (8.06) Spring 2008 Final Exam Solutions Quantum Physics III (8.6) Spring 8 Final Exam Solutions May 19, 8 1. Short answer questions (35 points) (a) ( points) α 4 mc (b) ( points) µ B B, where µ B = e m (c) (3 points) In the variational ansatz,

More information

Fractional exclusion statistics: A generalised Pauli principle

Fractional exclusion statistics: A generalised Pauli principle Fractional exclusion statistics: A generalised Pauli principle M.V.N. Murthy Institute of Mathematical Sciences, Chennai (murthy@imsc.res.in) work done with R. Shankar Indian Academy of Sciences, 27 Nov.

More information

Interference, vector spaces

Interference, vector spaces Interference, vector spaces Sourendu Gupta TIFR Graduate School Quantum Mechanics 1 August 6, 2008 Sourendu Gupta (TIFR Graduate School) Interference, vector spaces QM I 1 / 16 Outline 1 The double-slit

More information

Spinning strings and QED

Spinning strings and QED Spinning strings and QED James Edwards Oxford Particles and Fields Seminar January 2015 Based on arxiv:1409.4948 [hep-th] and arxiv:1410.3288 [hep-th] Outline Introduction Various relationships between

More information

Nonlinear stability of steady flow of Giesekus viscoelastic fluid

Nonlinear stability of steady flow of Giesekus viscoelastic fluid Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction

More information

NATURAL SCIENCES TRIPOS. Past questions. EXPERIMENTAL AND THEORETICAL PHYSICS Minor Topics. (27 February 2010)

NATURAL SCIENCES TRIPOS. Past questions. EXPERIMENTAL AND THEORETICAL PHYSICS Minor Topics. (27 February 2010) NATURAL SCIENCES TRIPOS Part III Past questions EXPERIMENTAL AND THEORETICAL PHYSICS Minor Topics (27 February 21) 1 In one-dimension, the q-state Potts model is defined by the lattice Hamiltonian βh =

More information

ACOUSTIC BLACK HOLES. MASSIMILIANO RINALDI Université de Genève

ACOUSTIC BLACK HOLES. MASSIMILIANO RINALDI Université de Genève ACOUSTIC BLACK HOLES MASSIMILIANO RINALDI Université de Genève OUTLINE Prelude: GR vs QM Hawking Radiation: a primer Acoustic Black Holes Hawking Radiation in Acoustic Black Holes Acoustic Black Holes

More information

Quantum Mechanics Exercises and solutions

Quantum Mechanics Exercises and solutions Quantum Mechanics Exercises and solutions P.J. Mulders Department of Physics and Astronomy, Faculty of Sciences, Vrije Universiteit Amsterdam De Boelelaan 181, 181 HV Amsterdam, the Netherlands email:

More information

Physics 3323, Fall 2016 Problem Set 2 due Sep 9, 2016

Physics 3323, Fall 2016 Problem Set 2 due Sep 9, 2016 Physics 3323, Fall 26 Problem Set 2 due Sep 9, 26. What s my charge? A spherical region of radius R is filled with a charge distribution that gives rise to an electric field inside of the form E E /R 2

More information

Free energy concept Free energy approach LBM implementation Parameters

Free energy concept Free energy approach LBM implementation Parameters BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH

More information

PHYS 5012 Radiation Physics and Dosimetry

PHYS 5012 Radiation Physics and Dosimetry Radiative PHYS 5012 Radiation Physics and Dosimetry Mean Tuesday 24 March 2009 Radiative Mean Radiative Mean Collisions between two particles involve a projectile and a target. Types of targets: whole

More information

12. MHD Approximation.

12. MHD Approximation. Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal

More information

Additional Mathematical Tools: Detail

Additional Mathematical Tools: Detail Additional Mathematical Tools: Detail September 9, 25 The material here is not required, but gives more detail on the additional mathmatical tools: coordinate systems, rotations, the Dirac delta function

More information

The Path Integral: Basics and Tricks (largely from Zee)

The Path Integral: Basics and Tricks (largely from Zee) The Path Integral: Basics and Tricks (largely from Zee) Yichen Shi Michaelmas 03 Path-Integral Derivation x f, t f x i, t i x f e H(t f t i) x i. If we chop the path into N intervals of length ɛ, then

More information

Quantum Mechanics I, Sheet 1, Spring 2015

Quantum Mechanics I, Sheet 1, Spring 2015 Quantum Mechanics I, Sheet 1, Spring 2015 February 18, 2015 (EP, Auditoire Stuckelberg) Prof. D. van der Marel (dirk.vandermarel@unige.ch) Tutorial: O. Peil (oleg.peil@unige.ch) I. PROBABILITY DISTRIBUTION

More information

Physics 486 Midterm Exam #1 Spring 2018 Thursday February 22, 9:30 am 10:50 am

Physics 486 Midterm Exam #1 Spring 2018 Thursday February 22, 9:30 am 10:50 am Physics 486 Midterm Exam #1 Spring 18 Thursday February, 9: am 1:5 am This is a closed book exam. No use of calculators or any other electronic devices is allowed. Work the problems only in your answer

More information

Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.

Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.

More information

Yangians In Deformed Super Yang-Mills Theories

Yangians In Deformed Super Yang-Mills Theories Yangians In Deformed Super Yang-Mills Theories arxiv:0802.3644 [hep-th] JHEP 04:051, 2008 Jay N. Ihry UNC-CH April 17, 2008 Jay N. Ihry UNC-CH Yangians In Deformed Super Yang-Mills Theories arxiv:0802.3644

More information

This is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or

This is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects

More information

Vectors. Three dimensions. (a) Cartesian coordinates ds is the distance from x to x + dx. ds 2 = dx 2 + dy 2 + dz 2 = g ij dx i dx j (1)

Vectors. Three dimensions. (a) Cartesian coordinates ds is the distance from x to x + dx. ds 2 = dx 2 + dy 2 + dz 2 = g ij dx i dx j (1) Vectors (Dated: September017 I. TENSORS Three dimensions (a Cartesian coordinates ds is the distance from x to x + dx ds dx + dy + dz g ij dx i dx j (1 Here dx 1 dx, dx dy, dx 3 dz, and tensor g ij is

More information

Interface Profiles in Field Theory

Interface Profiles in Field Theory Florian König Institut für Theoretische Physik Universität Münster January 10, 2011 / Forschungsseminar Quantenfeldtheorie Outline φ 4 -Theory in Statistical Physics Critical Phenomena and Order Parameter

More information

Solutions to Problems in Jackson, Classical Electrodynamics, Third Edition. Chapter 2: Problems 11-20

Solutions to Problems in Jackson, Classical Electrodynamics, Third Edition. Chapter 2: Problems 11-20 Solutions to Problems in Jackson, Classical Electrodynamics, Third Edition Homer Reid December 8, 999 Chapter : Problems - Problem A line charge with linear charge density τ is placed parallel to, and

More information

Noncollinear spins in QMC: spiral Spin Density Waves in the HEG

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

QFT Perturbation Theory

QFT Perturbation Theory QFT Perturbation Theory Ling-Fong Li Institute) Slide_04 1 / 44 Interaction Theory As an illustration, take electromagnetic interaction. Lagrangian density is The combination is the covariant derivative.

More information

Classical and Quantum Bianchi type I cosmology in K-essence theory

Classical and Quantum Bianchi type I cosmology in K-essence theory Classical and Quantum Bianchi type I cosmology in K-essence theory Luis O. Pimentel 1, J. Socorro 1,2, Abraham Espinoza-García 2 1 Departamento de Fisica de la Universidad Autonoma Metropolitana Iztapalapa,

More information

The particular integral :

The particular integral : Second order linear equation with constant coefficients The particular integral : d f df Lf = a + a + af = h( x) Solutions with combinations of driving functions d f df Lf = a + a + af = h( x) + h( x)

More information

The Ideal Gas. One particle in a box:

The Ideal Gas. One particle in a box: IDEAL GAS The Ideal Gas It is an important physical example that can be solved exactly. All real gases behave like ideal if the density is small enough. In order to derive the law, we have to do following:

More information

Modeling of Micro-Fluidics by a Dissipative Particle Dynamics Method. Justyna Czerwinska

Modeling of Micro-Fluidics by a Dissipative Particle Dynamics Method. Justyna Czerwinska Modeling of Micro-Fluidics by a Dissipative Particle Dynamics Method Justyna Czerwinska Scales and Physical Models years Time hours Engineering Design Limit Process Design minutes Continious Mechanics

More information

Curvilinear Coordinates

Curvilinear Coordinates University of Alabama Department of Physics and Astronomy PH 106-4 / LeClair Fall 2008 Curvilinear Coordinates Note that we use the convention that the cartesian unit vectors are ˆx, ŷ, and ẑ, rather than

More information

3 Chapter. Gauss s Law

3 Chapter. Gauss s Law 3 Chapter Gauss s Law 3.1 Electric Flux... 3-2 3.2 Gauss s Law (see also Gauss s Law Simulation in Section 3.10)... 3-4 Example 3.1: Infinitely Long Rod of Uniform Charge Density... 3-9 Example 3.2: Infinite

More information

Ginzburg-Landau Phenomenology

Ginzburg-Landau Phenomenology Chapter Ginzburg-Landau Phenomenology The divergence of the correlation length in the vicinity of a second-order phase transition indicates that the properties of the critical point are insensitive to

More information

Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension

Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey

More information

Path Integral methods for solving stochastic problems. Carson C. Chow, NIH

Path Integral methods for solving stochastic problems. Carson C. Chow, NIH Path Integral methods for solving stochastic problems Carson C. Chow, NIH Why? Often in neuroscience we run into stochastic ODEs of the form dx dt = f(x)+g(x)η(t) η(t) =0 η(t)η(t ) = δ(t t ) Examples Integrate-and-fire

More information

Physics Exam 1 Formulas

Physics Exam 1 Formulas INSTRUCTIONS: Write your NAME on the front of the blue exam booklet. The exam is closed book, and you may have only pens/pencils and a calculator (no stored equations or programs and no graphing). Show

More information

Analogous comments can be made for the regions where E < V, wherein the solution to the Schrödinger equation for constant V is

Analogous comments can be made for the regions where E < V, wherein the solution to the Schrödinger equation for constant V is 8. WKB Approximation The WKB approximation, named after Wentzel, Kramers, and Brillouin, is a method for obtaining an approximate solution to a time-independent one-dimensional differential equation, in

More information

Reminder Notes for the Course on Distribution Theory

Reminder Notes for the Course on Distribution Theory Reminder Notes for the Course on Distribution Theory T. C. Dorlas Dublin Institute for Advanced Studies School of Theoretical Physics 10 Burlington Road, Dublin 4, Ireland. Email: dorlas@stp.dias.ie March

More information

Jim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt

Jim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt Jim Lambers MAT 28 ummer emester 212-1 Practice Final Exam olution 1. Evaluate the line integral xy dx + e y dy + xz dz, where is given by r(t) t 4, t 2, t, t 1. olution From r (t) 4t, 2t, t 2, we obtain

More information

Chapter 4. Multiple Integrals. 1. Integrals on Rectangles

Chapter 4. Multiple Integrals. 1. Integrals on Rectangles Chapter 4. Multiple Integrals Let be a rectangle in I 2 given by 1. Integrals on ectangles = [a 1, b 1 [a 2, b 2 := {(x, y) I 2 : a 1 x b 1, a 2 y b 2 }. Let P 1 be a partition of [a 1, b 1 : P 1 : a 1

More information

PAPER 71 COSMOLOGY. Attempt THREE questions There are seven questions in total The questions carry equal weight

PAPER 71 COSMOLOGY. Attempt THREE questions There are seven questions in total The questions carry equal weight MATHEMATICAL TRIPOS Part III Friday 31 May 00 9 to 1 PAPER 71 COSMOLOGY Attempt THREE questions There are seven questions in total The questions carry equal weight You may make free use of the information

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent

More information

Basic Aspects of Discretization

Basic Aspects of Discretization Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)

More information

PHYSICS 234 HOMEWORK 2 SOLUTIONS. So the eigenvalues are 1, 2, 4. To find the eigenvector for ω = 1 we have b c.

PHYSICS 234 HOMEWORK 2 SOLUTIONS. So the eigenvalues are 1, 2, 4. To find the eigenvector for ω = 1 we have b c. PHYSICS 34 HOMEWORK SOUTIONS.8.. The matrix we have to diagonalize is Ω 3 ( 4 So the characteristic equation is ( ( ω(( ω(4 ω. So the eigenvalues are,, 4. To find the eigenvector for ω we have 3 a (3 b

More information

PHYS Handout 6

PHYS Handout 6 PHYS 060 Handout 6 Handout Contents Golden Equations for Lectures 8 to Answers to examples on Handout 5 Tricks of the Quantum trade (revision hints) Golden Equations (Lectures 8 to ) ψ Â φ ψ (x)âφ(x)dxn

More information