Overview
Neural interfaces represent an interdisciplinary convergence of microelectronics, clinical neurophysiology, materials science, and computational signal processing. To establish communication with, decode intention from, or deliver therapeutic electrical stimulation to the central and peripheral nervous systems, biomedical engineers have developed an array of distinct structural architectures. The foundational engineering trade-off across all brain-computer interfaces (BCIs) and bioelectronic medicines is the balance between surgical invasiveness and physical signal resolution. This sandbox simulator acts as an interactive biophysical testbed modeling volume conduction, biological low-pass filtering, thermal noise generation, and spatial resolution across leading clinical neurotechnology paradigms.
At the microscopic level, neural signaling is driven by ionic flux across neuronal cell membranes. During an action potential, transient opening of voltage-gated sodium ($\text{Na}^+$) and potassium ($\text{K}^+$) ion channels creates transmembrane current sinks and sources. Under the quasistatic approximation of Maxwell's equations ($\omega \epsilon \ll \sigma$), brain tissue behaves primarily as a purely resistive, volume-conducting medium. For a point current source $I$ in an isotropic, homogeneous medium with electrical conductivity $\sigma$, the extracellular potential $\Phi$ at radial distance $r$ is governed by:
$$\Phi(r) = \frac{I}{4 \pi \sigma r}$$
As micro-currents travel from deep cortical layers through successive anatomical boundaries—including gray matter ($\sigma_{\text{cortex}} \approx 0.33\ \text{S/m}$), cerebrospinal fluid ($\sigma_{\text{CSF}} \approx 1.79\ \text{S/m}$), the dura mater ($\sigma_{\text{dura}} \approx 0.18\ \text{S/m}$), cranial bone ($\sigma_{\text{bone}} \approx 0.006\ \text{S/m}$), and scalp skin ($\sigma_{\text{skin}} \approx 0.43\ \text{S/m}$)—the skull acts as a spatial low-pass filter with extreme attenuation. This high-resistivity boundary smears high-frequency, localized extracellular spikes ($\sim 300\text{--}3000\ \text{Hz}$) into aggregated, low-frequency local field potentials (LFPs) and electroencephalographic rhythms ($\sim 0.5\text{--}100\ \text{Hz}$). This simulation accurately portrays why intracortical microelectrodes can isolate single-unit action potentials, whereas endovascular, sub-scalp, and subcutaneous arrays capture spatially synchronized population dynamics.
How to Use
The workspace is divided into two synchronized visual viewports: the Anatomical Cross-Section (left pane), depicting tissue layers and physical electrode locations, and the real-time Oscilloscope & Spectral Density Analyzer (right pane).
Use the following controls to investigate the biophysical realities of each neural interface paradigm:
- Implant Paradigm Selection: Select any of the six clinical neurotechnology architectures from the dropdown menu to inspect its anatomical placement and signal output:
- Neuralink N1 (Intracortical): Polyimide microthreads inserted directly into cortical layers II–V ($d \approx 1.5\ \text{mm}$), capturing high-bandwidth extracellular action potentials (spikes) with minimal volume attenuation.
- Synchron Stentrode™ (Endovascular): An expandable nitinol stent deployed via the internal jugular vein into the superior sagittal sinus vessel lumen adjacent to the motor cortex, acquiring transvascular local field potentials without open craniotomy.
- Epiminder Minder iCEM® (Sub-scalp): A sub-pericranial ribbon electrode array placed beneath the scalp musculature outside the cranial bone, providing stable long-term regional EEG monitoring for ambulatory seizure forecasting.
- UNEEG SubQ (Subcutaneous): A discrete subcutaneous wire electrode implanted in the temporal region to track chronic regional EEG trends over months to years.
- NeuroPace RNS (Depth Closed-Loop): Deep parenchymal or cortical strip leads configured for real-time responsive neurostimulation; detects epileptiform spike-wave envelopes and triggers high-frequency burst suppression stimulation.
- Cochlear Implant (Auditory Spiral): A flexible silicone electrode array threaded into the scala tympani of the cochlea, utilizing tonotopic spatial mapping to stimulate auditory nerve fibers via frequency-modulated charge pulses.
- Depth and Impedance Diagnostics: Adjust the Surgical Depth and Tissue Impedance ($Z$) sliders. Increasing tissue impedance models glial encapsulation, fibrotic sheath formation, or skull resistance, directly increasing baseline Johnson-Nyquist thermal noise.
- Channel Scaling & Telemetry: Vary the Active Channels slider to observe changes in total uncompressed transmission bandwidth ($B = N_{\text{chan}} \cdot f_s \cdot \text{bit depth}$), illustrating the thermal and power constraints of wireless implants.
- Artifact Infiltration: Enable the Simulate Cranial/Ocular Artifact toggle to observe how high-amplitude physiological interference (e.g., ocular blinks and temporalis muscle contractions) severely corrupts sub-scalp and subcutaneous traces, while deep cortical arrays remain protected.
- Epileptiform Discharges & Closed-Loop Burst Stimulation: Check Simulate Epileptic Discharge / RNS to transition from resting baseline rhythms into high-voltage $3\ \text{Hz}$ spike-wave paroxysms. When active in NeuroPace mode, observe the closed-loop detector deliver high-frequency biphasic stimulation pulses that suppress the seizure envelope.
- Real-Time Audio Laboratory: Click Sound ON/OFF to activate the Web Audio API synthesis engine. Listen to discrete action potential crackles (multi-unit hash) for intracortical interfaces, low-frequency field oscillation hums, or pitch-modulated tonotopic carriers for the cochlear implant.
- Fullscreen Clinical Monitor & Slide-out HUD: Click 🖥️ FULLSCREEN to enter an immersive display mode. In fullscreen, click ⚙️ CONTROLS to open the animated slide-out HUD drawer for live adjustments.
Technical Details
The digital signal synthesizer executes a multi-stage physical simulation pipeline inside a synchronized 60 FPS requestAnimationFrame loop. Action potentials are synthesized using an asymmetric biphasic voltage function modeling fast $\text{Na}^+$ inward currents followed by delayed $\text{K}^+$ repolarization and hyperpolarization:
$$V_{\text{spike}}(t) = A \cdot \left(1 - \frac{t - t_0}{\tau_1}\right) \cdot \exp\left(-\frac{t - t_0}{\tau_2}\right)$$
where $A$ represents the peak extracellular potential ($\approx 40\text{--}120\ \mu\text{V}$ for intimate cortical contact), $t_0$ is the spike onset, $\tau_1 = 0.4\ \text{ms}$ governs depolarization, and $\tau_2 = 1.1\ \text{ms}$ governs the repolarization decay.
To simulate volume conduction through layered biological tissue, the signal undergoes frequency-dependent attenuation modeled by spatial transfer functions. The attenuation factor $\alpha(f, d)$ is calculated as:
$$\alpha(f, d) = \exp\left(-d \cdot \sqrt{\pi f \mu \sigma}\right) \cdot \left(\frac{\sigma_{\text{eff}}}{\sigma_{\text{eff}} + 2\pi f C_m}\right)$$
where $d$ is placement depth in millimeters, $f$ is frequency in Hertz, $\sigma_{\text{eff}}$ is the effective tissue layer conductivity, and $C_m$ accounts for cellular membrane capacitive shunting.
Thermal background noise is generated using the Johnson-Nyquist noise formulation across the electrode-electrolyte interface and biological tissue impedance $Z$:
$$V_{\text{RMS}} = \sqrt{4 k_B T \cdot Z \cdot \Delta f}$$
where $k_B = 1.380649 \times 10^{-23}\ \text{J/K}$, $T = 310.15\ \text{K}$ ($37^\circ\text{C}$ physiological body temperature), and $\Delta f$ is the recording analog bandwidth. The Signal-to-Noise Ratio (SNR) in decibels is calculated dynamically:
$$\text{SNR}_{\text{dB}} = 10 \log_{10}\left(\frac{P_{\text{signal}}}{P_{\text{noise}}}\right) = 20 \log_{10}\left(\frac{V_{\text{signal, RMS}}}{V_{\text{noise, RMS}}}\right)$$
The frequency domain analyzer computes a real-time 32-bin spectral power distribution across canonical neurological frequency bands: Delta ($\delta$: $0.5\text{--}4\ \text{Hz}$), Theta ($\theta$: $4\text{--}8\ \text{Hz}$), Alpha ($\alpha$: $8\text{--}13\ \text{Hz}$), Beta ($\beta$: $13\text{--}30\ \text{Hz}$), Gamma ($\gamma$: $30\text{--}80\ \text{Hz}$), and High-Frequency Multi-Unit Spiking ($300\text{--}3000\ \text{Hz}$). KaTeX mathematical formatting ensures clean, precise symbolic rendering across all devices.
Future Directions
The next frontier in neural interface engineering focuses on addressing chronic biological biocompatibility and micro-motion injury. Current silicon Utah arrays and rigid microwires experience mechanical compliance mismatch with brain tissue ($E_{\text{silicon}} \approx 170\ \text{GPa}$ versus $E_{\text{brain}} \approx 1\text{--}3\ \text{kPa}$). Over months in vivo, shear stresses trigger continuous foreign body reactions, activating microglia and reactive astrocytes to form a dense glial scar sheath. This scar increases tissue impedance $Z(t)$ and displaces viable neuronal somas outside the recording radius ($r > 100\ \mu\text{m}$), degrading high-frequency unit isolation.
Future releases of this simulation sandbox will incorporate dynamic glial encapsulation kinetics modeled over longitudinal timeframes:
$$Z(t) = Z_0 + \Delta Z_{\text{scar}} \cdot \left(1 - e^{-t / \tau_{\text{encap}}}\right)$$
Additionally, emerging paradigms—including ultra-flexible conductive polymer threads (PEDOT:PSS), carbon nanotube yarns, injectable mesh electronics, ultrasonic neural dust, and optogenetic wave-guided stimulation arrays—will be integrated to allow comparative evaluation of chronic signal stability and longevity profiles.
Explore Related Laboratory Simulators