Experimental predictions for strongly correlating liquids

Size: px
Start display at page:

Download "Experimental predictions for strongly correlating liquids"

Transcription

1 Experimental predictions for strongly correlating liquids A new look at the Prigogine-Defay ratio Ulf Rørbæk Pedersen Glass and Time - DNRF centre for viscous liquids dynamics, IMFUFA, Department of Science, Systems and Models, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark Collaborators: Nicholas Bailey, Nicoletta Gnan, Thomas B. Schrøder and Jeppe C. Dyre April 3, 2009 Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

2 Strong W-U Correlations in the Lennard-Jones liquid The well-studied Lennard-Jones liquid: [ ( U pair (r) = 4ε σ ) 12 ( r σ ) ] 6 r E = K + U pv = Nk B T + W Fluctuations of conf. parts: U(t) = U(t) U W (t) = W (t) W Constant V and T simulations. [Pedersen et al. (2008) PRL 100:015701] W γ U at constant VT Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

3 Relating to experiments Strongly correlating liquids (SCL) have W γ U at constant VT : ar n + br + c U pair LJ Two important numbers Degree of correlation: R W U VT ( W ) 2 VT ( U) 2 VT General problem Recall isomorphs: F (ρ γ /T ) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

4 Relating to experiments Strongly correlating liquids (SCL) have W γ U at constant VT : ar n + br + c U pair LJ Two important numbers Degree of correlation: R W U VT ( W ) 2 VT ( U) 2 VT Apparent exponent: γ W VT U VT = n/3 General problem Recall isomorphs: F (ρ γ /T ) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

5 Relating to experiments Strongly correlating liquids (SCL) have W γ U at constant VT : ar n + br + c U pair LJ Two important numbers Degree of correlation: R W U VT ( W ) 2 VT ( U) 2 VT Apparent exponent: γ W VT U VT = n/3 General problem Subtract off ideal gas -terms: W = pv Nk B T and U = E K. Recall isomorphs: F (ρ γ /T ) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

6 Relating to experiments Strongly correlating liquids (SCL) have W γ U at constant VT : ar n + br + c U pair LJ Two important numbers Degree of correlation: R W U VT ( W ) 2 VT ( U) 2 VT Apparent exponent: γ W VT U VT = n/3 General problem Subtract off ideal gas -terms: W = pv Nk B T and U = E K. Constant VT -ensemble (where we understand SCL) vs. constant pt -ensemble (experiments) Recall isomorphs: F (ρ γ /T ) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

7 γ related to crystallization γ = W a W b U a U b β V ( ) p T V = ( ) S V T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

8 γ related to crystallization γ = W a W b = Vr (pa p b ) U a U b E a E b β V ( ) p T V = ( ) S V T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

9 γ related to crystallization γ = W a W b U a U b = Vr (pa p b ) E a E b = Vr (pa p b ) T r (S a S b ) β V ( ) p T V = ( ) S V T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

10 γ related to crystallization β V ( ) p T V = ( ) S V T γ = W a W b U a U b = Vr (pa p b ) E a E b = Vr (pa p b ) T r (S a S b ) (1) Cons. pt melting (2) compress liq.: p a p b = 0 Vr + fus V (p b ) V r S a S b = fus S(p b )+ VK a T dv Vr + fus V (p b ) V r β b V dv Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

11 Pressure and energy fluctuations of viscous liquids. C EE (t) = C EE (t)/c EE (0), C EE (t) = E(0) E(t) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

12 Pressure and energy fluctuations of viscous liquids. C EE (t) = C EE (t)/c EE (0), C EE (t) = E(0) E(t) C EE (t 0) = C EE = ( E) 2 liquid ( E) 2 solid. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

13 Pressure and energy fluctuations of viscous liquids. C EE (t) = C EE (t)/c EE (0), C EE (t) = E(0) E(t) C EE (t 0) = C EE = ( E) 2 liquid ( E) 2 solid. Subtracting off ideal gas terms: C EE = C UU C pp = V 2 C WW C pe = VC WU. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

14 Pressure and energy fluctuations of viscous liquids. C EE (t) = C EE (t)/c EE (0), C EE (t) = E(0) E(t) C EE (t 0) = C EE = ( E) 2 liquid ( E) 2 solid. Subtracting off ideal gas terms: C EE = C UU C pp = V 2 C WW C pe = VC WU. R C pe C pp C EE Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

15 Pressure and energy fluctuations of viscous liquids. C EE (t) = C EE (t)/c EE (0), C EE (t) = E(0) E(t) C EE (t 0) = C EE = ( E) 2 liquid ( E) 2 solid. Subtracting off ideal gas terms: C EE = C UU C pp = V 2 C WW C pe = VC WU. R γ C pe C pp CEE V 2 Cpp CEE Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

16 R and γ from (static) response functions R = C pe C pp C EE = p E liquid q VT p E solid VT [ ( E) 2 liquid VT ( E)2 solid VT ][ ( p)2 liquid VT ( p)2 solid VT ] Fluctuation dissipation theorem [Allen & Tildesley (1987)] Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

17 R and γ from (static) response functions R = C pe C pp C EE = p E liquid q VT p E solid VT [ ( E) 2 liquid VT ( E)2 solid VT ][ ( p)2 liquid VT ( p)2 solid VT ] Fluctuation dissipation theorem [Allen & Tildesley (1987)] Definitions: c V 1 ( E ) ( ) ( ) V T V, K T V p V and β V p T T ( E) 2 VT = k B T 2 Vc V V. R W U ( W ) 2 ( U) 2 Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

18 R and γ from (static) response functions R = C pe C pp C EE = p E liquid q VT p E solid VT [ ( E) 2 liquid VT ( E)2 solid VT ][ ( p)2 liquid VT ( p)2 solid VT ] Fluctuation dissipation theorem [Allen & Tildesley (1987)] Definitions: c V 1 ( E ) ( ) ( ) V T V, K T V p V and β V p T T ( E) 2 VT = k B T 2 Vc V ( p) 2 VT = k BT V [ 2Nk BT 3V + p VT K T + X VT V ], V. R W U ( W ) 2 ( U) 2 Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

