Neural quantum states have a credibility problem. Every few months, a paper announces that a neural network has tackled some quantum many-body problem, and every few months, a closer look reveals the system sizes were too small, the timescales too short, or the model itself too tailored to the specific physics to be interesting. Markus Heyl's group at the University of Augsburg has published something different.
In a paper posted to arXiv on April 6, 2026, Heyl and collaborators show that a residual 3D convolutional neural network can simulate real-time quantum dynamics for systems of up to 1000 qubits in three spatial dimensions. The scale is real. The 3D geometry is real. And the physics being studied is genuinely hard to access by other means.
The central problem is one that has defined computational quantum physics for decades. Exact simulation of quantum many-body systems scales exponentially with system size, which is what makes quantum computers theoretically attractive and classical simulation maddeningly difficult. For static properties, DMRG and quantum Monte Carlo have carved out powerful niches. For real-time dynamics in three dimensions, particularly near critical points, the established toolkit has gaps.
The team used a residual 3D CNN with 3 by 3 by 3 convolutional layers, GELU activation, residual connections, and gated linear units. The time-dependent variational principle drove the real-time evolution. Architecture details matter here: the channel width was fixed at 4 and the depth was either 3 or 4 layers. Because the architecture uses weight sharing across all lattice sites, the variational parameters do not scale with system size. Train on a small cube, apply to a larger one.
Applied to the 3D transverse-field Ising model, the network produced what the authors call the first large-scale numerical demonstration of the quantum Kibble-Zurek mechanism in three dimensions at the upper critical dimension. The upper critical dimension is where mean-field theory becomes marginally correct and logarithmic corrections to power-law scaling appear. For the 3D Ising universality class, the relevant combination is d plus z equals 4, where d is spatial dimension and z is the dynamical critical exponent. Getting numerical confirmation of the scaling predictions in this regime has been difficult because the logarithmic corrections are subtle and the system sizes needed to see them are large.
The team validated against renormalization group predictions including corrections up to two-loop order. They found a critical field of 5.158136J. They tracked the KZ scaling of defect density and confirmed it matches the RG predictions within their numerical resolution. And they computed the quantum Fisher information density, which reached 29, implying genuine multipartite entanglement across at least 30 spins. That number is verifiable from the data in the paper.
There is a caveat that needs to be stated clearly: this is a classical neural network running on classical hardware. The 1000 qubits is a count of quantum degrees of freedom being simulated by a classical computer, not a quantum processor performing computations. The relevant comparison is not quantum hardware but other classical methods. For the specific problem of 3D quantum critical dynamics, the combination of system size and accessible timescale is beyond what straightforward exact diagonalization or standard quantum Monte Carlo can reach in that geometry. Monte Carlo struggles with real-time dynamics. Exact diagonalization hits exponential walls fast in 3D. The NQS approach, for this class of problems, occupies a regime that is otherwise nearly intractable.
The entanglement result deserves a similar caveat. A quantum Fisher information density of 29 across 30 spins means the network is capturing genuinely quantum correlations among a subset of the system. It does not mean the full 1000-qubit state is deeply entangled. For a problem this hard, the partial entanglement is the point, not a limitation.
For quantum computing itself, this paper does not change what is possible on actual quantum hardware. For quantum many-body physics, it shows that neural quantum states can produce verifiable, physically meaningful results for 3D nonequilibrium dynamics at a scale that matters for testing RG predictions. The question is not whether NQS replaces other approaches but which problems it makes newly tractable.
The answer in this case is narrow but genuine: 3D quantum critical dynamics near the upper critical dimension, accessible via real-time evolution at system sizes large enough to see the logarithmic corrections that mean-field theory misses. That is not a small thing. For a field that has spent decades trying to verify scaling predictions at the boundary of tractability, it is exactly the kind of numerical result that moves the conversation forward.
Heyl's group at University of Augsburg has shown that neural quantum states can do something non-trivial in three dimensions at a scale that counts. Whether that scale is large enough to matter broadly will depend on what problems physicists bring to the method next.