Chapter 1 Solubility Parameters - An Introduction
Chapter 1 Solubility Parameters - An Introduction
Chapter 1 Solubility Parameters - An Introduction
1 Solubility Parameters —
An Introduction
Charles M. Hansen
ABSTRACT
Solubility parameters have found their greatest use in the coatings industry to aid in the selection
of solvents. They are used in other industries, however, to predict compatibility of polymers,
chemical resistance, and permeation rates, and even to characterize the surfaces of pigments, fibers,
and fillers. Liquids with similar solubility parameters will be miscible, and polymers will dissolve
in solvents whose solubility parameters are not too different from their own. The basic principle
has been “like dissolves like.” More recently, this has been modified to “like seeks like,” as many
surface characterizations have also been made, and surfaces do not (usually) dissolve. Solubility
parameters help put numbers into this simple qualitative idea. This chapter describes the tools
commonly used in Hansen solubility parameter (HSP) studies. These include liquids used as energy
probes and computer programs to process data. The goal is to arrive at the HSP for interesting
materials either by calculation or, if necessary, by experiment and preferably with agreement
between the two.
INTRODUCTION
The solubility parameter has been used for many years to select solvents for coatings materials. A
lack of total success has stimulated further research. The skill with which solvents can be optimally
selected with respect to cost, solvency, workplace environment, external environment, evaporation
rate, flash point, etc., has improved over the years as a result of a series of improvements in the
solubility parameter concept and widespread use of computer techniques. Most commercial sup-
pliers of solvents have computer programs to help with solvent selection. One can now easily
predict how to dissolve a given polymer in a mixture of two solvents, neither of which can dissolve
the polymer by itself.
Unfortunately, this book cannot include discussion of all the significant efforts leading to our
present knowledge of the solubility parameters. An attempt is made to outline developments, provide
some background for a basic understanding, and give examples of uses in practice. The key factor
is to determine those affinities that the important components in a system have for each other. For
many products this means evaluating or estimating the relative affinities of solvents, polymers,
additives, pigment surfaces, filler surfaces, fiber surfaces, and substrates.
It is noteworthy that the concepts presented here have developed toward not just predicting
solubility that requires high affinity between solvent and solute, but for predicting affinities between
different polymers, leading to compatibility, and affinities to surfaces to improve dispersion and
adhesion. In these applications the solubility parameter has become a tool, using well-defined
liquids as energy probes, to measure the similarity, or lack of the same, of key components. Materials
with widely different chemical structures may be very close in affinities. Only those materials that
interact differently with different solvents can be characterized in this manner. It can be expected
that many inorganic materials, such as fillers, will not interact differently with these energy probes
as their energies are very much higher. An adsorbed layer of water on the high-energy surface can
also play an important role. Regardless of these concerns, it has been possible to characterize
pigments, both organic and inorganic, as well as fillers like barium sulfate, zinc oxide, etc., and
also inorganic fibers (see Chapter 7). Changing the surface energies by various treatments can lead
to a surface that can be characterized more readily and often interacts more strongly with given
organic solvents. When the same solvents that dissolve a polymeric binder are those which interact
most strongly with a surface, it can be expected that the binder and the surface have high affinity
for each other.
Solubility parameters are sometimes called cohesion energy parameters as they are derived
from the energy required to convert a liquid to a gas. The energy of vaporization is a direct measure
of the total (cohesive) energy holding the liquid’s molecules together. All types of bonds holding
the liquid together are broken by evaporation, and this has led to the concepts described in more
detail later. The term cohesion energy parameter is more appropriately used when referring to
surface phenomena.
δ = (E/V)1/2 (1.1)
Where V is the molar volume of the pure solvent, and E is its (measurable) energy of vaporization
(see Equation 1.15). The numerical value of the solubility parameter in MPa1/2 is 2.0455 times
larger than that in (cal/cm3)1/2. The solubility parameter is an important quantity for predicting
solubility relations, as can be seen from the following brief introduction.
Thermodynamics requires that the free energy of mixing must be zero or negative for the
solution process to occur spontaneously. The free energy change for the solution process is given
by the relation:
where ΔGM is the free energy of mixing, ΔHM is the heat of mixing, T is the absolute temperature,
and ΔSM is the entropy change in the mixing process.
Equation 1.3 gives the heat of mixing as proposed by Hildebrand and Scott:
The φ1 and φ2 are volume fractions of solvent and polymer, and VM is the volume of the mixture.
Equation 1.3 is not correct, and it has often been cited as a shortcoming of this theory in that only
positive heats of mixing are allowed. It has been shown by Patterson, Delmas, and coworkers that
ΔGMnoncomb is given by the right-hand side of Equation 1.3 and not ΔGM. This is discussed more in
Chapter 2. The correct relation is3–8:
The noncombinatorial free energy of solution, ΔGMnoncomb, includes all free energy effects other
than the combinatorial entropy of solution that results by simply mixing the components. Equation
1.4 is consistent with the Prigogine corresponding states theory (CST) of polymer solutions (see
Chapter 2) and can be differentiated to give expressions3,4 predicting both positive and negative
heats of mixing. Therefore, both positive and negative heats of mixing can be expected from
theoretical considerations and have been measured accordingly. It has been clearly shown that
solubility parameters can be used to predict both positive and negative heats of mixing. Previous
objections to the effect that only positive values are allowed in this theory are incorrect.
This discussion clearly demonstrates why the solubility parameter should be considered as a
free energy parameter. This is more in agreement with the use of the solubility parameter plots to
follow. These use solubility parameters as axes and have experimentally determined boundaries of
solubility defined by the fact that the free energy of mixing is zero. The combinatorial entropy
enters as a constant factor in the plots of solubility in different solvents, for example, as the
concentrations are usually constant for a given study.
It is important to note that the solubility parameter, or rather the difference in solubility
parameters for the solvent–solute combination, is important in determining the solubility of the
system. It is clear that a match in solubility parameters leads to a zero change in noncombinatorial
free energy, and the positive entropy change (the combinatorial entropy change), found on simple
mixing to result in a disordered mixture compared to the pure components, will ensure that a
solution is possible from a thermodynamic point of view. The maximum difference in solubility
parameters that can be tolerated where the solution still occurs is found by setting the noncombi-
natorial free energy change equal to the combinatorial entropy change:
This equation clearly shows that an alternate view of the solubility situation at the limit of
solubility is that it is the entropy change that dictates how closely the solubility parameters must
match each other for the solution to occur.
It will be seen in Chapter 2 that solvents with smaller molecular volumes will be thermody-
namically better than larger ones having identical solubility parameters. A practical aspect of this
effect is that solvents with relatively low molecular volumes, such as methanol and acetone, can
dissolve a polymer at larger solubility parameter differences than might be expected from compar-
isons with other solvents with larger molecular volumes. An average solvent molecular volume is
usually taken as about 100 cc/mol. The converse is also true. Larger molecular species may not
dissolve, even though solubility parameter considerations might predict they would. This can be a
difficulty in predicting the behavior of plasticizers solely based on data for lower molecular weight
solvents. These effects are also discussed elsewhere in this book, particularly in Chapter 2, Chapter
12, Chapter 13, and Chapter 16.
A shortcoming of the earlier solubility parameter work is that the approach was limited to
regular solutions, as defined by Hildebrand and Scott,2 and does not account for association between
molecules, such as those that polar and hydrogen-bonding interactions would require. The latter
problem seems to have been largely solved with the use of multicomponent solubility parameters;
however, the lack of accuracy with which the solubility parameters can be assigned will always
remain a problem. Using the difference between two large numbers to calculate a relatively small
heat of mixing, for example, will always be problematic.
A more detailed description of the theory presented by Hildebrand, and the succession of
research reports that have attempted to improve on it, can be found in Barton’s extensive handbook.9
The slightly older, excellent contribution of Gardon and Teas10 is also a good source of related
information, particularly for coatings and adhesion phenomena. The approach of Burrell,11 who
divided solvents into hydrogen bonding classes, has found numerous practical applications; the
approach of Blanks and Prausnitz12 divided the solubility parameter into two components, “non-
polar” and “polar.” Both are worthy of mention, however, in that the first has found wide use and
the second greatly influenced the author’s earlier activities. The Prausnitz article, in particular, was
farsighted in that a corresponding states procedure was introduced to calculate the dispersion energy
contribution to the cohesive energy. This is discussed in more detail in Chapter 2.
It can be seen from Equation 1.2 that the entropy change is beneficial to mixing. When
multiplied by the temperature, this will work in the direction of promoting a more negative free
energy of mixing. This is the usual case, although there are exceptions. Increasing temperature
does not always lead to improved solubility relations. Indeed, this was the basis of the pioneering
work of Patterson and coworkers,3–8 to show that subsequent increases in temperature can predict-
ably lead to insolubility. Their work was done in essentially nonpolar systems. Increasing temper-
ature can also lead to a nonsolvent becoming a solvent and, subsequently, a nonsolvent again with
still further increase in temperature. Polymer solubility parameters do not change much with
temperature, but those of a liquid frequently decrease rapidly with temperature. This situation allows
a nonsolvent, with a solubility parameter that is initially too high, to pass through a soluble condition
to once more become a nonsolvent as the temperature increases. These are usually “boundary”
solvents on solubility parameter plots.
The entropy changes associated with polymer solutions will be smaller than those associated
with liquid–liquid miscibility, for example, as the “monomers” are already bound into the config-
uration dictated by the polymer they make up. They are no longer free in the sense of a liquid
solvent and cannot mix freely to contribute to a larger entropy change. This is one reason poly-
mer–polymer miscibility is difficult to achieve. The free energy criterion dictates that polymer
solubility parameters match extremely well for mutual compatibility, as there is little to be gained
from the entropy contribution when progressively larger molecules are involved. However, poly-
mer–polymer miscibility can be promoted by the introduction of suitable copolymers or comono-
mers that interact specifically within the system. Further discussion of these phenomena is beyond
the scope of the present discussion; however, see Chapter 5.
polymers.13 This could never have been predicted by Hildebrand parameters, whereas the HSP
concept readily confirms the reason for this effect.
There are three major types of interactions in common organic materials. The most general are
the nonpolar interactions. These are derived from atomic forces and have also been called dispersion
interactions in the literature. As molecules are built up from atoms, all molecules contain those
types of attractive forces. For the saturated aliphatic hydrocarbons, for example, these are essentially
the only cohesive interactions, and the energy of vaporization is assumed to be the same as the
dispersion cohesive energy, ED. Finding the dispersion cohesive energy as the cohesion energy of
the homomorph, or hydrocarbon counterpart, is the starting point for calculating the three Hansen
parameters for a given liquid. As discussed in more detail later, this is based on a corresponding
states calculation.
The permanent dipole–permanent dipole interactions cause a second type of cohesion energy,
the polar cohesive energy, EP. These are inherently molecular interactions and are found in most
molecules to one extent or another. The dipole moment is the primary parameter used to calculate
these interactions. A molecule can be mainly polar in character without being water soluble, hence
there is a misuse of the term polar in the general literature. The polar solubility parameters referred
to here are well-defined, experimentally verified, and can be estimated from molecular parameters
as described later. As noted previously, the most polar of the solvents include those with relatively
high total solubility parameters that are not particularly water soluble, such as nitroparaffins,
propylene carbonate, and tri-n-butyl phosphate. Induced dipoles have not been treated specifically
in this approach but are recognized as a potentially important factor, particularly for solvents with
zero dipole moments (see the Calculation of the Polar Solubility Parameter section).
The third major cohesive energy source is hydrogen bonding, EH. This can be called more
generally an electron exchange parameter. Hydrogen bonding is a molecular interaction and resem-
bles the polar interactions in this respect. The basis of this type of cohesive energy is attraction
among molecules because of the hydrogen bonds. In this perhaps oversimplified approach, the
hydrogen bonding parameter has been used to more or less collect the energies from interactions
not included in the other two parameters. Alcohols, glycols, carboxylic acids, and other hydrophilic
materials have high-hydrogen-bonding parameters. Other researchers have divided this parameter
into separate parts — for example, acid and base cohesion parameters — to allow both positive
and negative heats of mixing. These approaches will not be dealt with here but are described in
Barton’s handbook9 and elsewhere.19–21 The most extensive division of the cohesive energy has
been done by Karger et al.,22 who developed a system with five parameters: dispersion, orientation,
induction, proton donor, and proton acceptor. As a single parameter, the Hansen hydrogen bonding
parameter has served remarkably well in the experience of the author and keeps the number of
parameters to a level that allows ready practical usage.
It is clear that there are other sources of cohesion energy arising in various types of molecules
from, for example, induced dipoles, metallic bonds, electrostatic interactions, or whatever type of
separate energy can be defined. The author stopped with the three major types found in organic
molecules. It has been recognized that additional parameters could be assigned to separate energy
types. For example, the description of organometallic compounds could be an intriguing study.
This would presumably parallel similar characterizations of surface-active materials, where each
part of the molecule requires separate characterization for completeness. The Hansen parameters
have mainly been used in connection with solubility relations, mostly, but not exclusively, in the
coatings and related industries.
Solubility and swelling have been used to confirm the solubility parameter assignments of many
of the liquids. Group contribution methods and suitable equations based on molecular properties
were then derived from these. They make possible estimates of the three parameters for additional
liquids. The goal of a prediction is to determine the similarity or difference of the cohesion energy
parameters. The strength of a particular type of hydrogen bond or any other bond is important only
to the extent that it influences the cohesive energy density.
HSPs do have direct applications in other scientific disciplines, such as surface science, where
they have been used to characterize the wettability of various surfaces and the adsorption properties
of pigment surfaces,10,14,16,23–26 and have even led to systematic surface treatment of inorganic fibers
so that they could be readily incorporated into polymers of low-solubility parameters such as
polypropylene27 (see also Chapter 7). Many widely different applications have been discussed by
Barton9 and Gardon.28 Surface characterizations have not been given the attention deserved in terms
of a unified similarity-of-energy approach. The author can certify that thinking in terms of similarity
of energy, whether surface or cohesive energies, can lead to rapid decisions and plans of action in
critical situations that lack data. In other words, the everyday industrial crisis situation often can
be reduced in scope by appropriate systematic approaches based on similarity of energy. The success
of the HSPs for surface applications are not surprising in view of the similarity of predictions
offered by these, and the Prigogine corresponding states theory of polymer solutions discussed in
Chapter 2. Flory also emphasized that it is the surface of molecules that interact to produce
solutions,29 so the interactions of molecules residing in surfaces should clearly be included in any
general approach to interactions among molecules. Surface mobility and surface rotation are
important factors in environmental stress cracking (Chapter 14), certain biological phenomena
(Chapter 15), the wetting of surfaces, and in other important phenomena relating to nanotechnology
(Chapter 18).
The basic equation governing the assignment of Hansen parameters is that the total cohesion
energy, E, must be the sum of the individual energies that make it up.
E = ED + EP + EH (1.6)
Dividing this by the molar volume gives the square of the total (or Hildebrand) solubility
parameter as the sum of the squares of the Hansen D, P, and H components.
δ2 = δ2 D + δ2 P + δ2 H (1.8)
To sum up this section, it is emphasized that HSPs quantitatively account for the cohesion
energy (density). Up to this point of time, an experimental latent heat of vaporization has been
considered a more reliable method to arrive at a cohesion energy rather than using molecular orbital
calculations or other calculations based on potential functions. Indeed, the goal of such extensive
calculations for polar and hydrogen bonding molecules should be to accurately arrive at the energy
of vaporization. The statistical thermodynamics approach of Panayiotou and coworkers reported in
Chapter 3 may have changed this. An alternative method of calculating the three parameters has
been presented, but full evaluation of this new information has not been possible as yet.
that mixtures of two nonsolvents could be systematically and synergistically (but predictably) found
to dissolve given polymers. This meant that these had parameters placing them on opposite sides
of the solubility region, a spheroid. Using a large number of such predictably synergistic systems
as a basis, reasonably accurate divisions into the three energy types were possible.
Using the experimentally established, approximate, δP and δH parameters, Hansen and Skaarup15
found that the Böttcher equation (Equation 10.25) could be used to calculate the polar parameter
quite well, and this led to a revision of the earlier values to those now accepted for the same liquids.
These values were also consistent with the experimental solubility data for 32 polymers available
at that time and with Equation 1.6. Furthermore, Skaarup developed the equation for the solubility
parameter “distance,” Ra, between two materials based on their respective partial solubility para-
meter components:
This equation was developed from plots of experimental data where the constant “4” was found
convenient and correctly represented the solubility data as a sphere encompassing the good solvents
(see Chapter 5). When the scale for the dispersion parameter is doubled, in comparison with the
other two parameters essentially spherical, rather than spheroidal, regions of solubility are found.
This greatly aids two-dimensional plotting and visualization. There are, of course, boundary regions
where deviations can occur. These are most frequently found to involve the larger molecular species
as being less effective solvents compared to their smaller counterparts that define the solubility
sphere. Likewise, smaller molecular species such as acetone, methanol, nitromethane, and others
often appear as outliers, in that they dissolve a polymer even though they have solubility parameters
placing them at a distance greater than the experimentally-determined radius of the solubility sphere,
Ro. This dependence on molar volume is inherent in the theory developed by Hildebrand, Scott,
and Scatchard discussed previously. Smaller molar volume favors lower ΔGM, as discussed in
Chapter 2. This in turn promotes solubility. Such smaller-molecular-volume species that dissolve
“better” than predicted by comparisons, based on solubility parameters alone, should not necessarily
be considered outliers.
The molar volume is frequently and successfully used as a fourth parameter to describe the
effects of molecular size. For example, these are especially important in correlating diffusional
phenomena with HSP (see Chapter 12, Chapter 13, and Chapter 16). The author has preferred to
retain the three, well-defined partial-solubility parameters with a fourth, separate, molar volume
parameter, rather than multiplying the solubility parameters by the molar volume raised to some
power to redefine them.
The reason for the experimentally determined constant 4 in Equation 1.9 will be discussed in
more detail in Chapter 2. It will be noted here, however, that the constant 4 is theoretically predicted
by the Prigogine corresponding states theory of polymer solutions when the geometric mean is
used to estimate the interaction in mixtures of dissimilar molecules.30 The constant 4 differentiates
between atomic and molecular intereactions. This is exceptionally strong evidence that dispersion,
permanent dipole–permanent dipole, and hydrogen bonding interactions all follow the geometric
mean rule. Patterson and coworkers have been especially instrumental in relating the Prigogine
theory to solubility parameters and to the Flory–Huggins theory of polymer solutions.3–8 The HSP
approach of dividing the cohesive energy into parts derived from different types of cohesive forces
has been confirmed both by experimental studies, as well as the Prigogine theory. The use of the
geometric mean is basic to this agreement between the HSP approach and that of Prigogine (see
Chapter 2).
The approach of optimizing solubility data to spheres is still very much in use. Plotting regions
of solubility based on experimental solubility data, or computer-optimizing boundaries of solubility
by locating the maximum difference in solubility parameters allowed by Equation 1.9 are both
used. The total free energy of mixing, ΔGM, is equal to zero on the boundary. It should be recognized
that using the solubility parameters relating to ΔGMnoncomb in Equation 1.4 differs from this by the
combinatorial entropy of mixing.
Another promising approach to arrive at the HSP for materials based on experimental data is
to use multivariable analysis of one type or another, as discussed in Chapter 5. This type of approach
has not been attempted by the author, but it clearly has advantages in some cases. The author’s
preferred approach of locating the polymer HSP as the center of a sphere has a problem in that it
is, in reality, the poor solvents or nonsolvents located near the boundary of the sphere that fix the
boundary (and center) rather than the best solvents in the middle. This may present problems for
smaller sets of data, but it is an advantage when extrapolating into regions of HSP higher than
those of any liquid that can be used in testing. This is discussed in more detail in Chapter 5 and
the definition of the limited segment of the boundary of the HSP sphere derivable from such
correlations is based on Equation 1.9.
Equation 1.9 is readily used on a computer (or on a hand calculator), and supplementary
relations allow easier scanning of large sets for data. It is obvious that solubility, or high affinity,
requires that Ra be less than Ro. The ratio Ra/Ro has been called the RED number, reflecting the
relative energy difference.
A RED number of 0 is found for no energy difference, RED numbers less than 1.0 indicate
high affinity; RED equal to or close to 1.0 is a boundary condition; and progressively higher RED
numbers indicate progressively lower affinities. Scanning a computer output for RED numbers less
than 1.0, for example, rapidly allows location of the most interesting liquids for a given application.
Parenthetically, it should be noted that the ratio of Ra to Ro is really a ratio of quantities having
the same units as the solubility parameter. The ratio (Ra/Ro)2 = (RED)2 is a ratio of cohesion
energies. The latter quantity is important for relating the HSP approach to that of Huggins and
Flory, as discussed in Chapter 2.
The revised set of parameters for the 90 original solvents was the basis for group contribution
procedures developed (most notably) by van Krevelen,31 Beerbower,32 and Hansen and Beerbower,17
who also used Fedors’ work.33 These various developments have been summarized by Barton,9
although Beerbower’s latest values have only appeared in the National Aeronautics and Space
Administration (NASA) document.32 Table 1.1 is an expanded table of Beerbower group contribu-
tions, which was distributed among those who were in contact with Beerbower in the late 1970s.
The majority of the data in this table, as well as Table 1.2, have also appeared in Reference 34.
Beerbower also developed a simple equation for the polar parameter,17 which involved only the
dipole moment and the square root of the molar volume. This is also given later (Equation 1.13)
and has been found quite reliable by Koenhen and Smolders.35 This equation has been found reliable
by the author as well, giving results generally consistent with Equation 1.6 to Equation 1.8, which,
again, is the basis of the whole approach. Koenhen and Smolders also give correlation coefficients
for other calculation procedures to arrive at the individual Hansen parameters.
The group contributions in Table 1.1 have been used extensively to arrive at the collection of
HSP data in Appendix Table A.1. Most of the chemicals of primary interest for which full data are
available are presumably already in this table. The trend has been to calculate HSP for larger and
still larger molecules. Many of these have multiple groups, and it becomes more and more difficult
to make decisions as to how to treat them best. At times the HSP for the larger molecules can be
estimated from the HSP of larger segments that make them up. Rather than expanding Table 1.1
with additional data, except as noted briefly later, the usual practice has been to locate chemicals
with similar groups and to use their HSP values in a group contribution-type calculation.
The procedure has developed to the point where its principle features can be identified in the
following table. If a boiling point is available, the procedures for calculating δD have been used.
If a boiling point is not available, the similarity with related molecules has been used. If a dipole
moment is available, the procedures given here were used in preference to group contributions. If
necessary, group contributions can be derived from similar molecules to the one in question, when
dipole moments are available for these, and not for the molecule in question. There is often a
change in the group contribution as a function of molecular size. This is the main reason for the
lack of expansion of Table 1.1. It is thought best that the uncertainty be clear to the user. For
example, it has been found that the group contribution for the polar component of aliphatic esters
should not be less than 300 cal/mol as given in Table 1.1. This is necessary to prevent δP for
materials like plasticizers from being clearly too low, based on their compatibility with, for example,
polyvinyl chloride. Sulfur containing compounds have also been somewhat difficult in this respect,
with major changes in estimated group contributions depending on molecular weight of the chemical
in question. Group contributions for sulfur, amides, and other groups not found in Table 1.1 can
be easily derived from the data on similar compounds reported in Appendix Table A.1. The same
is true of the δH component. The problem with this procedure is obvious: any error or distortion
of value for a class of compounds is perpetuated. This has been recognized and dealt with to the
extent possible, but there are limits to what can be done with limited data. The scope of this situation
has been beyond the resources available for its fully satisfactory resolution. The extensive list of
group contributions at the end of Chapter 3 provides what may be a partial replacement and/or a
supplement for Table 1.1. This requires some experience with the techniques involved.
A sizable number of materials have been assigned HSPs using the procedures described here.
Many of these have not been published. Exxon Chemical Corporation36,37 has indicated a computer
program with data for over 500 solvents and plasticizers, 450 resins and polymers, and 500
pesticides. The author’s files contain the three parameters for about 1200 chemicals (See Appendix
Table A.1), although several of them appear with two sets of possible values awaiting experimental
confirmation. In some cases, this is due to questionable physical data, for example, for latent heats
of vaporization, or wide variations in reported dipole moments. Another reason is that some liquids
are chameleonic,38 as defined by Hoy, in that they adopt configurations depending on their envi-
ronment. Hoy38 cites the formation of cyclic structures for glycol ethers with (nominally) linear
structure. The formation of hydrogen-bonded tetramers of alcohols in a fluoropolymer has also
been pointed out.39 The term compound formation can be found in the older literature, particularly
where mixtures with water were involved and structured species were postulated to explain phe-
nomena based on specific interactions among the components of the mixtures. Barton has discussed
some of the situations where cohesion parameters need a more careful use, and points out that
Hildebrand or Hansen parameters must be used with particular caution where the extent of
donor–acceptor interactions, especially, hydrogen bonding within a compound, is very different
from that between compounds.18 Amines, for example, are known to associate with each other.
Pure component data cannot be expected to predict the behavior in such cases.
Still another reason for difficulties is the large variation of dipole moments reported for the
same liquid. The dipole moment for some liquids depends on their environment, as discussed later.
A given solvent can be listed with different values in files to keep these phenomena in mind.
Large data sources greatly enhance a search for similar materials and the locating of new
solvents, as an example, for a polymer for which there are limited data. Unfortunately, different
authors have used different group contribution techniques, and there is a proliferation of different
“Hansen” parameters for the same chemicals in the literature. This would seem to be an unfortunate
situation, but may ultimately provide benefits. In particular, partial solubility parameter values
found in Hoy’s extensive tables9,40 are not compatible with the customary Hansen parameters
reported here. Hoy has provided an excellent source of total solubility parameters. He independently
arrived at the same type of division of cohesion energies as Hansen, although the methods of
calculation were quite different.
Many solvent suppliers have also presented tables of solvent properties and/or use computer
techniques with these tables in their technical service. Partial solubility parameters not taken directly
from earlier well-documented sources should be used with caution. The Hoy dispersion parameter,
Source: From Hansen, C. M., Paint Testing Manual, Manual 17, Koleske, J. V., Ed., American Society for Testing and Materials, Philadelphia, 1995, 388. Copyright ASTM. Reprinted with permission.
11
© 2007 by Taylor & Francis Group, LLC
7248_C001.fm Page 12 Monday, April 23, 2007 1:56 PM
TABLE 1.2
Lydersen Group Constants
Aliphatic, Cyclic, Aliphatic, Cyclic,
Group ΔT ΔT ΔPT ΔP ΔP
–NCO — — 0.0539 — —
HCON< — — 0.0546 — —
–CONH– — — 0.0843 — —
–CON< — — 0.0729 — —
–CONH2 — — 0.0897 — —
–OCONH– — — 0.0938 — —
Conjugation — — 0.0035 — —
cis double bond — — –0.0010 — —
trans double bond — — –0.0020 — —
Ortho — — 0.0015 — —
Meta — — 0.0010 — —
Para — — 0.0060 — —
Bicycloheptyl — — 0.0034 — —
Tricyclodecane — — 0.0095 — —
Source: Hansen, C.M., Solubility parameters, in Paint Testing Manual, Manual 17,
Koleske, J.V., Ed., American Society for Testing and Materials, Philadelphia, PA,
1995, 383–404. Reprinted with permission.
in particular, is consistently lower than that found by Hansen. Hoy subtracts estimated values of
the polar and hydrogen-bonding energies from the total energy to find the dispersion energy. This
allows for more calculational error and underestimates the dispersion energy, as the Hoy procedure
does not appear to fully separate the polar and hydrogen-bonding energies. The van Krevelen
dispersion parameters appear to be too low. The author has not attempted these calculations, being
completely dedicated to the full procedure based on corresponding states described here, but values
estimated independently using the van Krevelen dispersion parameters are clearly low. A compar-
ison with related compounds or the similarity principle gives better results than those found from
the van Krevelen dispersion group contributions.
In the following, calculation procedures and experience are presented according to the proce-
dures most reliable for the experimental and/or physical data available for a given liquid.
Δ ED kJ/mol
70
60
0
= 0.4
5
Tr 0.4
0
50 0.5
5
0.5
0.60
40 0.65
0.70
30
20
10
0
50 100 150 200 250 V, cm3/mol
FIGURE 1.1 Energy of vaporization for straight chain hydrocarbons as a function of molar volume and
reduced temperature. (From Hansen, C. M., Paint Testing Manual, Manual 17, Koleske, J. V., Ed., American
Society for Testing and Materials, Philadelphia, 1995, 389. Copyright ASTM. Reprinted with permission.)
400
350 Tr =
0.40
c, MPa
0.45
300
0.50
0.55
250 0.60
0.65
0.70
50 60 70 80 90 100 110 120 130
3
V, cm /mol
FIGURE 1.2 Cohesive energy density for cycloalkanes as a function of molar volume and reduced temper-
ature. (From Hansen, C. M., Paint Testing Manual, Manual 17, Koleske, J. V., Ed., American Society for
Testing and Materials, Philadelphia, 1995, 389. Copyright ASTM. Reprinted with permission.)
400
Tr =
0.40
350
c, MPa
0.45
0.50
300
0.55
0.60
250
0.65
0.70
80 90 100 110 120 130 140 150 160 170
V, cm3/mol
FIGURE 1.3 Cohesive energy density for aromatic hydrocarbons as a function of molar volume and reduced
temperature. (From Hansen, C. M., Paint Testing Manual, Manual 17, Koleske, J. V., Ed., American Society
for Testing and Materials, Philadelphia, 1995, 389. Copyright ASTM. Reprinted with permission.)
is inconsistent with the energy figure and in error. Its use is not recommended. When substituted
cycloaliphatics or substituted aromatics are considered, simultaneous consideration of the two
separate parts of the molecule is required. The dispersion energies are evaluated for each of the
molecules involved, and a weighted average is taken for the molecule of interest based on the
number of significant atoms. For example, hexyl benzene would be the arithmetic average of the
dispersion energies for an aliphatic and an aromatic liquid, each with the given molar volume of
hexyl benzene. Liquids such as chlorobenzene, toluene, and ring compounds with alkyl substitutions
that have only two or three carbon atoms have been considered only as cyclic compounds. Such
weighting has been found necessary to satisfy Equation 1.6.
The critical temperature, Tc, is required to use the dispersion energy figures. If the critical
temperature cannot be found, it must be estimated. A table of the Lydersen group contributions,41
ΔT , as given by Hoy40 for calculation of the critical temperature is included as Table 1.2. In some
cases, the desired groups may not be in the table, which requires some educated guessing. The end
result does not appear too sensitive to these situations. The normal boiling temperature, Tb, is also
required in this calculation. This is not always available and must be estimated by similarity, group
contribution, or some other technique. The Lydersen group contribution method involves the use
of Equation 1.11 and Equation 1.12 as follows:
and
Tr = T/Tc (1.12)
taken as 1650 J/mol for each of these atoms in the molecule. Dividing by the molar volume and
then taking the square root gives the (large atom corrected) dispersion solubility parameter.
The need for these corrections has been confirmed many times, both for interpretation of
experimental data and to allow Equation 1.6 to Equation 1.8 to balance. Research is definitely
needed in this area. The impact of these corrections is, of course, larger for the smaller molecular
species. Taking square roots of the larger numbers involved with the larger molecular species
reduces the errors involved in these cases, as the corrections are relatively small.
It can be seen from the dispersion parameters of the cyclic compounds that the ring has an
effect similar to increasing the effective size of the interacting species. The dispersion energies for
cycloaliphatic compounds are larger than their aliphatic counterparts, and they are higher for
aromatic compounds than their corresponding cycloaliphatics. Similar effects also appear with the
ester group. This group appears to act as if it were, in effect, an entity that is larger than the
corresponding compound containing only carbon (i.e., its homomorph), and it has a higher disper-
sion solubility parameter without any special need for corrections.
The careful evaluation of the dispersion cohesive energy may not have a major impact on the
value of the dispersion solubility parameter, because square roots of rather large numbers are taken.
Larger problems arise because of Equation 1.6. Energy assigned to the dispersion portion cannot
be reused when finding the other partial parameters using Equation 1.6 (or Equation 1.8). This is
one reason group contributions are recommended in some cases, as discussed later.
δP = 37.4(DM)/V1/2 (1.13)
have considerable polarity. The lower dipole moments seem to fit in other instances. Diethyl ether
has also presented problems as an outlier in terms of dissolving or not and permeating rapidly or
not. Here, the reported dipole moments42 vary from 0.74 to 2.0, with a preferred value of 1.17 and
1.79 in chloroform. Choosing a given value seems rather arbitrary. The chameleonic cyclic forms
of the linear glycol ethers would also seem to provide a basis for altered dipole moments in various
media.38
When Equation 1.13 cannot be used, the polar solubility parameter has been found using the
Beerbower table of group contributions, by similarity to related compounds and/or by subtraction
of the dispersion and hydrogen bonding cohesive energies from the total cohesive energy. The
question in each case is, “Which data are available and judged most reliable?” New group contri-
butions can also be developed from related compounds whose dipole moments are available. These
new polar group contributions then become supplementary to the Beerbower table.
For large molecules, especially those with long hydrocarbon chains, the accurate calculation
of the relatively small polar (and hydrogen bonding) contributions present special difficulties. The
latent heats are not generally available with sufficient accuracy to allow subtraction of two large
numbers from each other to find a very small one. In such cases, the similarity and group contribution
methods are thought best. Unfortunately, latent heats found in a widely used handbook45 are not
clearly reported as to the reference temperature. There is an indication that these are 25°C data,
but checking indicated many of the data to be identical with boiling point data reported elsewhere
in the literature. Subsequent editions of this handbook46 have a completely different section for the
latent heat of evaporation. Again, even moderate variations in reported heats of vaporization can
cause severe problems in calculating the polar (or hydrogen bonding) parameter when Equation
1.6 or Equation 1.8 are strictly adhered to.
This is done even if the melting point of the compound being considered is higher than 25°C.
The result is consistent with all the other parameters, and to date no problems with particularly
faulty predictions have been noted. It appears as if the predictions are not significantly in error
when experimental data are available for checking. When the latent heat at the boiling point is
given in cal/mol, Equation 1.14 is used to estimate the latent heat at 25°C. RT equal to 592 cal/mol
is then subtracted from this according to Equation 1.15, to find the total cohesion energy, E, in cgs
units at this temperature:
TEMPERATURE DEPENDENCE
Only very limited attempts have been made to calculate solubility parameters at a higher temperature
prior to the second edition of this handbook. The inclusion of Chapter 3 and Chapter 10 in this
handbook helps by providing a more accurate treatment of temperature dependence when the
situation warrants it. Solubility parameter correlations of phenomena at higher temperatures have
generally been found satisfactory when the established 25°C parameters have been used. Recalcu-
lation to higher temperatures is possible but has not generally been found necessary. In this direct
but approximate approach, it is assumed that the parameters all demonstrate the same temperature
dependence, which, of course, is not the case. It might be noted in this connection that the hydrogen-
bonding parameter, in particular, is the most sensitive to temperature. As the temperature increases,
more and more hydrogen bonds are progressively broken or weakened, and this parameter will
decrease more rapidly than the others.
The gas-phase dipole moment is not temperature dependent, although the volume of a fluid
does change with the temperature, which will also change its cohesive energy density. The change
of the δD, δP , and δH parameters for liquids with temperature, T, can be estimated by the following
equations, where α is the coefficient of thermal expansion17:
that solvent molecular size can be an important fourth parameter in solubility and, in some cases,
in chemical resistance. Specific examples are given in Chapter 5 and Chapter 12.
The size and shape of the solvent molecule are also very important for kinetic phenomena such
as diffusion, permeation, and attainment of equilibrium. Smaller and more linear molecules diffuse
more rapidly than larger and more bulky ones. The diffusion coefficient may be so low that
equilibrium is not attained for hundreds of years at room temperature. This was demonstrated in
common solvent exposures of rigid polymers like polyphenylene sulfide (PPS) with thicknesses of
several millimeters.49 Likewise, the second stage in the two-stage drying process in polymer film
formation by solvent evaporation can last for many years.16,50 Polymer samples used for solubility
parameter or other testing may well retain solvent or monomer for many years, and this may affect
the evaluations.
Attempts to include the molecular volume in new composite solubility and size parameters
have not been particularly successful.20,21 This may be because the size effect is most often not
caused due to the thermodynamic considerations on which the solubility parameters are based, but
rather through a kinetic effect of diffusion rates or other free volume considerations. The similarities
in the HSP approach and the Prigogine theory, discussed in Chapter 2, indicate a remarkably close,
if not identical, relation between the Prigogine ρ (segment size parameter) and the δD parameter,
suggesting that molecular size differences are at least partially accounted for in the δD parameter.
The Prigogine theory also has a parameter to describe “structural effects,” including size of polymer
molecules, but this has not been explored in relation to the present discussion. The increase of δD
with increasing molecular size among the aliphatic hydrocarbons, the higher δD values for the larger
units represented by cycloaliphatic and aromatic rings, and the need for corrections for larger atoms
discussed earlier tend to support the molecular size differences.
Sorting output data according to the molecular volume of the test solvents in a computer analysis
helps to discover whether solvent molecular size is indeed an additional significant factor in a given
correlation or testing program.
COMPUTER PROGRAMS
The author has used two computer programs extensively in his own studies and in collecting material
for this book. These are called SPHERE and SPHERE1. They are very similar, the only difference
being that SPHERE optimizes the polymer (or other material, of course) parameters based on all
the data, whereas SPHERE1 considers data for only those solvents considered as “good.” It neglects
the nonsolvent data. SPHERE1 has been most useful in correlations with pigments, fillers, and
fibers, as described in Chapter 7.
The data input is by solvent number followed by an indication of the quality of interaction
with that solvent. A “1” indicates a “good” solvent, whereas a “0” is used for a “bad” solvent.
What is considered good or bad varies according to the level of interaction being studied. This can
be solution or not, a given percentage of swelling or uptake, breakthrough time being less than a
given interval, permeation coefficients higher than a given value, long-time suspension of a pigment,
etc.
The program systematically evaluates the input data using a quality-of-fit function called the
desirability function.51 This suggestion was made by a reputed statistician many years ago as the
most appropriate statistical treatment for this type of problem. It has been in use since the late
1960s. The function has the form:
where n is the number of solvents for which there is experimental data in the correlation. The
DATA FIT approaches 1.0 as the fit improves during an optimization and reaches 1.0 when all the
good solvents are included within the sphere and all the bad ones are outside the sphere.
The Ai quotient for a given good solvent within the sphere and a bad solvent outside the sphere
will be 1.0. The error distance is the distance of the solvent in error to the sphere boundary. It
could either denote being good and outside the sphere or being bad and inside the sphere.
Ro is the radius of the sphere, and Ra is the distance from a given solvent point to the center
of the sphere. For a good solvent outside the sphere, an error enters the DATA FIT according to:
Such errors are often found for solvents having low molecular volumes.
For a bad solvent inside the sphere, the contribution to the DATA FIT is
Such errors can sometimes be found for larger molecular species such as plasticizers. This is
not unexpected for the reasons mentioned earlier.
Solvents with large and/or small molecules that give the “errors” can sometimes be (temporarily)
disregarded by generating a new correlation; this gives an excellent DATA FIT for an abbreviated
range of molecular volumes. There is a special printout with the solvents arranged in order of
molecular volume that helps to analyze such situations. The computer printouts all include a column
for the RED number.
The program assumes a starting point, based on an average of each of the HSP for the good
solvents only. The program then evaluates eight points at the corners of a cube, with the current
best values as center. Different radii are evaluated at each of these points in the optimization process.
When better fits are found among the eight points, the point with the best fit is taken as a new
center, and eight points around it are evaluated in a similar manner. This continues until the DATA
FIT cannot be improved upon. The length of the edge of the cube is then reduced in size to fine-
tune the fit. The initial length of the cube is 1 unit, which is reduced to 0.3 unit, and finally to 0.1
unit in the last optimization step.
Experimental data for the solvents are entered with solvent number (comma) and a “1” for a
good solvent, or a “0” for a bad one.
Errors in the correlations are indicated with an “*” in the SOLUB column where the experi-
mental input data are indicated. As stated above, systematic errors can sometimes be seen in the
molar volume printout. This may suggest a new analysis of the data. Nonsystematic errors may be
real, such as for reactions or some extraneous effect not predictable by the solubility parameter.
They may also be bad data, and rechecking data indicated with an “*” in the output has become
a routine practice. The output of this program is for the least radius allowing the maximum DATA
FIT. An example is found in Table 5.4.
Results from the SPHERE program reported in this book generally include the HSP, given as
D (δD), P (δP), and H (δH), respectively, and Ro for the correlation in question, as well as the DATA
FIT, the number of good solvents (G), and the total solvents (T) in the correlation. This latter
information has not always been recorded and may be lacking for some correlations, especially the
older ones.
TABLE 1.3
HSP Correlations Related to Water
Correlation δD δP δH Ro FIT G/T
is not recommended. This is true of water as a pure liquid or in mixtures. Table 1.3 includes data
from various HSP analyses of the behavior of water. The first set of data is derived from the energy
of vaporization of water at 25°C. The second set of data is based on a correlation of the solubility
of various solvents in water, where “good” solvents are soluble to more than 1% in water. “Bad”
solvents dissolve to a lesser extent. The third set of data is for a correlation of total miscibility of
the given solvents with water. The second and third entries in Table 1.3 are based on the SPHERE
program where both good and bad solvents affect the DATA FIT and hence the result of the
optimization. The last entry in Table 1.3 is for an analysis using the SPHERE1 program. The HSP
data are for the minimum sphere that encompasses only the good solvents. The bad solvents are
simply not considered in the data processing. This type of comparison usually results in some of
the parameters being lower than when all the data are included. A frequent problem is that a
considerable portion of the HSP spheres, such as in the case for water, covers such high energies
that no liquid can be found. The cohesion energy is so high as to require solids. The constant 4 in
the correlations (Equation 1.9) is still used for these correlations, primarily based on successes at
lower levels of cohesion energies, but this is also supported by the comparison with the Prigogine
cst of polymer solutions, discussed at some length in Chapter 2. The HSP for water as a single
molecule, based on the latent heat at 25°C is sometimes used in connection with mixtures with
water to estimate average HSP. More recently, it has been found in a study involving water, ethanol,
and 1,2-propanediol that the HSP for water indicated by the total water solubility correlation could
be used to explain the behavior of the mixtures involved. The averaged values are very questionable
as water can associate, and water has a very small molar volume as a single molecule. It almost
appears to have a dual character. The data for the 1% correlation,52 as well as for the total water
miscibility, suggest that about six water molecules associate into units.
Traditionally, solvents are considered as points. This is practical and almost necessary from an
experimental point of view as most solvents are so miscible as to not allow any experimental
characterization in terms of a solubility sphere. An exception to this is the data for water reported
in Table 1.3. The HSP reported here are the center points of HSP spheres where the good solvents
are either those that are completely miscible or those that are miscible to only 1% or more, as
discussed previously. It should also be mentioned that amines were a major source of outliers in
these correlations. No solids were included. Their use to predict solubility relations for amines and
for solids must therefore be done with caution.
CONCLUSION
This chapter has been dedicated to describing the tools with which different HSP characterizations
can be made and some of the pitfalls that may be encountered in the process. The justification for
Alcohols Esters
Glycols Ethers
Amides Amines
Phenols Methyl ethyl ketone
Chloroform Cyclohexanone
Chlorinated
Dichloromethane THF
Nitriles Toluene
Dimethylsulfoxide Isooctane
Nitromethane
24
20 Nitriles
N
A
S
16
δP, Polar Parameter
Amides
F Alcohols
B
12 M
A E
Ketones
M P Glycols
E
8
Chlorinated Phenols
C
E T P B
C
4 B Esters
Ethers O
E
P
Tol Amines
0
0 I 4 8 12 16 20 24 28
δH, Hydrogen Bonding Parameter
FIGURE 1.4 δp vs. δH plot showing the location of various common solvents. The glycols are ethylene glycol
and propylene glycol. The alcohols include methanol (M), ethanol (E), 1-butanol (B), and 1-octanol (O). The
amides include dimethyl formamide (F) and dimethyl acetamide (A). The nitriles are acetonitrile (A) and
butyronitrile (B). The esters are ethyl acetate (E) and n-butyl acetate (B). The amines are ethyl amine (E) and
propyl amine (P). The phenols are phenol (P) and m-cresol (C). The ethers are symbolized by diethyl ether.
the tools is further confirmed in Chapter 2 and Chapter 3, and their use is demonstrated in all the
subsequent chapters. Figure 1.4 is included to show where many common solvents are located on
a δp vs. δH plot.
REFERENCES
1. Hildebrand, J. and Scott, R.L., The Solubility of Nonelectrolytes, 3rd ed., Reinhold, New York, 1950.
2. Hildebrand, J. and Scott, R.L., Regular Solutions, Prentice-Hall, Englewood Cliffs, NJ, 1962.
3. Patterson, D. and Delmas, G., New aspects of polymer solution thermodynamics, Off. Dig. Fed. Soc.
Paint Technol., 34(450), 677–692, 1962.
4. Delmas, D., Patterson, D., and Somcynsky, T., Thermodynamics of polyisobutylene-n-alkane systems,
J. Polym. Sci., 57, 79–98, 1962.
5. Bhattacharyya, S.N., Patterson, D., and Somcynsky, T., The principle of corresponding states and the
excess functions of n-alkane mixtures, Physica, 30, 1276–1292, 1964.
6. Patterson, D., Role of free volume changes in polymer solution thermodynamics, J. Polym. Sci. Part
C, 16, 3379–3389, 1968.
7. Patterson, D.D., Introduction to thermodynamics of polymer solubility, J. Paint Technol., 41(536),
489–493, 1969.
8. Biros, J., Zeman, L., and Patterson, D., Prediction of the C parameter by the solubility parameter and
corresponding states theories, Macromolecules, 4(1), 30–35, 1971.
9. Barton, A.F.M., Handbook of Solubility Parameters and Other Cohesion Parameters, CRC Press,
Boca Raton, FL, 1983; 2nd ed., 1991.
10. Gardon, J.L. and Teas, J.P., Solubility parameters, in Treatise on Coatings, Vol. 2, Characterization
of Coatings: Physical Techniques, Part II, Myers, R.R. and Long, J.S., Eds., Marcel Dekker, New
York, 1976, chap. 8.
11. Burrell, H., Solubility parameters for film formers, Off. Dig. Fed. Soc. Paint Technol., 27(369),
726–758, 1972; Burrell, H., A solvent formulating chart, Off. Dig. Fed. Soc. Paint Technol., 29(394),
1159–1173, 1957; Burrell, H., The use of the solubility parameter concept in the United States, VI
Federation d’Associations de Techniciens des Industries des Peintures, Vernis, Emaux et Encres
d’Imprimerie de l’Europe Continentale, Congress Book, 21–30, 1962.
12. Blanks, R.F. and Prausnitz, J.M., Thermodynamics of polymer solubility in polar and nonpolar
systems, Ind. Eng. Chem. Fundam., 3(1), 1–8, 1964.
13. Hansen, C.M., The three dimensional solubility parameter — key to paint component affinities I, J.
Paint Technol., 39(505), 104–117, 1967.
14. Hansen, C.M., The three dimensional solubility parameter — key to paint component affinities II, J.
Paint Technol., 39(511), 505–510, 1967.
15. Hansen, C.M. and Skaarup, K., The three dimensional solubility parameter — key to paint component
affinities III, J. Paint Technol., 39(511), 511–514, 1967.
16. Hansen, C.M., The Three Dimensional Solubility Parameter and Solvent Diffusion Coefficient, Doc-
toral dissertation, Danish Technical Press, Copenhagen, 1967.
17. Hansen, C.M. and Beerbower, A., Solubility parameters, in Kirk-Othmer Encyclopedia of Chemical
Technology, Suppl. Vol., 2nd ed., Standen, A., Ed., Interscience, New York, 1971, pp. 889–910.
18. Barton, A.F.M., Applications of solubility parameters and other cohesion energy parameters, Polym.
Sci. Technol. Pure Appl. Chem., 57(7), 905–912, 1985.
19. Sørensen, P., Application of the acid/base concept describing the interaction between pigments,
binders, and solvents, J. Paint Technol., 47(602), 31–39, 1975.
20. Van Dyk, J.W., Paper presented at the Fourth Chemical Congress of America, New York, August
25–30, 1991.
21. Anonymous [Note: This was, in fact, Van Dyk, J.W., but this does not appear on the bulletin], Using
Dimethyl Sulfoxide (DMSO) in Industrial Formulations, Bulletin No. 102, Gaylord Chemical Corp.,
Slidell, LA, 1992.
22. Karger, B.L., Snyder, L.R., and Eon, C., Expanded solubility parameter treatment for classification
and use of chromatographic solvents and adsorbents, Anal. Chem., 50(14), 2126–2136, 1978.
23. Hansen, C.M. and Wallström, E., On the use of cohesion parameters to characterize surfaces, J. Adhes.,
15, 275–286, 1983.
24. Hansen, C.M., Characterization of surfaces by spreading liquids, J. Paint Technol., 42(550), 660–664,
1970.
25. Hansen, C.M., Surface dewetting and coatings performance, J. Paint Technol., 44(570), 57–60, 1972.
26. Hansen, C.M. and Pierce, P.E., Surface effects in coatings processes, Ind. Eng. Chem. Prod. Res. Dev.,
13(4), 218–225, 1974.
27. Hennissen, L., Systematic Modification of Filler/Fiber Surfaces to Achieve Maximum Compatibility
with Matrix Polymers, Lecture for the Danish Society for Polymer Technology, Copenhagen, February
10, 1996.
28. Gardon, J.L., Critical review of concepts common to cohesive energy density, surface tension, tensile
strength, heat of mixing, interfacial tension and butt joint strength, J. Colloid Interface Sci., 59(3),
582–596, 1977.
29. Flory, P.J., Principles of Polymer Chemistry, Cornell University Press, New York, 1953.
30. Prigogine, I. (with the collaboration of Bellemans, A. and Mathot, A.), The Molecular Theory of
Solutions, North-Holland, Amsterdam, 1957, chap. 16, 17.
31. van Krevelen, D.W. and Hoftyzer, P.J., Properties of Polymers: Their Estimation and Correlation with
Chemical Structure, 2nd ed., Elsevier, Amsterdam, 1976.
32. Beerbower, A., Environmental Capability of Liquids, in Interdisciplinary Approach to Liquid Lubri-
cant Technology, NASA Publication SP-318, 1973, 365–431.
33. Fedors, R.F., A method for estimating both the solubility parameters and molar volumes of liquids,
Polym. Eng. Sci., 14(2), 147–154, 472, 1974.
34. Hansen, C.M., Solubility parameters, in Paint Testing Manual, Manual 17, Koleske, J.V., Ed., American
Society for Testing and Materials, Philadelphia, PA, 1995, pp. 383–404.
35. Koenhen, D.N. and Smolders, C.A., The determination of solubility parameters of solvents and
polymers by means of correlation with other physical quantities, J. Appl. Polym. Sci., 19, 1163–1179,
1975.
36. Anonymous, Brochure: Co-Act — A Dynamic Program for Solvent Selection, Exxon Chemical
International Inc., 1989.
37. Dante, M.F., Bittar, A.D., and Caillault, J.J., Program calculates solvent properties and solubility
parameters, Mod. Paint Coat., 79(9), 46–51, 1989.
38. Hoy, K.L., New values of the solubility parameters from vapor pressure data, J. Paint Technol., 42(541),
76–118, 1970.
39. Myers, M.M. and Abu-Isa, I.A., Elastomer solvent interactions III — effects of methanol mixtures
on fluorocarbon elastomers, J. Appl. Polym. Sci., 32, 3515–3539, 1986.
40. Hoy, K.L., Tables of Solubility Parameters, Union Carbide Corp., Research and Development Dept.,
South Charleston, WV, 1985; 1st ed. 1969.
41. Reid, R.C. and Sherwood, T.K., Properties of Gases and Liquids, McGraw-Hill, New York, 1958
(Lydersen Method — see also Reference 31).
42. McClellan, A.L., Tables of Experimental Dipole Moments, W.H. Freeman, San Francisco, 1963.
43. Tables of Physical and Thermodynamic Properties of Pure Compounds, American Institute of Chem-
ical Engineers Design Institute for Physical Property Research, Project 801, Data Compilation, Danner,
R.P. and Daubert, T.E., Project Supervisors, DIPPR Data Compilation Project, Department of Chemical
Engineering, Pennsylvania State University, University Park.
44. Hansen, C.M., Selection of Chemicals for Permeation Testing Based on New Solubility Parameter
Models for Challenge 5100 and Challenge 5200, under contract DTCG50-89-P-0333 for the U.S.
Coast Guard, June 1989, Danish Isotope Centre, Copenhagen.
45. Weast, R.C., (Editor-in-Chief), CRC Handbook of Chemistry and Physics, 65th ed., CRC Press, Boca
Raton, FL, 1988–1989, pp. C-672–C-683.
46. Majer, V., Enthalpy of vaporization of organic compounds, in Handbook of Chemistry and Physics,
72nd ed., Lide, D.R., (Editor-in-Chief), CRC Press, Boca Raton, FL, 1991–1992, pp. 6-100–6-107.
47. Fishtine, S.H., Reliable latent heats of vaporization, Ind. Eng. Chem., 55(4), 20–28, 1963; Ind. Eng.
Chem., 55(5), 55–60; Ind. Eng. Chem., 55(6), 47–56.
48. Hansen, C.M., New developments in corrosion and blister formation in coatings, Prog. Org. Coat.,
26, 113–120, 1995.
49. Hansen, C.M., Solvent Resistance of Polymer Composites — Glass Fiber Reinforced Polyphenylene
Sulfide, Centre for Polymer Composites (Denmark), Danish Technological Institute, Taastrup, 1993,
1–62, ISBN 87-7756-286-0.
50. Hansen, C.M., A mathematical description of film drying by solvent evaporation, J. Oil Color Chem.
Assoc., 51(1), 27–43, 1968.
51. Harrington, E.C., Jr., The desirability function, Ind. Qual. Control, 21(10), 494–498, April 1965.
52. Hansen, C.M. and Andersen, B.H., The affinities of organic solvents in biological systems, Am. Ind.
Hyg. Assoc. J., 49(6), 301–308, 1988.