Mathematical Derivation of Subsurface Light Scattering Factors in Highly Filled Uncoated Mechanical Furnishes

Flux

Light transport within a turbid cellulose matrix operates through continuous absorption and isotropic redistribution across optical phase boundaries. A photon entering an uncoated sheet encounters high-yield wood fibres rich in native lignin, mineral filler aggregates, and microscopic air voids created by incomplete conformability during web consolidation. Classical radiative transfer describes this transport across an infinitesimal thickness element dz of basis weight dw, measured in grams per square metre under ISO 536.

The balance of downward diffuse radiant energy i and upward diffuse radiant energy j inside the sheet follows two coupled linear differential equations:

-di / dw = -(s + k) i + s j

dj / dw = -(s + k) j + s i

The parameter s denotes the specific Kubelka-Munk scattering coefficient in square metres per kilogram, and k denotes the specific absorption coefficient in the same units. Integrating these equations across an optical depth from the illuminated top surface where w = 0 to the backing interface where w = W requires defining the physical reflectance R as the quotient j / i. Dividing the differential relations yields the Riccati differential equation for local reflectance:

dR / dw = s – 2(s + k)R + s R² = s

Setting the parameter a = 1 + k/s and b = (a² – 1)^1/2 = ^1/2 allows direct separation of variables. When sheet basis weight approaches optical infinity, reflectance reaches R_∞, where the differential gradient dR / dw equals zero:

s = 0

Solving this quadratic relation yields the standard expression for infinite sheet reflectance:

R_∞ = a – b = 1 + (k/s) – ^1/2

Inversion of this equation produces the classic relationship between light absorption, light dispersion, and observed infinite sheet reflectance:

k / s = (1 – R_∞)² / (2 R_∞)

Mechanical pulps contain between twenty-four and thirty percent native lignin by dry weight, giving thermomechanical pulp an unbleached absorption factor k exceeding 6.5 m²/kg at 457 nm under ISO 2470 test conditions. The resulting reflectance values reflect strong blue absorption bands. Caliper gauges read thirty micrometres.

Light enters the sheet surface.

Uncoated mechanical furnishes conditioned at 23 degrees Celsius and 50 percent relative humidity lose four percent absolute opacity when moisture expands unbonded fibre interfaces.

Determining the true specific scattering coefficient of a finite sheet of basis weight W over an ideal black backing where background reflectance R₀ = 0 requires integrating the Riccati expression from w = 0 to w = W:

1. Boundary integration step sets the definite integral limits from zero reflectance over the zero-reflecting trap to the measured single-sheet reflectance factor R₀.
2. Hyperbolic transformation step converts algebraic roots into hyperbolic functions through the substitution b s W = Arsh .
3. Scattering isolation step extracts the absolute factor s by isolating basis weight into the standard formulation s = (1 / b W) ln{ / }.
4. Background correction step incorporates finite backing reflectance R_g under ISO 2471 when test pieces run over white calibration tiles.

The hyperbolic solutions express reflectance R as a direct function of total basis weight W, intrinsic infinite reflectance R_∞, and underlying background reflectance R_g:

R = /

When measuring an unbacked sheet where R_g = 0, this relationship condenses to:

R₀ = sinh(b s W) /

Solving directly for the subsurface factor yields the operational Kubelka-Munk equation used on mill spectrophotometers:

s = (1 / b W) Arctanh

These classical formulations treat the furnish as an optical continuum without discrete structural phase boundaries. Real mechanical sheets contain dense fibrillar bundles, ray cells, mineral particles, and micro-voids that alter actual photon paths through refraction and diffraction.

Ash

Mineral pigmentation alters light distribution by multiplying internal refractive discontinuities per unit volume. In an uncoated mechanical sheet running between twenty and thirty-eight percent mineral filler, the total specific scattering coefficient deviates from simple linear mixture rules. The additive model proposed by standard furnish calculations assumes each component scatters independently according to its mass fraction:

s_furnish = x_f s_f + x_m s_m + x_p s_p

Here x_f, x_m, and x_p denote mass fractions of chemical pulp, mechanical pulp, and pigment ash determined by combustion at 525 degrees Celsius under ISO 1762. Ash content reaches thirty percent. Clay platelets orient along machine direction.

Pigment interfaces diffuse the photons.

