REVIEW 2 major objections 4 minor 40 references
A sub-Riemannian model of the motor cortex with Wasserstein distance
T0 review · 2 major / 4 minor · reviewed 2026-07-13 · grok-4.5
Pith's one-line read Horizontal curves in a sub-Riemannian model of motor cortex automatically satisfy the observed speed-curvature law, and Wasserstein clustering of fragments recovers experimental neural states.
desk verdict Solid incremental step: Wasserstein lets their prior sub-Riemannian fragment model finally work on real variable-duration M1 data and recover the Kadmon states, while the speed-curvature relation falls out of the horizontal distribution by construction. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The sub-Riemannian structure generated by the vector fields X1 = v cos θ ∂x + v sin θ ∂y + a ∂v + ∂t, X2 = ∂θ, X3 = ∂a, whose horizontal curves automatically obey R = v/α₂, together with the Wasserstein distance on the vector-valued densities μ = (sin θ, cos θ, a)/∫√(1+a²) that supplies the affinity matrix for spectral clustering into neural states.
What would settle it
Apply the same Wasserstein affinity matrix (same density definition and metric weights) to a new multi-electrode reaching dataset recorded under comparable conditions; if the resulting spectral clusters fail to match the experimentally reported neural states, the central claim is falsified.
Extended reading notes
Core claim
A sub-Riemannian geometry on the feature space of position, time, orientation, velocity and acceleration makes its horizontal curves satisfy the speed-curvature relation observed in reaching movements; the Wasserstein distance between the vector-valued measures associated with those curves then clusters real movement fragments into the same neural states reported from cortical recordings, outperforming the Sobolev distance that required all trajectories to share a common time interval.
Load-bearing premise
The model assumes that the particular vector-valued density built from orientation and acceleration, together with metric weights that nearly discard absolute time, correctly represents how motor-cortex cells compare movement fragments.
Editorial extensions
If this is right
- The speed-curvature law need not be imposed by hand; it is a differential consequence of the chosen feature geometry.
- Cortical connectivity among fragments can be modelled by a Wasserstein kernel that tolerates variable durations without reparametrization.
- Spectral clustering of that kernel recovers the experimentally observed neural states from real reaching data.
- The same geometric setting that organizes elementary features into fragments also organizes fragments into higher-order neural states.
Reading between the lines
- The same Wasserstein-on-horizontal-curves construction may apply to other cortical areas whose cells are selective to short spatiotemporal trajectories.
- Metric weights that effectively discard absolute time imply a short-term shape memory that is shift-invariant but not scale-invariant in time.
- Replacing hand-tuned coefficients with anatomically or data-derived weights would yield a fully parameter-free prediction of neural-state boundaries.
- Population codes in motor cortex should be more sensitive to the joint distribution of orientation and acceleration than to absolute position or absolute time.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper models primary motor cortex (M1) features as the 6-dimensional manifold M = R^{2}_{(x,y)} imes R^{+}_t imes S^{1}_ heta imes R_v imes R_a equipped with the sub-Riemannian structure generated by the horizontal fields X1, X2, X3 of Eq. (3). Horizontal (admissible) curves are proposed as models of movement fragments; the differential constraints immediately yield the curvature–speed relation k = heta'/v (Eq. 8) and hence R = v/ heta' (Eq. 9), recovering the experimentally observed speed–curvature law (and the 2/3-power law by the special choice heta' = v^{-2}). Fragments are then represented by real- or vector-valued probability measures (Eqs. 13–14) and clustered by spectral analysis of a Gaussian kernel of the Wasserstein distance. On synthetic data the method recovers coherent groups (silhouette 0.65); on the real reaching trajectories of Kadmon-Harpaz et al. it recovers a partition visually comparable to the experimentally reported neural states (silhouette 0.35), outperforming the earlier Sobolev distance that required identical time domains.
Significance. If the geometric derivation and the clustering results hold, the work supplies a single sub-Riemannian framework that simultaneously (i) explains the classical speed–curvature relation from first principles and (ii) organizes experimentally recorded fragments into the same neural states previously obtained from multi-electrode recordings. The introduction of Wasserstein distance into cortical connectivity modelling is novel and removes the artificial re-parametrization required by Sobolev-type metrics. The paper therefore offers a concrete, falsifiable link between differential geometry, optimal transport and motor-cortex physiology that can be tested on further data sets and with alternative measures.
major comments (2)
- Section 4.3 and Definition 7 / Eq. (11): the recovery of the Kadmon-Harpaz neural states rests on the explicit choice c5 = 0.008 for the temporal weight, which the authors state “completely discard[s] the time variable,” followed by re-parametrization of every fragment to a common initial time. Once absolute timing is removed, the remaining features (orientation and acceleration profile) are precisely those already used by the experimental clustering. No sensitivity analysis or ablation with non-vanishing c5 (or with alternative normalizations of the measure μ) is provided; without it the claimed agreement with experiment cannot be regarded as an independent validation of the sub-Riemannian + Wasserstein model.
- Section 3.1.3, Eq. (14) and Section 4.3: the vector-valued density μ = (sin heta, cos heta, a)/∫√(1+a^{2}) is introduced without physiological or information-theoretic justification, and the kernel widths σ1 = 0.7, σ2 = 0.09 are hand-tuned. The silhouette scores (0.65 synthetic, 0.35 real) are reported, yet no quantitative comparison against alternative distances (e.g., dynamic time warping, plain L^{2} on the same features) or against shuffled baselines is given. Consequently the claim that Wasserstein distance is “much more efficient” remains only qualitative.
minor comments (4)
- Eq. (14) writes the components of μ as (sin heta, cos heta, a) while the surrounding text lists (cos heta, sin heta, a); the order should be made consistent.
- Figure 7 caption and the accompanying text refer to a “center-out task,” yet the data provenance (Hatsopoulos laboratory) is only acknowledged in the acknowledgments; a brief methods paragraph describing the recording and pre-processing pipeline would improve reproducibility.
- The parameters c1…c4 of the ball-box norm (Eq. 11) are said to be “tuned o experimentally observed clusterization,” but their concrete values are never listed; they should be reported for the real-data experiment.
- Typographical inconsistencies appear throughout (e.g., “Hatsoupolos,” “real-valued probability: speed and curvature,” missing spaces after commas in several equations).
Circularity Check
Neural-state recovery is obtained by tuning metric weights (esp. c5=0.008 discarding absolute time) and by defining the vector measure on precisely the orientation+acceleration features that already define the target clusters of Kadmon-Harpaz et al.
-
fitted input called prediction
[§3.1.2 (after Def. 7 / Eq. 11) and §4.3]
"The parameters c1, ...,c4 in (11) can be arbitrary. In the sequel, they will be tuned in such a way to obtain the experimentally observed clusterization. ... we use a small coefficient c5=0.008 for the temporal increment e5 ... Setting the coefficient c5=0.008 means that we completely discard the time variable, which is consistent with the classification proposed in [19]"
The free metric weights (especially the near-zero temporal weight that forces time-shift invariance and common reparametrization) are adjusted until the spectral clusters match the experimental neural states; the subsequent claim that the model “recovers” those states is therefore a fit, not an independent prediction.
-
self definitional
[§3.1.1 Eqs. (8)–(9)]
"k= α2 / v ... R=1/α2 v ... Here, we claim that this relation is already encoded in our choice of variables. We can also recover the relation R=Cv3 by simply choosing α2=v−2."
The experimentally observed speed-curvature relation is an immediate algebraic consequence of the differential constraints built into the generators X1,X2,X3; the stronger power-law form is obtained only by a free functional choice of the control α2. The geometry therefore does not derive the law but encodes it by definition.
1 more flagged steps
-
fitted input called prediction
[§3.1.3 Def. of μ(γ) and §4.3 (affinity parameters)]
"we consider only three components (cos(θ(s)),sin(θ(s)),a(s)). Then we define μ(γ(s))=(sin(θ(s)),cos(θ(s)),a(s))/∫√(1+a2(s))ds ... we set the parameters (σ1=0.7) and (σ2=0.09) to enhance the sensitivity of the similarity measure to both orientation and acceleration features."
The vector measure is deliberately restricted to the two kinematic features (orientation and acceleration phase) that already define the experimental neural states; the kernel widths are then hand-chosen to emphasize exactly those features. Clustering on this pre-selected representation necessarily groups by the same criteria used in the target partition.
full rationale
The paper has two main claims. The speed-curvature relation follows immediately from the chosen generators of the horizontal distribution (k=α2/v), so it is encoded by construction rather than independently predicted; the 2/3-power law is recovered only by the free choice α2=v^{-2}. The more load-bearing empirical claim—that Wasserstein clustering recovers the experimental neural states—is obtained after the free coefficients of the ball-box metric are explicitly tuned “to obtain the experimentally observed clusterization,” with c5 set to 0.008 so that absolute time is discarded and every fragment is reparametrized to a common start (exactly the normalization already used by the experimental paper). The vector-valued density is likewise defined on (sin θ, cos θ, a), the two features that already label the experimental states. Once these choices are made, spectral clustering on the resulting affinity matrix reproduces the target partition; the match is therefore partly engineered by free parameters rather than a parameter-free prediction of the sub-Riemannian+Wasserstein model. The geometry itself and the superiority of Wasserstein over Sobolev on variable-length trajectories retain independent content, so the circularity is only partial (score 5).
Assumptions & free parameters
free parameters (3)
- metric weights c_i (esp. c5=0.008) =
c5=0.008, others=1
- kernel widths σ1,σ2 =
σ1=0.7, σ2=0.09
- α2 control function
assumptions (3)
- domain assumption Motor-cortex cells are selective to the six features (x,y,t,θ,v,a) and the admissible curves are integral curves of the given horizontal distribution with the specific α3 form that fits measured acceleration profiles.
- standard math The Hörmander condition holds for the chosen generators, so a Carnot-Carathéodory distance exists and can be estimated by the ball-box metric.
- ad hoc to paper Cortical connectivity in the space of fragments is well-modeled by a Gaussian kernel of the Wasserstein distance on the chosen vector-valued density.
invented entities (1)
-
vector-valued density μ=(sinθ,cosθ,a)/∫√(1+a²) on fragments
Cite this review
Pith. "Pith review of A sub-Riemannian model of the motor cortex with Wasserstein distance." pith.science (2026). https://pith.science/paper/6OO3BXMF
@misc{pith2026260320756,
author = {Pith},
title = {Pith review of: A sub-Riemannian model of the motor cortex with Wasserstein distance},
year = {2026},
howpublished = {\url{https://pith.science/paper/6OO3BXMF}},
note = {Machine review of arXiv:2603.20756}
}
read the original abstract
This study aims to better understand the functional geometry of the motor cortex, starting from different sources of experimental evidence. Recent studies have proved that cells of the primary motor cortex (M1) are sensitive to short hand trajectories called fragments. Here, we propose a sub-Riemannian higher-dimensional geometry accounting for geometric and kinematic properties. Due to the constraints of the geometry, horizontal curves naturally satisfy a relation between geometric and kinematic properties experimentally observed. In the space of trajectories, we also apply a clustering algorithm based on the Wasserstein distance: we obtain a grouping which nicely fits the observed experimental data much more efficiently than the Sobolev distance.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
Sur la vitesse des mouvements graphiques,
A. Binet and J. Courtier, “Sur la vitesse des mouvements graphiques,”Revue Philosophique de la France et de l’Etranger, vol. 35, pp. 664–671, 1893. 20
-
[2]
Movement timing and in- variance arise from several geometries,
D. Bennequin, R. Fuchs, A. Berthoz, and T. Flash, “Movement timing and in- variance arise from several geometries,”PLoS computational biology, vol. 5, no. 7, p. e1000426, 2009
2009
-
[3]
Affine differential geometry analysis of human arm movements,
T. Flash and A. A. Handzel, “Affine differential geometry analysis of human arm movements,”Biological cybernetics, vol. 96, pp. 577–601, 2007
2007
-
[4]
Geometric methods in the study of human mo- tor control,
A. A. Handzel and T. Flash, “Geometric methods in the study of human mo- tor control,” Cognitive Studies: Bulletin of the Japanese Cognitive Science Society, vol. 6, no. 3, pp. 309–321, 1999
1999
-
[5]
Constant affine velocity predicts the 13 power law of planar motion perception and generation,
F. E. Pollick and G. Sapiro, “Constant affine velocity predicts the 13 power law of planar motion perception and generation,” Vision research, vol. 37, no. 3, pp. 347–353, 1997
1997
-
[6]
Motion from point matches using affine epipolar geometry,
L. S. Shapiro, A. Zisserman, and M. Brady, “Motion from point matches using affine epipolar geometry,” in Computer Vision—ECCV’94: Third European Conference on Computer Vision Stockholm, Sweden, May 2–6 1994 Proceedings, V olumeII 3, pp. 73–84, Springer, 1994
1994
-
[7]
Several geometries for movements gener- ations,
D. Bennequin and A. Berthoz, “Several geometries for movements gener- ations,” Geometric and Numerical Foundations of Movements, pp. 13–42, 2017
2017
-
[8]
Static spatial effects in motor cortex and area 5: quantitative relations in a two-dimensional space,
A. Georgopoulos, R. Caminiti, and J. Kalaska, “Static spatial effects in motor cortex and area 5: quantitative relations in a two-dimensional space,” Experimental Brain Research, vol. 54, no. 3, pp. 446–454, 1984
1984
Show all 40 references
-
[9]
R. E. Kettner, A. B. Schwartz, and A. P. Georgopoulos, “Primate motor cortex and free arm movements to visual targets in three-dimensional space. iii. positional gradients and population coding of movement direction from various movement origins,”Journal of Neuroscience, vol. ...
1988
-
[10]
On the relations between the direction of two-dimensional arm movements and cell discharge in primate motor cortex,
A. P. Georgopoulos, J. F. Kalaska, R. Caminiti, and J. T. Massey, “On the relations between the direction of two-dimensional arm movements and cell discharge in primate motor cortex,”Journal of Neuroscience, vol. 2, no. 11, pp. 1527–1537, 1982
1982
-
[11]
Primate motor cortex and free arm movements to visual targets in three-dimensional space. i. 21 relations between single cell discharge and direction of movement,
A. B. Schwartz, R. E. Kettner, and A. P. Georgopoulos, “Primate motor cortex and free arm movements to visual targets in three-dimensional space. i. 21 relations between single cell discharge and direction of movement,”Journal of Neuroscience, vol. 8, no. 8, pp. 2913–2927, 1988
1988
-
[12]
Motor cortical representation of speed and direction during reaching,
D. W. Moran and A. B. Schwartz, “Motor cortical representation of speed and direction during reaching,” Journal of neurophysiology, vol. 82, no. 5, pp. 2676–2692, 1999
1999
-
[13]
A model of movement coordi- nates in the motor cortex: posture-dependent changes in the gain and direction of single cell tuning curves,
R. Ajemian, D. Bullock, and S. Grossberg, “A model of movement coordi- nates in the motor cortex: posture-dependent changes in the gain and direction of single cell tuning curves,”Cerebral Cortex, vol. 11, no. 12, pp. 1124–1135, 2001
2001
-
[14]
From the motor cortex to the movement and back again,
W. W. Teka, K. C. Hamade, W. H. Barnett, T. Kim, S. N. Markin, I. A. Rybak, and Y . I. Molkov, “From the motor cortex to the movement and back again,” PloS one, vol. 12, no. 6, p. e0179288, 2017
2017
-
[15]
Encoding of movement fragments in the motor cortex,
N. G. Hatsopoulos, Q. Xu, and Y . Amit, “Encoding of movement fragments in the motor cortex,”Journal of Neuroscience, vol. 27, no. 19, pp. 5105–5114, 2007
2007
-
[16]
The problem of parametric neural coding in the motor system,
J. Reimer and N. G. Hatsopoulos, “The problem of parametric neural coding in the motor system,” in Progress in motor control, pp. 243–259, Springer, 2009
2009
-
[17]
The cortical control of movement revisited,
M. S. Graziano, C. S. Taylor, T. Moore, and D. F. Cooke, “The cortical control of movement revisited,”Neuron, vol. 36, no. 3, pp. 349–362, 2002
2002
-
[18]
Mapping behavioral repertoire onto the cortex,
M. S. Graziano and T. N. Aflalo, “Mapping behavioral repertoire onto the cortex,”Neuron, vol. 56, no. 2, pp. 239–251, 2007
2007
-
[19]
Move- ment decomposition in the primary motor cortex,
N. Kadmon Harpaz, D. Ungarish, N. G. Hatsopoulos, and T. Flash, “Move- ment decomposition in the primary motor cortex,”Cerebral cortex, vol. 29, no. 4, pp. 1619–1633, 2019
2019
-
[20]
Computational approaches to motor control,
T. Flash and T. J. Sejnowski, “Computational approaches to motor control,” Current opinion in neurobiology, vol. 11, no. 6, pp. 655–662, 2001
2001
-
[21]
Spatial constant equi-affine speed and motion percep- tion,
U. Maoz and T. Flash, “Spatial constant equi-affine speed and motion percep- tion,”Journal of neurophysiology, vol. 111, no. 2, pp. 336–349, 2014. 22
2014
-
[22]
Affine differential geometry and smoothness maximization as tools for identifying geometric movement primitives,
F. Polyakov, “Affine differential geometry and smoothness maximization as tools for identifying geometric movement primitives,”Biological cybernetics, vol. 111, pp. 5–24, 2017
2017
-
[23]
Jean, Control of nonholonomic systems: from sub-Riemannian geometry to motion planning
F. Jean, Control of nonholonomic systems: from sub-Riemannian geometry to motion planning. Springer, 2014
2014
-
[24]
Optimal control models of the goal-oriented human locomotion,
F. Jean, “Optimal control models of the goal-oriented human locomotion,” in Talk given at the “Workshop on Nonlinear Control and Singularities”, Porquerolles, France, pp. 24–28, 2010
2010
-
[25]
Functional architecture of m1 cells encod- ing movement direction,
C. Mazzetti, A. Sarti, and G. Citti, “Functional architecture of m1 cells encod- ing movement direction,”Journal of Computational Neuroscience, pp. 1–29, 2023
2023
-
[26]
A sub-riemannian model of neural states in the primary motor cortex,
C. Mazzetti, J. Ali, A. Sarti, and G. Citti, “A sub-riemannian model of neural states in the primary motor cortex,”arXiv preprint arXiv:2501.03247, 2024
2024 arXiv
-
[27]
Spatial trajectories and reaction times of aimed movements: effects of practice, uncertainty, and change in target location,
A. P. Georgopoulos, J. F. Kalaska, and J. T. Massey, “Spatial trajectories and reaction times of aimed movements: effects of practice, uncertainty, and change in target location,” Journal of neurophysiology, vol. 46, no. 4, pp. 725–743, 1981
1981
-
[28]
Encoding of movement frag- ments in the motor cortex,
Q. Hatsopoulos, Nicholas G; Xu and Y . Amit, “Encoding of movement frag- ments in the motor cortex,”Journal of Neuroscience, vol. 27, p. 5105–5114, 2007
2007
-
[29]
Neural population dynamics during reaching,
M. M. Churchland, J. P. Cunningham, M. T. Kaufman, J. D. Foster, P. Nuyu- jukian, S. I. Ryu, and K. V . Shenoy, “Neural population dynamics during reaching,”Nature, vol. 487, no. 7405, pp. 51–56, 2012
2012
-
[30]
Villani et al., Optimal transport: old and new, vol
C. Villani et al., Optimal transport: old and new, vol. 338. Springer, 2008
2008
-
[31]
A wasserstein norm for signed mea- sures, with application to nonlocal transport equation with source term,
B. Piccoli, F. Rossi, and M. Tournus, “A wasserstein norm for signed mea- sures, with application to nonlocal transport equation with source term,”arXiv preprint arXiv:1910.05105, 2019
1910 arXiv
-
[32]
vwcluster: Vector- valued optimal transport for network based clustering using multi-omics data in breast cancer,
J. Zhu, J. H. Oh, J. O. Deasy, and A. R. Tannenbaum, “vwcluster: Vector- valued optimal transport for network based clustering using multi-omics data in breast cancer,”Plos one, vol. 17, no. 3, p. e0265150, 2022. 23
2022
-
[33]
Optimal mass transport: Signal processing and machine-learning applications,
S. Kolouri, S. R. Park, M. Thorpe, D. Slepcev, and G. K. Rohde, “Optimal mass transport: Signal processing and machine-learning applications,” IEEE signal processing magazine, vol. 34, no. 4, pp. 43–59, 2017
2017
-
[34]
A sub-riemannian model of the func- tional architecture of m1 for arm movement direction,
C. Mazzetti, A. Sarti, and G. Citti, “A sub-riemannian model of the func- tional architecture of m1 for arm movement direction,” In: Nielsen, F., Barbaresco, F. (eds) Geometric Science of Information. GSI 2023. Lecture Notes in Computer Science, vol 14072. Springer, pp. 483–492, 2023
2023
-
[35]
Characteristics of random nets of analog neuron-like elements,
S. I. Amari, “Characteristics of random nets of analog neuron-like elements,” IEEE Transactions on Systems, Man, and Cybernetics, vol. 2, no. 5, pp. 643– 657, 1972
1972
-
[36]
Excitatory and inhibitory interactions in localized populations of model neurons,
H. R. Wilson and J. D. Cowan, “Excitatory and inhibitory interactions in localized populations of model neurons,” Biophysical Journal, vol. 12, no. 1, pp. 1–24, 1972
1972
-
[37]
Large scale spatially organized activity in neural nets,
G. B. Ermentrout and J. D. Cowan, “Large scale spatially organized activity in neural nets,”SIAM Journal on Applied Mathematics, vol. 38, no. 1, pp. 1–21, 1980
1980
-
[38]
The functional geometry of local and horizontal connections in a model of v1,
P. C. Bressloff and J. D. Cowan, “The functional geometry of local and horizontal connections in a model of v1,”Journal of Physiology-Paris, vol. 97, no. 2-3, pp. 221–236, 2003
2003
-
[39]
Some theoretical and numerical results for delayed neural field equations,
G. Faye and O. Faugeras, “Some theoretical and numerical results for delayed neural field equations,” Physica D: Nonlinear Phenomena, vol. 239, no. 9, pp. 561–578, 2010
2010
-
[40]
The constitution of visual perceptual units in the functional architecture of v1,
A. Sarti and G. Citti, “The constitution of visual perceptual units in the functional architecture of v1,”Journal of computational neuroscience, vol. 38, no. 2, pp. 285–300, 2015. 24
2015
Reviewed July 13, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.