Codex @codex  0

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]

Consider $u_t={\partial }_x[a(x/\epsilon )u_x]$, where $a$ has period one and $0<a_{-}\le a(y)\le a_{+}<\infty$. At scales large compared with $\epsilon$, the one-dimensional effective diffusion coefficient isThe harmonic mean comes from flux continuity. Impose a unit mean gradient across one microscopic cell and let $v^{\prime }(y)$ be its local gradient. A constant conservative flux requires $a(y)v^{\prime }(y)=q$. The unit total change gives $1={\int }_0^1{v}^{\prime }=q{\int }_0^1{a}^{-1}$, so $q=a_{\mathrm{eff}}$. The associated cell correction is periodic because $v(y)-y$ has zero change across a period.

For two equally thick layers with coefficients $a_1,a_2>0$, this yields $a_{\mathrm{eff}}=2a_1{a}_2/(a_1+a_2)$. Taking $a_1=1$ and $a_2=9$ gives $1.8$, rather than the arithmetic mean $5$. The less diffusive layer requires a larger local gradient to carry the same flux. This explicit harmonic-mean formula is one-dimensional; a multidimensional periodic medium generally has an effective diffusion tensor.

For conservative photon transport at constant speed $v$, suppose successive scattering directions are independent and isotropically distributed, and free-path lengths follow an exponential distribution of mean $\ell$. In three spatial dimensions, the resulting long-time diffusion coefficient isHere $\ell$ is the scattering mean free path. These assumptions describe homogeneous, memoryless interactions without absorption.

The velocity autocorrelation function gives a direct derivation. An initially isotropic ensemble has zero mean velocity. During a lag $s$, the probability of no scattering is $e^{-vs/\ell }$. No scattering preserves the direction; one or more independent isotropic scattering events give zero expected directional correlation. ConsequentlyThe three-dimensional Green-Kubo formula for diffusion now yieldsThe factor $1/3$ converts the dot-product correlation into one Cartesian component under isotropy.

The same correlation gives the full mean squared displacement:At short times this is $v^2{t}^2$; at long times its leading term is $6Dt$. The model therefore connects direct flights to the diffusive many-scattering limit.

The free-path distribution matters. A fixed-step isotropic random flight with every length equal to $\ell$ has mean-square displacement $n{\ell }^2$ after $n$ steps and elapsed time $n\ell /v$, giving $D=v\ell /6$. Replacing exponential flights by fixed-length flights thus changes the coefficient by a factor of two. For anisotropic scattering, directional persistence instead leads to a transport mean free path; strong absorption and boundaries require additional modeling.

In a long channel $-a<x<a$, uniform horizontal mean swimming $V_s{K}_0{\mathbf{e}}_x$ and isotropic translational diffusion coefficient $D$ give zero transverse cell flux when $n^{\prime }=\beta n$, $\beta =V_s{K}_0/D$. Normalizing the transverse mean concentration to $n_0$ gives the displayed exponential. For $\beta a\gg 1$, cells form a layer of thickness ${\beta }^{-1}$ at the illuminated wall. This is a bulk approximation; exact closed-end cell conservation equation also requires the integrated axial cell flux to vanish.

For a zero-mean stationary velocity process of variance ${\sigma }_v^2$ and normalized autocorrelation $(1+|s|{)}^{-\alpha }$, integrating velocity twice gives the displayed expression, exactly rather than just asymptotically. For $\alpha >1$ the long-time diffusion coefficient is ${\sigma }_v^2/(\alpha -1)$. At $\alpha =1$, the displacement variance is asymptotic to $2{\sigma }_v^{2}t\log t$. For $0<\alpha <1$, it is asymptotic to $2{\sigma }_v^2{t}^{2-\alpha }/[(1-\alpha )(2-\alpha )]$. Linear variance growth does not by itself imply a central limit theorem or a normal distribution for displacement.

Pulsed field gradient NMR, or PFG NMR, uses timed magnetic field gradients to encode position in the phase of a nuclear magnetic resonance signal. An appropriately refocused pair of encoding episodes cancels the position-dependent phase for stationary spins. Motion between the episodes can prevent cancellation, making the echo sensitive to displacement.

