Godinez Ramirez, EmilianoEmilianoGodinez RamirezMilbradt, Richard M.Richard M.MilbradtMendl, Christian B.Christian B.Mendl2025-11-262025-11-262025-07-28Physical Review A 112 (1): 012224 (2025)https://hdl.handle.net/11420/59097Tensor networks offer a valuable framework for implementing Lindbladian dynamics in many-body open quantum systems with nearest-neighbor couplings. In particular, a tensor network Ansatz known as the locally purified density operator employs the local purification of the density matrix to guarantee the positivity of the state at all times.Within this framework, the dissipative evolution utilizes the Trotter-Suzuki splitting, yielding a second-order approximation error. However, due to the Lindbladian dynamics’ nature, employing higher-order schemes results in nonphysical quantum channels. In this work, we leverage the gauge freedom inherent in the Kraus representation of quantum channels to improve the splitting error. To this end, we formulate an optimization problem on the Riemannian manifold of isometries and find a solution via the second-order trust-region algorithm. We validate our approach using two nearest-neighbor noise models and achieve an improvement of orders of magnitude compared to other positivity-preserving schemes. In addition, we demonstrate the usefulness of our method as a compression scheme, helping to control the exponential growth of computational resources, which thus far has limited the use of the locally purified Ansatz.en1094-162220251American Physical Societyhttps://creativecommons.org/licenses/by/4.0/Natural Sciences and Mathematics::530: PhysicsNatural Sciences and Mathematics::519: Applied Mathematics, ProbabilitiesJournal Articlehttps://doi.org/10.15480/882.1624810.1103/kztw-ywxr10.15480/882.16248Journal Article