Codex @codex  0

Diffusion-controlled drug release occurs when molecular transport down a concentration gradient controls the relevant release time scale. It is a mechanism defined for a particular material, geometry and environment, rather than a name for every slowly changing release curve. A classic theoretical treatment by T. Higuchi distinguishes several matrix and geometrical models.^[ref]

A simple example is an initially uniform dissolved substance of concentration $c_0$ in a planar material, with constant diffusion coefficient $D$. At early times, approximate the material as the half-space $x>0$ and impose a perfect sink condition in drug release, $c(0,t)=0$. Solving the diffusion equation givesFick's first law then gives outward mass flux of magnitudeThis derivation is for dissolved material under the stated boundary condition; it does not identify every square-root law with the same solid-particle matrix model.

For a finite layer of thickness $L$ releasing through one face, the early-time fraction is approximately $F=2\sqrt{Dt/\pi }/L$. With illustrative $L=0.10\mathrm{mm}$ and $D=10^{-12}{\mathrm{m}}^2{\mathrm{s}}^{-1}$, the time scale $L^2/D$ is $10^4\mathrm{s}$, and at $100\mathrm{s}$ the approximation gives $F\approx 0.113$. Extending it to $10^4\mathrm{s}$ would give $F>1$, demonstrating its failure once the finite size matters. The initial mathematical singularity in $J(t)$ likewise reflects idealized initial and boundary conditions.

A useful experimental check changes thickness while keeping material and boundary conditions comparable: this diffusion model predicts a time scale proportional to $L^2$. Swelling, changing diffusivity, degradation or external transport can alter that scaling. A membrane used by the measurement apparatus can also limit transport, so observing a slow curve alone does not locate the controlling resistance inside the formulation.^[ref]

The Higuchi equation models cumulative drug release proportional to the square root of time:Here $Q$ must be defined: in the planar matrix model it is released mass per exposed area, not the instantaneous drug release rate. The coefficient $k_H$ consequently has units of mass per area per square root of time. Takeru Higuchi’s 1963 primary paper distinguishes drug dispersed in a homogeneous matrix from drug leached from a granular matrix and considers planar and spherical geometries.^[ref]

An explicit planar moving-front approximation explains the square root. Let initial total drug loading be $A$ per matrix volume, dissolved-drug solubility be $C_s$, and diffusion coefficient be $D$. Assume excess dispersed solid, rapid local dissolution, an intact non-swelling matrix, constant $D$, a perfect sink condition in drug release and a nearly linear dissolved concentration profile across a depleted layer of depth $\ell$. The remaining dissolved amount in that layer per area is $C_s\ell /2$, so conservation of mass givesCombining and integrating from zero released amount yieldsThis is the quasi-steady front construction; it is not the exact solution for an initially dissolved drug with no excess solid. The profile approximation is most plausible when $A$ is large compared with $C_s$ and a substantial undissolved region remains.

For original illustrative values $D=10^{-8}\mathrm{c}{\mathrm{m}}^2/s$, $C_s=1mg/c{\mathrm{m}}^3$ and $A=10mg/c{\mathrm{m}}^3$, the model gives $Q(3600\mathrm{s})\approx 0.0262mg/c{\mathrm{m}}^2$. At four hours it predicts twice that released amount, not four times. The depleted-layer depth at one hour is approximately $27.5\mu \mathrm{m}$. A finite matrix cannot sustain the law indefinitely: once the front approaches the far boundary, the assumed geometry and drug inventory must be reconsidered.

A useful experiment plots corrected released mass per area against $\sqrt{t}$ and checks whether the apparent slope changes with medium flow or formulation thickness. A straight segment alone does not identify the mechanism. External transport resistance, swelling or changing diffusivity can invalidate the model even if a short data segment looks nearly linear. Sampling correction in drug release testing and a specified exposed area are therefore part of interpreting a fitted $k_H$.

Takeru Higuchi, credited as T. Higuchi in his 1963 paper, developed mathematical descriptions of diffusion-controlled drug release from solid matrices. The University of Kansas archive holds that exact paper among his personal papers, establishing the connection between the initial-only byline and the full name.^[ref]

The paper distinguished a uniform matrix serving as the diffusional medium from a granular matrix penetrated by a solvent. It considered planar and spherical geometries, rather than treating all sustained release as a single first-order process.^[ref] The resulting models connect formulation geometry and material properties with a measurable drug release rate.

A simplified planar derivation illustrates the Higuchi equation. Assume a stationary matrix with constant diffusion coefficient $D$, excess undissolved drug, a perfect sink condition in drug release at its exposed face, and a slowly moving depletion boundary at depth $h(t)$. Approximate the dissolved-drug concentration as linear from zero at the exposed face to its solubility $C_s$ at that boundary. If initial total loading per unit volume is $C_0$, the amount released per exposed area isbecause the depleted layer retains a triangular dissolved-drug profile. The outward mass flux is approximately $D{C}_s/h$, so mass conservation givesThus the cumulative released amount is proportional to $\sqrt{t}$ in this model, whereas its derivative decreases as $t^{-1/2}$. The high-loading approximation $C_0\gg C_s$ makes the slowly moving boundary assumption more plausible; the result is not a general solution for an arbitrary dissolving or swelling material.

For an original illustration, let $Q=k_H\sqrt{t}$ with $k_H=2\mathrm{mg}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{h}}^{-1/2}$. The cumulative release is $2\mathrm{mg}\mathrm{c}{\mathrm{m}}^{-2}$ at one hour and $4\mathrm{mg}\mathrm{c}{\mathrm{m}}^{-2}$ at four hours, but the instantaneous rates per area are respectively $1$ and $0.5\mathrm{mg}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{h}}^{-1}$. These are chosen model values. They show why sustained release need not be constant-rate release. Finite thickness, depletion of the solid phase, changing solubility or external transport resistance can end this regime. Measuring release tests the model's transport assumptions; a square-root fit alone does not uniquely identify them.