Real-time Dynamics in 3D for up to 1000 Qubits with Neural Quantum States: Quenches and the Quantum Kibble--Zurek Mechanism

Exponential complexity of many-body wave functions limits accurate numerical simulations of real-time dynamics, especially beyond 1D, where rapid entanglement growth poses severe challenges. Neural Quantum States (NQS) have emerged as a powerful approach for real-time dynamics in 2D, but their scalability and accuracy in 3D have remained an open challenge. Here, we establish NQS as a scalable framework for 3D quantum dynamics by introducing a residual-based convolutional architecture tailored to cubic spin lattices. Focusing on the 3D transverse-field Ising model, we demonstrate that NQS reliably capture distinct quench regimes, including collapse-and-revival dynamics and, most challengingly, the dynamics following a sudden quench to the quantum critical point. We perform finite-rate quenches to the critical point on lattices containing up to 1000 qubits, an unprecedented system size for numerical simulations of real-time dynamics beyond 1D. This enables the first large-scale numerical demonstration of the 3D quantum Kibble--Zurek mechanism. The QKZM in 3D is particularly intriguing because it lies at the upper critical dimension of the Ising universality class, where the standard power laws are modified by logarithmic factors together with prominent sub-leading logarithmic corrections. By deriving these corrections from renormalization-group flow equations up to two-loop order, we obtain a robust data collapse across all simulated system sizes for the correlation function, the excess energy, and the quantum Fisher information, the latter revealing universal multipartite-entanglement dynamics. In all cases, we find compelling agreement with the expected scaling dimensions. Our findings establish NQS as a scalable and reliable tool for exploring nonequilibrium phenomena in 3D quantum matter and for providing numerical benchmarks for 3D quantum simulators.