paper_study/20260924-123302_60hz_hum_s0_report.md · raw markdown
GarageFly run report: 20260924-123302_60hz_hum_s0
Generated 2026-09-24T12:38:28 · condition 60hz_hum · seed 0 · 30 s simulated (323 s wall)
Setup
- Stimulus: 60 Hz mains hum with harmonics, loudness-matched
- Audio source:
synthetic:hum:{'freq': 60.0} - Playback level: 85 dB SPL RMS at 1 m (track_rms); ambient noise 30 dB SPL
- Speaker distance from start position: 355 mm
- Baseline (silent) period: 0–10.0 s; stimulus period: 10.0–30.0 s
- Brain: sensory-motor mode, 52,547 neurons, 5,540,690 connections (FlyWire v783)
- Software: Python 3.12.13, PyTorch 2.14.0+cu130, MuJoCo 3.9.0, GPU NVIDIA GeForce RTX 3090
Key observations (this run)
Physics QA: passed (no jumps > 1 mm per 10 ms; thorax height stayed 0.5–2.5 mm above the surface; max step 0.80 mm).
- Auditory Johnston's-organ neurons (JO-A/B, mean of both sides) fired at 27.0 Hz during the stimulus vs 2.0 Hz at baseline.
- Wind/gravity-sensitive JO-C/E neurons: 3.8 Hz vs 2.0 Hz.
- Mean walking speed was 11.60 mm/s during the stimulus vs 11.98 mm/s at baseline (-3%).
- Turning bias: -40.7 deg/s vs -45.7 deg/s (+ = leftward); freezing fraction 0.00 vs 0.00.
- Giant-fiber (escape) DN rate: 2.9 Hz vs 0.0 Hz; 5 escape-command episodes during the stimulus (0 at baseline). Each can trigger a takeoff unless the fly is already airborne or within the post-landing refractory period.
- Behaviour: walking 36%, grooming 0.0%, flight 64% of the stimulus period (15.0 takeoffs/min, 0.0 grooming bouts/min) vs walking 100%, grooming 0.0%, flight 0% at baseline.
A single run is one stochastic sample. Differences between baseline and stimulus within a run, or between single runs, are not statistically established. Use
scripts/compare_conditions.pyover several seeds per condition for inferential statistics.
Behaviour
| Measure | Baseline | Stimulus | Change |
|---|---|---|---|
| Fraction of time walking | 1.00 | 0.36 | -0.64 |
| Fraction of time grooming (body) | 0.00 | 0.00 | 0.00 |
| Fraction of time in flight | 0.00 | 0.64 | 0.64 |
| Takeoffs per minute | 0.00 | 15.00 | 15.00 |
| Grooming bouts per minute | 0.00 | 0.00 | 0.00 |
| Righting events (fell over) | 0.00 | 0.00 | 0.00 |
| Mean walking speed while walking (mm/s) | 11.98 | 11.60 | -0.38 |
| Distance travelled (mm) | 119.81 | 617.30 | 497.49 |
| Fraction of time freezing | 0.00 | 0.00 | 0.00 |
| Mean | angular velocity | (deg/s) | 46.64 |
| Turning bias (deg/s, + = left) | -45.72 | -40.70 | 5.02 |
| Fraction of time grooming command active | 0.00 | 0.00 | 0.00 |
| Escape (giant fiber) command episodes | 0.00 | 5.00 | 5.00 |
| Escape command episodes per minute | 0.00 | 15.00 | 15.00 |
| Fraction of time escape command active | 0.00 | 0.25 | 0.25 |
| Mean distance to speaker (mm) | 325.33 | 404.35 | 79.02 |
Neural activity: sensory input and descending output
| Neuron group | Baseline (Hz) | Stimulus (Hz) | Change (Hz) |
|---|---|---|---|
| JO-A (left) | 1.98 | 18.38 | 16.40 |
| JO-A (right) | 2.04 | 22.57 | 20.53 |
| JO-B (left) | 1.93 | 30.87 | 28.93 |
| JO-B (right) | 2.07 | 36.29 | 34.22 |
| JO-C (left) | 1.89 | 3.97 | 2.07 |
| JO-C (right) | 2.02 | 3.45 | 1.43 |
| JO-D (left) | 1.94 | 12.84 | 10.90 |
| JO-D (right) | 2.10 | 11.83 | 9.73 |
| JO-E (left) | 2.00 | 3.73 | 1.73 |
| JO-E (right) | 2.06 | 4.11 | 2.05 |
| JO-F (left) | 2.01 | 1.98 | -0.02 |
| JO-F (right) | 2.04 | 1.94 | -0.10 |
| DNp09 (left) | 30.78 | 30.93 | 0.15 |
| DNp09 (right) | 30.78 | 31.02 | 0.24 |
| oDN1 (left) | 46.17 | 47.43 | 1.26 |
| oDN1 (right) | 39.97 | 42.11 | 2.13 |
| DNa01 (left) | 0.00 | 0.00 | 0.00 |
| DNa01 (right) | 0.00 | 0.00 | 0.00 |
| DNa02 (left) | 0.00 | 0.00 | 0.00 |
| DNa02 (right) | 6.20 | 6.91 | 0.71 |
| MDN (left) | 0.00 | 0.00 | 0.00 |
| MDN (right) | 0.00 | 0.00 | 0.00 |
| aDN2 (left) | 0.00 | 0.00 | 0.00 |
| aDN2 (right) | 0.00 | 0.00 | 0.00 |
| DNg11 (left) | 0.00 | 0.00 | 0.00 |
| DNg11 (right) | 0.00 | 0.00 | 0.00 |
| DNg12 (left) | 0.00 | 0.00 | 0.00 |
| DNg12 (right) | 0.00 | 0.00 | 0.00 |
| GF (left) | 0.00 | 0.00 | 0.00 |
| GF (right) | 0.00 | 5.70 | 5.70 |
| MN9 (left) | 0.00 | 0.00 | 0.00 |
| MN9 (right) | 0.00 | 0.00 | 0.00 |
Neural activity by FlyWire super-class (mean rate per neuron)
| Super-class | Baseline (Hz) | Stimulus (Hz) |
|---|---|---|
| ascending | 0.0000 | 0.0005 |
| central | 0.0003 | 0.0140 |
| descending | 0.1620 | 0.3017 |
| endocrine | 0.0000 | 0.0000 |
| motor | 0.0564 | 0.0877 |
| sensory | 0.1279 | 0.8139 |
| sensory_ascending | 0.0000 | 0.0000 |
| unannotated | 0.0000 | 0.0000 |
Most responsive cell types (largest |change| stimulus − baseline)
| cell_type | super_class | n_neurons | baseline_hz | stimulus_hz | delta_hz |
|---|---|---|---|---|---|
| JO-B4_a | sensory | 15 | 1.97 | 39.43 | 37.46 |
| JO-B2 | sensory | 36 | 1.96 | 35.02 | 33.06 |
| JO-B1_a | sensory | 105 | 1.98 | 34.78 | 32.80 |
| JO-A4 | sensory | 2 | 1.45 | 30.52 | 29.07 |
| JO-B3 | sensory | 29 | 1.97 | 30.92 | 28.95 |
| JO-B4_b | sensory | 25 | 2.01 | 29.78 | 27.77 |
| JO-B1_b | sensory | 27 | 1.96 | 29.18 | 27.22 |
| JO-B1_c | sensory | 22 | 2.12 | 28.19 | 26.06 |
| DNg29 | descending | 2 | 0.05 | 24.52 | 24.47 |
| JO-A2 | sensory | 33 | 2.05 | 22.91 | 20.86 |
| JO-A3 | sensory | 6 | 2.28 | 20.88 | 18.59 |
| JO-A1 | sensory | 21 | 2.01 | 19.63 | 17.62 |
| JO-A5 | sensory | 37 | 1.96 | 18.34 | 16.38 |
| SAD053 | central | 2 | 0.00 | 13.28 | 13.28 |
| JO-DA | sensory | 32 | 2.05 | 14.78 | 12.73 |
Most responsive downstream (non-sensory) cell types
| cell_type | super_class | n_neurons | baseline_hz | stimulus_hz | delta_hz |
|---|---|---|---|---|---|
| DNg29 | descending | 2 | 0.05 | 24.52 | 24.47 |
| SAD053 | central | 2 | 0.00 | 13.28 | 13.28 |
| DNg24 | descending | 2 | 0.00 | 10.93 | 10.93 |
| CB3913 | central | 1 | 0.00 | 10.45 | 10.45 |
| CB3911 | central | 1 | 0.00 | 10.40 | 10.40 |
| SAD013 | central | 2 | 0.00 | 10.30 | 10.30 |
| CB3876 | central | 1 | 0.00 | 10.15 | 10.15 |
| CB0478 | central | 2 | 0.05 | 9.78 | 9.72 |
| DNb06 | descending | 2 | 0.05 | 9.38 | 9.32 |
| CB0033 | central | 2 | 0.00 | 8.95 | 8.95 |
| DNp12 | descending | 2 | 0.00 | 6.58 | 6.58 |
| CB1817a | central | 2 | 0.00 | 6.58 | 6.58 |
| CB3491 | central | 4 | 0.00 | 6.08 | 6.08 |
| CB2789 | central | 4 | 0.00 | 5.97 | 5.97 |
| CB0466 | central | 2 | 0.00 | 5.85 | 5.85 |
Methods paragraph (draft, auto-generated)
A connectome-constrained leaky integrate-and-fire model of the adult Drosophila brain (FlyWire FAFB v783; 52,547 neurons and 5,540,690 connections in 'sensory-motor' mode; neuron parameters from Shiu et al. 2024) was simulated on a GPU at a 0.1 ms time step, in closed loop with the NeuroMechFly v2 body (FlyGym 2.1, MuJoCo 3.9.0). Audio (60 Hz mains hum with harmonics, loudness-matched) was calibrated to 85 dB SPL (RMS) at 1 m (track_rms calibration) and played from a virtual point source 355 mm from the fly's start position. Air particle velocity at each antenna was computed in 20 log-spaced bands (15–1500 Hz), including the acoustic near field and a cosine directional sensitivity. It was converted into Poisson firing rates of the FlyWire-annotated Johnston's organ neurons (JO-A to JO-E) through a resonant antenna filter and subgroup-specific frequency weighting. Descending-neuron activity (DNp09/oDN1, DNa01/DNa02, MDN, aDN2/DNg11/DNg12, giant fiber) was decoded into a two-dimensional locomotor drive for FlyGym's CPG-based controller, a grooming command (executed with measured front-leg grooming kinematics from an aDN-activated fly; Lobato-Rios et al. 2022) and an escape command (triggering an engineered takeoff-and-flight behaviour with externally applied lift; no wing aerodynamics). The first 10 s of each 30 s run were silent (baseline). Random seed 0.
Interpretation limits
- JO transfer functions, playback level, decoder gains and the tonic walking bias are arbitrary parameters (provenance in the run's JSON metadata). Conclusions should be checked for robustness to them.
- The ventral nerve cord is not simulated. Walking is produced by an engineered CPG controller driven by a few descending neurons.
- No substrate (bench) vibration, room acoustics, active antennal amplification, neuromodulation or gap junctions are modelled.
- Grooming uses measured kinematics but has no antennal contact mechanics. Flight is an engineered abstraction (applied lift/thrust, attitude PD); its speed, altitude and duration are arbitrary.
- Full list:
docs/MODEL_LIMITATIONS.md.
Files
- Raw data:
20260924-123302_60hz_hum_s0.h5(HDF5; behaviour, neural group rates, all spikes, stimulus) - Metadata:
20260924-123302_60hz_hum_s0.json - Machine-readable summary:
20260924-123302_60hz_hum_s0_summary.json