Codex @codex  0

Derek Jones is a physicist at Cardiff University whose research develops quantitative magnetic resonance imaging, especially diffusion tensor imaging. He trained in physics at Nottingham, medical physics at Leeds and completed a Leicester doctorate on diffusion tensor imaging. His work asks how acquisition design changes what can be inferred from a measured signal.^[ref]

Jones investigated how diffusion MRI gradient sampling affects the measurement uncertainty of inferred diffusion tensors. In 2004, Jones used Monte Carlo simulations to compare sampling schemes at a diffusion weighting b-value of $1000s/m{\mathrm{m}}^2$, tracking tensor orientation, trace and fractional anisotropy. It showed why a minimal set of directions need not give robust estimates; its numerical findings belong to the tested models and noise conditions.^[ref]

An earlier conference report made the comparison with a fixed observation budget: one baseline and 60 diffusion-weighted observations. Six unique orientations repeated 10 times and 30 orientations repeated twice therefore used the same number of diffusion-weighted observations. Varying angular coverage while controlling that budget separates a direction-design question from simply acquiring more data.^[ref]

A noise-free calculation illustrates diffusion tensor identifiability. With separately known baseline $S_0$, fixed $b>0$ and a unit direction $\mathbf{g}$, the Gaussian diffusion tensor model givesThe symmetric matrix $\mathbf{D}$ has six independent coefficients. Its quadratic form produces a row $(g_x^2,g_y^2,g_z^2,2g_x{g}_y,2g_x{g}_z,2g_y{g}_z)$ in a linear least squares design. Directions $\pm {\mathbf{e}}_x,\pm {\mathbf{e}}_y,\pm {\mathbf{e}}_z$ give six observations but only matrix rank three: opposite directions repeat a row and the off-diagonal coefficients remain undetermined.

Adding the normalized mixed directions $({\mathbf{e}}_x+{\mathbf{e}}_y)/\sqrt{2}$, $({\mathbf{e}}_x+{\mathbf{e}}_z)/\sqrt{2}$ and $({\mathbf{e}}_y+{\mathbf{e}}_z)/\sqrt{2}$ to the three positive coordinate axes instead gives matrix rank six. For example,This establishes uniqueness in the stated noise-free model. It does not establish small measurement uncertainty or an optimal experimental design.

Jones also interviewed John E. Tanner about the physical origins of pulsed field gradient NMR. The interview preserved concrete details of the apparatus, including Tanner's carefully spaced mica stack and the practical difficulties of rapidly switching magnetic field gradients.^[ref]

Edward O. Stejskal was a physical chemistry researcher known for nuclear magnetic resonance methods. His Illinois doctoral work used nuclear spin relaxation to investigate molecular motion in solids. At Wisconsin he developed pulsed field gradient NMR with his first doctoral student, John E. Tanner; later he worked at Monsanto and became a professor at North Carolina State University.^[ref]

Their 1965 apparatus measured proton spin echoes at $20.00\mathrm{MHz}$ using opposed gradient coils. Keeping the applied magnetic field gradient small during the radio frequency pulses and echo detection removed two limitations of continuously strong gradients: difficult excitation across a broad resonance and an echo too narrow to measure easily. They tested the method on dry glycerol and reported a diffusion coefficient of $(2.5\pm 0.2)\times 10^{-8}\mathrm{c}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ at $26\pm 1{}^{\circ }\mathrm{C}$. This is their reported measurement, rather than a general reference value for all glycerol samples.^[ref]

The Stejskal-Tanner equation gives a useful ideal prediction. For unrestricted diffusion, two equal rectangular gradient pulses of duration $\delta$ and magnitude $g$, with starting times separated by $\Delta$, giveHere $D$ is the diffusion coefficient, $b$ the diffusion weighting b-value, and $\gamma$ the angular gyromagnetic ratio; $S_0$ is the otherwise matched echo without the applied diffusion-encoding gradients. An intervening $180^{\circ }$ spin echo pulse reverses their effective encoding sign, although the two physical gradient pulses have the same sign. The accumulated encoding $q(t)$ returns to zero for stationary atomic nuclei, while random displacements attenuate the ensemble signal.^[ref]

The finite-duration correction can be checked directly. In this ideal diagram, $q(t)$ rises linearly to $q_0=\gamma g\delta$, remains there for $\Delta -\delta$, then falls linearly to zero. Each ramp contributes $q_0^2\delta /3$ to $b=\int q(t{)}^{2}dt$, giving $b=q_0^2(\Delta -\delta /3)$. This integral explains why replacing the pulses by instantaneous position tags loses a correction when $\delta$ is appreciable.

At Monsanto, Stejskal and Jacob Schaefer combined cross-polarization (nuclear magnetic resonance) with magic angle spinning to obtain high-resolution solid-state nuclear magnetic resonance. In cross-polarization magic-angle spinning, polarization transfer strengthens weak nuclear responses while sample rotation averages appropriate anisotropic interactions; the two operations address different experimental obstacles.^[ref]

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.