Finite-sample deviations and convergence in the statistics of Bohmian trajectory ensembles

We analyze finite-sample statistics of Bohmian trajectories for single spinless and spin-1/2 particles. Equivariance ensures agreement with |ψ|² in the quantum equilibrium limit, yet experiments and simulations necessarily use finite ensembles. We show that in regular flows (e.g., wavepackets or low-mode superpositions of eigenstates of harmonic oscillators) sample means and/or variances over modest N are consistent with Born-rule moments. In contrast, degenerate superpositions of 3D oscillators with nodal barriers and chaotic Bohmian dynamics exhibit sensitive dependence on initial conditions and complex flow partitioning, which can yield noticeable finite-sample deviations in the mean and variance. For the spin-1/2 particle, both convective and Pauli currents conserve |ψ|², but they are associated with different velocity fields and thus might yield different finite-sample trajectory statistics. These findings calibrate the interpretation of trajectory-based uncertainty and provide practical guidance for numerical Bohmian simulations of spin and transport, without challenging the equivalence to orthodox quantum mechanics in the quantum equilibrium ensemble.