We introduce a framework for fermionic variational Monte Carlo in which a continuous normalizing flow refines a fixed antisymmetric base wavefunction. The flow is a permutation-equivariant neural ordinary differential equation (ODE), a smooth, topology-preserving map; equivariance preserves the antisymmetry of the base, so the flow can in principle improve any antisymmetric ansatz that can be sampled efficiently. The flow deforms the nodal surface of the base but cannot change its topology, so the base must supply the correct nodal topology. We use Slater and Jastrow-Slater bases in this work. Exact samples from the flow's Born distribution are obtained by pushing pre-cached base samples through the forward ODE, requiring no Markov chain Monte Carlo at training time. Base samples are generated offline and reused across batches and runs, enabling embarrassingly parallel training across multiple GPUs. We introduce three permutation-equivariant vector field architectures, Pairwise Deep Sets, FermiNet Vector Fields, and Pairwise Deep Sets Gradient, each balancing expressivity and computational cost differently. We further introduce an augmented dynamics formulation that co-evolves the derivatives required for the kinetic energy as ODE state variables, eliminating differentiation through the ODE trajectory and significantly reducing wall-clock time and memory. On harmonically trapped spinless electrons, the flow reaches ground-state energies below configuration interaction with single and double excitations, on par with a FermiNet baseline trained with the same optimizer. Training shows near-ideal strong scaling to 128 NVIDIA A100 GPUs on NERSC's Perlmutter for the largest systems tested, up to N = 48 particles in three dimensions.