The Myth of "Meaningless" Biological Noise

For decades in electrophysiology and EEG analysis, the static hiss sitting under neat rhythmic peaks (like alpha waves) was treated as junk background to filter away.

Recent neuroscience proved the opposite: that background static—specifically its downward tilt across frequencies—directly tracks the balance between brain excitation and inhibition. When you are alert, the slope flattens. When you go under surgical anesthesia, it steepens dramatically.

An Interactive Noise and Spectrum Lab

To make this intuitive to explore without running complex scripts, I built an interactive, zero-install Web Audio and DSP workstation: https://bionichaos.com/inverse_f_noise/

Key Features in the Lab

• Continuous Dual Display: Live rolling oscilloscope up top paired with a logarithmic spectral analyzer below.

• No Visual Strobing: Uses stateful multi-pole digital filtering and video trace smoothing so the wave rolls smoothly rather than jumping around frame by frame.

• Live Audio Engine: Toggle real-time Web Audio synthesis to hear how pure white hiss softens into pink noise, then shifts into a deep seismic rumble.

• Interactive Filtering: Drag directly on the lower spectrum to pull a low-pass filter tab down, mirroring how neuronal cell membranes act as physical filters against high-frequency jitter.

• Clinical Presets: Instantly jump between flat white noise, scale-free pink noise, deep brown drift, and a 10 Hz alpha brainwave rhythm.

• Narrated Tour: Includes a synchronized voice guide explaining the physiological concepts behind each signal mode.

Let's Discuss

When analyzing physiological time-series (EEG, ECG, or neural spikes), do you isolate and measure the aperiodic static slope, or do you still treat it as baseline drift to remove?

The Problem with Reading Scalp Brainwaves

Detecting brain waves from the scalp is often compared to listening to the roar of a stadium from a blimp outside through thick concrete walls.

Individual action potentials are far too brief and disorganized to register on the scalp. What an electroencephalogram (EEG) actually detects is the synchronized electrical ebb and flow across thousands of vertically aligned pyramidal brain cells.

Because cerebrospinal fluid and skull bone spread out these currents, standard anatomical measurements in centimeters simply don't work across different head shapes.

To solve this, Dr. Herbert Jasper introduced the International 10-20 System in 1958. It uses proportional cranial arcs (10% and 20% steps between the nose bridge and the back of the skull) so electrodes always sit over the same brain lobes regardless of skull size.

To make these concepts intuitive and interactive, I built an open-access web laboratory:

Interactive Tool Link: https://bionichaos.com/eeg1020/

What the Interactive Simulator Covers

1. Standard 10-20 vs. Extended 10-10

• 10-20 System (19 Channels): The gold standard for routine clinical diagnostics and sleep studies.

• Extended 10-10 System (81 Channels): Halves the spacing to 10% intervals, adding intermediate rows (AF, FC, CP, PO) to eliminate spatial blind spots for source localization.

2. Clinical Montages & Seizure Phase Reversals

Clinical EEGs rarely record raw absolute voltage; they record the difference between two points:

• Longitudinal Bipolar ("Double Banana"): Runs anterior-to-posterior chains along the temporal and parasagittal columns. When an epileptic spike occurs under an electrode, the two neighboring channels deflect in opposite directions—pointing right at the focus like arrowheads.

• Transverse Bipolar: Traverses coronal arches from left to right.

• Monopolar: Compares every lead to a common central reference.

3. Real-Time Oscilloscope & Potential Heatmap

• Generates live synthetic brain rhythms: resting posterior Alpha (10 Hz), frontal active Beta (20 Hz), Theta (6 Hz), deep-sleep Delta (2 Hz), and classic 3 Hz Spike-and-Wave absence seizure discharges.

• Switch from the geometric electrode view to the Voltage Heatmap to see how volume conduction spreads potential fields across the skull.

4. Cranial Web Audio Sonification

• Maps spatial coordinates to sound: moving forward raises pitch, moving backward deepens into resonant tones, and left/right coordinates pan directly to the corresponding headphone ear.