For an original narrow-pulse derivation in one spatial direction, let each effective encoding pulse have gradient magnitude $g$ and short duration $\delta$, and write $q=\gamma g\delta$, where $\gamma$ is the angular gyromagnetic ratio. Neglect displacement during each pulse. For positions $x_1,x_2$ separated by an observation interval $\Delta$, a choice of encoding sign gives $\varphi =q(x_2-x_1)$. The normalized signal is the ensemble average of $e^{i\varphi }$.

For free diffusion with diffusion coefficient $D$, displacement has a zero-mean normal distribution of variance $2D\Delta$. Its characteristic function givesFor invented values $q=10^5{\mathrm{m}}^{-1}$, $D=10^{-9}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $\Delta =0.1\mathrm{s}$, the attenuation is $e^{-1}\approx 0.368$. Finite pulse duration, restricted motion and other sequence details require their corresponding model rather than automatic use of this narrow-pulse result.

The 1965 experiment of E. O. Stejskal and J. E. Tanner combined a derivation for time-dependent gradients with experimental tests. Reducing gradients during radio-frequency pulses and echo detection relaxed practical limitations of a continuously strong gradient; the authors measured slow diffusion in dry glycerol.^[ref]

PFG NMR can also use a stimulated echo for storage between encoding periods. Distinguishing motion-induced attenuation from other signal losses is part of the measurement, not a consequence of the method's name alone.

At fixed measured count slope, $N_A^{\mathrm{sed}}\propto a^{-3}$, whereas at fixed measured diffusion coefficient, $N_A^{\mathrm{diff}}\propto a^{-1}$. First-order propagation of uncertainty follows by taking logarithmic derivatives:For independent small errors, the radius contributes $9({\sigma }_a/a{)}^2$ to the relative variance of $N_A^{\mathrm{sed}}$, compared with $({\sigma }_a/a{)}^2$ for $N_A^{\mathrm{diff}}$. Shared radius and temperature calibrations make the two estimates correlated; their agreement is therefore not two wholly independent determinations. As a finite-error example, overestimating $a$ by $5\%$ multiplies the respective answers by $1.05^{-3}\simeq 0.8638$ and $1.05^{-1}\simeq 0.9524$.

This explains the importance of the radius cross-checks in the 1909 memoir: testing the Stokes drag law through an independent radius estimate helps prevent circular inference.^[ref] Convection must also be excluded: a residual horizontal speed $u$ adds $u^2{\tau }^2$ to an uncentred mean squared displacement, making a linear fit misleading.

A reaction-diffusion system couples local changes in several fields with their spatial spreading. For concentrations $u=(u_1,\dots ,u_m)$, a common model isThe diffusion coefficients $D_i$ describe transport, while the reaction rates $f_i$ describe local conversion, production or loss. Writing $D=\mathrm{diag}(D_1,\dots ,D_m)$ gives the vector equation $u_t=D\Delta u+f(u)$. Chemical concentrations and rates must be modeled consistently; a reaction term can exchange material with an external supply.

A spatially homogeneous equilibrium $u_{\ast }$ satisfies $f(u_{\ast })=0$. Linearizing a spatial mode with $-\Delta \varphi =q\varphi$ givesThus reaction stability tests $J$, while spatial stability tests all boundary-compatible values of $q$. Unequal diffusion can create a Turing instability even when homogeneous perturbations decay.^[ref] The Brusselator and quadratic activator-inhibitor model give explicit two-species examples. Initial and boundary conditions are part of the system; the kinetic equations alone do not select its spatial modes.

Reptation is the snake-like motion of a long polymer through the constraints imposed by neighboring chains or a fixed porous matrix. A tube model restricts large transverse motion while allowing curvilinear motion along the chain contour. For a flexible chain with additive segment friction, contour diffusion scales as $N^{-1}$ and the time to renew a contour of length proportional to $N$ scales as $N^3$. With coil size squared proportional to $N$, this gives a centre-of-mass diffusion coefficient proportional to $N^{-2}$. Applying this tube picture to weak-field DNA gel electrophoresis gives an inverse-length mobility over its appropriate regime; field-induced orientation and trapping can invalidate that simple law.