Linear models fail at commercial filler additions because mineral particles induce two conflicting physical mechanisms within the sheet structure. First, mineral grains prevent adjacent cellulose fibrils from consolidating during wet pressing, creating unbonded fibre surface area that scatters additional light. Second, when filler loading exceeds eighteen percent by mass, pigment particles agglomerate, causing optical crowding.

Inter-particle distances fall below half the wavelength of visible light, reducing scattering cross-sections below theoretical values predicted by Mie theory for isolated particles. Groundwood fibres resist chemical swelling.

Accounting for crowding and debonding demands an expanded derivation where the composite factor incorporates an interaction parameter s_interaction:

s_sheet = (1 – x_ash) s_pulp + x_ash s_ash + s_debond(x_ash) – s_crowd(x_ash)

The debonding term s_debond relates to the unbonded fibre area generated per gram of retained filler. Nitrogen adsorption measurements under the Brunauer-Emmett-Teller method demonstrate that unbonded specific surface area A_free increases proportionally with filler content up to a threshold x_crit:

s_debond = σ_fibre ΔA_free = σ_fibre

Here σ_fibre represents the specific scattering factor of pure cellulose-lignin surfaces, roughly 0.038 m²/g for spruce groundwood, and x_max represents maximum structural packing where the matrix loses tensile cohesion. Mechanical fibres carry residual lignin.

Particle crowding follows modified van de Hulst approximations for dependent scattering in concentrated suspensions. As the volume fraction of filler φ rises, spatial correlation functions create destructive phase interference between forward-scattered wavelets:

s_crowd = x_ash s_ash

Mineral suppliers claim that high-shear slurry dispersion eliminates dependent scattering losses across all commercial retention thresholds.

Interface

Boundary conditions at external and internal sheet surfaces dictate how much diffuse illumination penetrates the fibrous core. Spectrophotometers measure external reflectance factors, designated as R_measured, under diffuse or directed geometries such as d/0 or 45/0 specified in ISO 2469. Raw paper surfaces exhibit an abrupt change in refractive index between ambient air (n₀ = 1.000) and the wet-pressed cellulose and filler web (n_sheet ≈ 1.540).

Photons crossing this boundary undergo external Fresnel reflection before entering the bulk furnish, while scattered light traveling outward encounters total internal reflection at angles exceeding the critical angle θ_c = arcsin(1 / n_sheet) ≈ 40.5°.

The Saunderson relationship reconciles theoretical Kubelka-Munk reflectance R_KM with observed reflectance values R_measured:

R_measured = k₁ + /

The constant k₁ represents external surface reflectance for incident light, while k₂ represents internal surface reflectance for diffuse backscattered light. For smooth polymeric surfaces, Fresnel equations yield theoretical values of k₁ ≈ 0.040 for collimated normal incidence and k₂ ≈ 0.600 for isotropic internal radiation. Uncoated mechanical papers, however, possess irregular micro-topography.

Surface roughness, measured via Parker Print-Surf under ISO 8791-4 at 1000 kPa clamping pressure, distorts these reflection factors. Wet pressing packs the furnish tighter.

Standard supply agreements specifying ISO brightness under ISO 2470 permit batch rejections when delivered sheets show a two-unit drop in directional reflectance.

On supercalendered mechanical grades like SC-A and SC-B, surface fibres and mineral platelets experience intense shear and compressive plastic deformation inside multi-roll calender nips. Calendering collapses light scattering voids. Sheet opacity drops after calender nips.

These surface distortions complicate optical verification across high-ash mechanical deliveries:

• Micro-facet tilt divergence scatters specular gloss beams into diffuse collection apertures, falsifying absorption values calculated from low-angle spectrophotometers.
• Pore collapse during densification eliminates sub-micron air cavities within three micrometres of the sheet surface, driving local internal reflectance k₂ above 0.680.
• Two-sided mineral distribution produces an asymmetric scattering profile between the wire side and top side of fourdrinier-formed webs.
• Refractive matching by processing oils occurs when offset fountain solutions or calender lubricants penetrate surface-bound filler aggregates, suppressing refractive index mismatches.

Evaluating the internal reflectance constant k₂ requires solving the Saunderson relation inversely against empirical measurements over certified standard backings:

R_KM = (R_measured – k₁) /

Substituting this corrected R_KM into the Kubelka-Munk derivation yields true volumetric scattering factors unaffected by boundary gloss. The operational dispute centers on whether k₁ and k₂ remain constant across varying filler ratios or shift dynamically with localized surface mineral concentration.