19 R and γ from (static) response functions R = C pe C pp C EE = p E liquid q VT p E solid VT [ ( E) 2 liquid VT ( E)2 solid VT ][ ( p)2 liquid VT ( p)2 solid VT ] Fluctuation dissipation theorem [Allen & Tildesley (1987)] Definitions: c V 1 ( E ) ( ) ( ) V T V, K T V p V and β V p T T ( E) 2 VT = k B T 2 Vc V R ( p) 2 VT = k BT V [ 2Nk BT 3V + p VT K T + X VT V ], U p = k B T 2 (β V N V k B), W U ( W ) 2 ( U) 2 V. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

20 R and γ from (static) response functions R = C pe C pp C EE = p E liquid q VT p E solid VT [ ( E) 2 liquid VT ( E)2 solid VT ][ ( p)2 liquid VT ( p)2 solid VT ] Fluctuation dissipation theorem [Allen & Tildesley (1987)] Definitions: c V 1 ( E ) ( ) ( ) V T V, K T V p V and β V p T T ( E) 2 VT = k B T 2 Vc V R ( p) 2 VT = k BT V [ 2Nk BT 3V + p VT K T + X VT V ], U p = k B T 2 (β V N V k B), W U ( W ) 2 ( U) 2 V. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

21 R and γ from (static) response functions R = C pe C pp C EE = p E liquid q VT p E solid VT [ ( E) 2 liquid VT ( E)2 solid VT ][ ( p)2 liquid VT ( p)2 solid VT ] Fluctuation dissipation theorem [Allen & Tildesley (1987)] Definitions: c V 1 ( E ) ( ) ( ) V T V, K T V p V and β V p T T ( E) 2 VT = k B T 2 Vc V R γ ( p) 2 VT = k BT V [ 2Nk BT 3V + p VT K T + X VT V ], U p = k B T 2 (β V N V k B), W U ( W ) 2 ( U) 2 = ( W ) 2 ( U) 2 (β liquid V βv solid ) q T KT solid )(c liquid V cv solid )/T = V. β V KT c V /T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

22 R and γ from (static) response functions R = C pe C pp C EE = p E liquid q VT p E solid VT [ ( E) 2 liquid VT ( E)2 solid VT ][ ( p)2 liquid VT ( p)2 solid VT ] Fluctuation dissipation theorem [Allen & Tildesley (1987)] Definitions: c V 1 ( E ) ( ) ( ) V T V, K T V p V and β V p T T ( E) 2 VT = k B T 2 Vc V R γ ( p) 2 VT = k BT V [ 2Nk BT 3V + p VT K T + X VT V ], U p = k B T 2 (β V N V k B), W U ( W ) 2 ( U) 2 = ( W ) 2 ( U) 2 = (β liquid V βv solid ) q T KT solid )(c liquid V cv solid )/T T KT solid ) T (cv liquid c solid = KT v ) T c v = V. β V KT c V /T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

23 A standard Roskilde liquid - DC704 silicon oil (a SCL) K S (ω): [Tina Hecksher, unpub.] T ref = 214K. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

24 A standard Roskilde liquid - DC704 silicon oil (a SCL) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] T ref = 214K. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

25 A standard Roskilde liquid - DC704 silicon oil (a SCL) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

26 A standard Roskilde liquid - DC704 silicon oil (a SCL) γ ( W ) 2 ( U) 2 = T KT solid ) T (cv liquid cv solid ) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

27 A standard Roskilde liquid - DC704 silicon oil (a SCL) γ ( W ) 2 ( U) 2 = c V = c p T α 2 pk T T KT solid ) T (cv liquid cv solid ) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

28 A standard Roskilde liquid - DC704 silicon oil (a SCL) γ ( W ) 2 ( U) 2 = c V = c p T α 2 pk T K T = K S c V cp T KT solid ) T (cv liquid cv solid ) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

29 A standard Roskilde liquid - DC704 silicon oil (a SCL) γ ( W ) 2 ( U) 2 = c V = c p T α 2 pk T c K T = K V S cp ( c V = c p 1 + T α2 p K S c p T KT solid ) T (cv liquid cv solid ) ) 1 K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

30 A standard Roskilde liquid - DC704 silicon oil (a SCL) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) γ ( W ) 2 ( U) 2 = c V = c p T α 2 pk T c K T = K V S cp ( c V = c p 1 + T α2 p K S c p ( K T = K S 1 + T α2 p K S c p T KT solid ) T (cv liquid cv solid ) ) 1 ) 1 Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

31 A standard Roskilde liquid - DC704 silicon oil (a SCL) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) γ ( W ) 2 ( U) 2 = c V = c p T α 2 pk T c K T = K V S cp ( c V = c p 1 + T α2 p K S c p ( K T = K S 1 + T α2 p K S T (α liquid p ) 2 K liquid S T (α solid p c liquid p ) 2 KS solid c solid p c p = 0.13 = 0.01 T KT solid ) T (cv liquid cv solid ) ) 1 ) 1 Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

32 A standard Roskilde liquid - DC704 silicon oil (a SCL) K S (ω): [Tina Hecksher, unpub.] C p (ω): [Bo Jacobsen, unpub.] α liquid P αp solid T ref = 214K. = K 1 = K 1 (guess) γ ( W ) 2 ( U) 2 = c V = c p T α 2 pk T c K T = K V S cp ( c V = c p 1 + T α2 p K S c p ( K T = K S 1 + T α2 p K S T (α liquid p ) 2 K liquid S T (α solid p c liquid p ) 2 KS solid c solid p γ = 6.6 c p = 0.13 = 0.01 T KT solid ) T (cv liquid cv solid ) ) 1 ) 1 Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

33 A method of extrapolation [Takahara et al. (1999)] OTP-OPP data. R = q γ = T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

34 A method of extrapolation [Takahara et al. (1999)] OTP-OPP data. R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

35 A method of extrapolation [Takahara et al. (1999)] OTP-OPP data. R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) c V = c p T α 2 pk T )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

36 A method of extrapolation [Takahara et al. (1999)] OTP-OPP data. R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) c V = c p T α 2 pk T β V = α p K T )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

37 A method of extrapolation [Takahara et al. (1999)] OTP-OPP data. R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) c V = c p T α 2 pk T β V = α p K T R = 0.79 )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

38 A method of extrapolation [Takahara et al. (1999)] OTP-OPP data. R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) c V = c p T α 2 pk T β V = α p K T R = 0.79 γ = 6.0 )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