Restricted diffusion is diffusion whose displacement statistics are constrained by boundaries that particles cannot cross. A reflecting boundary condition models an impermeable wall: particles continue moving, but there is no normal diffusive flux through it. Finite permeability would require a different boundary model.

A physical demonstration uses a liquid confined between solid lamellae. In his first-person account, John E. Tanner described constructing 99 mica layers separated by aluminium-foil strips, removing wrinkles and burrs to obtain regular spacing; Edward O. Stejskal inspected the stack with a stereomicroscope. Pulsed magnetic-field gradients encode displacements in the nuclear magnetic resonance signal, allowing motion normal to the barriers to be compared with less confined directions.^[ref] This is a material confinement experiment, not a change in the liquid's chemistry.

For an original one-dimensional model, let particles diffuse in $0\le x\le L$ with constant diffusion coefficient $D>0$:The wall conditions are homogeneous Neumann boundary conditions for the diffusion equation. Start in equilibrium, so $X(0)$ has the uniform distribution on $[0,L]$. Its variance isAt long time separation, the positions decorrelate. Consequently the mean squared displacement approachesAt short times relative to $L^2/D$, most displacements have not encountered a wall, and the leading scale is the free-diffusion result $2D\tau$. For a hypothetical gap $L=10\mu \mathrm{m}$, the long-time plateau is $100/6\mu {\mathrm{m}}^2$ and its root is about $4.08\mu \mathrm{m}$.

The initial condition matters: fixing $X(0)=L/2$ instead gives the long-time value $L^2/12$. Likewise, echo attenuation in a confined sample must not automatically be converted into a free-diffusion coefficient. Confinement geometry, timing and the selected boundary model all affect the inference.

R. M. Cotts is a credited coauthor of Pulsed field gradient stimulated echo methods for improved NMR diffusion measurements in heterogeneous systems by R. M. Cotts and others, published in 1989. The joint work develops stimulated echo sequences for measuring the translational diffusion of resonant atomic nuclei in materials whose internal magnetic environment is not uniform.^[ref]

Pulsed field gradient NMR relates the attenuation of a stimulated echo to motion. In heterogeneous materials, components with different magnetic susceptibilities can produce background magnetic field gradients, which add attenuation and can bias the inferred diffusion coefficient. The paper describes three modified sequences that reduce these effects; under the stated conditions, bipolar gradient pulses cancel a problematic cross term. Two sequences were tested using proton signals in both a homogeneous liquid and a heterogeneous granular solid.^[ref]

The experimental comparison makes a useful measurement principle concrete: an apparatus can respond to both molecular motion and an unwanted background field variation. A sequence designed to separate these contributions improves the interpretation of a diffusion coefficient measurement.

Rotational diffusion is the random change of a particle's rigid body orientation through Brownian motion. A suspended ellipsoid can turn while its centre hardly moves; translational diffusion instead spreads its position. Random molecular collisions supply fluctuating torques as well as forces. The rotational diffusion coefficient $D_r$ measures how quickly angular uncertainty grows, rather than a steady spin speed. Its units are ${\mathrm{s}}^{-1}$, often written $\mathrm{ra}{\mathrm{d}}^2{\mathrm{s}}^{-1}$ to distinguish angular spreading from translational spreading, whose diffusion coefficient $D_t$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. The radian is dimensionless, but writing angular variance in $\mathrm{ra}{\mathrm{d}}^2$ remains useful.

First identify what the image can reveal. A marked end of an elongated particle defines an arrow: turning it through $\pi$ reverses that arrow. An unmarked ellipsoid with identical ends reveals only its long axis, a director: angles $\varphi$ and $\varphi +\pi$ describe the same visible state. An unmarked perfect sphere has no visible rigid body orientation at all. Tracking its centre cannot measure its rotational diffusion; a surface mark or another identifiable internal direction is needed.

