⚛️ quantum physics

Finite Systems in the Space of Projective Measurements

This paper investigates the dynamics of two-level systems under permutation-invariant Hamiltonians by constructing an approximate star product for projected $Q$ functions to rewrite the Liouville equation, demonstrating that the method accurately models evolution for permutation-invariant states while showing deviations for arbitrary states and noninvariant Hamiltonians.

Original authors: Miguel Flores R., Carlos A. Muñoz

Published 2026-10-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Miguel Flores R., Carlos A. Muñoz

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine trying to understand the behavior of a crowd by tracking every single person's heartbeat, step, and thought. In the quantum world, where particles like electrons or atoms behave according to the strange rules of the very small, this task becomes impossible. When scientists study a system made of many such particles, the amount of information required to describe it grows so fast that it quickly overwhelms even the most powerful computers. Furthermore, in many real-world situations, the particles are identical and indistinguishable; you cannot label one "Alice" and another "Bob" to follow them separately. This creates a massive hurdle for understanding how large groups of quantum particles move and change over time. To get around this, researchers often look for ways to describe the system as a whole, focusing on the collective properties that emerge when many particles act together, rather than getting lost in the details of each individual.

A team of physicists in Mexico has developed a new way to map this collective behavior, offering a clearer window into how complex quantum systems evolve. They focused on a specific mathematical tool used to visualize quantum states, known as a quasidistribution function. Think of this function as a topographical map that shows where the "weight" or probability of a quantum system is located. For systems with many particles, this map usually becomes a chaotic, high-resolution grid that is too crowded to read. The researchers' breakthrough was to project this messy, high-dimensional map onto a much simpler, three-dimensional space that captures only the essential, collective features of the system. By doing this, they could watch the system's evolution without being blinded by the sheer volume of data.

The core of their work involved creating a new rule for how these simplified maps change when the system is acted upon by energy or forces. In quantum mechanics, the way a system changes is governed by a specific equation involving a density matrix, which describes the state of the system, and a Hamiltonian, which describes the energy or forces acting on it. To predict the future state, one must multiply these two mathematical objects together. In the traditional, complex grid, this multiplication is handled by a specific operation called a star product. The researchers constructed an approximate version of this star product that works within their simplified three-dimensional space. This allowed them to rewrite the fundamental equation of motion, known as the Liouville equation, in a way that could be solved much more easily.

To test their method, the team simulated the evolution of several different types of quantum states under various conditions. They looked at states that are perfectly symmetrical, where swapping any two particles leaves the system unchanged, as well as states that are not symmetrical. They also tested their method against Hamiltonians that respect particle symmetry and those that do not. When they compared the results of their simplified, approximate method against the exact, computationally heavy calculations, they found the two matched almost perfectly for systems where the particles and the forces were symmetrical. The difference between the approximate map and the exact map was so small it was barely measurable, indicating that their simplified view retained all the necessary information for these cases.

Surprisingly, the method also worked very well for systems that were not symmetrical, provided the forces acting on them were symmetrical. Even when the initial arrangement of particles was messy or uneven, the simplified three-dimensional map tracked the system's evolution with high accuracy. This suggests that the collective view is robust enough to handle disorder in the initial state, as long as the rules governing the system's change are orderly. However, the researchers found that when the forces themselves were not symmetrical, the approximation began to drift away from the exact truth, especially as time went on. In these cases, the simplified map lost some details that were present in the full, complex description, leading to larger errors over longer periods.

The study also examined a specific model used in nuclear physics to understand phase transitions, which are sudden changes in a system's behavior, similar to water freezing into ice. The researchers tested their method as the system moved through this critical transition point. Even in this challenging scenario, where the system's behavior changes drastically, their simplified approach remained accurate over long periods of time. This resilience suggests that the method could be a powerful tool for studying large quantum systems where exact calculations are currently impossible. The researchers noted that while their current simulations were done with a specific number of particles, the nature of their approximation suggests it might become even more accurate as the number of particles increases, potentially making it ideal for studying truly macroscopic quantum systems.

Ultimately, this work provides a practical way to navigate the complexity of many-particle quantum systems. By projecting the full, overwhelming information of a quantum state onto a manageable three-dimensional space, the researchers have created a tool that preserves the essential dynamics of the system. While the method has limits, particularly when the forces acting on the system are irregular, it offers a highly accurate and efficient way to understand how quantum states evolve when the underlying rules are symmetrical. This approach opens the door to analyzing larger and more complex systems than previously possible, bringing the behavior of the quantum many-body world into sharper focus.

Technical Summary: Finite Systems in the Space of Projective Measurements

Problem Statement
Describing quantum systems with a large number of particles ( $N$ ) approaching the macroscopic limit presents significant theoretical and experimental challenges. The complexity of representing such systems grows exponentially with particle number, making standard discrete phase-space representations (e.g., $2^N \times 2^N$ grids) difficult to interpret due to excessive information density. Furthermore, in many macroscopic settings involving indistinguishable particles, individual particle measurements are inaccessible; only collective observables can be measured. Existing methods often fail to efficiently characterize global properties or visualize state evolution in this regime.

