Not yet reviewed by Pith; the record is open.

We extend correlated sampling from classical auxiliary-field quantum Monte Carlo to the quantum-classical (QC-AFQMC) framework, enabling accurate nuclear force computations crucial for geometry optimization and reaction dynamics. Stochastic electronic structure methods typically encounter prohibitive statistical noise when computing gradients via finite differences. To address this, our approach maximizes correlation between nearby geometries by synchronizing random number streams, aligning orbitals, using deterministic integral decompositions, and employing a consistent set of classical shadow measurements defined at a single reference geometry. Crucially, reusing this single, reference-defined shadow ensemble eliminates the need for additional quantum measurements at displaced geometries. Together, these methodological choices substantially reduce statistical variance in computed forces. We validate the method across hydrogen chains, confirming accuracy throughout varying correlation regimes, and demonstrate significant improvements over single-reference methods in force evaluations for N$_2$ and stretched linear H$_4$, particularly in strongly correlated regions where conventional coupled cluster approaches qualitatively fail. Orbital-optimized trial wave functions further boost accuracy for demanding cases such as stretched CO$_2$, without increasing quantum resource requirements. Finally, we apply our methodology to the MEA-CO$_2$ carbon capture reaction, employing quantum information metrics for active space selection and matchgate shadows for efficient overlap evaluations, establishing QC-AFQMC as a robust framework for exploring complex reaction pathways.