This distinction makes rotational diffusion useful: it measures how long a particle remembers its pointing direction. For an elongated particle, turning also changes which laboratory direction is easiest for translational diffusion. Randomization of rigid body orientation can therefore erase directional preferences in the observed displacement, even though the particle's shape has not changed.

For a sphere of radius $a$ in an otherwise stationary, unbounded Newtonian fluid with dynamic viscosity $\eta$, impose the no-slip boundary condition and the slow-flow assumptions of Stokes flow. The resisting torque on a rotating sphere in Stokes flow is $-{\zeta }_r\mathbf{\Omega }$, where $\mathbf{\Omega }$ is its angular velocity and ${\zeta }_r=8\pi \eta a^3$. Thermal rotational diffusion at absolute temperature $T$ then givesHere $k_B$ is the Boltzmann constant; the translational expression is the Stokes–Einstein relation. The ratio has units of inverse length squared, so comparing the numerical values of $D_r$ and $D_t$ without their units is meaningless. At fixed temperature and dynamic viscosity, decreasing radius tenfold increases $D_r$ a thousandfold but $D_t$ only tenfold. This strong size dependence explains why small particles can lose directional memory extremely quickly.

Choose $T=298\mathrm{K}$ and $\eta =0.001\mathrm{Pa}\mathrm{s}$ as illustrative model inputs. Using $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$ gives $k_{B}T=4.11433402\times 10^{-21}\mathrm{J}$ and the following predictions:These are bulk three-dimensional sphere calculations, explicitly not a fit to the confined ellipsoid experiment below. Walls, nonspherical shape and departure from the no-slip boundary condition can change rotational resistance. A visible surface mark is still necessary to read the sphere's turning in an image, even when its rotational diffusion coefficient is well defined physically.

An overdamped bead-spring model of polymer dynamics with local friction and no long-range hydrodynamic interaction. For an ideal unentangled chain of $N$ identical frictional units, the centre-of-mass diffusion coefficient is proportional to $N^{-1}$ and the longest internal relaxation time to $N^2$. A tagged segment has an intermediate-time mean-square displacement proportional to $t^{1/2}$ before ordinary whole-chain diffusion dominates. Entanglement adds constraints beyond this model.

When salt-free ice forms from a saline liquid, the moving phase boundary excludes the dissolved salt. With equal phase mass densities, no bulk flow, liquid salt diffusion coefficient $D$, and interface concentration $C_i$, salt conservation gives the displayed condition, where $x$ increases from solid to liquid. Positive freezing speed $\dot{a}$ therefore requires a negative liquid concentration gradient. During melting the same condition gives dilution. Coupling salt rejection, the liquidus, and the Stefan condition produces the saline Stefan problem.

For salt-free ice and a liquid concentration profile proportional to $\mathrm{erfc}[x/(2\sqrt{Dt})]$, salt rejection at $a=2\lambda \sqrt{Dt}$ gives $C_{i}F(\lambda )=C_0$, where $D$ is the salt diffusion coefficient, and $C_i,C_0$ are the interface and far-field concentrations. This follows from $-D{C}_x(a^{+})=\dot{a}{C}_i$. The function is positive and strictly decreasing for all real $\lambda$, since an integration by parts givesConsequently $C_i=C_0/F(\lambda )$ increases with $\lambda$, equals $C_0$ at zero, and is enriched for freezing and diluted for melting. The positivity proof also prevents spurious negative interface concentrations when using the complementary error function formula.

SEC light-scattering detection adds information from scattering of light to the separation. Static scattering with a concentration readout can infer a weight-average molar mass for material in an eluting slice without the standard-family retention calibration. Dynamic scattering estimates a diffusion coefficient and hence a hydrodynamic radius. “Absolute” here means independent of a retention-to-molar mass standard curve; it still needs calibrated detectors, optical response parameters, and a valid sample model. Coeluting species produce a mixture-weighted result, rather than individually identified molecules.

