Introduction & Context

Heterogeneous nucleation describes the formation of a new solid phase (e.g., ice) on a pre‑existing surface or impurity within a supersaturated liquid. In process engineering—particularly in frozen‑food production, cryopreservation, and atmospheric science—predicting the nucleation rate is essential for controlling crystal size distribution, product texture, and safety. The calculation presented here quantifies the energy barrier that must be overcome for a critical ice nucleus to form on a substrate characterized by a contact angle θ, and it converts that barrier into a nucleation rate J (units m⁻³·s⁻¹).

Methodology & Formulas

The procedure follows the classical nucleation theory (CNT) adapted for heterogeneous conditions. All quantities are expressed in SI units unless otherwise noted.

Temperature conversion – Convert the supplied temperature in degrees Celsius (T_C) to absolute Kelvin: \[ T_K = T_C + 273.15 \]
Contact‑angle conversion – Transform the contact angle from degrees (θ_deg) to radians for trigonometric functions: \[ \theta = \theta_{\text{deg}} \frac{\pi}{180} \]
Volumetric driving force – The thermodynamic driving force per unit volume for phase change is: \[ \Delta G_v = -\frac{R\,T_K}{V_m}\,\ln(\beta) \] where R is the universal gas constant, V_m the molar volume of the solid phase, and β the supersaturation ratio (dimensionless). The logarithm must be evaluated with a positive argument; a safeguard against non‑physical values is implied.
Homogeneous critical free‑energy barrier – For nucleation without a substrate: \[ \Delta G_{\text{hom}} = \frac{16\pi\,\sigma^{3}}{3\,\Delta G_v^{2}} \] where σ is the interfacial tension between the solid and liquid.
Geometric reduction factor – The presence of a surface reduces the barrier by a factor that depends only on the contact angle: \[ f(\theta) = \tfrac{1}{4}\,\bigl(2+\cos\theta\bigr)\,\bigl(1-\cos\theta\bigr)^{2} \]
Heterogeneous energy barrier – Multiply the homogeneous barrier by the reduction factor: \[ \Delta G_{\text{het}} = \Delta G_{\text{hom}}\;f(\theta) \]
Nucleation rate – Apply an Arrhenius‑type expression with a pre‑exponential factor J₀: \[ J = J_0 \,\exp\!\left(-\frac{\Delta G_{\text{het}}}{k_B\,T_K}\right) \] where k_B is Boltzmann’s constant.

Validity Checks & Applicability Limits

Parameter	Recommended Range	Purpose of Check
Temperature (T_C)	\(-5 \le T_C \le 0\) °C	Ensures the liquid is in the typical frozen‑food regime where ice nucleation is relevant.
Supersaturation ratio (β)	\(1 < \beta \le 1.5\)	Maintains the validity of the logarithmic driving‑force expression for modest supersaturation.
Contact angle (θ_deg)	\(0^\circ < \theta_{\text{deg}} < 180^\circ\)	Prevents undefined trigonometric values and ensures a physically meaningful reduction factor.
Interfacial tension (σ)	\(0.025 \le \sigma \le 0.035\) J·m⁻²	Reflects experimentally observed ice–water interfacial energies.
Volumetric driving force (ΔG_v)	\(\Delta G_v < 0\)	Negative sign indicates a thermodynamically favorable transition; a non‑negative value signals inconsistent input.

Result Presentation (Typical Output)

The script prints the computed quantities rounded to three decimal places: temperature in Kelvin, supersaturation ratio, contact angle in degrees, the volumetric driving force ΔG_v, the homogeneous barrier ΔG_hom, the geometric factor f(θ), the heterogeneous barrier ΔG_het, and finally the nucleation rate J.

In heterogeneous nucleation, the new phase forms on a pre‑existing surface (solid particle, wall, or impurity). The presence of this surface reduces the interfacial area that must be created, lowering the free‑energy cost. Consequently, the energy barrier ΔG*_het is a fraction of the homogeneous barrier ΔG*_hom, typically expressed as ΔG*_het = f(θ)·ΔG*_hom, where f(θ) is a geometric factor that depends on the contact angle θ between the nucleus and the substrate.

The contact angle can be obtained by:

