Gauge-Fixing Quantum Density Operators At Scale

We provide theory, algorithms, and simulations of non-equilibrium quantum systems using a one-dimensional (1D) completely-positive (CP), matrix-product (MP) density-operator (ρ) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MPρ factorization naturally emerges. In this work we analytically and numerically examine the virtual freedoms associated with the representation of quantum density operators. Using this augmented perspective, we simplify algorithms in certain limits to integrate the canonical form's master equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale-up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, due to loss from environmental interactions, non-local complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional non-equilibrium settings, beyond qubit engineering.