Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Supersaturation Ratios and Thermodynamic Driving Force

Dissolution and Solubility Modeling · Tutorial 250 of 799

Supersaturation Ratios and Thermodynamic Driving Force

Learn to calculate supersaturation and interpret its thermodynamic meaning while accounting for changing solubility and concentration across a dissolution profile.

Advanced 10 min read

What You'll Learn

  • Define the supersaturation ratio and degree relative to the equilibrium solubility of a specified solid form.
  • Relate activity-based supersaturation to the chemical-potential difference from equilibrium.
  • Classify undersaturated, saturated, and metastably supersaturated concentration states.
  • Calculate and interpret supersaturation along a hypothetical concentration profile.
  • Recognize when a concentration ratio is misleading because the reference solubility or measured concentration is inappropriate.
  • Compare ratio-based, absolute-excess, and activity-based descriptions of the driving force.

Why a Concentration Must Be Compared With Its Equilibrium Reference

A dissolved concentration does not, by itself, say whether a solution is undersaturated or supersaturated. A concentration of 0.20 mg/mL could be below equilibrium solubility in one medium and above it in another. The relevant comparison is with the equilibrium solubility for the specified compound form and conditions: temperature, pH, medium composition, and the solid phase that is in equilibrium with the solution.

This comparison is useful when interpreting dissolution profiles. As the dissolved concentration changes over time or between gastrointestinal segments, the equilibrium reference may change as well. A profile can therefore move between undersaturated and supersaturated states even if concentration changes only modestly. Supersaturation indicates a thermodynamic tendency for dissolved drug to leave solution and approach equilibrium; it does not, by itself, specify how quickly that will occur.

Definition: Let \(C\) be the concentration of dissolved drug in the solution of interest and \(C_{\mathrm{eq}}\) its equilibrium solubility under the same conditions, relative to a specified solid phase. The concentration-based supersaturation ratio is \(S_C=C/C_{\mathrm{eq}}\). The degree of supersaturation is \(S_C-1=(C-C_{\mathrm{eq}})/C_{\mathrm{eq}}\); as a percentage, it is \(100(S_C-1)\%\). Thus \(S_C<1\) is undersaturation, \(S_C=1\) is saturation, and \(S_C>1\) is supersaturation.

Both concentrations in the ratio must use the same units, and the denominator must be an equilibrium solubility appropriate to the conditions at which \(C\) is defined. For example, using a solubility measured at a different pH or for a different solid form can give the wrong classification. The ratio is dimensionless; the absolute excess \(C-C_{\mathrm{eq}}\) retains concentration units.

The concentration-based ratio is convenient, but thermodynamics is expressed more precisely in terms of activity. If \(a\) is the activity of dissolved drug and \(a_{\mathrm{eq}}\) is its activity at equilibrium with the specified solid phase, define the activity-based ratio \(S_a=a/a_{\mathrm{eq}}\). In a sufficiently dilute solution where activity coefficients are effectively constant, the activity ratio is approximated by the concentration ratio. When non-ideality matters, the activity-based comparison is the more direct thermodynamic measure.

From Supersaturation to Thermodynamic Driving Force

For the same molecular species and reference solid phase at fixed temperature and pressure, the chemical-potential difference between dissolved drug at activity \(a\) and at equilibrium activity \(a_{\mathrm{eq}}\) is

$$ \mu-\mu_{\mathrm{eq}}=RT\ln\left(\frac{a}{a_{\mathrm{eq}}}\right)=RT\ln S_a, $$

where \(R\) is the gas constant and \(T\) is absolute temperature. This equation states how far the solution’s chemical potential is from its equilibrium value. If \(S_a>1\), the difference is positive: moving drug from solution toward the equilibrium solid state is thermodynamically favorable. If \(S_a<1\), the solution chemical potential is below its equilibrium value, and transfer from solid into solution is favored. At \(S_a=1\), there is no chemical-potential difference driving net transfer.

The result follows from the standard relation \(d\mu=RT\,d\ln a\) at fixed temperature and pressure. Integrating between equilibrium activity and the current activity gives

$$ \mu-\mu_{\mathrm{eq}} =\int_{a_{\mathrm{eq}}}^{a} RT\,d\ln a =RT\ln(a)-RT\ln(a_{\mathrm{eq}}) =RT\ln\left(\frac{a}{a_{\mathrm{eq}}}\right). $$