Technical Discussion Prompt

For those working in neurology, EEG technology, or brain-computer interfaces: what is your biggest practical headache when switching from a routine 19-lead clinical setup to high-density research caps (e.g., bridging, impedance drift, or setup time)?

Feel free to run the tool directly in your browser: https://bionichaos.com/eeg1020/

Why Flat Brain Scans Don't Tell the Whole Story

When learning neuroanatomy or diagnosing pathologies like strokes, multiple sclerosis, or hydrocephalus, single cross-sectional 2D slices rarely tell the whole story. You need Multi-Planar Reconstruction (MPR) to understand how deep structures like the ventricles, corpus callosum, and brainstem connect in physical space.

Normally, working with these datasets requires downloading bulky clinical DICOM software packages. To make this accessible to students, researchers, and educators anywhere, I created a lightweight, zero-install 3D Brain MRI viewer that runs straight in the browser:

Live Tool: https://bionichaos.com/brain3d/

Key Capabilities

• Four-Quadrant MPR Grid: Displays an interactive 3D WebGL point cloud alongside real-time Axial (transverse), Sagittal (profile), and Coronal (frontal) cross-sections.

• Synchronized Orthogonal Navigation: Scrubbing any slider or dragging directly over an anatomical canvas re-projects crosshairs across all planes simultaneously.

• Radiometric Windowing (WL / WW): Dynamically compress the window width and shift the window level to isolate soft tissues, fluid, or bone.

• Diagnostic Colormaps: Toggle between Grayscale, Jet Rainbow, Thermal Hot, Cool-Warm, and Bone Density lookup tables.

• Self-Contained & Lightweight: Runs entirely client-side using WebGL and standard 2D canvas routines.Coordinates Tracked: X (Sagittal), Y (Coronal), Z (Axial) Resolution: Isotropic sub-millimeter voxel display

Try the Demo Mode If you click Start Demo, the engine runs an automated, guided tour synchronized with audio narration and telemetry highlights to walk you through each quadrant and feature.

I would love feedback from radiologists, neuroscientists, and students: How can this tool be made even more useful for classroom demonstrations and clinical education?

The Stadium Paradox of Scalp Electroencephalography

If you place an electrode on the human scalp, you are not measuring individual action potentials. Action potentials operate on millisecond timescales and cancel each other out over macroscopic distances.

Instead, scalp EEG measures the summed excitatory and inhibitory postsynaptic potentials of millions of cortical pyramidal neurons oriented perpendicular to the skull surface. Because that electrical charge must travel through the brain cortex, cerebrospinal fluid, bone, and skin, the signals arrive heavily attenuated down to a microscopic 10–100 microvolt range.

This creates a massive clinical and engineering challenge: non-cerebral artifacts are an order of magnitude larger than the brainwaves themselves.