Methodology
The authors propose a framework to analyze the dynamics of $N$ -qubit systems using discrete quasidistributions (specifically the Husimi $Q$ and Glauber $P$ functions) projected onto a three-dimensional "space of collective measurements."

Discrete Phase Space: The system is mapped to a discrete phase space defined over the finite field $\mathbb{F}_{2^N}$ . States are associated with functions on a $2^N \times 2^N$ grid using generalized Pauli operators ( $Z_\alpha, X_\beta$ ) and discrete coherent states $|\alpha, \beta\rangle$ .
Projection to Collective Space: To reduce dimensionality, the authors project the full $Q$ function onto a space defined by the Hamming weights of the phase-space coordinates: $m=h(\alpha)$ , $n=h(\beta)$ , and $k=h(\alpha+\beta)$ . This projection, denoted $\tilde{Q}(m,n,k)$ , retains all information necessary to calculate expectation values for permutation-invariant observables. The projection space is a discrete tetrahedron with $M = (N+1)(N+2)(N+3)/6$ points.
Projected Star Product: The core methodological contribution is the construction of an approximate analog to the star product for these projected functions. The authors derive an expression for the projected $Q$ $Q$ function of a product of two operators (specifically a density matrix $\hat{\rho}$ $\overset{ρ}{^}$ and a Hamiltonian $\hat{H}$ $\hat{H}$ ) in terms of their individual projected $Q$ $Q$ and $P$ $P$ functions.
- They introduce a kernel $\Psi_{m,n,k,m',n',k'}^{p,q,r}$ which acts as a convolution operator in the projected space.
- An approximation is made: $Q(\alpha, \beta) \approx R_{m,n,k}^{-1} \tilde{Q}(m,n,k)$ , where $R_{m,n,k}$ is the multiplicity factor. This relation is exact for permutation-invariant states but approximate for others.
Dynamics: The Liouville equation is rewritten in terms of these projected functions. The time evolution is solved numerically by expressing the system as a matrix differential equation involving the kernel $\Psi$ and the projected Hamiltonian function $\tilde{P}_H$ .

Key Contributions

Dimensionality Reduction: The paper successfully maps the high-dimensional discrete phase space to a manageable 3D tetrahedral space of collective measurements, facilitating the visualization of information density and global properties.
Approximate Star Product: The authors derive a specific mathematical formulation (Eq. 3.13–3.17) that approximates the product of operators in the projected space. This allows the rewriting of the Liouville equation without requiring the full $2^N$ -dimensional Hilbert space.
Kernel Construction: A numerical method is provided to calculate the kernel $\Psi$ once for a given $N$ , which can then be reused for any state and Hamiltonian.

Results
The performance of the approximation was evaluated using the Hellinger distance ( $D$ ) between the exact evolution and the approximate evolution in the projected space.

Permutation-Invariant Cases: For states and Hamiltonians that are invariant under particle permutations (e.g., GHZ, W, and coherent states evolving under Harper, Kerr, or LMG Hamiltonians), the method performs exceptionally well. The Hellinger distance is consistently of the order $10^{-5}$ or smaller, indicating negligible information loss.
Non-Invariant States with Invariant Hamiltonians: Surprisingly, the approximation remains highly accurate for non-invariant states (e.g., domain-wall, NOON, singlet states) when evolved under permutation-invariant Hamiltonians. The error remains low ( $D \sim 10^{-5}$ ), suggesting the method is robust even when the state symmetry is broken, provided the dynamics are symmetric.
Non-Invariant Hamiltonians: When the Hamiltonian itself is not permutation-invariant (e.g., Ising models with external fields), deviations arise.
- At short times, the approximation remains close to the exact evolution for various states.
- At longer times, information loss becomes apparent for certain states (e.g., coherent states under specific Ising Hamiltonians), with Hellinger distances increasing significantly.
- However, some non-invariant states (like the NOON state) still perform better than others under non-invariant Hamiltonians, indicating a complex relationship between initial state structure and Hamiltonian symmetry.
Phase Transitions: In the Lipkin-Meshkov-Glick (LMG) model, the approximation performed well across the quantum phase transition point ( $\tilde{g}=1$ ), maintaining low error over long evolution times.

Significance and Claims
The paper claims that the proposed method provides a viable tool for analyzing quantum dynamics in the macroscopic limit by focusing on collective properties. The authors emphasize that:

The method is highly effective for permutation-invariant systems, offering a clear visualization of state evolution that is otherwise obscured in full phase space.
The approximation is unexpectedly robust for non-invariant states when the Hamiltonian is symmetric.
While the method shows deviations for non-invariant Hamiltonians at long times, it remains useful for short-time dynamics and specific state-Hamiltonian combinations.
The error appears to decrease as the particle number $N$ increases, suggesting the approximation may become exact in the thermodynamic limit ( $N \gg 1$ ).
Future work is identified as simplifying the computational construction of the kernel $\Psi$ , which currently scales exponentially with $N$ , to enable studies of larger systems.

More like this