This is a thermodynamic driving force per mole, with units of energy per mole. It is not a dissolution rate. A kinetic model is needed to translate the driving force into a rate, and the rate may also depend on transport, surface area, or other processes. At low supersaturation, \(\ln S_a\) is approximately \(S_a-1\), but this approximation becomes less accurate as the ratio moves farther from one. The degree of supersaturation is therefore a convenient relative excess, while the logarithm of the activity ratio is the quantity in the chemical-potential expression.

Worked Example: Converting a Supersaturation Ratio to a Driving Force

In a hypothetical dissolution experiment at 310 K, the free dissolved concentration is 0.300 mg/mL and the equilibrium solubility for the reference solid phase is 0.200 mg/mL. Assume an ideal dilute solution so the activity ratio is approximated by the concentration ratio. Calculate the ratio, degree of supersaturation, and chemical-potential difference.

Ratio and degree.

$$ S_C=\frac{0.300\ \mathrm{mg/mL}}{0.200\ \mathrm{mg/mL}}=1.50, \qquad S_C-1=1.50-1=0.500=50.0\%. $$

The solution is 50.0% above its equilibrium solubility, relative to that solubility. For the thermodynamic calculation, \(S_a\approx1.50\). With \(R=8.314\ \mathrm{J\,mol^{-1}\,K^{-1}}\),

$$ \mu-\mu_{\mathrm{eq}} \approx (8.314\ \mathrm{J\,mol^{-1}\,K^{-1}})(310\ \mathrm{K})\ln(1.50) =1.05\times10^3\ \mathrm{J/mol} =1.05\ \mathrm{kJ/mol}. $$

The positive difference indicates a thermodynamic tendency toward the equilibrium state. It does not establish that precipitation will occur during the observation period; the ratio alone contains no information about the rate of that process.

Reading Supersaturation Along a Concentration Profile

A concentration-time curve becomes more informative when each point is compared with the relevant local equilibrium solubility. For a spatially varying system, such as a sequence of gut segments, use the concentration and equilibrium reference for the same location and time. Classify each point using \(S_C=C/C_{\mathrm{eq}}\), rather than comparing every location with one global solubility value.

Worked Example: Classifying a Hypothetical Segmental Profile

Suppose a hypothetical model reports free dissolved concentrations and equilibrium solubilities at three positions. Both quantities are in mg/mL and are assumed to refer to the same drug species and local conditions.

Position\(C\) (mg/mL)\(C_{\mathrm{eq}}\) (mg/mL)\(S_C=C/C_{\mathrm{eq}}\)State
10.1200.2000.600Undersaturated
20.2400.2001.20Supersaturated
30.1800.1501.20Supersaturated

At position 1, \(0.120/0.200=0.600\), so the degree is \(0.600-1=-0.400\), or 40.0% below equilibrium solubility. At position 2, \(0.240/0.200=1.20\), giving a degree of \(0.200\), or 20.0% above equilibrium. At position 3, \(0.180/0.150=1.20\) as well, so it is also 20.0% above its local equilibrium solubility.

Positions 2 and 3 have the same ratio even though their absolute concentrations differ. Their relative excesses above the respective local equilibrium values are the same; their absolute excesses are not: position 2 is \(0.040\ \mathrm{mg/mL}\) above equilibrium, while position 3 is \(0.030\ \mathrm{mg/mL}\) above. Reporting both the ratio and, when useful, the absolute excess preserves these different interpretations.

The label metastably supersaturated describes a solution with \(C>C_{\mathrm{eq}}\) that remains dissolved for a finite period rather than immediately returning to equilibrium. It is still thermodynamically supersaturated: “metastable” does not mean that the driving force has disappeared. Conversely, an observed value at or below equilibrium solubility is not evidence of a supersaturated state merely because the concentration is high in absolute terms.

Choosing the Right Comparison

Three related quantities answer different questions. The ratio \(S_C\) measures concentration relative to equilibrium solubility. The degree \(S_C-1\) expresses the fractional excess or deficit relative to that solubility. The absolute difference \(C-C_{\mathrm{eq}}\) gives the concentration excess or deficit in the original units. For thermodynamic driving force, the activity-based logarithm \(RT\ln S_a\) is the relevant expression.

QuantityExpressionUseful interpretationLimitation
Supersaturation ratio\(C/C_{\mathrm{eq}}\)Relative position above or below equilibriumDepends on a valid, condition-matched reference
Degree of supersaturation\(C/C_{\mathrm{eq}}-1\)Fractional excess or deficit; multiply by 100 for percentNot itself a rate or a chemical-potential difference
Absolute concentration difference\(C-C_{\mathrm{eq}}\)Amount above or below equilibrium in concentration unitsDoes not show relative excess across different solubilities
Activity-based driving force\(RT\ln(a/a_{\mathrm{eq}})\)Chemical-potential difference per moleRequires activities or a justified approximation

