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]

Drug release kinetics describes how the amount leaving a formulation changes with time. Let $M_0$ be the initially releasable amount and $M_{\mathrm{rel}}(t)$ the cumulative amount released. The released fraction and instantaneous drug release rate areA fraction has no units; a rate has units of amount per time. These definitions also require stating the releasable amount used for normalization, rather than silently identifying it with every molecule initially loaded.

Several simple models give distinct curve shapes. Constant-rate release gives $M_{\mathrm{rel}}=r_{0}t$ until the available amount is exhausted. First-order depletion of a reservoir gives $F=1-e^{-kt}$, with an amount released per time that declines as the reservoir empties. In a suitable early-time diffusion-controlled drug release regime, a model may instead give $F=a\sqrt{t}$. Thenso equal increments of released fraction do not occur in equal time intervals. The square-root expression is a limited-time model, not a law that applies from arbitrarily small times through complete release.

For an illustrative curve $F(t)=0.10\sqrt{t/\mathrm{h}}$, the fractions at one, four and nine hours are $0.10$, $0.20$ and $0.30$. The mean fraction released per hour in the intervals $0$–$1$, $1$–$4$ and $4$–$9$ hours is respectively $0.10$, $0.0333$ and $0.020$. A plot with equally spaced vertical increments can therefore conceal a strongly changing release rate.

The mechanism must be tested, not inferred solely from a convenient curve fit. Diffusion, material erosion, changing geometry and the measurement apparatus can produce overlapping features. The historical polymer sustained-release work and the modern FDA assay account illustrate the separate tasks of designing a carrier and measuring what it releases.^[ref][ref] A defensible comparison states the time range, normalization, sampling correction and physical assumptions, then checks whether the proposed mechanism is consistent with additional observations.

A perfect sink condition in drug release is an idealized receiving boundary that removes released drug fast enough for its concentration there to be treated as zero. In diffusion-controlled drug release, this preserves the largest concentration drop available at the exposed boundary. It is a condition on transport at the release surface, rather than a statement that the entire apparatus contains no drug.

For a planar matrix occupying $x>0$, a common diffusion equation model useswhere $c$ is dissolved-drug concentration and $D$ is its diffusion coefficient; the outward flux magnitude is used when reporting a positive drug release rate. A finite mass-transfer boundary in drug release instead allows a nonzero surface concentration. If the same concentration scale applies on both sides of the interface, an external transport model iswith transfer coefficient $h$, surface concentration $c_s$ and well-mixed bulk concentration $c_b$. A small $c_b$ does not force $c_s$ to be small when a stagnant layer or membrane imposes substantial resistance.

An original finite-reservoir example separates capacity from transport. Suppose a formulation can release $10\mathrm{mg}$ into $100\mathrm{mL}$, and measured solubility in the chosen medium is $1mg/mL$. Even after complete release the well-mixed bulk concentration is only $0.10mg/mL$. In a simple saturated-interface dissolution model $J=h(c_{\mathrm{sat}}-c_b)$, that concentration still reduces flux to $90\%$ of its zero-bulk-concentration value. Large capacity makes a zero-concentration approximation plausible to a chosen tolerance; it does not make the approximation exact.

A useful experimental test compares in vitro drug release test curves while increasing receiving volume or flow, with temperature and medium chemistry controlled. Sampling and medium replacement require sampling correction in drug release testing. If increasing external transport changes the curve, apparatus resistance or accumulation may still matter. The FDA’s primary account of adaptive perfusion describes pressure-driven filtration designed to avoid rate limitation by a passive dialysis membrane.^[ref]

The ideal boundary should be distinguished from a method’s practical sink criterion. The FDA’s 1997 dissolution guidance discusses medium, apparatus and agitation together and explicitly allows justified non-sink tests; a zero boundary is a modelling choice to be checked against the actual experiment.^[ref]

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.