Coupling

Extending two-flux models to multi-phase mechanical papers requires linking structural morphology to local optical cross-sections through radiative transfer theory. Classical Kubelka-Munk theory assumes isotropic diffuse flux distributions. In dense mechanical sheets, forward scattering predominates at individual fibre-pigment boundaries.

Chandrasekhar’s radiative transfer equation models this directional angular distribution across polar angle θ and azimuthal angle ψ:

μ (dI(τ, μ, ψ) / dτ) = I(τ, μ, ψ) – (ω₀ / 4π) ∫_-1¹ ∫₀^2π p(μ, ψ; μ’, ψ’) I(τ, μ’, ψ’) dψ’ dμ’

Here μ = cos θ, τ denotes optical thickness, I represents specific radiance, ω₀ is the single-scattering albedo defined as σ_s / (σ_s + σ_a), and p(μ, ψ; μ’, ψ’) represents the scattering phase function. For isotropic scattering, p equals unity, restoring two-flux symmetry. In mechanical sheets filled with scalenohedral precipitated calcium carbonate, the phase function elongates along the forward vector.

The asymmetry factor g quantifies this distortion:

g = (1/2) ∫_-1¹ p(cos Θ) cos Θ d(cos Θ)

The Henyey-Greenstein analytical phase function describes forward-peaked angular scattering in terms of this asymmetry factor:

p(cos Θ) = (1 – g²) / ^3/2

Pure mechanical groundwood displays an asymmetry factor g ≈ 0.72, while precipitated calcium carbonate particles exhibit g ≈ 0.85. The similarity principle links radiative transfer transport parameters directly to Kubelka-Munk coefficients:

s = (3/4) σ_s (1 – g)

k = 2 σ_a

This derivation establishes that the Kubelka-Munk factor s represents a reduced transport scattering coefficient governed entirely by the backscattering fraction (1 – g). Fines migrate toward the wire. Sizing alters surface pore geometry.

Lignin absorbs blue wavelengths heavily.

Deriving the subsurface scattering parameters of filled mechanical papers demands tracking the reduction of backscattering efficiency as filler loading rises. Increasing mineral mass packs filler particles into agglomerated clusters, driving the asymmetry factor g upward from 0.74 toward 0.83. This optical forward-peaking cancels out gains in specific surface area.

Highly calendered furnishes show divergence between calculated theoretical scattering and measured sheet performance. High ash fractions consistently depress sheet resilience during long storage runs.

Balance

Executing grade substitutions between low-ash standard mechanical sheets and super-filled stocks demands balancing optical performance against tensile integrity and commercial yield. In an uncoated 48 g/m² mechanical web, substituting unbleached thermomechanical pulp with scalenohedral precipitated calcium carbonate shifts sheet physics across three distinct operational fronts: specific light scattering, web thickness, and tensile index under ISO 1924-2. Mineral filler costs less per dry tonne than virgin wood furnish, tempting mill managers to maximize ash content.

Filler additions above twenty-five percent introduce severe structural penalties.

Every additional percentage point of mineral ash raises the specific scattering coefficient by roughly 1.2 m²/kg, but degrades cross-machine tensile energy absorption by nearly three percent. Mineral particles disrupt hydrogen bonding between mechanical pulp fibrils. In rotogravure and high-speed coldset offset printing, this bond loss triggers web breaks across press folders.

Groundwood sheets calendered to high density display catastrophic linting when filler particles detach from the sheet surface under tacky ink splitting forces. Mill certificates report ISO brightness.

Calendering pressure must be limited to prevent structural pin-holing in lightweight mechanical webs.

The commercial equation links optical opacity targets to finished sheet yield. Purchasing paper by the metric tonne and converting it into printed area establishes an inverse yield relationship governed by apparent density:

Yield (m² / tonne) = 1,000,000 / (Basis Weight in g/m²)

Increasing filler content from twelve to thirty percent raises apparent sheet density from 0.55 g/cm³ to 1.10 g/cm³. To preserve bending stiffness and prevent visual show-through under high ink coverage, converters cannot simply downgrade basis weight. If a buyer reduces nominal basis weight to recover yield while maintaining elevated ash fractions, opacity targets can be met via enhanced scattering, but sheet caliper drops precipitously.

Caliper losses compromise structural runnability in folder superstructures, leading to web snaps, press downtime, and pallet rejection across printing plants.