39 Density scaling OTP-OPP τ = f (ρ γ /T ) From [Roland et al., Rep. Prog. Phys. 68, 1405 (2005)] Reanalyzing OTP-OPP glass transition data from [Takanhara et al., 1999]: K T T c v = 6.0 γ scaling Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

40 A method of extrapolation. What are the uncertainties? R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) c V = c p T α 2 pk T β V = α p K T R = 0.79 γ = 6.0 (c liquid p c liquid V )/c liquid V = 26% )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

41 A method of extrapolation. What are the uncertainties? R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) c V = c p T α 2 pk T β V = α p K T R = 0.79 γ = 6.0 (c liquid p (c solid p c liquid V )/c liquid c solid V )/c solid V = 26% V = 3% )/T Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

42 A method of extrapolation. What are the uncertainties? R = q γ = T K T = 1/κ T (β liquid V βv solid ) KT solid )(cv liquid cv solid T KT solid ) T (cv liquid cv solid ) c V = c p T α 2 pk T β V = α p K T R = 0.79 γ = 6.0 (c liquid p (c solid p (c liquid V c liquid V )/c liquid c solid V )/c solid V = 26% V = 3% cv solid )/c liquid V = 19% Ulf R. Pedersen (Roskilde Unv.) April 3, / 14 )/T

43 Literature values of the (pt) Prigogine-Defay ratio Material Π pt Π 1/2 pt Pure SiO 2 glass Glycerol, C 3 H 5 (OH) GeO 2 6.9± ± Na 2 O10B 2 O 3 74SiO B 2 O Ca NO 3 0.6KNO Glucose, OC 6 H 7 (OH) Se ZrTiCuNiBe Na 2 O74SiO Polyvinylacetate 2.2(1.7) 0.67(0.77) n-propanol, C 3 H 7 OH Polyvinylchlorid, PVC Polystyrene, PS 1.3± ±0.03 OTP(75%)-TPCM(25%) OTP(67%)-OPP(33%) 1.20± ±0.02 Polystyrene, PS Phenoxy Polycarbonate Polysulfone Polyarylate Polyisobutene Π VT = K T c V T ( β V ) 2 R = Π 1/2 VT Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

44 Literature values of the (pt) Prigogine-Defay ratio Material Π pt Π 1/2 pt Pure SiO 2 glass Glycerol, C 3 H 5 (OH) GeO 2 6.9± ± Na 2 O10B 2 O 3 74SiO B 2 O Ca NO 3 0.6KNO Glucose, OC 6 H 7 (OH) Se ZrTiCuNiBe Na 2 O74SiO Polyvinylacetate 2.2(1.7) 0.67(0.77) n-propanol, C 3 H 7 OH Polyvinylchlorid, PVC Polystyrene, PS 1.3± ±0.03 OTP(75%)-TPCM(25%) OTP(67%)-OPP(33%) 1.20± ±0.02 Polystyrene, PS Phenoxy Polycarbonate Polysulfone Polyarylate Polyisobutene Π VT = K T c V T ( β V ) 2 R = Π 1/2 VT Π pt = κ T c p T ( α p) 2 1 = Π pt = Π VT [JCP 126, (2007)] Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

45 Literature values of the (pt) Prigogine-Defay ratio Material Π pt Π 1/2 pt Pure SiO 2 glass Glycerol, C 3 H 5 (OH) GeO 2 6.9± ± Na 2 O10B 2 O 3 74SiO B 2 O Ca NO 3 0.6KNO Glucose, OC 6 H 7 (OH) Se ZrTiCuNiBe Na 2 O74SiO Polyvinylacetate 2.2(1.7) 0.67(0.77) n-propanol, C 3 H 7 OH Polyvinylchlorid, PVC Polystyrene, PS 1.3± ±0.03 OTP(75%)-TPCM(25%) OTP(67%)-OPP(33%) 1.20± ±0.02 Polystyrene, PS Phenoxy Polycarbonate Polysulfone Polyarylate Polyisobutene Π VT = K T c V T ( β V ) 2 R = Π 1/2 VT Π pt = κ T c p T ( α p) 2 1 = Π pt = Π VT [JCP 126, (2007)] Π pt 1 for van der Waals bonded glass formers. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

46 Literature values of the (pt) Prigogine-Defay ratio Material Π pt Π 1/2 pt Pure SiO 2 glass Glycerol, C 3 H 5 (OH) GeO 2 6.9± ± Na 2 O10B 2 O 3 74SiO B 2 O Ca NO 3 0.6KNO Glucose, OC 6 H 7 (OH) Se ZrTiCuNiBe Na 2 O74SiO Polyvinylacetate 2.2(1.7) 0.67(0.77) n-propanol, C 3 H 7 OH Polyvinylchlorid, PVC Polystyrene, PS 1.3± ±0.03 OTP(75%)-TPCM(25%) OTP(67%)-OPP(33%) 1.20± ±0.02 Polystyrene, PS Phenoxy Polycarbonate Polysulfone Polyarylate Polyisobutene Π VT = K T c V T ( β V ) 2 R = Π 1/2 VT Π pt = κ T c p T ( α p) 2 1 = Π pt = Π VT [JCP 126, (2007)] Π pt 1 for van der Waals bonded glass formers.... not network forming. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

47 R and γ from imaginary parts of response c K β ω V (ω) = C Vk B T 2 ω { E(0) E(t) VT } T (ω) = ω Vk B T C ω{ p(0) p(t) VT } V (ω) = ω C Vk B T 2 ω { E(0) p(t) VT } Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

48 R and γ from imaginary parts of response c K β ω V (ω) = C Vk B T 2 ω { E(0) E(t) VT } T (ω) = ω Vk B T C ω{ p(0) p(t) VT } V (ω) = ω C Vk B T 2 ω { E(0) p(t) VT } R C ω{ E(0) p(t) VT } Cω{ E(0) E(t) VT }C ω{ p(0) p(t) VT } Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

49 R and γ from imaginary parts of response c K β ω V (ω) = C Vk B T 2 ω { E(0) E(t) VT } T (ω) = ω Vk B T C ω{ p(0) p(t) VT } V (ω) = ω C Vk B T 2 ω { E(0) p(t) VT } R C ω{ E(0) p(t) VT } Cω{ E(0) E(t) VT }C ω{ p(0) p(t) VT } R(ω) = β V (ω) K T (ω)c V (ω)/t Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