Worked Example: Why the Reference Solubility Changes the Classification

In a hypothetical sample, \(C=0.150\ \mathrm{mg/mL}\). Suppose the applicable equilibrium solubility is \(0.100\ \mathrm{mg/mL}\) at the sample’s conditions, while a value of \(0.200\ \mathrm{mg/mL}\) was measured under different conditions. Show how using each reference changes the result.

With the condition-matched reference, \(S_C=0.150/0.100=1.50\), and the degree is \(1.50-1=0.500\), or 50.0%. The sample is supersaturated relative to the applicable equilibrium state. With the mismatched reference, the calculated ratio would be \(0.150/0.200=0.750\), suggesting undersaturation. That classification is not valid for the sample because the denominator does not describe its conditions. The arithmetic is correct, but the comparison is scientifically inappropriate.

Assumptions, Diagnostics, and Common Pitfalls

Before interpreting a calculated ratio, verify that the concentration and equilibrium solubility refer to the same dissolved species, medium, temperature, and relevant solid phase. Also check that the concentration is the concentration that belongs in the comparison. For example, a measured total concentration that includes drug bound or solubilized in another form need not equal the free concentration used to define equilibrium with the solid. The introductory material on equilibrium and apparent solubility distinguishes equilibrium references from apparent measurements; the key requirement here is to use a reference consistent with the quantity being compared.

For a profile, inspect the reference values as well as the measured concentrations. A sudden jump in \(S_C\) can arise from a real concentration change, a change in \(C_{\mathrm{eq}}\), or an inconsistent definition or unit. Ratios very near one also require caution: measurement error in either quantity can change the classification. Report uncertainty or avoid over-interpreting small differences when the available data do not resolve whether \(S_C\) is above or below one.

  • Pitfall: treating high concentration as proof of supersaturation. A high \(C\) can still be below a higher \(C_{\mathrm{eq}}\). Calculate the ratio using the relevant equilibrium solubility.
  • Pitfall: calling \(S_C\) the percentage excess. If \(S_C=1.30\), the degree is \(0.30\), or 30% above equilibrium—not 130% above equilibrium.
  • Pitfall: assuming supersaturation guarantees immediate precipitation. The ratio describes thermodynamic status, not a timescale or a rate. A metastably supersaturated state may persist.
  • Pitfall: using one solubility reference across changing conditions. If equilibrium solubility varies along a profile, calculate the ratio with the local reference at each point.

The concentration-based ratio is often the practical choice when activities are unavailable and dilute-solution behavior is a reasonable approximation. The absolute difference is useful when the amount of excess concentration matters. An activity-based ratio is preferable for a thermodynamic interpretation when non-ideal solution behavior is important. None of these descriptions replaces a kinetic model when the question concerns how quickly concentration changes or how long a supersaturated state persists.

Key takeaway: Supersaturation is a relative state: compare dissolved concentration with the equilibrium solubility for the same conditions and reference solid. The ratio classifies the state, while \(RT\ln S_a\) expresses the activity-based chemical-potential difference. Neither quantity alone specifies a kinetic rate.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. A hypothetical solution has \(C=0.090\ \mathrm{mg/mL}\) and applicable \(C_{\mathrm{eq}}=0.120\ \mathrm{mg/mL}\). Calculate \(S_C\), the degree as a percentage, and classify the state.
  2. At \(C=0.250\ \mathrm{mg/mL}\), two locations have equilibrium solubilities of 0.200 and 0.300 mg/mL. Calculate and compare their ratios.
  3. Why does \(S_C>1\) not establish that precipitation occurs immediately?
  4. State the activity-based chemical-potential difference for a solution with \(S_a=0.80\), and give its sign and interpretation.

Solutions. (1) \(S_C=0.090/0.120=0.750\). The degree is \(0.750-1=-0.250\), or 25.0% below equilibrium; the state is undersaturated. (2) The first ratio is \(0.250/0.200=1.25\), so it is supersaturated by 25.0%. The second is \(0.250/0.300=0.833\), so it is undersaturated by about 16.7%. The same concentration has different classifications because the equilibrium references differ. (3) \(S_C\) indicates thermodynamic status, not a kinetic rate or timescale; a solution can remain metastably supersaturated. (4) \(\mu-\mu_{\mathrm{eq}}=RT\ln(0.80)\), which is negative because \(\ln(0.80)<0\). The solution’s chemical potential is below its equilibrium value, so transfer from the solid into solution is thermodynamically favored.