Neural Engineering

Derivation 1: Channelrhodopsin Photocurrent Model

The ChR2 photocurrent is modeled by a 3-state kinetic scheme: closed (C), open (O), and desensitized (D). The transition rates depend on light irradiance $\Phi$:

$$C \xrightarrow{\epsilon \Phi \sigma} O \xrightarrow{G_d} D \xrightarrow{G_r} C$$

The fraction of channels in the open state evolves as:

$$\frac{dO}{dt} = \epsilon \Phi \sigma \cdot C - G_d \cdot O$$

$$\frac{dD}{dt} = G_d \cdot O - G_r \cdot D$$

where $\epsilon$ is quantum efficiency (~0.5), $\sigma$ is cross-section (~1.2 $\times$ 10$^{-20}$ m$^2$), $G_d \approx 0.1$ ms$^{-1}$ (open-to-desensitized), and $G_r \approx 0.01$ ms$^{-1}$ (recovery). The photocurrent is:

$$I_{\text{ChR2}} = g_{\max} \cdot O \cdot (V - E_{\text{rev}})$$

where $g_{\max}$ is the maximal conductance and $E_{\text{rev}} \approx 0$ mV. The peak photocurrent occurs at light onset, followed by adaptation to a steady state as channels enter the desensitized state. This model predicts the characteristic transient-to-sustained response ratio and the frequency-dependent fidelity of optogenetic stimulation.

Derivation 2: Volume of Tissue Activated (VTA)

The electric field around a DBS electrode in a homogeneous medium with conductivity$\sigma$ is:

$$\mathbf{E}(\mathbf{r}) = -\nabla V = \frac{I}{4\pi\sigma |\mathbf{r} - \mathbf{r}_0|^2} \hat{\mathbf{r}}$$

Neural activation occurs where the second spatial derivative of the potential (activating function) exceeds the threshold:

$$f_n = \frac{\partial^2 V_e}{\partial x^2} > \frac{V_{\text{th}}}{R_a \Delta x}$$

The VTA radius for a point source is approximately:

$$r_{\text{VTA}} \approx \sqrt{\frac{I}{4\pi\sigma \cdot f_{\text{th}}}}$$

For typical DBS parameters (I = 3 mA, $\sigma$ = 0.2 S/m), the VTA radius is ~2–3 mm. Patient-specific finite element models incorporating brain anatomy and diffusion-tensor-derived conductivity tensors predict the actual VTA shape, which guides programming of directional DBS leads with segmented electrodes.

Derivation 3: Informational Lesion Hypothesis

High-frequency DBS (~130 Hz) is therapeutic despite opposing effects (excitation vs. inhibition) being debated. The informational lesion hypothesis proposes that DBS replaces pathological firing patterns with regular, high-frequency activity that carries less pathological information. The entropy of the stimulus-driven firing pattern is:

$$H_{\text{DBS}} = -\sum_i p_i \log_2 p_i \approx 0$$

since the ISI distribution becomes concentrated at $1/f_{\text{DBS}}$. The pathological information (e.g., beta-band oscillations in PD) is quantified by the mutual information between STN activity and downstream targets:

$$I_{\text{path}} = H_{\text{GPi}} - H_{\text{GPi}|\text{STN}} > 0$$

DBS reduces $I_{\text{path}}$ to near zero by replacing the pathological STN signal with a regular pattern that downstream nuclei cannot use to propagate abnormal oscillations. This explains why DBS frequency must exceed the pathological frequency (beta ~20 Hz) and why lower frequencies can worsen symptoms.

Derivation 4: Cochlear Implant Channel Interaction

A cochlear implant has $N_e$ electrodes (typically 12–22) that stimulate different tonotopic locations along the cochlea. The spread of excitation from electrode$j$ to auditory nerve fiber $i$ follows an exponential decay:

$$w_{ij} = \exp\left(-\frac{|x_i - x_j|}{\lambda_e}\right)$$

where $\lambda_e$ is the spread constant (~3–8 mm). The effective spectral resolution is limited by channel interaction. The number of independent information channels is:

$$N_{\text{eff}} = \frac{N_e}{1 + (N_e - 1)\bar{\rho}} \leq N_e$$

where $\bar{\rho}$ is the mean inter-electrode correlation. Typical cochlear implants achieve $N_{\text{eff}} \approx 4\text{--}8$ independent channels despite having 12–22 physical electrodes. This limits spectral resolution and explains why music perception remains challenging for CI users. Current-focusing strategies (tripolar stimulation, current steering) reduce $\lambda_e$ and increase $N_{\text{eff}}$.

Derivation 5: Organoid Growth and Nutrient Diffusion Limit

The maximum size of an organoid is limited by oxygen diffusion. The steady-state oxygen concentration in a spherical organoid of radius $R$ follows:

$$D_O \nabla^2 c = q_0, \quad c(R) = c_0$$

where $D_O$ is the oxygen diffusion coefficient (~2 $\times$ 10$^{-5}$ cm$^2$/s),$q_0$ is the metabolic consumption rate, and $c_0$ is the surface concentration. The solution in spherical coordinates is:

$$c(r) = c_0 - \frac{q_0}{6 D_O}(R^2 - r^2)$$

The center becomes anoxic ($c(0) = 0$) when the radius exceeds:

$$R_{\max} = \sqrt{\frac{6 D_O c_0}{q_0}} \approx 0.5\text{--}1 \text{ mm}$$

This limit constrains organoid size without vascularization. Current approaches include slicing organoids for culture, using bioreactors with enhanced oxygen delivery, and transplanting organoids into vascularized host tissue. Vascularized organoids grown with endothelial cells can exceed this size limit and develop more mature neural circuits.