50 R and γ from imaginary parts of response c K β ω V (ω) = C Vk B T 2 ω { E(0) E(t) VT } T (ω) = ω Vk B T C ω{ p(0) p(t) VT } V (ω) = ω C Vk B T 2 ω { E(0) p(t) VT } R C ω{ E(0) p(t) VT } Cω{ E(0) E(t) VT }C ω{ p(0) p(t) VT } R(ω) = γ(ω) = β V (ω) K T (ω)c V (ω)/t K T (ω) Tc V (ω). Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

51 Conclusions Crystallization of SCL: γ = Vr (pa p b ) T r (S a S b ) at constant VT Strongly correlating liquid (SCL) W γ U at constant VT R W U VT ( W ) 2 VT ( U) 2 VT γ W VT U VT = n/3. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

52 Thank you for your attention. Ulf R. Pedersen (Roskilde Unv.) April 3, / 14

Experimental predictions for strongly correlating liquids

Experimental predictions for strongly correlating liquids Experimental predictions for strongly correlating liquids A new look at the Prigogine-Defay ratio, and how to estimate γ Ulf Rørbæk Pedersen Glass and Time - DNRF centre for viscous liquids dynamics, IMFUFA,

More information

Feasibility of a single-parameter description of equilibrium viscous liquid dynamics

Feasibility of a single-parameter description of equilibrium viscous liquid dynamics PHYSICAL REVIEW E 77, 2 28 Feasibility of a single-parameter description of equilibrium viscous liquid dynamics Ulf R. Pedersen, Tage Christensen, Thomas B. Schrøder, and Jeppe C. Dyre DNRF centre Glass

More information

arxiv: v1 [cond-mat.soft] 7 Apr 2010

arxiv: v1 [cond-mat.soft] 7 Apr 2010 Strongly correlating liquids and their isomorphs Ulf R. Pedersen Department of Chemistry, University of California, Berkeley, California 94720, USA arxiv:1004.1182v1 [cond-mat.soft] 7 Apr 2010 Nicoletta

More information

New from Ulf in Berzerkeley

New from Ulf in Berzerkeley New from Ulf in Berzerkeley from crystallization to statistics of density fluctuations Ulf Rørbæk Pedersen Department of Chemistry, University of California, Berkeley, USA Roskilde, December 16th, 21 Ulf

More information

Intermolecular distance and density scaling of dynamics in molecular liquids

Intermolecular distance and density scaling of dynamics in molecular liquids Intermolecular distance and density scaling of dynamics in molecular liquids D. Fragiadakis and C. M. Roland Naval Research Laboratory, Chemistry Division, Washington, DC 20375-5342, USA (April 2, 2019)

More information

arxiv: v1 [physics.chem-ph] 4 Feb 2014

arxiv: v1 [physics.chem-ph] 4 Feb 2014 Discrete dynamics versus analytic dynamics Søren Toxvaerd DNRF centre Glass and Time, IMFUFA, Department of Sciences, Roskilde University, arxiv:1402.0679v1 [physics.chem-ph] 4 Feb 2014 Postbox 260, DK-4000

More information

Linear-response theory and the fluctuation-dissipation theorem: An executive summary

Linear-response theory and the fluctuation-dissipation theorem: An executive summary Notes prepared for the 27 Summer School at Søminestationen, Holbæk, July 1-8 Linear-response theory and the fluctuation-dissipation theorem: An executive summary Jeppe C. Dyre DNRF centre Glass and Time,

More information

arxiv:cond-mat/ v2 [cond-mat.soft] 10 Mar 2006

arxiv:cond-mat/ v2 [cond-mat.soft] 10 Mar 2006 An energy landscape model for glass-forming liquids in three dimensions arxiv:cond-mat/v [cond-mat.soft] Mar Ulf R. Pedersen, Tina Hecksher, Jeppe C. Dyre, and Thomas B. Schrøder Department of Mathematics

More information

arxiv: v1 [cond-mat.soft] 28 Jul 2018

arxiv: v1 [cond-mat.soft] 28 Jul 2018 The pair-potential system. I. Fluid phase isotherms, isochores, and quasiuniversality Andreas Kvist Bacher, Thomas B. Schrøder, and Jeppe C. Dyre Glass and Time, IMFUFA, Department of Science and Environment,

More information

Computing phase diagrams of model liquids & self-assembly of (bio)membranes

Computing phase diagrams of model liquids & self-assembly of (bio)membranes Computing phase diagrams of model liquids & self-assembly of (bio)membranes Ulf Rørbæk Pedersen (ulf@urp.dk) Crystal structure of a simplistic model of a molecule. What is the melting temperature? A typical

More information

arxiv: v1 [cond-mat.stat-mech] 26 Sep 2013

arxiv: v1 [cond-mat.stat-mech] 26 Sep 2013 EPJ manuscript No. (will be inserted by the editor) Density-temperature scaling of the fragility in a model glass-former arxiv:1309.6891v1 [cond-mat.stat-mech] 26 Sep 2013 Shiladitya Sengupta 1,a homas

More information

Density-temperature scaling of the fragility in a model glass-former

Density-temperature scaling of the fragility in a model glass-former Eur. Phys. J. E (2013) 36: 141 DOI 10.1140/epje/i20133141-9 Regular Article HE EUROPEAN PHYSICAL JOURNAL E Density-temperature scaling of the fragility in a model glass-former Shiladitya Sengupta 1,a,

More information

arxiv: v2 [cond-mat.soft] 6 Jun 2016

arxiv: v2 [cond-mat.soft] 6 Jun 2016 Studies of the Lennard-Jones fluid in 2, 3, and 4 dimensions highlight the need for a liquid-state 1/d expansion arxiv:1605.06007v2 [cond-mat.soft] 6 Jun 2016 Lorenzo Costigliola, Thomas B. Schrøder, and

More information

Chapter 4 Phase Transitions. 4.1 Phenomenology Basic ideas. Partition function?!?! Thermodynamic limit Statistical Mechanics 1 Week 4

Chapter 4 Phase Transitions. 4.1 Phenomenology Basic ideas. Partition function?!?! Thermodynamic limit Statistical Mechanics 1 Week 4 Chapter 4 Phase Transitions 4.1 Phenomenology 4.1.1 Basic ideas Partition function?!?! Thermodynamic limit 4211 Statistical Mechanics 1 Week 4 4.1.2 Phase diagrams p S S+L S+G L S+G L+G G G T p solid triple

More information

Imperfect Gases. NC State University

Imperfect Gases. NC State University Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is

More information

Effect of density on the physical aging of pressure-densified polymethylmethacrylate

Effect of density on the physical aging of pressure-densified polymethylmethacrylate ABSTRACT Effect of density on the physical aging of pressure-densified polymethylmethacrylate R. Casalini and C.M. Roland Naval Research Laboratory, Chemistry Division, Washington DC 20375-5342 (July 6,

More information

Part 2: Molecular Dynamics. Literature History Practical Issues

Part 2: Molecular Dynamics. Literature History Practical Issues Part 2: Molecular Dynamics Literature History Practical Issues 1 Literature A. K. Hartmann Big Practical Guide to Computer Simulations (World Scientific, 2015) M. P. Allen and D.J. Tildesley Computer Simulations

More information

Microscopic Picture of Aging in SiO 2 : A Computer Simulation

Microscopic Picture of Aging in SiO 2 : A Computer Simulation Microscopic Picture of Aging in SiO 2 : A Computer Simulation Katharina Vollmayr-Lee, Robin Bjorkquist, Landon M. Chambers Bucknell University & Göttingen 7 td tb r n (t) 6 R 5 4 3 2 ti t [ns] waiting

More information

Outline Review Example Problem 1 Example Problem 2. Thermodynamics. Review and Example Problems. X Bai. SDSMT, Physics. Fall 2013

Outline Review Example Problem 1 Example Problem 2. Thermodynamics. Review and Example Problems. X Bai. SDSMT, Physics. Fall 2013 Review and Example Problems SDSMT, Physics Fall 013 1 Review Example Problem 1 Exponents of phase transformation 3 Example Problem Application of Thermodynamic Identity : contents 1 Basic Concepts: Temperature,

More information

Section B: Electrical Engineering Answer ALL questions

Section B: Electrical Engineering Answer ALL questions Section B: Electrical Engineering Answer ALL questions EE8) (0 Marks) The transient current flowing in an RL electrical circuit is a solution of the first order differential equation: di R + i = 0. dt

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

Pre-yield non-affine fluctuations and a hidden critical point in strained crystals

Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Supplementary Information for: Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Tamoghna Das, a,b Saswati Ganguly, b Surajit Sengupta c and Madan Rao d a Collective Interactions

More information

arxiv: v2 [physics.chem-ph] 4 Sep 2014

arxiv: v2 [physics.chem-ph] 4 Sep 2014 Cooee bitumen II: Stability of linear asphaltene nanoaggregates arxiv:1406.7592v2 [physics.chem-ph] 4 Sep 2014 Claire A. Lemarchand DNRF Centre Glass and Time, IMFUFA, Department of Sciences, Roskilde

More information

An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation

