Room BI 115 Lecture: Monday, Wednesday Friday 2-3pm. Ch121a Atomic Level Simulations of Materials and Molecules. Lecture 6 and 7, April 16 and 21, 2013 MD3: vibrations. Lecture 6 Presented by Jason Crowley. William A. Goddard III, [email protected] - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Homework and Research ProjectFirst 5 weeks: The homework each week uses generally available computer software implementing the basic methods on applications aimed at exposing the students to understanding how to use atomistic simulations to solve problems. Each calculation requires making decisions on the specific approaches and parameters relevant and how to analyze the results. Midterm: each student submits proposal for a project using the methods of Ch121a to solve a research problem that can be completed in the final 5 weeks.The homework for the last 5 weeks is to turn in a one page report on progress with the projectThe final is a research report describing the calculations and conclusions
Vibration of moleculesClassical and quantum harmonic oscillatorsInternal vibrations and normal modesRotations and selection rulesExperimentally probing the vibrationsDipoles and polarizabilitiesIR and Raman spectraSelection rulesThermodynamics of moleculesDefinition of functionsRelationship to normal modesDeviations from ideal classical behavior
Starting with an atom inside a molecule at equilibrium, we can expand its potential energy as a power series. The second order term gives the local spring constantWe conceptualize molecular vibrations as coupled quantum mechanical harmonic oscillators (which have constant differences between energy levels)Including Anharmonicity in the interactions, the energy levels become closer with higher energySome (but not all) of the vibrational modes of molecules interact with or emit photons This provides a spectroscopic fingerprint to characterize the molecule
Consider a one dimensional spring with equilibrium length xe which is fixed at one end with a mass M at the other. If we extend the spring to some new distance x and let go, it will oscillate with some frequency, w, which is related to the M and spring constant k.To determine the relation we solve Newton’s equation M (d2x/dt2) = F = -k (x-xe)Assume x-x0=d = A cos(wt) then –Mw2 Acos(wt) = -k A cos(wt) Hence –Mw2 = -k or w = Sqrt(k/M). Stiffer force constant k higher w and higher M lower w
CM = Center of mass Fix Rcm = (M1R1 + M2R2)/(M1+ M2) = 0Relative coordinate R=(R2-R1)Then Pcm = (M1+ M2)*Vcm = 0 And P2 = - P1
Thus KE = ½ P12/M1 + ½ P2
2/M2 = ½ P12/m
Where 1/m = (1/M1 + 1/M2) or m = M1M2/(M1+ M2) Is the reduced mass.Get w = sqrt(k/m). Thus we can treat the diatomic molecule as a simple mass on a spring but with a reduced mass, m
Fk = -(∂E(Rnew)/∂Rk) = -(∂E/∂Rk)0 - Sm (∂2E/∂Rk∂Rm) (dR)m
Where we have neglected terms of order d2. Writing the 2nd derivatives as a matrix (the Hessian) Hkm = (∂2E/∂Rk∂Rm) and setting (∂E/∂Rk)0 = 0, we get
Fk = - Sm Hkm (dR)m = Mk (∂2Rk/∂t2) To find the normal modes we write (dR)m = Am cos wt leading to
Mk(∂2Rk/∂t2) = Mk w2 (Ak cos wt) = Sm Hkm (Amcos wt)
Here the coefficient of cos wt must be {Mk w2 Ak - Sm Hkm Am}=0
There are 3N degrees of freedom (dof) which we collect together into the 3N vector, Rk where k=1,2..3NThe interactions then lead to 3N net forces, Fk = -(∂E(Rnew)/∂Rk) all of which are zero at equilibrium, R0
Now consider that every particle is moved a small amount leading to a 3N distortion vector, (dR)m = Rnew – R0
Solving for the Vibrational modesThe normal modes satisfy {Mk w2 Ak - Sm Hkm Am}=0To solve this we mass weight the coordinates as Bk = sqrt(Mk)Ak
leading toSqrt(Mk) w2 Bk - Sm Hkm [1/sqrt(Mm)]Bm}=0 leading to
Sm Gkm Bm = wk2 Bk where Gkm = Hkm/sqrt(MkMm)
G is referred to as the reduced HessianFor M degrees of freedom this has M eigenstatesSm Gkm Bmp = dkp Bk (w2)p
where the M eigenvalues (w2)p are the squares of the vibrational energies.If the Hessian includes the 6 translation and rotation modes then there will be 6 zero frequency modes, wp = 0
If the point of interest were a saddle point rather than a minimum, G would have one negative eigenvalue, (w2)p = - A2 where A is a positive number This leads to an imaginary frequency, wp = iA ,
For practical simulationsWe can obtain reasonably accurate vibrational modes from just the classical harmonic oscillators, usually within a few %N atoms => 3N degrees of freedomHowever, there are 3 degrees for translation, n = 0 3 degrees for rotation for non-linear molecules, n = 0 2 degrees if linearThe rest are vibrational modes
•EM energy absorbed by interatomic bonds in organic compounds•frequencies between 4000 and 400 cm-1 (wavenumbers) •Useful for resolving molecular vibrations
4 independent CH bonds 4 CH stretch modes, by symmetry one is triply degenerate6 possible angle terms HCH 5 HCH modes, one doubly degenate, on triply deg.Reason only 5 linearly independent HCH
IRVibrations at same frequency as radiationTo be observable, there must be a finite dipole derivativeThus homonuclear diatomic molecule (O2 , N2 ,etc.) does not lead to IR absorption or emission. Raman spectroscopy is complementary to IR spectroscopy.radiation at some frequency, n, is scattered by the molecule to frequency, n’, shifted observed frequency shifts are related to vibrational modes in the molecule IR and Raman have symmetry based selection rules that specify active or inactive modes
Dx=0 Dy=0 Dz=Dz •center of mass rotation (nonlinear molecules) Dx=0 Dy=-cDqx Dz=bDqx Dx= cDqy Dy=0 Dz=-aDqy Dx= -bDqz Dy=aDqx Dz=0•linear molecules have only 2 rotational degrees of freedom•The translational and rotational degrees of freedom can be removed beforehand by using internal coordinates or by transforming to a new coordinate system in which these 6 modes are separated out
Translation and Rotation Modes
E is a constant dE/dx = 0 d2E/dx2 = 0Thus the eigenmode l=0
E is a constant dE/dx = 0 d2E/dx2 = 0Thus the eigenmode l=0
The moment of inertia about an axis q is defined as
)(2 qxmIk
kkqq = xk(q) is the perpendicular distance to the axis q
Can also define a moment of inertia tensor where (just replace the mass density with point masses and the integral with a summation. Diagonalization of this matrix gives the principle moments of inertia!
For rotationsWavefunctions are spherical harmonicsProject the dipole and polarizability due to rotationIt can be shown that for IRDelta J changes by +/- 1Delta MJ changes by 0 or +/-1Delta K does not changeFor RamanDelta J could be 1 or 2Delta K = 0But for K=0, delta J cannot be +/- 1
• Phonons are the normal modes of lattice vibrations (thermal + zero point energy)
• When a photon absorbs/emits a single phonon, momentum and energy conservation the photon gains/loses the energy and the crystal momentum of the phonon. – q ~ q` => K = 0– The process is called anti-Stokes for absorption and
Stokes for emission.– Alternatively, one could look at the process as a
Doppler shift in the incident photon caused by a first order Bragg reflection off the phonon with group velocity v = (ω/ k)*k
•The external EM field is monochromatic•Dipole moment of the system•Interaction between the field and the molecules•Probability for a transition from the state i to the state f (the Golden Rule)
•Rate of energy loss from the radiation to the system
•The flux of the incident radiation
)ωcos(ε)( 0 tEtE =
==
==N
j
ij
ij
n
ii qr
1i
1
μ DipoleMolecular μμ Dipole Total
)(μ)( tEt=
)]ωω(δ)ωω(δ[|με|2
π)ω(2
20 = fififi ifEP
fffi ωωω =
)ω(ωρ)ω( fifii f
irad PE =
20π8
EcnS =c: speed of lightn: index of refraction of the medium
Anharmonicity – bonds do eventually dissociateCoriolis forces
Interaction between vibration and rotationInversion doublingIdentical atoms on rotation – need to obey the Pauli PrincipleTotal wavefunction symmetric for Boson and antisymmetric for Fermion ),()1(),( ,, qq
How does a Molecule response to an oscillating external electric field (of frequency w)? Absorption of radiation via exciting to a higher energy state ħw ~ (Ef - Ei)
Using the vibrational modes: thermodynamicsIn QM and MM the Energy at minima = motionless state at 0KBUT, experiments are made at finite T, hence corrections are required to allow for rotational, translational and vibrational motion.The internal energy of the system: U(T)=Urot(T)+Utran(T)+Uvib(T)+Uvib(0)
32-AJB
Uvib T )= hn i
2 hn i
exp hn i KBT 1
i=1
Nmod
The vibrational frequencies ni) of the normal modes (Nmod) calculated from the eigenvalues λi) of the force-constant equivalent of Hessian matrix of second derivatives
Vibrational frequencies can be used to calculate entropies and free energies, or to compare with results of spectroscopic experiments
n i =λi
2
From equipartition theorem: Urot(T) = (3/2)KBT , Utran(T) = (3/2)KBT per molecule (except Urot(T)=KBT for linear molecules)BUT, vibrational energy levels are often only partially excited at room T, thus Uvib(T) requires knowledge of vibrational frequenciesUvib(T) = vibrational enthalpy @ T - vibrational enthalpy @ 0K
Validation of 2PT method for entropy and free energy for amino acid side chains
Tod A PascalRavi AbrolWilliam A Goddard IIIMSC 2013 Research Conference
TOD A. PASCAL MSC, CALTECH
General OutlineOverview of the 2PT method for
calculating accurate entropies and free energies
Application to the solvation thermodynamics of Amino Acid sidechains
TOD A. PASCAL MSC, CALTECH
Entropy
“Any method involving the notion of entropy, the very existence of which depends on the second law of thermodynamics, will doubtless seem to many far-fetched, and may repel beginners as obscure and difficult of comprehension.”
--Willard Gibbs, Graphical Methods in the Thermodynamics of Fluids (1873)
Rudolf Clausius - originator of the concept of "entropy".
Entropy: How to calculate it?• Test Particle Methods (insertion or deletion)
• Good for low density systems • Perturbation Methods • Thermodynamic integration, Thermodynamic perturbation• Applicable to most problems • Require long simulations to maintain “reversibility”
• Nonequilibrium Methods (Jarzinski’s equality)• Obtaining differential equilibrium properties from irreversible processes • Require multiple samplings to ensure good statistics
• Normal mode Methods• Good for gas and solids• Fast• Not applicable for liquids
References: Frenkel, D.; Smit, B. Understanding Molecular Simulation from Algorithms to Applications. Academic press: Ed., New York, 2002. McQuarrie, A. A. Statistical Mechanics. Harper & Row: Ed., New York, 1976. Jarzynski, C. Nonequilibrium Equality for Free Energy Differences. Phys. Rev. Lett. 1997, 78, 2690.
49
Remaining issue: experimental energies are free energies, need to calculate entropy
Free Energy, F = U – TS = − kBT ln Q(N,V,T)Kirkwood thermodynamic integration
J. G. Kirkwood. Statistical mechanics of fluid mixtures, J. Chem. Phys., 3:300-313,1935
Enormous computational cost required for complete sampling of the thermally relevant configurations of the system often makes this impractical for realistic systems. Additional complexities, choice of the appropriate approximation formalism or somewhat ad-hoc parameterization of the “reaction coordinate”
General approach to predict Entropy, S, and Free Energy
50
Solvation free energies amino acid side chainsPande and Shirts (JCP 122 134508 (2005) Thermodynamic
integration leads to accurate differential free energies
51
Solvation free energies amino acid side chainsPande and Shirts (JCP 122 134508 (2005) Thermodynamic
integration leads to accurate differential free energies
But costs 8.4 CPU-years on 2.8 GHz processor
52
Starting point: Get Density of states from MD (Velocity autocorrelation) Partition Function entropy
DoS(n) is the vibrational density of States
Problem: as n 0 get DoS ∞ unless DoS(0) = 0
zero
zero zero
DoS(n)
Calculate entropy from DoS(n)
Velocity autocorrelation function
53
S ( )Finite density of states at n =0Proportional to diffusion coefficient
Also strong anharmonicity at low frequencies
Problem with Liquids: S(0)≠0
DoS(n)
n
where D is the diffusion coefficientN=number of particlesM = mass
zero
zero zero
54
•Decompose liquid DoS(n) to a gas and a solid contribution•DoS(n) total = DoS(n) gas + DoS(n) solid •S(0) attributed to gas phase diffusion, fit to hard sphere theory•Gas component contains small n anharmonic effects•Solid component contains quantum effects
S ( )
Gasexponentialdecay
S ( )
SolidDebye crystalS(v) ~v2
S ( )
solid-likegas-like
Two-Phase Thermodynamics Model (2PT)
Total
=
)()()()(00
gP
gHOP
s WSdWSdP
=Property =The two-phase model for calculating thermodynamic properties of liquids from molecular dynamics: Validation for the phase diagram of Lennard-Jones fluids; Lin, Blanco, Goddard; JCP, 119:11792(2003) wag536
Computational SchemeRun a MD simulation(trajectory information saved)
Calculate VAC
Calculate DoS (FFT of VAC)Apply HO
approximationTo S()
1PT thermodynamic
predictions
Calculate S(0) and D
Solve for f
Determine Sgas(), Ssolid()
Apply HO statisticsTo Ssolid()
Apply HS statistics to Sgas()
2PT thermodynamic
predictions
TOD A. PASCAL MSC, CALTECH
Test case 1: Lennard Jones Fluid
●stable
●metastable
●unstable
T - diagram for Lennard Jones Fluid
0.6
1.0
1.4
1.8
0.0 0.4 0.8 1.2*
T*
Solid
LiquidGas
Supercritical Fluid
• Intermolecular potential
• Phase diagram
= 612 )()(4)(
rrrV e
Lennard-Jones Potential
r = V(r)
r
0-e
(T*=kT/e *=3)
– critical point
– triple point006.0304.0
006.0316.1*
*
=
=
c
cT
69.0* tpT
TOD A. PASCAL MSC, CALTECH
LJ VAC and DoS
0
10
20
30
40
50
60
70
80
90
100
0 20 40 60 80 100
frequency v(cm-1)
DoS
S(c
m)
=====
Density of StatesVelocity Autocorrelation
-800
-400
0
400
800
1200
1600
0 0.5 1 1.5 2 2.5 3time (ps)
C (t)
=====
gasliquidsolid
gasliquidsolid
)()0()( tvvtC = dtetCkT
S ti
=
πυ2)(lim 2)υ(
TOD A. PASCAL MSC, CALTECH
2PT DoS Decomposition
0
200
400
600
800
1000
1200
0 5 10 [cm-1]
S(
) [cm
]
0
5
10
15
20
25
30
0 50 100 150 [cm-1]
S(
) [cm
]
0
5
10
15
20
25
30
35
0 50 100 150 [cm-1]
S(
) [cm
]
• Examples
LJ gas liquid FCC solid
solid-likegas-likegas-like
solid-like
solid-likegas-like
TOD A. PASCAL MSC, CALTECH
Pressure and Energy for LJ GasPressure
-2
0
2
4
6
8
10
12
14
16
18
0 0.2 0.4 0.6 0.8 1 1.2 *
P *
T*=1.8T*=1.4T*=1.1T*=0.9MD
-7
-6
-5
-4
-3
-2
-1
0
1
2
3
0 0.2 0.4 0.6 0.8 1 1.2 *
E *
T*=1.8T*=1.4T*=1.1T*=0.9MD2PT(Q)
Pressures and MD Energies agree with EOS valuesQuantum Effect (ZPE) most significant for crystals (~2%)
Total Energy
TOD A. PASCAL MSC, CALTECH
Entropy
4
6
8
10
12
14
16
18
20
0 0.2 0.4 0.6 0.8 1 1.2 *
S *
T*=1.8
T*=1.4
T*=1.1
T*=0.9
1PT(Q)
4
6
8
10
12
14
16
18
20
0 0.2 0.4 0.6 0.8 1 1.2 *
S *
T*=1.8T*=1.4T*=1.1T*=0.92PT(Q)2PT(C)
1PT 2PT model
•Overestimate entropy for low density gases
•Underestimate entropy for liquids•Accurate for crystals
•Accurate for gas, liquid, and crystal•Accurate in metastable regime•Quantum Effects most important for crystals (~1.5%)
gas
liquid
crystal
TOD A. PASCAL MSC, CALTECH
Gibbs Free Energy
-30
-25
-20
-15
-10
-5
0
5
0 0.2 0.4 0.6 0.8 1 1.2 *
G *
T*=1.8T*=1.4T*=1.1T*=0.91PT(Q)
1PT 2PT model
•Underestimate free energy for low density gases
•overestimate entropy for liquids•Accurate for crystals
•Accurate for gas, liquid, and crystal•Accurate in metastable regime
gas
liquid
crystal
-30
-25
-20
-15
-10
-5
0
5
0 0.2 0.4 0.6 0.8 1 1.2 *
G *
T*=1.8T*=1.4T*=1.1T*=0.92PT(Q)2PT(C)
TOD A. PASCAL MSC, CALTECH
Why does 2PT work?
0
1
2
3
4
5
0 50 100 150 [cm-1]
WS(
)
0
5
10
15
20
25
30
0 50 100 150 [cm-1]
S(
) [cm
]
0
200
400
600
800
1000
1200
0 2 4 6 8 10 [cm-1]
S(
) [cm
]
HS fy = 0.036
QHO
HS fy = 0.309
CHO
)()()()(00
gP
gHOP
s WSdWSdP
=
gas
liquid • 1PT overestimates Wsgas for gas for
modes < 5 cm-1
• 1PT underestimates Wsgas for liquid for
modes between 5 and 100 cm-1
• 2PT properly corrects these errors
TOD A. PASCAL MSC, CALTECH
Convergence of 2PT
6.5
7.5
8.5
9.5
10.5
11.5
12.5
13.5
14.5
15.5
100 1000 10000 100000 1000000
MD steps
S*
2PT(Q)2PT(C) MBWR EOS
gas (*=0.05 T*=1.8)
liquid (*=0.85 T*=0.9)
• For gas, the entropy
converges to within 0.2% with
2500 MD steps (20 ps)
• For liquid, the entropy
converges to within 1.5% with
2500 MD steps (20 ps).
TOD A. PASCAL MSC, CALTECH
Melting and Solidification
0.6
0.8
1.0
1.2
1.4
1.6
1.8
2.0
0.00 0.40 0.80 1.20
T*
Simulation conditions
solid
supercritical fluid
3
4
5
6
7
8
0.80 1.20 1.60 2.00T*
S*liquid (EOS)solid (EOS)heatingcoolingclassical
Ent
ropy
solidmetatstable
unstablesupercritical
fluid
starting withfcc crystal
starting with amorphous liquid
• Initial amorphous structure is used in the cooling process
•The fluid remains amorphous in simulation even down to T*=0.8 (supercooled)
•The predicted entropy for the fluid and supercooled fluid agree well with EOS for LJ fluids
• Initial fcc crystal is used in the heating process
•The crystal appears stable in simulation even up to T*=1.8 (superheated)
•The predicted entropies for the crystal and superheated crystal agree well with EOS for LJ solids
TOD A. PASCAL MSC, CALTECH
0
2
4
6
8
10
12
14
16
0 500 1000 1500 2000 2500 3000
S_hs(v)[cm]S_s(v)[cm]Stot(v)[cm]
Water Power Spectrum (DoS) (25 ps, 1fs steps)
w (cm-1)
w)
Power spectrum for water at 300 K. The power spectrum is decomposed into a gas (diffusive) and a solid (fixed) spectra and their contributions added to yield the free energy of the liquid state.
0
2
4
6
8
10
12
14
16
1 10 100 1000 10000
S_hs(v)[cm]S_s(v)[cm]Stot(v)[cm]
Power spectrum for water at 300 K. The power spectrum is decomposed into a gas (diffusive) and a solid (fixed) spectra and their contributions added to yield the free energy of the liquid state.
TOD A. PASCAL MSC, CALTECHAccurate entropies of liquid systems Fourier Transform VAC
to get Density of Vibrational StatesGet Enthalpy and Free Energy using quantum partition function
Do 25 picosec MD, Extract Velocity Autocorrelation (VAC) Function
0
10
20
30
40
50
60
70
80
90
100
0 20 40 60 80 100
frequency v(cm-1)
DoS
S(c
m)
=====
gasliquidsolid
-800
-400
0
400
800
1200
1600
0 0.5 1 1.5 2 2.5 3time (ps)
C (t)
=====
gasliquidsolid
Works well for solids
TOD A. PASCAL MSC, CALTECH
0
2
4
6
8
10
12
14
16
0 500 1000 1500 2000 2500 3000
S_hs(v)[cm]S_s(v)[cm]Stot(v)[cm]
Water Power Spectrum (DoS) (25 ps, 1fs steps)
w (cm-1)
w)
Power spectrum for water at 300 K. The power spectrum is decomposed into a gas (diffusive) and a solid (fixed) spectra and their contributions added to yield the free energy of the liquid state.
0
2
4
6
8
10
12
14
16
1 10 100 1000 10000
S_hs(v)[cm]S_s(v)[cm]Stot(v)[cm]
Power spectrum for water at 300 K. The power spectrum is decomposed into a gas (diffusive) and a solid (fixed) spectra and their contributions added to yield the free energy of the liquid state.
Accuracy predicted entropy only limited by accuracy of force field
Thermodynamics of liquids: standard molar entropies and heat capacities of common solvents from 2PT molecular dynamics; Pascal; Lin; Goddard; Phys. Chem. Chem. Phys., 13: 169-181 (2011) wag897
TOD A. PASCAL MSC, CALTECH
How accurate is it? water model F3C SPC SPC/E TIP3P TIP4P-
TOD A. PASCAL MSC, CALTECHCalculating Interfacial Surface TensionDirect evaluation of surface free energy
Requires evaluation of the surface entropy Requires extensive simulation time for
convergence Can be approximated from potential of mean
force calculations – large uncertainties
𝛾=( 𝜕𝐺𝜕 𝐴 )𝑁 ,𝑇 ,𝑃
=(𝐺𝑠𝑢𝑟𝑓𝑎𝑐𝑒−𝐺𝑏𝑢𝑙𝑘
𝜕 𝐴 )𝑁 ,𝑇 ,𝑃
𝐺=𝐻− 𝑃𝑉¿𝑈 −𝑇𝑆−𝑃𝑉
surface
bulk
TOD A. PASCAL MSC, CALTECH
The Recipe: the 2PT Method
92
S ( )
Finite density of states at n =0Proportional to diffusion coefficient
Harmonic Approximation Entropy= ∞•Also strong anharmonicity at low frequencies
Liquid
S ( )
Gasexponentialdecay
S ( )
SolidDebye crystalS(v) ~v2New Model
2 phase theory (2PT )Liquid Solid + Gas •Two-Phase Thermodynamics
Model (2PT)•Decompose liquid S(v) to a gas and a solid contribution
•S(0) attributed to gas phase diffusion•Gas component contains anharmonic effects
•Solid component contains quantum effects
S ( )
solid-likegas-like
The two-phase model for calculating thermodynamic properties of liquids from molecular dynamics: Validation for the phase diagram of Lennard-Jones fluids; Lin, Blanco, Goddard; JCP, 119:11792(2003)
)/2βexp(-1)/2βexp()βexp()( n
e
hhq
n
QHO
==
TOD A. PASCAL MSC, CALTECH
Comparison of Surface Tensions calculated with both methods
Excellent agreement Systematic underestimation of experimental value: Deficiency of forcefield (SPC-Ew) used
TOD A. PASCAL MSC, CALTECH
The thermodynamic factors that determine the IFT
The decrease in surface tension with increasing temperature is an entropically driven process
TOD A. PASCAL MSC, CALTECH
Where does the increased entropy with temperature originate?
Both rotational and translational entropy increases with increasing temperature
TOD A. PASCAL MSC, CALTECH
What happens to the enthalpy with increasing temperature?
Hydrogen bonding in supercooled water is different from ambient water Subsurface water molecules are enthalpically stabilized Effect is reduced with increasing temperature
TOD A. PASCAL MSC, CALTECH
The thermodynamics stability of interfacial water
Surface tension effects propagate into subsurface Implications for propensity of ions at the interface
TOD A. PASCAL MSC, CALTECH
ConclusionsCan evaluate the surface tension of
liquids from direct evaluation of the surface energy
Reduction in surface tension with temperature is entropically driven
Sub-surface water molecules are preferentially stabilized enthalpically, especially for super-cooled water
TOD A. PASCAL MSC, CALTECH
Hydrogen bonding and entropy of confined water molecules under high pressure
Collaborators:Yousung Jung (KAIST – Korea)William A. Goddard III (Caltech)
TOD A. PASCAL MSC, CALTECHThe Hydrophobic effect: Importance/Consequences
Driving force in formation of Membranes Micelles Globular proteins
Hydrophobic DNA bases stack so as to exclude water molecules
"The antipathy of the paraffin chain for water is, however, frequently misunderstood. There is no question of actual repulsion between individual water molecules and paraffin chains, nor is there any very strong attraction of paraffin chains for one another. There is, however, a very strong attraction of water molecules for one another in comparison with which the paraffin-paraffin or paraffin-water attractions are slight." - G. S. Hartley 1936
Nonpolar groups (alkane chains) are hydrophobic
Polar groups (alcohols) are hydrophilic
What are the microscopic thermodynamic forces involved? Enthalpy and/or Entropy?
How does structure of water molecules change as you move away from the surface?
How do these microscopic differences impact the macroscopic bulk properties?
(11,11 and onwards) to entropic gain Translational entropy (due to reduced density
near hydrophobic interface) responsible for rest
TOD A. PASCAL MSC, CALTECH
Finally a complete picture of free energy!
water molecules inside the CNTs have lower free energies than bulk water
entropy dominates for tube diameters less than 1.0 nm (gas phase), the enthalpy dominates for tubes between 1.1 and 1.2 nm (ice phase), and both energies compensate for tubes larger than 1.4 nm (liquid phase)
TOD A. PASCAL MSC, CALTECH
The origin of the trends in rotational entropy
Water inside (6,6) and (7,7) resemble as gas – increase rotational entropy Water inside (8,8) and (9,9) resemble ice/water – decrease rotational entropy Water inside larger CNTs resemble bulk water – same rotational entropy
TOD A. PASCAL MSC, CALTECHThere is something special about water…
LJ liquid with same interactions with CNT as water (M3B) has unfavorable free energies
Thermodynamics recovered by including 3-body H-bond (mW)
SPC-E
qH: +0.4238 e-
M3B
mWStillinger-Weber 3-body
TOD A. PASCAL MSC, CALTECHSo waters flow into CNTs because…
Favorable local chemical potential inside CNT Lower free energy due to lower
enthalpy for ice-like CNTs but higher entropy for all others
Loss of hydrogen bonding inside tube overcome by increased entropy due to confinement