Direct optical measurement on a representative substrate using a goniometer.
Inferring θ from known surface energies via Young’s equation: γ_sv = γ_sl + γ_lvcosθ.
Using literature values for similar material combinations when direct measurement is impractical.

For process engineers, the most common approach is to use published data for the specific catalyst or wall material and the fluid of interest.

The barrier is affected by:

Interfacial tensions (γ_lv, γ_sl, γ_sv).
Supersaturation or supercooling level, which determines the driving force Δμ.
Surface roughness and chemistry, which modify the effective contact angle.
Particle size and curvature of the nucleating surface.
Presence of contaminants that can either promote or inhibit nucleation.

To lower ΔG*_het and encourage nucleation:

Introduce finely divided seed particles with a low contact angle for the target phase.
Increase supersaturation by adjusting temperature or concentration gradients.
Modify surface chemistry (e.g., coating reactor walls) to make them more wettable.
Employ agitation or ultrasonic fields to create additional nucleation sites.

Conversely, to suppress unwanted nucleation, you can increase the contact angle by using inert coatings or reducing supersaturation.

Worked Example – Estimating the Heterogeneous Nucleation Barrier for Ice Formation on a Chilled Stainless-Steel Valve

A process engineer is validating the start-up procedure for a cryogenic water line that will be cooled to −5 °C. To avoid blockage, the team must confirm that ice crystals will nucleate readily on the 30°-contact-angle stainless-steel valve seat rather than in the bulk. The key design check is the reduction in nucleation energy compared with homogeneous freezing.

Knowns

Temperature: \( T = -5\,^\circ\mathrm{C} = 268.15\,\mathrm{K} \)
Water–ice interfacial energy: \( \sigma = 0.030\,\mathrm{J\,m^{-2}} \)
Molar volume of ice: \( V_\mathrm{m} = 1.96\times10^{-5}\,\mathrm{m^{3}\,mol^{-1}} \)
Contact angle on steel: \( \theta = 30^\circ \)
Pre-exponential factor: \( J_{0} = 1\times10^{36}\,\mathrm{m^{-3}\,s^{-1}} \)
Thermodynamic beta: \( \beta = 1.2 \)

Step-by-Step Calculation

Compute the volumetric free-energy driving force: \[ \Delta G_\mathrm{v} = -\frac{\beta\,R\,T}{V_\mathrm{m}} = -\frac{1.2\times8.314\times268.15}{1.96\times10^{-5}} = -2.074\times10^{7}\,\mathrm{J\,m^{-3}} \]
Evaluate the homogeneous nucleation barrier: \[ \Delta G_\mathrm{hom} = \frac{16\pi}{3}\frac{\sigma^{3}}{(\Delta G_\mathrm{v})^{2}} = \frac{16\pi}{3}\frac{(0.030)^{3}}{(-2.074\times10^{7})^{2}} = 1.052\times10^{-18}\,\mathrm{J} \]
Calculate the geometrical correction factor for a 30° contact angle: \[ f(\theta) = \frac{1}{4}(2+\cos\theta)(1-\cos\theta)^{2} = \frac{1}{4}(2+\cos30^\circ)(1-\cos30^\circ)^{2} = 0.0129 \]
Determine the heterogeneous barrier: \[ \Delta G_\mathrm{het} = f(\theta)\,\Delta G_\mathrm{hom} = 0.0129\times1.052\times10^{-18} = 1.353\times10^{-20}\,\mathrm{J} \]
Estimate the nucleation rate on the valve seat: \[ J = J_{0}\exp\!\left(-\frac{\Delta G_\mathrm{het}}{k_\mathrm{B}T}\right) = 1\times10^{36}\exp\!\left(-\frac{1.353\times10^{-20}}{1.38\times10^{-23}\times268.15}\right) = 2.58\times10^{34}\,\mathrm{m^{-3}\,s^{-1}} \]

Final Answer

The heterogeneous nucleation energy barrier on the 30°-contact-angle steel surface is 1.35 × 10⁻²⁰ J per critical cluster, two orders of magnitude lower than the homogeneous barrier. The corresponding nucleation rate exceeds 10³⁴ m⁻³ s⁻¹, confirming rapid ice formation on the valve seat at −5 °C.

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle