Interactions Chemistry
Interactions Chemistry
1. Motivation
As technology moves more towards miniaturization in novel product
developments, it is imperative to integrate interfacial interactions into design strategies.
Consequently, interfacial forces have to be explored. Interfacial forces are on the order of
10-6 to 10-10 N, strong enough, for instance, to freeze gears in micro-electrical mechanical
systems (MEMS), to affect the stability of colloidal system, or to wipe out magnetically
stored data information in hard drives. There are multiple ways of exploring the strength
of interfacial interactions, one of which is by force spectroscopy, also known as force-
displacement (FD) analysis. The FD analysis involves a nanometer sharp scanning force
microscopy (SFM) tip that is moved relative to the sample surface in nanometer to
micrometer per second, as illustrated at end of this document in Figure 10. Before we
discuss FD analysis, we first discuss interaction forces, particularly weak interactions
between molecules and solids.
1
(VdW) interactions. They arise from dipole-dipole interactions. Both hydrogen and VdW
interactions can be responsible for cooperation and structuring in fluidic systems, but are
also strong enough to build up condensed phases. Following is a description of these
short range forces:
A. Ionic Bonds: These are simple Coulombic forces, which are a result of electron
transfer. For example in lithium fluoride, lithium transfers its 2s electron to the
fluorine 2p state. Consequently the shells of the atoms are filled up, but the
lithium has a net positive charge and the Flourine has a net negative charge.
These ions attract each other by Coulombic interaction which stabilizes the ionic
crystal in the rock-salt structure.
B. Covalent Bond: The standard example for a covalent bond is the hydrogen
molecule. When the wave-function overlap is considerable, the electrons of the
hydrogen atoms will be indistinguishable. The total energy will be decreased by
the “exchange energy”, which causes the attractive force. The characteristic
property of covalent bonds is a concentration of the electron charge density
between two nuclei. The force is strongly directed and falls off within a few
Ǻngstroms.
C. Metallic Bonds and Interaction: The strong metallic bonds are only observed
when the atoms are condensed in a crystal. They originate from the free valence
electron sea which holds together the ionic core. A similar effect is observed
when two metallic surfaces approach each other. The electron clouds have the
tendency to spread out in order to minimize the surface energy. Thus a strong
exponentially decreasing, attractive interaction is observed.
D. Dipole Interactions:
D.1. Hydrogen Bond Interaction: Strong type of directional dipole-dipole
interaction
D.2. Van der Waals Interaction: The relevance of VdW interactions goes beyond
of building up matter (e.g., Van der Waals organic crystals (Naphthalene)).
Because of their “medium” range interaction length of a few Ǻngstroms to
hundreds of Ǻngstroms, VdW forces are significant in fluidic systems (e.g,
colloidal fluids), and for adhesion between microscopic bodies. VdW forces
can be divided into three groups:
o Dipole-dipole force: Molecules having permanent dipoles will interact by
dipole-dipole interaction.
o Dipole-induced dipole forces: The field of a permanent dipole induces a
dipole in a non-polar atom or molecule.
o Dispersion force: Due to charge fluctuations of the atoms there is an
instantaneous displacement of the center of positive charge against the
center of the negative charge. Thus, at a certain moment, a dipole exists
and induces a dipole in another atom. Therefore non-polar atoms (e.g.
neon) or molecules attract each other.
2
Force Spectroscopy Analysis
3
between two atoms (or molecules) for φ(r) = 0. The LJ potential is depicted in Figure 4..
4. Surface Forces
The integral form of interaction forces between surfaces of macroscopic bodies through
a third medium (e.g., vacuum and vapor) are called surfaces forces. To apply the VdW
formalism to macroscopic bodies, one has to integrate the point interaction form
presented above. Consequently, the dipole-dipole interaction strength C but also the
exponent of the distance dependence become geometry dependent. For instance, while
for point-point particles the exponent is -6, it is -1 and -3 for macroscopic sphere-sphere
and sphere-plane interactions, respectively. Thus, while, VdW point particle interactions
are very short ranged (~1/r6), macroscopic VdW interactions are long ranged (e.g.,
sphere-sphere: ~1/D, where D represents the shortest distance between the two
macroscopic objects). Table 3 provides a list of geometry dependent non-retarded VdW
interaction strengths and exponents.
In vacuum, the main contributors to long-range surface interactions are the Van
der Waals and electromagnetic interactions. At separation distance < 2 nm one might also
have to consider short range retardation due to covalent or metallic bonding forces. Van
der Waals and electromagnetic interactions can be both attractive or repulsive. In the case
of a vapor environment as the third medium (e.g., atmospheric air containing water and
organic molecules), one also has to consider modifications by the vapor due to surface
adsorption or interaction shielding. This can lead to force modification or additional
forces such as the strong attractive capillary forces.
The SFM tip-sample interaction potential W are typically modeled as a sphere-
plane interaction, i.e.,
− AR
W( D ) = (2a)
6D
with the force
/W − AR
− = F( D ) = (2b)
/D 6D 2
where R is the radius of curvature of the tip, and D is the distance between the tip and
the plane. The interaction constant A, is called the Hamaker constant, defined as
A = π2Cρ1ρ2, with the interaction parameter of the point-point interaction C, and the
number density of the molecules in both solids ρi (i = 1,2). The Hamaker constant is
based on the mean-field Lifshitz theory. If known, A provides the means to deduce the
material specific (i.e., geometry independent) interaction parameter C. Typical values
for A, C and ρ are provided in Table 2. Table 3 summarizes the Van der Waals
interaction potential for various geometries.
4
Force Spectroscopy Analysis
− πCρ
Atom-Surface
6D 3
− A R1 R2
Sphere-Sphere
6 D ( R1 + R2 )
Plane-Sphere
− AR
6D
Body Interaction
12
AL ⎛ R1 R2 ⎞
Two Cylinders ⎜⎜ ⎟⎟
12 2 D 3 2 ⎝ ( R1 + R2 ) ⎠
− A R1 R2
Two Crossed Cylinders
6D
−A
Plane-Plane
12πD 2
5. Hamaker Constant
Originally the Hamaker constant was determined based on a purely additive
method in which polarization was ignored. The Lifshitz theory has overcome the problem
of additivity. It is a continuum theory which neglects the atomic structure. The input
parameters are the dielectric constants, ε, and refractive indices, n. The Hamaker constant
for two macroscopic phases 1 and 2 interacting across a medium 3 is approximated as:
3 ⎛ ε 1 − ε 3 ⎞⎛ ε 2 − ε 3 ⎞ 3hυ e
2 2 2
(
n1 − n3 n2 − n3 )(
2
)
){ ( )}
A ≈ kT ⎜⎜ ⎟⎜ ⎟+ (3)
1 (3 2 )(
4 ⎝ ε 1 + ε 3 ⎟⎠⎜⎝ ε 2 + ε 3 ⎟⎠ 8 2 n 2 + n 2 n 2 + n 2
3 n1
2
+ n 3
2
+ n 2
2
+ n 3) (
2
15
where νe is the absorption frequency (e.g., for H2O: νe = 3 x 10 Hz). Table 4 provides
non-retarded Hamaker constants determined with the Lifshitz theory (eq. 3).
In general, there is an attractive VDW interaction for A > 0, and the two
macroscopic phases are attracted to each other. In cases where it is desired to have
repulsive forces, the medium must have dielectric properties which are intermediate to
the macroscopic phases.
5
Table 4: Non-retarded Hamaker constants for two interacting media across a vacuum
(air) (Source: intermolecular & Surface Forces, J. Israelachvili, Academic Press) 3
a
Dielectric constant Refractive Index Absorption frequency Hamaker Constant
ε n ν Amedium/air/medium
15 -1 -20
Medium (10 s ) (10 )
6
Force Spectroscopy Analysis
which was obtained by pairwise summation of energies between all the atoms of medium
1 with medium 2. The summation of atom interactions within the same medium have
been neglected, which yields additional energy terms, i.e.,
A
W2 = − const . + 2
(6)
12πDo
consisting of a bulk cohesive energy term (assumed to be constant), and an energy term
related to unsaturated "bonds" at the two surfaces in contact (i.e., D = Do). Notice that
contact cannot be defined as D = 0 due to molecular repulsive forces. Do is called the
"cutoff distance". Hence the total energy of two planar surfaces at a distance D ≥ Do apart
is (neglecting the bulk cohesive energy)
A ⎛ 1 1 ⎞ A ⎛ Do ⎞
2
W = W1 + W2 = − ⎜⎜ 2 − 2 ⎟⎟ = ⎜ 1 − ⎟. (7)
12π ⎝ Do D ⎠ 12πDo 2 ⎜⎝ D 2 ⎟⎠
In contact (i.e., D=Do) W = 0. In the case of isolated surfaces, i.e., D = ∞,
A
W= . (8)
12πDo 2
Thus, in order to separate the two surfaces one has to overcome the energy difference
A
ΔW=W(Do)- W(D=∞)=- 2
, (9)
12πDo
which corresponds to the adhesive energy per unit area of W''=2γ12. Hence, the interfacial
energy can expressed as function of the Hamaker constant and the cutoff distance:
A
γ 12 = , (10)
24πDo2
24π ⎜ ⎟
⎝ 2.5 ⎠
For σ = 0.4 nm and γ12 = A/(24πDo2) it follows that Do = 0.16 nm. Do = 0.16 nm is a
remarkable "universal constant" yielding values for surface energies γ that are in good
agreement with experiments as shown in the Table 5.
7
Table 5: Surface energies based on Lifshitz theory and experimental values.(Source:
intermolecular & Surface Forces, J. Israelachvili, Academic Press) 3
Lifshiz Theory
Material A Experimental*
A/24 π Do2
-20
(10 ) {Do=0.165nm} (20oC)
8
Force Spectroscopy Analysis
γ LV
rK = , (13)
⎛ p⎞
RT log ⎜⎜ ⎟⎟
⎝ ps ⎠
where γL is the surface tension, R the gas constant, T the temperature, V the mol volume
and p/ps the relative vapor pressure (relative humidity for water). The surface tension γL
of water is 0.074N/m (T=20°C) leading to a critical Van der Waals distance of water of
γLV/RT = 5.4 Å. Consequentially, we obtain for p/ps=0.9 a Kelvin radius of 100 Å. At
small vapor pressures, the Kelvin radius gets comparable to the dimensions of the
molecules, and thus, the Kelvin equation breaks down.
The meniscus forces between two objects of spherical and planar geometry can be
approximated, for D « R, as:
4πRγ L cos Θ
F R >> D = , (14)
(1 + D / d )
where R is the radius of the sphere, d the length of PQ , see Figure 5, D the distance
between the sphere and the plate, and θ the meniscus contact angle.
Figure 5: Capillary meniscus between two objects of spherical and planar geometry
R >> d
The maximum force, found at at D = 0 (contact), is Fmax = 4πRγ cos θ . While this
expression estimates the capillary forces of relatively large spheres fairly accurately, the
capillary forces of highly wetted nanoscale spheres requires a geometrical factor K.
(1 + cos φ ) 2
K= (15)
4 ⋅ cos φ
where φ is the filling angle.
9
minimum thickness requirement of a liquid precursor film for spreading. The height of
the precursor film can not drop below a certain minimum, e, which is
1/ 2 1/ 2
⎛γ ⎞ ⎛ A ⎞
e = a0 ⎜ ⎟ ; a0 = ⎜⎜ ⎟⎟ ; S = γ SO − γ SL − γ , (15)
⎝S⎠ ⎝ 6πγ ⎠
where ao is a molecular length, S the spreading coefficient, A the Hamaker constant, γSO
the solid-vacuum interfacial energy, and γSL the solid-liquid interfacial energy. As the
water vapor film thickness depends on the RH (i.e., p/ps), a relative humidity smaller than
40 % does not provide a minimum thickness for the formation of a capillary neck. Once a
capillary neck forms between the SFM tip and the substrate surfaces, the pull-off force
increases suddenly, and provides over regime II a pull-off force that contains both, VdW
and capillary forces. VdW forces from SFM FD analysis as determined, for instance,
from regime I, see Figure 6(b), are on the order of 1-10 nN. The capillary force, on the
other hand, is on the order of up to 100 nN, and thus, dominates VdW interactions in
regime II.
Pull-off Force
I II III
Relative Humidity
Figure 6(a): Generic sketch of the functional Figure 6(b): Pull-off force vs. RH measured
relationship between the pull-off force and the between a hydrophilic silicon oxide SFM tip
relative humidity (RH). Regimes I, II and III and a ultra-smooth silicon oxide wafer. ●
represent the van der Waals regime, mixed van measured for increasing RH, ▼ measured for
der Waals – capillary regime, and capillary decreasing RH. .1
regime decreased by repulsive forces,
respectively.
In the high RH regime (III) the pull-off force decreases with increasing RH for
hydrophilic counter-surfaces. At such high humidity, the water vapor film thickness
dimensions exceeds the contact size ( asperity flooding), and the effect of the capillary
interface decreases.
10
Force Spectroscopy Analysis
As shown in Figure 7, the forces, Fstv, Fstw, and Fcap, are components of the measured
pull-off force Fmea. When the relative humidity is below the transition regime, i.e., ϕ <
ϕ0, the Fmea consists Fstv only, represents the lower limit of the sigmoidal fit. Above the
transition regime, Fmea is the sum of Fstw and Fcap, represents the upper limit of the
sigmoidal fit.
The Fstw and Fstv can be expressed by assuming that the contact is between an
incompressible sphere and a hard flat surface, i.e. Bradley’s model (see page 24),
Fstv = 2π ⋅ R ⋅ Wstv or Fstw = 2π ⋅ R ⋅ Wstw (17a or 17b)
where R is the sphere radius (i.e. SFM tip radius), W is the work of adhesion which is
expressed as,
Wijm = γ im + γ jm + γ ij (18)
where γ is the interfacial energies of the two materials, and i, j, m represents solid i, solid
j, and the medium m in which the contact take place, respectively. If the contact is
between two solids with the same material, i.e., i = j, Eq. (18) reduces to Wijm = 2γ im .
In order to determine the tip radius R, Eq. (17a) is solved for R using experimentally
determined Fstv. The R value is then used to determine Fstw through Eq. (17b), and Fcap is
deduced. Employing the geometric coefficient K for the capillary force equation,
(1 + cos φ ) 2
Fcap = 4π ⋅ R ⋅ γ water cosθ ⋅ (19)
4 ⋅ cos φ
11
the filling angle φ can be deduced. For example, the result obtained by He et. al.1 on the
silicon wafer surface, was analyzed using this model. Using the value of γ SiO/air 100
mJ/m25, γSiO/water 24.5 mJ/m25, γ water 72.8 mJ/m2, and the contact angle Θ of 0 °, the tip
radius R and the filling angle φ was determined to be 8.7 nm and 85.6 ° respectively.
Figure 8: The contact angle measurement. The contact angle Θ is the measure of
hydrophobicity (wettability). Left: hydrophilic surface. Right: hydrophobic surface.
12
Force Spectroscopy Analysis
F(r)
A0 D D0
A
C
B0
C0
B
(a) (b)
Figure 10: (a) The actual path taken by a SFM cantilever. The inset illustrates the snap-in
instability at A0 where the second derivative of the interaction potential exceeds the spring constant
of the cantilever. (b) Typical force distance curve. D = displacement, F(D) = force.
1. Line 1-A0: The probe and sample are not in contact but the tip is moving toward
the sample.
2. Line A0-B0: Jump into contact caused by the attractive van der Waals forces
outweighing the force of the cantilever spring between the tip and the sample
causing the cantilever to bend.
3. Line B0-2: Shows upward deflection of the cantilever in response to the sample
motion after they are in contact. The shape of the segment indicates whether the
sample is deforming in response to the force from the cantilever. (may not always
be straight) If the sample is assumed to be a hard surface, the slope of this line is
the sensitivity (springiness) of the cantilever.
13
4. Line 2-C0: As the tip moves away, the slope follows the slope of line B0-2 closely.
If line 2-C0 is parallel to line B0-2, no additional information can be determined.
However, if there is a difference in the in and out-going curves (hysteresis) gives
information on the plastic deformation of the sample. Once it passes point 2’, the
cantilever begins to deflect downward due to adhesive forces..
5. Line C0-D0: A jump out of contact occurs when the cantilever force exceeds the
adhesive forces.
The jump out of contact distance will always be greater than the jump into contact
distance because of few possible causes are:
a. During contact, some adhesive bonds are created.
b. During contact, the sample buckles and “wraps” around the tip, increasing the
contact area.
c. Hysteresis contributions
d. Capillary forces exerted by contaminants such as water.
FD analysis is widely used for adhesion and force interaction studies. Recently
biological materials have been studied by force spectroscopy, such as adsorption strength
of proteins on a substrate and folding/unfolding energy of DNAs.
References
1
M. He, A. Blum, D. E. Aston, C. Buenviaje, and R. M. Overney, Journal of chemical
physics 114 (3), 1355 (2001).
2
H. Hertz, J. Reine und Angewandte Mathematik 92, 156 (1882).
3
J. N. Israelachvili, Intermolecular and surface forces. (Academic Press, London,
1992).
4
D. L. Sedin and K. L. Rowlen, Analytical Chemistry 72, 2183 (2000).
5
C. Ziebert and K. H. Zum Gahr, Tribology Letters 17 (4), 901 (2004).
6
J. A. Greenwood, Proc. R. Soc. Lond. A 453, 1277 (1997).
Recommended Reading
Contact Mechanics by K. L. Johnson, Cambridge University Press, Cambridge, 1985.
Gases, liquids and solids and other states of matter by D. Tabor, Cambridge Univ.
Press, Cambridge, 3rd ed. 2000.
Intermolecular and surface forces by J. N. Israelachvili, Academic Press, London, 1992.
Nanoscience: Friction and Rheology on the Nanometer Scale by E. Meyer, R. M.
Overney, K. Dransfeld, and T. Gyalog, World Scientific Publ., Singapore, 1998.
14