An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation J. Comput. Chem. Jpn., Vol. 8, o. 3, pp. 97 14 (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA

More information

c i c c T T h I i I k L m u u v v V

c i c c T T h I i I k L m u u v v V A a 2 c i c c E = F = f = f g h I i I k L m p N n p = P = P = Q = r R Ω = r s t T T u u v v V = V W = x λ ρ 1 1000 1 10 6 1 10 9 1 0.01 1 0.001 1 10 6 1 10 9 1 2 1 1 0.01 0.01 10 4 2 s = vt v u = at F

More information

arxiv: v4 [cond-mat.soft] 18 Jan 2012

arxiv: v4 [cond-mat.soft] 18 Jan 2012 The dynamic thermal expansivity of liquids near the glass transition. Kristine Niss, Ditte Gundermann, Tage Christensen, and Jeppe C. Dyre DNRF Centre Glass and Time, IMFUFA, Department of Sciences, Roskilde

More information

Important practical questions:

Important practical questions: Colloidal stability Important practical questions: - Does dispersion change when T, P or... is changed? - If T, P or... is changed and the dispersion phase separates, what are then the final products?

More information

Density scaling of the diffusion coefficient at various pressures in viscous liquids [accepted for publication in Phys. Rev.

Density scaling of the diffusion coefficient at various pressures in viscous liquids [accepted for publication in Phys. Rev. Density scaling of the diffusion coefficient at various pressures in viscous liquids [accepted for publication in Phys. Rev. E (2009)] Anthony N. Papathanassiou University of Athens, Physics Department,

More information

arxiv: v1 [cond-mat.other] 5 Oct 2015

arxiv: v1 [cond-mat.other] 5 Oct 2015 a) st@ruc.dk Dynamics of homogeneous nucleation Søren Toxvaerd a) arxiv:1510.01065v1 [cond-mat.other] 5 Oct 2015 DNRF centre Glass and Time, IMFUFA, Department of Sciences, Roskilde University, Postbox

More information

arxiv: v1 [physics.chem-ph] 9 Jul 2013

arxiv: v1 [physics.chem-ph] 9 Jul 2013 Four-component united-atom model of bitumen J. S. Hansen DNRF Centre Glass and Time, IMFUFA, Department of Science, Systems and Models, Roskilde University, Postbox 26, DK-4 Roskilde, Denmark Claire A.

More information

Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses

Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses A normal liquid will crystallize at a melting temperature T m as it is cooled via a first-order phase transition (see figure above).

More information

T. Egami. Model System of Dense Random Packing (DRP)

T. Egami. Model System of Dense Random Packing (DRP) Introduction to Metallic Glasses: How they are different/similar to other glasses T. Egami Model System of Dense Random Packing (DRP) Hard Sphere vs. Soft Sphere Glass transition Universal behavior History:

More information

Spatially heterogeneous dynamics in supercooled organic liquids

Spatially heterogeneous dynamics in supercooled organic liquids Spatially heterogeneous dynamics in supercooled organic liquids Stephen Swallen, Marcus Cicerone, Marie Mapes, Mark Ediger, Robert McMahon, Lian Yu UW-Madison NSF Chemistry 1 Image from Weeks and Weitz,

More information

Nanoscale simulation lectures Statistical Mechanics

Nanoscale simulation lectures Statistical Mechanics Nanoscale simulation lectures 2008 Lectures: Thursdays 4 to 6 PM Course contents: - Thermodynamics and statistical mechanics - Structure and scattering - Mean-field approaches - Inhomogeneous systems -

More information

1.1 Atoms. 1.1 Atoms

1.1 Atoms. 1.1 Atoms 1. Chemical bonding and crystal structure 19 21 Hydrogen atom Scanning electron microscopy Ni surface Cleaved surface ZnO, TiO 2, NiO, NaCl, Si, Ge, GaAs, InP Crystals are build by small repeating units

More information

arxiv: v2 [physics.chem-ph] 3 Sep 2013

arxiv: v2 [physics.chem-ph] 3 Sep 2013 Cooee bitumen. Chemical aging arxiv:1307.2460v2 [physics.chem-ph] 3 Sep 2013 Claire A. Lemarchand DNRF Centre Glass and Time, IMFUFA, Department of Sciences, Roskilde University, Postbox 260, DK-4000 Roskilde,

More information

The correlation between fragility, density and atomic interaction in glassforming

The correlation between fragility, density and atomic interaction in glassforming The correlation between fragility, density and atomic interaction in glassforming liquids Lijin Wang 1, Pengfei Guan 1*, and W. H. Wang 2 1 Beijing Computational Science Research Center, Beijing, 100193,

More information

5.62 Physical Chemistry II Spring 2008

5.62 Physical Chemistry II Spring 2008 MIT OpenCourseWare http://ocw.mit.edu 5.62 Physical Chemistry II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.62 Lecture #20: Virial

More information

Scientific Computing II

Scientific Computing II Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the

More information

Speeding up path integral simulations

Speeding up path integral simulations Speeding up path integral simulations Thomas Markland and David Manolopoulos Department of Chemistry University of Oxford Funded by the US Office of Naval Research and the UK EPSRC Outline 1. Ring polymer

More information

Electrostatic Charged Two-Phase Turbulent Flow Model

Electrostatic Charged Two-Phase Turbulent Flow Model ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.5(2008) No.1,pp.65-70 Electrostatic Charged Two-Phase Turbulent Flow Model Jianlong Wen, Jing Wang, Zhentao Wang,

More information

Supporting Information

Supporting Information Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos

More information

Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State

Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State 23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke

More information

Phases of matter and phase diagrams

Phases of matter and phase diagrams Phases of matter and phase diagrams Transition to Supercritical CO2 Water Ice Vapor Pressure and Boiling Point Liquids boil when the external pressure equals the vapor pressure. Temperature of boiling

More information

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

Thermodynamic behaviour of mixtures containing CO 2. A molecular simulation study

Thermodynamic behaviour of mixtures containing CO 2. A molecular simulation study Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues

More information

What is Classical Molecular Dynamics?

What is Classical Molecular Dynamics? What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated

More information

Coarse-grained Models for Oligomer-grafted Silica Nanoparticles

Coarse-grained Models for Oligomer-grafted Silica Nanoparticles Modeling So+ Ma-er: Linking Mul3ple Length and Time Scales KITP Conference, Santa Barbara, June 4-8, 2012 Coarse-grained Models for Oligomer-grafted Silica Nanoparticles Bingbing Hong Alexandros Chremos

More information

kg C 10 C = J J = J kg C 20 C = J J = J J

kg C 10 C = J J = J kg C 20 C = J J = J J Seat: PHYS 1500 (Spring 2007) Exam #3, V1 Name: 5 pts 1. A pendulum is made with a length of string of negligible mass with a 0.25 kg mass at the end. A 2nd pendulum is identical except the mass is 0.50

More information

1. A Covalent bonding B Polar covalent bonding C Metallic bonding D Hydrogen bonding E Ionic bonding F London dispersion forces

1. A Covalent bonding B Polar covalent bonding C Metallic bonding D Hydrogen bonding E Ionic bonding F London dispersion forces Higher (National 6) Unit 1: Chemical changes and structure 1c Bonding and structure Problem sheets 1. A Covalent bonding B Polar covalent bonding C Metallic bonding D Hydrogen bonding E Ionic bonding F

More information

MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS. Vinu Unnithan and J. N. Reddy

MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS. Vinu Unnithan and J. N. Reddy MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS Vinu Unnithan and J. N. Reddy US-South American Workshop: Mechanics and Advanced Materials Research and Education

More information

Superfluidity in bosonic systems

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

Outline Review Example Problem 1. Thermodynamics. Review and Example Problems: Part-2. X Bai. SDSMT, Physics. Fall 2014

Outline Review Example Problem 1. Thermodynamics. Review and Example Problems: Part-2. X Bai. SDSMT, Physics. Fall 2014 Review and Example Problems: Part- SDSMT, Physics Fall 014 1 Review Example Problem 1 Exponents of phase transformation : contents 1 Basic Concepts: Temperature, Work, Energy, Thermal systems, Ideal Gas,

More information

G : Statistical Mechanics

G : Statistical Mechanics G5.651: Statistical Mechanics Notes for Lecture 8 I. DISTRIBUTION FUNCTIONS AND PERTURBATION THEORY A. General formulation Recall the expression for the configurational partition function: Z N = dr 1 dr

More information

The Temperature Dependence of the Relaxation Time in ultraviscous liquids

The Temperature Dependence of the Relaxation Time in ultraviscous liquids The Temperature Dependence of the Relaxation Time in ultraviscous liquids Is there evidence for a dynamic divergence in data? Tina Hecksher, Albena I. Nielsen, Niels B. Olsen, and Jeppe C. Dyre DNRF Centre

More information

Sub -T g Relaxation in Thin Glass

Sub -T g Relaxation in Thin Glass Sub -T g Relaxation in Thin Glass Prabhat Gupta The Ohio State University ( Go Bucks! ) Kyoto (January 7, 2008) 2008/01/07 PK Gupta(Kyoto) 1 Outline 1. Phenomenology (Review). A. Liquid to Glass Transition

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 46 Outline 1 Lecture 6: Electromechanical Systems

More information

Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions. Let's get together.

Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions. Let's get together. Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions Let's get together. Most gases are NOT ideal except at very low pressures: Z=1 for ideal gases Intermolecular interactions come

More information

Monte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects

Monte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Lecture 1 2 Lecture 1&2 3 Lecture 1&3 4 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT

More information

3.Which of the following has the highest melting temperature? A) H 2 O B) CO 2 C) S 8 D) MgF 2 E) P 4