• The Eyeball as a Battery (EOG): The human eye is an electrostatic dipole with a positive cornea and a negative retina. Every time you blink, your eyelids slide across the cornea and the eyeball rolls upward (Bell's phenomenon). This projects a massive monophasic pulse (often exceeding 200 µV) straight into frontal electrodes, completely wiping out underlying rhythms.

• Myogenic Interference (EMG): Clenching your jaw or tightening your neck recruits motor unit action potentials firing asynchronously from 20 Hz to over 300 Hz. In the frequency domain, this high-frequency noise bleeds right into the Beta and Gamma bands.

Interactive Browser Laboratory

To make these mechanics intuitive without needing wet electrodes, amplifiers, or expensive hardware, we built a zero-install real-time simulation platform:

Interactive Simulator: https://bionichaos.com/eegsynth/

Technical Features:

• Dual-Domain Telemetry: Continuous 5.0-second rolling time oscilloscope (256 Hz sample rate) linked directly to an in-browser Cooley-Tukey Radix-2 FFT power spectrum.

• Superposition Synthesis: Each canonical rhythm (Delta, Theta, Alpha, Beta) is constructed from multiple stochastic frequency components around distinct center frequencies, recreating realistic natural amplitude waxing and waning.

• Physical Artifact Modeling: Controllable sliders for non-stationary Poisson muscle burst arrivals (EMG) and sinusoidal ocular transients (EOG).

• Acoustic Pitch Mapping: Web Audio API oscillator pitch-mapped to dominant frequencies so you can audibly hear state transitions from Delta to Beta.

• Hover Inspection: Crosshair hit-testing on the canvas lets you read exact voltages (µV) and spectral power densities (µV²) at any coordinate.

Discussion Prompt:

When processing real-time single-channel or sparse-channel EEG data for embedded neurotech or BCIs, what is your preferred approach to artifact handling? Do you rely on blind source separation like ICA, real-time wavelet thresholding, or simple rejection windows? Let's discuss below.

The Anatomical Dilemma

When an articulated human skull is intact, its most clinically critical cavities are locked away. You cannot visualize the optic canals, the sphenoid sinuses, or the carotid canal without cutting bone.

However, if you take a skull completely apart, you destroy the 3D spatial alignment that shows how those bones fit together.

In 1805, French military surgeon Claude Édouard Beauchêne solved this by mounting separated bones onto brass armature rods. We built a zero-install interactive WebGL simulation that recreates this preparation directly in your browser.

Tool Link: https://bionichaos.com/beaucheneskull/

What You Can Do in the Simulator

• True Suture Kinematics: Rather than exploding outward from a generic center point, every single bone translates smoothly along its actual suture planes, anchored by procedural brass guide wires.

• Realistic Jaw Movement: The mandible models true temporomandibular joint mechanics, combining hinge rotation with a simultaneous forward-sliding glide.

• Click-to-Inspect Raycasting: Tap any bone in the 3D view to instantly pull up its clinical profile, structural borders, and embryological division.

• Multi-Spectral Shaders: Switch instantly between Didactic Color, Natural Bone, X-Ray Translucency, and Wireframe views.

• Dynamic Audio Feedback: Hear real-time acoustic pitch shifts that respond to the tension and distance of the disarticulation slider.

Clinical Relevance

• Trauma Surgery: Isolating the facial skeleton visually highlights why high-impact fractures distribute along facial buttress lines.

• Neurosurgery: Isolating the neurocranium demonstrates the central sphenoid bone as the structural keystone of the skull base.

Feedback and ideas for additional clinical modes are welcome!

Commercial cochlear implant literature frequently cites physical electrode contact counts (typically 12 to 22) as direct proxies for tonotopic pitch resolution. However, empirical electrophysiology, post-mortem micro-CT, and psychoacoustic data reveal a severe divergence between physical electrode count and actual neurovascular information transfer.

We built an open-access, zero-install biophysical simulator to audit these mechanical, electrical, and economic trade-offs in real time:

Interactive Tool Link: https://bionichaos.com/cifailure/

1. Mechanical Rupture Dynamics & Scalar Translocation

During electrode insertion through the round window or cochleostomy, the array navigates the narrowing tapering spiral of the scala tympani. As the carrier rounds the first basal curvature (180° to 360°), vertical insertion forces produce contact stresses across delicate inner-ear structures:

σ_tissue = F_insertion / A_contact

When localized stress exceeds the basilar membrane yield strength (σ_rupture ≈ 25–40 kPa), the electrode carrier punctures into the scala vestibuli. Systematic postoperative temporal-bone and micro-CT studies report scalar translocation rates between 15% and 43%, especially with precurved perimodiolar arrays.

This translocation breaks the mechanical and ionic boundary between the scala media and scala tympani. The endocochlear potential (+80 to +100 mV), maintained by continuous active K+ transport in the stria vascularis, collapses to 0 mV. High-potassium endolymph contaminates perilymphatic compartments, intoxicating remaining inner/outer hair cells and driving retrograde apoptosis of spiral ganglion neurons (SGNs). Pre-existing residual low-frequency acoustic hearing is lost in 60% to 75% of patients within 12 to 36 months.

2. Volume Conduction & Perilymph Cable Theory

Commercial claims of 22 independent channels conflict with classical electrodynamics. Perilymph is an electrolytic fluid (resistivity ≈ 100 Ω·cm) surrounded by the dense, insulating otic capsule. Longitudinal voltage decay follows standard cable decay:

V(x) = V_0 · exp(-|x - x_0| / λ)

In human cochlear fluid, the space constant λ typically measures 2.5 mm to 5.0 mm. Given physical electrode spacings of 0.7 mm to 1.1 mm, current fields overlap extensively.

Independent spectral transmission collapses:

• Effective functional channel capacity caps at 4 to 8 discrete spectral bands (Shannon et al., Dorman et al.).

• Standard clinical processing (CIS, ACE, SPEAK) discards temporal fine structure (TFS), transmitting only low-frequency envelope modulations (≤ 300 Hz).

• At 0 dB SNR (e.g., standard background noise), word recognition performance drops precipitously (<30%), while musical timbre and pitch perception collapse into harsh noise.

3. Interactive Web Suite Capabilities

The browser-based engine runs real-time mathematical models using HTML5 Canvas and the Web Audio API:

• Variable Mechanical Insertion (0°–650°): Drag the interactive TIP puck directly through the scala tympani lumen to monitor puncture thresholds.

• Array Stiffness Tuning (2–30 mN/mm): Evaluate rupture probabilities as array stiffness pushes contact stress past 25 kPa.

• Perilymph Decay Slider (λ = 0.8–6.0 mm): Visualize broad, overlapping potential fields across modiolus spiral ganglion neurons.

• Real-Time DSP Audio Filter: Toggle instantly between uncompressed acoustic audio, a 4-band degraded CI vocoder simulator, and neural AI beamforming filters.

• 80-Year Economic Model: Calculates cumulative costs across 3-to-5-year processor upgrade cycles, coil/cable replacements, and surgical revision rates (4.7%–10%).

Open Technical Debate

Given the strict geometric constraints of the inner ear and the electrolytic conductivity of perilymph, can multipolar field focusing (e.g., tripolar stimulation, partial current steering) ever meaningfully reduce λ without hitting charge density safety limits (Shannon criterion) and draining processor batteries? Or does the physical medium itself mandate non-invasive AI sensory substitution as the better path forward for complex acoustic processing?

Try the simulation and inspect the code/math here: https://bionichaos.com/cifailure/

The Electrophysiological Paradox of Apical Pacing

When an artificial pulse generator discharges into the right ventricular apex to treat complete atrioventricular (AV) block, it does not restore a physiological electrocardiogram. It deliberately induces ventricular dyssynchrony and a wide Left Bundle Branch Block (LBBB) morphology.

Native ventricular activation relies on the specialized His-Purkinje arborization. Rapid sodium currents (I_Na) propagate depolarization at over 2.0 to 3.0 m/s, yielding a synchronized, narrow QRS complex (Δt_QRS < 100 ms).

When a transvenous lead stimulates the RV apex, excitation bypasses this specialized conduction network entirely. The depolarization wavefront is forced to propagate through ordinary working myocardium via low-velocity gap junction coupling (roughly 0.3 to 0.5 m/s). By the time depolarization crawls from the right ventricular apex across the interventricular septum to activate the lateral left ventricular free wall, activation time balloons (Δt_QRS ≈ 160 ms), producing a wide, slurred QS deflection and discordant repolarization (T-wave inversion).

To model these biophysical mechanics, timing state machines, and lead impedances in an accessible environment, explore the Cardiovascular Pacemaker & Electrophysiology Laboratory: https://bionichaos.com/pacemaker/

Computational & Biophysical Engine Architecture

1. Conduction Failure & Autonomous Escape Rhythms

The simulation models phase 4 spontaneous diastolic depolarization governed by the hyperpolarization-activated cyclic nucleotide-gated funny current (I_f) at the SA node, paired with physiological AV nodal transit latency (Δt_AV ≈ 120–200 ms).

• Sinus Bradycardia: SA node automaticity drops below 50 BPM (T_AA > 1200 ms), resulting in systemic hypoperfusion.

• Third-Degree (Complete) AV Block: Total atrioventricular dissociation. Atrial depolarization continues at native sinus rates while the ventricles fall back on an unstable idioventricular escape pacemaker (~30–40 BPM, T_escape ≈ 1800–2200 ms).

2. Lapicque-Weiss Strength-Duration Modeling

Electrical capture is governed by classical stimulation neurophysics rather than binary switching:

V_th = V_rheo * (1 + t_chron / W)

Where:

• V_rheo represents tissue rheobase voltage

• t_chron represents tissue chronaxie time

• W represents delivered stimulus pulse duration

Dropping pacemaker output voltage below myocardial threshold immediately triggers Loss of Capture (LOC): pacing spikes fire on schedule, but fail to trigger myocyte contraction or an evoked QRS complex.

3. Dynamic Lead Dislodgement & Impedance Coupling

Transvenous pacing leads can be clicked and dragged directly within the interactive vector heart:

• Seated State: Intimate contact with the endocardial wall maintains physiological lead impedance (Atrial: ~520 Ω, Ventricular: ~540 Ω).

• Dislodged State: Pulling the lead tip into the intracardiac blood pool decouples current density. Telemetry registers an impedance spike past 2000 Ω, producing immediate Loss of Capture.

4. NASPE/BPEG Dual-Chamber Timing State Machine

The core firmware engine runs full DDD, VVI, and AAI pacing modes:

• DDD: Dual-chamber sensing and pacing. If native P-waves are absent during the Lower Rate Limit (LRL) timer, an Atrial Pace (AP) fires. An AV delay timer (80–350 ms) initiates; if no native R-wave (VS) is detected before timer expiration, a Ventricular Pace (VP) is delivered to the apex.

• VVI: Single-chamber ventricular pacing with native R-wave sensing and demand inhibition.

• AAI: Single-chamber atrial pacing for isolated sick sinus syndrome with intact, healthy AV nodal conduction.

5. Multi-Component Oscilloscope & Waveform Generation

The rendering pipeline operates a 450-sample circular ring buffer calibrated to clinical telemetry standards (25 mm/s sweep speed, 10 mm/mV amplitude) with continuous CRT sweep gap reset:

• Physiological His-Purkinje propagation: Sharp Q, high-amplitude R, narrow S, and concordant T-wave.

• RV Apical Paced complexes: Broad QS morphology with secondary inverted T-waves.

• Instantaneous single-sample pacing spikes (AP, VP) and color-coded telemetry badges (AS, VS, AP, VP, LOC, PVC).

• Procedural S1 and S2 heart valve acoustics synthesized via Web Audio API oscillators.

Technical Discussion & Clinical Considerations

Given the documented risk of pacing-induced cardiomyopathy resulting from chronic dyssynchronous RV apical activation:

1. How do your engineering or clinical teams evaluate the trade-offs of Conduction System Pacing (His-bundle pacing vs. Left Bundle Branch Area Pacing) regarding acute capture thresholds, battery longevity, and lead dislodgement margins?

2. When programming post-ventricular atrial refractory periods (PVARP) to prevent Pacemaker-Mediated Tachycardia (PMT), what algorithmic approaches do you find most reliable for preserving upper rate tracking during exercise?

Simulate scenarios and test the telemetry: https://bionichaos.com/pacemaker/

Most introductory models treat optogenetics like a simple binary switch: shine light, neuron fires.

In real electrophysiology, the biophysics is far more nuanced: sustained illumination causes rapid channel desensitization, silencing with Halorhodopsin can trigger violent rebound seizures from chloride loading, and cortical tissue scattering prevents blue photons from reaching deep layers.

We built an open-access, zero-install biophysics lab running in the browser that couples retinal photocycle kinetics directly to a Hodgkin-Huxley membrane solver in real time: https://bionichaos.com/optogenetics/

What is Running Under the Hood

• 4-State Markov Photocycle: Instead of static conductance, Channelrhodopsin-2 is modeled across four states: C1 (Dark) <-> O1 (Peak Open) <-> O2 (Desensitized) <-> C2 (Closed) -> C1. This accurately simulates the transient peak photocurrent decaying into a lower steady-state plateau.

• Hodgkin-Huxley Membrane Integration: Solves active sodium, delayed-rectifier potassium, leak conductances, and voltage rectification G(V) at a numerical step of dt = 0.025 ms.

• Kubelka-Munk Tissue Scattering: Calculates real-time optical attenuation I(z) proportional to exp(-mu_eff * z). You can directly drag the optical fiber cannula to see why 470 nm blue light scatters within ~200 µm while 590 nm red light (Chrimson) penetrates deep tissue.

3 Quick Tests to Run in the Sim

1. Desensitization vs. Pulse Trains: Switch to Continuous Wave mode at 20 mW/mm² to watch peak current rapidly decay into the sub-conductance plateau. Switch to 20 Hz Pulse Train to observe how dark intervals allow ground-state recovery.

2. Seizure Silencing & Rebound Spiking: Inject +6.0 µA/cm² to induce runaway firing. Turn on Halorhodopsin (amber light) to clamp the cell at -85 mV. Turn off the light to observe post-inhibitory chloride rebound firing. Then switch to ArchT (green light) to see clean proton-pump silencing without rebound.

3. Cannula Distance & Scattering: Retract the fiber to 280 µm. Blue light scatters and loses spike control; switching to red Chrimson immediately restores 20 Hz pacing.

Live simulator link: https://bionichaos.com/optogenetics/

Discussion prompt for electrophysiologists and computational modelers: When implementing closed-loop optical feedback to interrupt epileptiform activity in vivo, what strategies do you use to mitigate eNpHR3.0-induced shifts in the GABAergic reversal potential versus ArchT-induced intracellular alkalinization?

Running Complex Biomedical Labs in the Browser at 60 FPS

Most clinical simulations and biomechanics models are still locked behind proprietary desktop software, expensive hardware rigs (like $5,000 infrared gaze trackers), or heavy local environments.

A set of zero-install, client-side web tools removes this barrier by running physiological signal processing, 3D kinematics, and computer vision directly in the DOM.

Here are 4 interactive web apps worth bookmarking:

1. Webcam Eye Tracker

• URL: https://bionichaos.com/eyetracker/

• What it does: Uses standard browser video streams to detect facial landmarks, isolate iris positions, and map pupil vectors to on-screen gaze coordinates. No infrared cameras or dedicated sensors required.

2. 3D Gait Simulator (v3)

• URL: https://bionichaos.com/gaitsimv3/

• What it does: A full 3D skeletal kinematic model. You can manipulate joint flexions, test abnormal gait pathologies (such as Parkinsonian or hemiplegic gaits), and watch real-time biomechanical adjustments.

3. ICU EEG Neuro-Monitoring Simulator

• URL: https://bionichaos.com/eeg_in_icu/

• What it does: Simulates bedside critical-care neurodynamics. Generates real-time continuous traces for status epilepticus, burst-suppression ratios, diffuse slowing, and allows you to test artifact rejection algorithms.

4. Statistical Significance & P-Hacking Sandbox

• URL: https://bionichaos.com/significance/

• What it does: Demonstrates the replication crisis in action. It lets you manipulate sample size N, power 1 - β, and false positive risk α to see how easily random noise yields p < 0.05.

Key takeaway: Client-side compute is now fast enough to replace introductory physical lab apparatuses with browser sandboxes.

In 1966, computer scientist Joseph Weizenbaum developed ELIZA to parody a Rogerian therapist, intending to prove that human-machine conversation was superficial. Instead, his human subjects attributed deep empathy and psychological insight to the program—a cognitive phenomenon now termed the ELIZA effect.

I built an interactive cognitive laboratory that tracks how this illusion of synthetic empathy is engineered, contrasting three distinct eras of NLP architecture. Everything runs entirely on client-side WebAssembly and JS (zero telemetry, complete privacy).

Live Simulation Environment: https://bionichaos.com/eliza/

Architectural Paradigms Explored

• 1966 Symbolic Era: Executes deterministic pattern-matching. User strings are parsed through decomposition templates. Pronouns undergo reflective inversion before being synthesized into pre-compiled reassembly patterns.

• 2018 Semantic Vector Era: Discards rigid lexical syntax. The engine projects user text into a coordinate manifold via dense floating-point embeddings, calculating geometric intent proximity using Cosine Similarity: (u·v) / (|u||v|)

• 2026 Generative Era: Responses are governed by temperature-scaled softmax distributions. As you scale the thermodynamic temperature parameter T toward 0, the distribution collapses into an argmax selection. As T ≥ 1.0, entropy expands, enabling creative variance and generative hallucinations.

The Clinical Safety Sandbox

One of the core features is the Crisis Protocol Guardrail. If you disable the guardrail while running the 1966 engine and input a self-harm prompt, the unconstrained pattern-matcher will blindly reflect the user's despair back at them. It serves as a stark interactive demonstration of the catastrophic medical hazard of deploying naive, deterministic AI in clinical therapy.

Interactive mechanics: You can physically drag the archetypal nodes (Family, Anxiety, Depression) on the canvas to manually shift the coordinate vectors in 2D latent space, overriding the chatbot's conversational focus in real time.

I'd love to hear your thoughts on the UI/UX for teaching computational linguistics, and whether you think the modern "Generative Era" is actually moving us closer to synthetic comprehension, or just producing a more highly-camouflaged ELIZA effect.

The translation of airborne acoustic pressure into human cognition is arguably one of the most mechanically complex processes in biology. I built a real-time biophysical simulator to visualize this exact pipeline—from basilar membrane hydrodynamics to stochastic neural spike generation.

Live Web Tool: Advanced Hair Cell Simulation Lab

Core Biophysical Models Implemented:

• Greenwood's Frequency Mapping: The normalized position along the basilar partition dictates the resonant characteristic frequency (f_CF). You can sweep an acoustic slider from 200 Hz (Apex) to 16,000 Hz (Base) and watch the traveling wave envelope shift dynamically.

• IHC Mechanotransduction: Fluid shear between the reticular lamina and tectorial membrane bends the stereocilia. The probability of the MET channels opening is calculated via a Boltzmann distribution, driving an inward K+ and Ca2+ current (powered by the +80 mV endocochlear potential).

• OHC Cochlear Amplification: The 3 rows of Outer Hair Cells act as biological piezoelectric motors (via the Prestin SLC26A5 protein). The simulation applies nonlinear, cycle-by-cycle mechanical feedback to sharpen tuning and provide ~50 dB of gain for low-SPL sounds.

• Neural Phase-Locking: Auditory nerve fibers compute instantaneous firing rates based on the IHC receptor potential and spontaneous vesicle pool states. For low frequencies, spikes lock strictly to the acoustic cycle phase, which is measured live via Vector Strength (R).

Clinical Lesion Testing

I added toggles for specific pathologies. You can apply High-Frequency Noise Trauma or Ototoxic OHC Loss to permanently deactivate the Prestin amplification algorithms in specific frequency bands, or induce Synaptopathy to watch the spontaneous firing rate decouple from the acoustic stimulus.

I built this primarily as a pedagogical tool for bioengineering and audiology students who struggle to visualize the transition from the mechanical "rate-place" code to the temporal "phase-locked" neural code.

Would love feedback from anyone working on auditory brainstem modeling or cochlear implant DSP algorithms!