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]

John E. Tanner was a physical chemistry researcher whose work with Edward O. Stejskal established pulsed field gradient NMR as a method for measuring self-diffusion. He studied physical chemistry at Oberlin, Indiana and Wisconsin, and later worked at Idaho National Laboratory. His early NMR work received renewed recognition with an honorary fellowship in 2016.^[ref]

The appeal of pulsed field gradient NMR was practical as well as theoretical: strong magnetic field gradients could encode motion without remaining strong during excitation or detection. In a 2017 interview with Derek Jones, Tanner recalled choosing Stejskal's laboratory because it had a magnet and because he wanted to measure slow diffusion. He also described rebuilding the probe to deal with eddy currents generated by switching the magnetic field gradients.^[ref]

Tanner made restricted diffusion experimentally concrete. He described assembling 99 layers of mica with aluminium foil spacers, carefully removing wrinkles and burrs to obtain regular gaps. Such a controlled geometry lets an echo probe displacement between known barriers, rather than treating every loss of NMR signal as free diffusion with a smaller diffusion coefficient.^[ref]

His 1970 work used a stimulated echo to extend the observation interval when spin-lattice relaxation is slower than transverse relaxation. The interval can store the relevant pattern as longitudinal magnetization, rather than keeping it transverse throughout.^[ref] For an illustrative ideal comparison, take $T_1=1\mathrm{s}$, $T_2=0.1\mathrm{s}$ and storage time $T=0.5\mathrm{s}$. The storage-dependent nuclear spin relaxation factors are $e^{-T/T_1}\simeq 0.607$ and $e^{-T/T_2}\simeq 0.00674$. Even allowing the ideal three-pulse pathway's initial factor of $1/2$, the former factor is about 45 times larger. This compares only that storage contribution: pulse losses, the other transverse intervals and diffusion attenuation still matter.^[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.