3.Which of the following has the highest melting temperature? A) H 2 O B) CO 2 C) S 8 D) MgF 2 E) P 4 2. Which if the following is the correct order of boiling points for KNO 3, CH 3 OH, C 2 H 6, Ne? A) Ne < CH 3 OH < C 2 H 6 < KNO 3 B) KNO 3 < CH 3 OH < C 2 H 6 < Ne C) Ne < C 2 H 6 < KNO 3 < CH 3 OH D)

More information

Nanotube AFM Probe Resolution

Nanotube AFM Probe Resolution Influence of Elastic Deformation on Single-Wall Carbon Nanotube AFM Probe Resolution Ian R. Shapiro, Santiago D. Solares, Maria J. Esplandiu, Lawrence A. Wade, William A. Goddard,* and C. Patrick Collier*

More information

The Second Virial Coefficient & van der Waals Equation

The Second Virial Coefficient & van der Waals Equation V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by

More information

V.E Mean Field Theory of Condensation

V.E Mean Field Theory of Condensation V.E Mean Field heory of Condensation In principle, all properties of the interacting system, including phase separation, are contained within the thermodynamic potentials that can be obtained by evaluating

More information

arxiv: v1 [cond-mat.stat-mech] 6 Jan 2014

arxiv: v1 [cond-mat.stat-mech] 6 Jan 2014 arxiv:1401.1181v1 [cond-mat.stat-mech] 6 Jan 2014 Determination of Forces from a Potential in Molecular Dynamics (note) Bernard Monasse Mines-ParisTech, Cemef bernard.monasse@mines-paristech.fr January

More information

Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies

Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies Supporting Information for Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies Shuntaro Chiba, 1* Tadaomi Furuta, 2 and Seishi

More information

Damped harmonic motion

Damped harmonic motion Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,

More information

PROPERTIES OF SOLIDS SCH4U1

PROPERTIES OF SOLIDS SCH4U1 PROPERTIES OF SOLIDS SCH4U1 Intra vs. Intermolecular Bonds The properties of a substance are influenced by the force of attraction within and between the molecules. Intra vs. Intermolecular Bonds Intramolecular

More information

! 94

