We introduce a new theoretical formalism to compute solid-state vibrational circular dichroism (VCD) spectra from molecular dynamics simulations. Having solved the origin-dependence problem of the periodic magnetic gauge, we present IR and VCD spectra of (1S,2S )-trans-1,2-cyclohexanediol obtained from first-principles molecular dynamics calculations and nuclear velocity perturbation theory, along with the experimental results. As the structure model imposes periodic boundary conditions, the common origin of the rotational strength has hitherto been ill-defined and was approximated by means of averaging multiple origins. The new formalism relies on the velocity representation of VCD and the time-correlation function, whose symmetry properties are exploited to reconnect the periodic model with the finite physical system. It restores the gauge freedom of finite models, but still allows fully accounting for non-local spatial coupling terms from the gauge transport term. We show that even for small simulation cells the rich nature of solid-state VCD spectra found in experiments can be reproduced to a very satisfactory level. While the general workflow to compute solid-state VCD spectra relies on the interplay of experimental data and theoretical simulation, expressions of VCD in liquid state and for isolated systems can be derived as simplification of the general equations.