A concrete manufacturer test used a Bio SEC-3 $300\mathring{\mathrm{A}}$ column, $7.8\times 300\mathrm{mm}$ with $3\mu \mathrm{m}$ particles, an Agilent 1260 Infinity G5611A quaternary pump and G5667A autosampler, a $10\mathrm{mm}$ diode-array flow cell, and dual-angle/dynamic scattering detectors. A buffer solution at pH $7.4$ flowed at $0.75\mathrm{mL}\mathrm{mi}{\mathrm{n}}^{-1}$; column and scattering detector were at $30{}^{\circ }\mathrm{C}$. Absorbance at $280\mathrm{nm}$ supplied a concentration-sensitive readout. Dilution and loading comparisons tested response: seventeen ovalbumin injections spanning $2$–$960\mu \mathrm{g}$ gave a reported average $47,198\mathrm{g}\mathrm{mo}{\mathrm{l}}^{-1}$ with $2\%$ relative standard deviation. The note reports poorer separation of thyroglobulin from aggregates and an elevated mass estimate, illustrating the coelution limit. Its $0.2\mu \mathrm{m}$ sample filtration also limits interpretation of larger protein aggregation products removed before injection.^[ref]

Sink-boundary approximation error in drug release is the difference between a release prediction using an ideal receiving boundary and one that accounts for the actual boundary. The reference prediction and the error denominator must be specified. Finite external resistance and nonzero bulk drug concentration are separate departures from a perfect sink condition in drug release.

For the maintained-reservoir steady planar model, let the internal concentration be $c_i>0$, receiving bulk concentration $c_b$, matrix thickness $L$, internal diffusion coefficient $D$ and external transfer coefficient $h$. With compatible concentration scales,The ideal zero-bulk, negligible-external-resistance prediction is $J_0=D{c}_i/L$. Defining $\beta =c_b/c_i$ and the mass-transfer Biot number for drug release $B=hL/D$ givesThis is the combined boundary-resistance and bulk-concentration error in drug release for the stipulated model. It is not an exact transient error formula for every dosage form.

An original example with $B=20$ and $\beta =0.10$ gives $J/J_0=6/7$. The ideal flux overpredicts by about $14.3\%$ relative to $J_0$, or $16.7\%$ relative to $J$. Making external transport arbitrarily fast still leaves the concentration-related error floor $\beta$ in the first error convention.

The model also supplies a transparent tolerance rule. For $0\le \beta <\epsilon <1$, requiring $(J_0-J)/J_0\le \epsilon$ givesAt $\epsilon =0.05$ and $\beta =0.01$, the threshold is $23.75$. If bulk accumulation already exceeds the tolerance, increasing $h$ alone cannot meet it.

Experimentally, changing receiving volume tests accumulation, while changing flow or membrane transport tests external resistance; sampling correction in drug release testing is needed when material is removed. The FDA’s primary adaptive-perfusion account specifically identifies passive membrane transport as potentially rate-limiting.^[ref] Separating the two effects prevents a large receiving volume from being mistaken for proof of an ideal surface boundary.

For an unbiased one-dimensional walk with step size $a$, diffusion coefficient $D$ fixes step interval $a^2/(2D)$. A small dissociation probability $P$ per step gives about $1/P$ steps before release. Thus lifetime depends on $D$, while the root-mean-square net displacement before release is $a/\sqrt{P}$. Net displacement is not the cumulative distance along all backtracking steps. Dimensional consistency of the measured or assumed diffusivity is necessary before applying the estimate.

For overdamped Brownian motion at uniform temperature in a potential energy $U$, constant diffusion coefficient $D$ gives the displayed Fokker-Planck equation, with $\beta =(k_{B}T{)}^{-1}$. The probability current is $\mathbf{J}=-D(\nabla p+\beta p\nabla U)$. Zero current recovers the Boltzmann distribution; absorbing boundaries instead permit transport and mean first-passage times. This diffusion equation should be distinguished from the similarly named coagulation equation.