! 94 ! 94 4 : - : : / : : : : ( :) : : : - : / : / : : - 4 : -4 : : : : : -5 () ( ) : -6 : - - : : : () : : : :4 : -7. : : -8. (. : ( : -9 : ( ( ( (5 (4 4 : -0! : ( : ( :. : (. (. (. (4. ( ( ( : ( 4 : - : :

More information

Inelastic X ray Scattering

Inelastic X ray Scattering Inelastic X ray Scattering with mev energy resolution Tullio Scopigno University of Rome La Sapienza INFM - Center for Complex Dynamics in Structured Systems Theoretical background: the scattering cross

More information

Vacuum Technology and film growth. Diffusion Resistor

Vacuum Technology and film growth. Diffusion Resistor Vacuum Technology and film growth Poly Gate pmos Polycrystaline Silicon Source Gate p-channel Metal-Oxide-Semiconductor (MOSFET) Drain polysilicon n-si ion-implanted Diffusion Resistor Poly Si Resistor

More information

Liquid to Glass Transition

Liquid to Glass Transition Landscape Approach to Glass Transition and Relaxation (Lecture # 3, March 30) Liquid to Glass Transition Instructor: Prabhat Gupta The Ohio State University (gupta.3@osu.edu) PKGupta(OSU) Landscape #3

More information

Dynamics. capacities of. (A.I.1) (E b ), bond. E total stretching (A.I.2) (A.I.3) A.I.1. Ebond force (A.I.4) (A.I.5)

Dynamics. capacities of. (A.I.1) (E b ), bond. E total stretching (A.I.2) (A.I.3) A.I.1. Ebond force (A.I.4) (A.I.5) Thermodynamics of liquids: Standard d molar entropies and heat capacities of common solvents from 2PT Molecular Dynamics Tod A Pascal, Shiang-Tai Lin and William A Goddard III Supplementary Materials Appendix

More information

Letting p shows that {B t } t 0. Definition 0.5. For λ R let δ λ : A (V ) A (V ) be defined by. 1 = g (symmetric), and. 3. g

Letting p shows that {B t } t 0. Definition 0.5. For λ R let δ λ : A (V ) A (V ) be defined by. 1 = g (symmetric), and. 3. g 4 Contents.1 Lie group p variation results Suppose G, d) is a group equipped with a left invariant metric, i.e. Let a := d e, a), then d ca, cb) = d a, b) for all a, b, c G. d a, b) = d e, a 1 b ) = a

More information

The Periodic Table and Chemical Reactivity

The Periodic Table and Chemical Reactivity The and Chemical Reactivity Noble gases Less electronegative elements More electronegative elements Then what is electronegativity? The tendency of an atom to attract an electron (or electron density)

More information

arxiv: v1 [cond-mat.soft] 8 May 2013

arxiv: v1 [cond-mat.soft] 8 May 2013 The Rosenfeld-Tarazona expression for liquids specific heat: A numerical investigation of eighteen systems arxiv:1305.1776v1 [cond-mat.soft] 8 May 2013 Trond S. Ingebrigtsen, 1, Arno A. eldhorst, 1 Thomas

More information

MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing

MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,

More information

The simplified model now consists only of Eq. 5. Degrees of freedom for the simplified model: 2-1

The simplified model now consists only of Eq. 5. Degrees of freedom for the simplified model: 2-1 . a) Overall mass balance: d( ρv ) Energy balance: = w + w w () d V T Tref C = wc ( T Tref ) + wc( T Tref ) w C T Because ρ = constant and ( Tref ) V = V = constant, Eq. becomes: () w = + () w w b) From

More information

Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases

Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song

More information

Dimensionless Numbers

Dimensionless Numbers 1 06.10.2017, 09:49 Dimensionless Numbers A. Salih Dept. of Aerospace Engineering IIST, Thiruvananthapuram The nondimensionalization of the governing equations of fluid flow is important for both theoretical

More information

Chapter 2 Experimental sources of intermolecular potentials

Chapter 2 Experimental sources of intermolecular potentials Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical

More information

CH352 Assignment 3: Due Thursday, 27 April 2017

CH352 Assignment 3: Due Thursday, 27 April 2017 CH352 Assignment 3: Due Thursday, 27 April 2017 Prof. David Coker Thursday, 20 April 2017 Q1 Adiabatic quasi-static volume and temperature changes in ideal gases In the last assignment you showed that

More information

Long-time simulations of viscous liquids from strong correlations to crystallization

Long-time simulations of viscous liquids from strong correlations to crystallization phd thesis Long-time simulations of viscous liquids from strong correlations to crystallization Ulf Rørbæk Pedersen Supervisor: Thomas B. Schrøder Danish National Research Foundation Centre Glass and Time,

More information

Chem Hughbanks Exam 3, April 19, 2012

Chem Hughbanks Exam 3, April 19, 2012 Chem 107 - Hughbanks Exam 3, April 19, 2012 Name (Print) UIN # Section 503 Exam 3, Version # A On the last page of this exam, you ve been given a periodic table and some physical constants. You ll probably

More information

Wiener Measure and Brownian Motion

Wiener Measure and Brownian Motion Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u

More information

Pressure densification of a simple liquid

Pressure densification of a simple liquid 1 Pressure densification of a simple liquid R. Casalini and C.M. Roland Naval Research Laboratory, Chemistry Division, Washington, DC 20375-53342 (Aug 9, 2017) ABSTRACT The magnitude of the high frequency,

More information

Quiz 1. Introduction to Polymers

Quiz 1. Introduction to Polymers 100406 Quiz 1. Introduction to Polymers 1) Polymers are different than low-molecular weight oligomers. For example an oligomeric polyethylene is wax, oligomeric polystyrene is similar to naphthalene (moth

More information

mean-square displacement in some fixed axis direction i, (b.x; (t)). In terins of the canonical equilibrium probability,

mean-square displacement in some fixed axis direction i, (b.x; (t)). In terins of the canonical equilibrium probability, PHYSICAL REVIEW E VOLUME 5, NUMBER 3 SEPTEMBER 1995 Universal time dependence of the mean-square displacement in extremely rugged energy landscapes with equal minima Jeppe C. Dyre and Jacob M. Jacobsen

More information

Chem Hughbanks Exam 3A, Solutions

Chem Hughbanks Exam 3A, Solutions Chem 107 - Hughbanks Exam 3A, Solutions Name (Print) UIN # Section 503 Exam 3, Version # A On the last page of this exam, you ve been given a periodic table and some physical constants. You ll probably

More information

MD simulation of methane in nanochannels

MD simulation of methane in nanochannels MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall

More information

Thermal & Statistical Physics Study Questions for the Spring 2018 Department Exam December 6, 2017

Thermal & Statistical Physics Study Questions for the Spring 2018 Department Exam December 6, 2017 Thermal & Statistical Physics Study Questions for the Spring 018 Department Exam December 6, 017 1. a. Define the chemical potential. Show that two systems are in diffusive equilibrium if 1. You may start

More information

States of matter. Chapter 11. Kinetic Molecular Theory of Liquids and Solids. Kinetic Molecular Theory of Solids Intermolecular Forces

States of matter. Chapter 11. Kinetic Molecular Theory of Liquids and Solids. Kinetic Molecular Theory of Solids Intermolecular Forces States of matter Chapter 11 Intermolecular Forces Liquids and Solids By changing the T and P, any matter can exist as solid, liquid or gas. Forces of attraction determine physical state Phase homogeneous

More information

Stoichiometry. Mole Concept. Balancing Chemical Equations

Stoichiometry. Mole Concept. Balancing Chemical Equations Stoichiometry The story so far The structure of an atom protons, neutrons & electrons Electron structure & the Periodic Table Shapes of electron orbitals (Quantum Numbers) Essential and toxic elements

More information

Ionic Bonding. Example: Atomic Radius: Na (r = 0.192nm) Cl (r = 0.099nm) Ionic Radius : Na (r = 0.095nm) Cl (r = 0.181nm)

Ionic Bonding. Example: Atomic Radius: Na (r = 0.192nm) Cl (r = 0.099nm) Ionic Radius : Na (r = 0.095nm) Cl (r = 0.181nm) Ionic Bonding Ion: an atom or molecule that gains or loses electrons (acquires an electrical charge). Atoms form cations (+charge), when they lose electrons, or anions (- charge), when they gain electrons.

More information

a) 1.3 x 10 3 atm b) 2.44 atm c) 8.35 atm d) 4.21 x 10-3 atm e) 86.5 atm

a) 1.3 x 10 3 atm b) 2.44 atm c) 8.35 atm d) 4.21 x 10-3 atm e) 86.5 atm 1. (6 pts) A sample of gas with a volume of 750 ml exerts a pressure of 756 mm Hg at 30.0 0 C. What pressure (atm) will the sample exert when it is compressed to 250 ml and cooled to -25.0 0 C? a) 1.3

More information

Thermal impedance calculations for micro-electronic structures. Bjorn Vermeersch Seminar Fysische Elektronica March 9, 2006

Thermal impedance calculations for micro-electronic structures. Bjorn Vermeersch Seminar Fysische Elektronica March 9, 2006 Thermal impedance calculations for micro-electronic structures Bjorn Vermeersch Seminar Fysische Elektronica March 9, 2006 Overview Introduction Thermal characterization Numerical calculations Influence

More information

Fast and slow dynamics of hydrogen bonds in liquid water. Abstract

Fast and slow dynamics of hydrogen bonds in liquid water. Abstract Fast and slow dynamics of hydrogen bonds in liquid water Francis W. Starr 1, Johannes K. Nielsen 1,2 & H. Eugene Stanley 1 1 Center for Polymer Studies, Center for Computational Science, and Department

More information

From Atoms to Materials: Predictive Theory and Simulations

From Atoms to Materials: Predictive Theory and Simulations From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology

More information

Plastics Testing and Characterization Industrial Applications

Plastics Testing and Characterization Industrial Applications Alberto Naranjo C., Maria del Pilar Noriega E., Tim A. Osswald, Alejandro Rojan, Juan Diego Sierra M. Plastics Testing and Characterization Industrial Applications ISBN-10: 3-446-41315-4 ISBN-13: 978-3-446-41315-3

More information

A Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential

More information

Nonequilibrium transitions in glassy flows II. Peter Schall University of Amsterdam

Nonequilibrium transitions in glassy flows II. Peter Schall University of Amsterdam Nonequilibrium transitions in glassy flows II Peter Schall University of Amsterdam Flow of glasses 1st Lecture Glass Phenomenology Basic concepts: Free volume, Dynamic correlations Apply Shear Glass? Soft

More information